* wool.sas, clean wool experiment, Table 17.12 (page 634),
* in text code (page 632);
* (see also Figure 17.3 (page 633) and Figure 17.4 (page 634));
options ls=72;
;
DATA WOOL;
  INPUT BALE CONTENT;
  LINES;
   1 52.33
   1 56.26
   1 62.86
   1 50.46
   2 56.99
   2 58.69
   2 58.20
   2 57.35
   3 54.64
   3 57.48
   3 59.29
   3 57.51
   4 54.90
   4 60.08
   4 58.72
   4 55.61
   5 59.89
   5 57.76
   5 60.26
   5 57.53
   6 57.76
   6 59.68
   6 59.58
   6 58.08
   7 60.27
   7 60.30
   7 61.09
   7 61.45
;
PROC PRINT;
;
* Generate Figure 17.3;
PROC GLM;  
  CLASS BALE;  
  MODEL CONTENT = BALE;  
  OUTPUT OUT=STATS P=PRED R=Z;  
PROC STANDARD STD=1.0 MEAN=0.0;  
  VAR Z;  
PROC RANK NORMAL=BLOM;  
  VAR Z;  
  RANKS NSCORE;  
PROC PLOT;  
  PLOT Z*PRED=BALE / VPOS=15 HPOS=50 VREF=0; 
;
* Generate Figure 17.4;
DATA WOOL2; SET WOOL;
  IF CONTENT NE 62.86;
PROC GLM;  
  CLASS BALE;  
  MODEL CONTENT = BALE;  
  OUTPUT OUT=STATS P=PRED R=Z;  
PROC STANDARD STD=1.0 MEAN=0.0;  
  VAR Z;  
PROC RANK NORMAL=BLOM;  
  VAR Z;  
  RANKS NSCORE;  
PROC PLOT;  
  PLOT Z*PRED=BALE / VPOS=15 HPOS=50 VREF=0; 
;
* Table 17.12 program, generating Figure 17.5;
PROC GLM DATA=WOOL; 
  CLASS BALE; 
  MODEL CONTENT = BALE; 
  OUTPUT OUT=STATS P=PRED; 
PROC STANDARD STD=1.0 MEAN=0.0; 
  VAR PRED; 
PROC RANK NORMAL=BLOM; 
  VAR PRED; 
  RANKS NSCORE; 
PROC PLOT; 
  PLOT PRED*NSCORE=BALE / VPOS=15 HPOS=40 HREF=0 VREF=0;