全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
1619 6
2014-12-21
数据分析(第2版)_范金城2010,最后一页有一个计算二元回归的proc iml例子,怎么运行没有任何反应呀?
我对编程一无所知,只是出于好奇,录入了全部句子,结果发现运行没结果。提示如下:
=====================
8    x={1 274 2450,1 180 3254,1 375 3802,1 205 2838,1 86 2347,1 265 3782,1 98 3008,1 330 2450,18  ! 195 2137,1 53 2560,1 430 4020,1 372 4427,1 236 2660,1 157 2088,1 370 2605};9    y={162,120,323,131,67,169,81,192,116,55,252,232,144,103,212};10   a=inv(x'            -            22            7610 ! *x);ERROR 22-322: Syntax error, expecting one of the following: #, ##, &, (, (|, *, **, +, ',', -,              /, //, :, <, <=, <>, =, >, ><, >=, @, [, ^=, |, ||.ERROR 76-322: Syntax error, statement will be ignored.=========================



proc iml;
x={1 274 2450,1 180 3254,1 375 3802,1 205 2838,1 86 2347,1 265 3782,1 98 3008,1 330 2450,1 195 2137,1 53 2560,1 430 4020,1 372 4427,1 236 2660,1 157 2088,1 370 2605};
y={162,120,323,131,67,169,81,192,116,55,252,232,144,103,212};
a=inv(x'*x);
beta=a*x'*y;
print beta;
h=x*a*x';
yhat=h*y;
sse=(y-yhat)'*(y-yhat);
sigma=sse/12;
print yhat sse sigma;
mssr=y*(h-(1/15)#j(15,15,1))*y;
f=mssr/sigma;
fpvalue=1-probf(1,2,12);
print f fpvalue;
varbeta=sigma#a;
do k=1 to3;
t=beta[k,1]/sqrt(varbeta[k,k]);
abst=abs(t);
tpvalue=2#probt(-abst,12);
print t,tpvalue;
end;
quit;


二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2014-12-21 12:39:44
这逗号是小数点?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-12-21 15:17:26
teqel 发表于 2014-12-21 12:39
这逗号是小数点?
我不清楚你说的是哪个地方?
附件列表
1.jpg

原图尺寸 38.83 KB

1.jpg

2.jpg

原图尺寸 70.02 KB

2.jpg

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-12-21 15:33:32
是我看错了。也有可能你从书上copy过来,其中的'是中文符号里面的,程序不认。试试看删掉粘帖过来的  ’ ,自己输一个。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-12-21 15:36:13
teqel 发表于 2014-12-21 15:33
是我看错了。也有可能你从书上copy过来,其中的'是中文符号里面的,程序不认。试试看删掉粘帖过来的  ’ , ...
我贴的那些命令本就可以复制的啊。是自己在英文状态下自己输入的,不是拷贝来的。
谢谢你。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2014-12-21 23:05:17
proc iml;
x={1 274 2450,1 180 3254,1 375 3802,1 205 2838,1 86 2347,1 265 3782,1 98 3008,1 330 2450,1 195 2137,1 53 2560,1 430 4020,1 372 4427,1 236 2660,1 157 2088,1 370 2605};
y={162,120,323,131,67,169,81,192,116,55,252,232,144,103,212};
a=inv(x`*x);
beta=a*x`*y;
print a beta;
h=x*a*x`;
yhat=h*y;
sse=(y-yhat)`*(y-yhat);
sigma=sse/12;
print yhat sse sigma;
mssr=y`*(h-(1/15)#j(15,15,1))*y;
f=mssr/sigma;
fpvalue=1-probf(1,2,12);
print f fpvalue;
varbeta=sigma#a;
do k=1 to 3;
t=beta[k,1]/sqrt(varbeta[k,k]);
abst=abs(t);
tpvalue=2#probt(-abst,12);
print t,tpvalue;
end;
quit;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群