作者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