options ls=75 nodate;

data one;
   infile "sleeping.dat" firstobs=5;
   input subject drinking $ 25-34 reading $ 39-48 sleepqual;
run;
proc print;run;

proc gplot;
   plot sleepqual*drinking=reading;
run;


proc glm data=one;
   class subject drinking reading;
   model sleepqual=drinking|reading subject(drinking);
   *lsmeans drinking|reading; 
   lsmeans drinking*reading/slice=drinking; 
   *random subject(drinking)/test;
run;
proc mixed method=type3;
   class subject drinking reading;
   model sleepqual=drinking|reading; 
   random subject(drinking);
   lsmeans drinking*reading/slice=drinking;
   lsmeans drinking|reading/diff;
run;