求解sas中的iml模块问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

题目是SAS/IML编写一个程序，计算岭参数k从0到1.000（每次变化量为0.001）时岭估计，一直运行不出来
代码proc iml;
start estimate(x,y,z);
   b=(inv(x`*x+z))*x`*y;
   return(b);
finish;
I={1,1,1,1,1,1,1,1,1,1,1};
X={1 149.3 4.2 108.1,
   1 161.2 4.1 114.8,
   1 171.5 3.1 123.2,
   1 175.5 3.1 126.9,
   1 180.8 1.1 132.1,
   1 190.7 2.2 137.7,
   1 202.1 2.1 146.0,
   1 212.4 5.6 154.1,
   1 226.1 5.0 162.3,
   1 231.9 5.1 164.3,
   1 239.0 0.7 167.6};
Y={15.9,16.4,19,19.1,18.8,20.4,22.7,26.5,28.1,27.6,26.3};
do k=0 to 1.000 by 0.001;
Z=k*I;

beta(k)=estimate(X,Y,Z);
end;
print beta(k);
run;
quit;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：estimate RETURN finish print Start SAS iml

qiaoluo 发表于 2010-4-1 09:25
题目是SAS/IML编写一个程序，计算岭参数k从0到1.000（每次变化量为0.001）时岭估计，一直运行不出来
代码proc iml;
start estimate(x,y,z);
   b=(inv(x`*x+z))*x`*y;
   return(b);
finish;
I={1,1,1,1,1,1,1,1,1,1,1};
X={1 149.3 4.2 108.1,
   1 161.2 4.1 114.8,
   1 171.5 3.1 123.2,
   1 175.5 3.1 126.9,
   1 180.8 1.1 132.1,
   1 190.7 2.2 137.7,
   1 202.1 2.1 146.0,
   1 212.4 5.6 154.1,
   1 226.1 5.0 162.3,
   1 231.9 5.1 164.3,
   1 239.0 0.7 167.6};
Y={15.9,16.4,19,19.1,18.8,20.4,22.7,26.5,28.1,27.6,26.3};
do k=0 to 1.000 by 0.001;
Z=k*I;

beta(k)=estimate(X,Y,Z);
end;
print beta(k);
run;
quit;

A quick fix and let you you figure it the rest.

proc iml;
start estimate(x,y,z);
   b=(inv(x`*x+z))*x`*y;
   return(b);
finish;
I={1,1,1,1,1,1,1,1,1,1,1};
X={1 149.3 4.2 108.1,
   1 161.2 4.1 114.8,
   1 171.5 3.1 123.2,
   1 175.5 3.1 126.9,
   1 180.8 1.1 132.1,
   1 190.7 2.2 137.7,
   1 202.1 2.1 146.0,
   1 212.4 5.6 154.1,
   1 226.1 5.0 162.3,
   1 231.9 5.1 164.3,
   1 239.0 0.7 167.6};
Y={15.9,16.4,19,19.1,18.8,20.4,22.7,26.5,28.1,27.6,26.3};
do k=0 to 1.000 by 0.01;
Z=k*I(4);
beta=estimate(X,Y,Z);
print k beta;
end;

run;
quit;

谢谢，后来解决了