/* ch 6 one obs per cell -- orthogonal contrasts */
options ls=72;
data beads;
  INPUT sizelet $ sizecat finglet $ finglnth fingcat time a1 a2 b1 b2;
  htime=log10(tan((time*3.14159/180)/2)-.3);
  code=sizecat*10+fingcat;
  letcode=trim(sizelet)||trim(finglet);
   a1b1=a1*b1;
   a1b2=a1*b2;
   a2b1=a2*b1;
   a2b2=a2*b2;
  LINES;
   s 1 t 64 1 100 -1 1 -1 1
   s 1 a 72 2 140  -1 1 0 -2
   s 1 b 76 3 93  -1 1 1 1
   m 2 t 70 1 85 0 -2 -1 1
   m 2 a 71 2 45  0 -2 0 -2
   m 2 b 79 3 105  0 -2 1 1
   l 3 t 70 1 40  1 1 -1 1
   l 3 a 75 2 37  1 1 0 -2
   l 3 b 77 3 35  1 1 1 1
  ;
proc print;
proc glm;
model htime= a1 a2 b1 b2 a1b1 a2b2;
run; quit;