* tempr.sas, temperature experiment, Table 17.13 (page 636),
* in text code (page 636);
options ls=72;
;
DATA TEMPR; 
  INPUT THERM SITE SUBJ TIME; 
  LINES;
   1 1 1  62.16
   1 2 1  61.53
   2 1 1 154.42
   2 2 1 310.46
   3 1 1  95.98
   3 2 1 225.65
   1 1 2  65.63
   1 2 2  63.70
   2 1 2 132.30
   2 2 2 284.64
   3 1 2  98.50
   3 2 2 241.63
   1 1 3  63.12
   1 2 3  61.34
   2 1 3 105.52
   2 2 3 315.61
   3 1 3 110.05
   3 2 3 364.07
   1 1 4  61.51
   1 2 4  61.54
   2 1 4  94.88
   2 2 4 294.16
   3 1 4 107.93
   3 2 4 304.58
;
* Generate Table 17.14; 
* Note that the option TEST gives correct denominators; 
PROC PRINT; 
PROC GLM; 
  CLASSES THERM SITE SUBJ; 
  MODEL TIME = SUBJ THERM SITE THERM*SITE SUBJ*SITE; 
  RANDOM SUBJ SUBJ*SITE / TEST; 
  OUTPUT OUT=RESIDS PREDICTED=PTIME RESIDUAL=Z; 
  CONTRAST 'SITE1-SITE2' SITE  1 -1 / E=SUBJ*SITE;   
;
* Other SAS procedures;
PROC VARCOMP  METHOD=TYPE1; 
  CLASSES THERM SITE SUBJ; 
  MODEL TIME = THERM SITE THERM*SITE SUBJ SUBJ*SITE / FIXED=3;