/* ch 6 2-way cell means model... note additional data */
options ls=72;
DATA beads;
  INPUT sizelet $ sizecat finglet $ finglnth fingcat time;
  htime=log10(tan((time*3.14159/180)/2)-.3);
  code=sizecat*10+fingcat;
  letcode=trim(sizelet)||trim(finglet);
  LINES;
   s 1 t 64 1 100
   s 1 a 72 2 140
   s 1 b 76 3 93
   m 2 t 70 1 85
   m 2 a 71 2 45
   m 2 b 79 3 105
   l 3 t 70 1 40
   l 3 a 75 2 37
   l 3 b 77 3 35
   s 1 t 63 1 119
   s 1 a 72 2 127
   s 1 b 79 3  70
   m 2 t 69 1  49
   m 2 a 71 2 133
   m 2 b 80 3  86
   l 3 t 65 1  44
   l 3 a 75 2  35
   l 3 b 77 3  53
;
proc glm;
  classes code;
  model htime = code;
 means code/bon cldiff alpha=.10;
run; quit;
/* means ok here no need for lsmeans since like a 1-way */