options nocenter ps=60 ls=78;
goptions colors=(none);
symbol1 v=circle; axis1 offset=(5); 

data wash;
 input stain soap y @@;
 if y=37 then y=.;
 cards;
1 1 45 1 2 47 1 3 48 1 4 42
2 1 43 2 2 46 2 3 50 2 4 37
3 1 51 3 2 52 3 3 55 3 4 49
;

proc glm;
 class stain soap;
 model y=soap stain;
 output out=diag p=pred r=res;
 means soap / lsd lines;
 lsmeans soap / stderr;

data new1;
 set wash;
 if y=. then y=42.166666666;

proc glm;
 class stain soap;
 model y=soap stain;
 output out=diag p=pred r=res;
 means soap / lsd lines;
run;