options nocenter ls=75 ps=60;

data new;
infile "u:\.www\datasets1\sunscreen.dat";
input subject lotion resp;

proc mixed covtest cl maxiter=20;
 class subject lotion;
 model resp=lotion /ddfm=kr;
 random subject subject*lotion;
 lsmeans lotion / diff cl;
run;
quit;