OPTIONS LS=72;
TITLE 'NUTRITION DEFICIENCY P.112, GENERALIZED LINEAR MODEL';
DATA NUTRIDEF;
   INPUT WGT HGT AGE;
   LABEL WGT='WEIGHT'
         HGT='HEIGHT';
      CARDS;
64 57 8
71 59 10
53 49 6
67 62 11
55 51 8
58 50 7
77 55 10
57 48 9
56 42 10
51 42 6
76 61 12
68 57 9
;
PROC PRINT;

DATA NUTRIDF1;
   SET NUTRIDEF;
   AGE2=AGE*AGE;

PROC CORR DATA=NUTRIDF1;

PROC GLM DATA=NUTRIDF1;
      MODEL WGT=HGT AGE AGE2 /P CLM;
      OUTPUT OUT=REGOUT P=PRED_WGT R=RESID ;
PROC PLOT;
     PLOT WGT*HGT='A' PRED_WGT*HGT='P'/OVERLAY;
     PLOT RESID*PRED_WGT;