/* ch 5 with untransformed */
options ls=72;
data beads;
  input size finglnth fingcat time;
  LINES;
   1 64 1 100
   1 72 2 140
   1 76 3 93
   2 70 1 85
   2 71 2 45
   2 79 3 105
   3 70 1 40
   3 75 2 37
   3 77 3 35
;
proc print;
PROC GLM;
  CLASSES size   ;
  MODEL time   = size  ;
output out=beads2 predicted=ypred residual=z;
  Means size;
proc standard std=1.0 data=beads2;
     var z;
proc rank normal=blom;
 var z;
 ranks nscore;
proc print;
proc plot;
 plot z*size z*ypred z*fingcat/vref=0 vpos=19 hpos=50;
 plot z*nscore/vref=0 href=0 hpos=50;
run; quit;
/* ch 5 untransformed */