/* ch 6 2-way additive model....may be inappropriate */
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
;
proc glm;
classes sizecat fingcat;
model htime=sizecat fingcat;
means sizecat fingcat/bon cldiff alpha=.10;
lsmeans sizecat fingcat/pdiff=all cl adjust=bon alpha=.10;
contrast 's_vs_m' sizecat 1 -1 0;
estimate 't_vs_b' fingcat 1 0 -1;
run; quit;
/* ch6 2-way additive model ... may be inappropriate */
/* next : in class include one new line and must use lsmeans*/
/*  l 3 b 78 3 44     */
/* lsmeans sizecat fingcat/pdiff=all cl adjust=bon alpha=.10; */
/* the means statement is wrong .  use lsmeans */