OPTIONS LS=72;
TITLE 'EXAMPLE 9. POLYNOMIAL REGRESSION: FISH LENGTH DATA';
DATA FISH;
INPUT X Y1 Y2 Y3 Y4;
X2=X*X;
X3=X*X*X;
CARDS;
 14    620    625    590    590
 21    910    820    910    910
 28   1315   1215   1305   1205
 35   1635   1515   1730   1605
 42   2120   2110   2140   1915
 49   2300   2320   2725   2035
 56   2600   2805   2890   2140
 63   2925   2940   3685   2520
 70   3110   3255   3920   2710
 77   3315   3620   4325   2870
 84   3535   4015   4410   3020
 91   3710   4235   4485   3025
 98   3935   4315   4515   3030
105   4145   4435   4480   3025
112   4465   4495   4520   3040
119   4510   4475   4545   3177
126   4530   4535   4525   3180
133   4545   4520   4560   3180
140   4570   4600   4565   3257
147   4605   4600   4626   3166
154   4600   4600   4566   3214
;
PROC PRINT DATA=FISH;
RUN;

PROC REG DATA=FISH;
MODEL Y3=X/R P SS1 SS2;
OUTPUT OUT=OUT1 P=YPRED R=YRES;
RUN;
PROC PLOT DATA=OUT1;
PLOT Y3*X;
PLOT YRES*YPRED;
PLOT YRES*X;
RUN;
PROC UNIVARIATE DATA=OUT1 PLOT FREQ NORMAL;
VAR YRES;
RUN;

PROC REG DATA=FISH;
MODEL Y3=X X2/R P SS1 SS2;
OUTPUT OUT=OUT2 P=YPRED R=YRES;
RUN;
PROC PLOT DATA=OUT2;
PLOT YRES*YPRED;
PLOT YRES*X;
RUN;
PROC UNIVARIATE DATA=OUT2 PLOT FREQ NORMAL;
VAR YRES;
RUN;

PROC REG DATA=FISH;
MODEL Y3=X X2 X3/R P SS1 SS2;
OUTPUT OUT=OUT3 P=YPRED R=YRES;
RUN;
PROC PLOT DATA=OUT3;
PLOT YRES*YPRED;
PLOT YRES*X;
RUN;
PROC UNIVARIATE DATA=OUT3 PLOT FREQ NORMAL;
VAR YRES;
RUN;