/* this example taken from Kuehl */

options ls=75;
data one;
   input method @;
   do day=1 to 4;
      do rep=1 to 2;
         input y @;
         output;
      end;
   end;
cards;
1 142.3 144.0 134.9 146.3 148.6 156.5 152.0 151.4
2 142.9 147.4 125.9 127.6 135.5 138.9 142.9 142.3 
;
run;

proc print;run;

/*proc glm;
   class day method;
   model y=day|method;
   random day day*method;
   test h=day e=day*method;
   test h=method e=day*method;
   means method;
   means day;
   means day*method;
run;

proc varcomp;
   class day method;
   model y=method|day/fixed=1;
run;*/

proc mixed cl method=type3;  /* will symmetric c.i.'s Type 3 ss*/
*proc mixed cl;  /* will report Wald c.i.'s using Satterthwaite */
   class day method;
   model y=method/s cl ddfm=satterth;
   random day day*method;
   *lsmeans method/stderr tdiff cl;  /*stderr is not a valid option */
   /* a fixed statement is given below */
   lsmeans method/diff cl;
   estimate "avg. method effect" method -1 1; /* redundant */
run;