options nocenter ls=78; goptions colors=(none);

data wash;
 input stain soap 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;
 random stain;
 lsmeans soap / stderr tdiff;
 
proc mixed;
 class stain soap;
 model y = soap;
 random stain;
 lsmeans soap / diff;
run;