data mod92;
input y1 xn1 xp1 xk1 @@;
n=_n_;
xn11=xn1**2; xp11=xp1**2; xk11=xk1**2; y11=log(y1);
cards;
15.18 0 196 372
21.36 34 196 372
25.72 67 196 372
32.29 101 196 372
34.03 135 196 372
39.45 202 196 372
43.15 259 196 372
43.46 336 196 372
40.83 404 196 372
30.75 471 196 372
33.46 259 0 372
32.47 259 24 372
36.06 259 49 372
37.96 259 73 372
41.04 259 98 372
40.09 259 147 372
41.26 259 196 372
42.17 259 245 372
40.36 259 294 372
42.73 259 342 372
18.98 259 196 0
27.35 259 196 47
34.86 259 196 93
38.52 259 196 140
38.44 259 196 186
37.73 259 196 279
38.43 259 196 372
43.87 259 196 465
42.77 259 196 558
65.22 259 196 651
;
proc reg data=mod92 outest=a press;
model y1=xn1 xp1 xk1 / r;
output out=om predicted=pred residual=resid
student=student h=leverage;
run;
proc sort data=om;
by y1;
run;
proc print data=om;
var n y1 xn1 xp1 xk1 pred resid;
run;
proc gplot data=om;
plot student*pred;
run;
proc sgplot data=om;
scatter x=pred y=student / markerattrs=(symbol=asterisk)
datalabel=n;
run;
好多命令不知道什什么意思,希望大家多多帮忙