请问这段程序怎么运行没有任何反应呀

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

数据分析(第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;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：何反应 following statement Statemen pvalue

这逗号是小数点？

teqel 发表于 2014-12-21 12:39
这逗号是小数点？

我不清楚你说的是哪个地方？

是我看错了。也有可能你从书上copy过来，其中的'是中文符号里面的，程序不认。试试看删掉粘帖过来的  ’ ，自己输一个。

teqel 发表于 2014-12-21 15:33
是我看错了。也有可能你从书上copy过来，其中的'是中文符号里面的，程序不认。试试看删掉粘帖过来的  ’ ， ...

我贴的那些命令本就可以复制的啊。是自己在英文状态下自己输入的，不是拷贝来的。
谢谢你。

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;

jl60156 发表于 2014-12-21 23:05
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, ...

真心感谢！经您这么细心查看并修改后，就一切正常了！
向您学习！