作者waterage (熱情洋溢~YA!)
看板PH-91
標題高統作業上週SAS程式
時間Fri Apr 30 02:19:24 2004
同學以下是上週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
1F:→ ago:謝謝助教!!!!!!!!!!!!!!!!!! 推 210.85.66.70 04/30
2F:→ lyr19:謝謝~ 昨天生統網頁打不開 只好自己打T_T 推 218.166.28.252 04/30
3F:→ banturtle:天阿.助教真是天使@@ 推 140.112.7.59 04/30