/* ch 11 balanced incomplete block*/
options ls=72;
DATA beads;
  INPUT blocktod size time;
  LINES;
1.00 1.00 60.00
1.00 2.00 50.00
1.00 3.00 40.00
2.00 1.00 65.00
2.00 2.00 54.00
2.00 4.00 35.00
3.00 1.00 58.00
3.00 3.00 38.00
3.00 4.00 31.00
4.00 2.00 52.00
4.00 3.00 44.00
4.00 4.00 33.00
;
proc glm;
 classes blocktod size;
 model time=blocktod size ;
 lmeans blocktod;
 lmeans size   / pDIFF=all cl adjust=scheffe ALPHA=0.10;
 /* means only ok with equal ss */
run; quit;