* ice.sas, ice cream experiment, in text code (page 638),
* generating Figure 17.6 (page 638);
* (Section 17.10.2, on covariates);
options ls=72;
;
DATA ICE; 
  INPUT FLAVOR MELTTIME ORDER; 
  X=ORDER-16.5; 
  X2=X*X; 
  LINES;
   1  924  1
   1  876  2
   1 1150  5
   1 1053  7
   1 1041 10
   1 1037 12
   1 1125 15
   1 1075 16
   1 1066 20
   1  977 22
   1  886 25
   2  891  3
   2  982  4
   2 1041  8
   2 1135 13
   2 1019 14
   2 1093 18
   2  994 27
   2  960 30
   2  889 31
   2  967 32
   2  838 33
   3  817  6
   3 1032  9
   3  844 11
   3  841 17
   3  785 19
   3  823 21
   3  846 23
   3  840 24
   3  848 26
   3  848 28
   3  832 29
;
PROC GLM; 
  CLASS FLAVOR; 
  MODEL MELTTIME = X X2 FLAVOR; 
  RANDOM FLAVOR / TEST;  
;
* Generate Figure 17.6;
PROC GLM;
  CLASS FLAVOR;
  MODEL MELTTIME = X FLAVOR;
  OUTPUT OUT=STATS RESIDUAL=Z;
PROC STANDARD STD=1.0;
  VAR Z;
PROC PLOT;
  PLOT Z*ORDER / VREF=0 VPOS=15 HPOS=40; *44;