高統作業上週SAS程式

[waterage (熱情洋溢~YA!)]

同學以下是上週SAS程式
直接copy至SAS就可以跑了
我有稍微對程式作一下註解,
希望你們可以知道程式的意思
y1 y2 y3 y4 y5的值是分別對Y做轉換後的值
如果程式有什麼不懂的地方再問我吧

                                        高統助教


**********************************************************

data a;
input x y y1 y2 y3 y4 y5;
cards;
0       98      3.957   6.259   9.899   15.658  24.766
1       135     4.356   7.114   11.619  18.976  30.991
2       162     4.601   7.653   12.728  21.169  35.210
3       178     4.733   7.946   13.342  22.400  37.609
4       221     5.050   8.665   14.866  25.506  43.760
5       232     5.124   8.835   15.232  26.260  45.273
6       283     5.439   9.566   16.823  29.585  52.030
7       300     5.535   9.791   17.321  30.639  54.198
8       374     5.914   10.694  19.339  34.972  63.243
9       395     6.011   10.930  19.875  36.138  65.708
;
run;
title 'λ=1';
proc reg;
/*對X、Y作迴歸 */
model y=x;
/*對residual與fitted的Y值(Yhat)畫圖；畫出殘差的常態機率*/
plot r.*pred. r.*npp.;
/*印出predicted的Y值與residual值(Y-Yhat)*/
print p r
run;




title 'λ=0.3';
proc reg;
model y1=x;
plot r.*pred. r.*nqq.;
print p r;
run;

title 'λ=0.4';
proc reg;
model y2=x;
plot r.*pred. r.*nqq.;
print p r;
run;

title 'λ=0.5';
proc reg;
model y3=x;
plot r.*pred. r.*nqq.;
print p r;
run;

title 'λ=0.6';
proc reg;
model y4=x;
plot r.*pred. r.*nqq.;
print p r;
run;

title 'λ=0.7';
proc reg;
model y5=x;
plot r.*pred. r.*nqq.;
print p r;
run;




--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 61.59.172.52
→ ago:謝謝助教!!!!!!!!!!!!!!!!!!
→ lyr19:謝謝～ 昨天生統網頁打不開 只好自己打T_T
→ banturtle:天阿.助教真是天使@@