IML中循环使用矩阵的赋值失败问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

向各位大侠求救：
其实我只是想用循环来做筛选A矩阵中每一列大于30的元素的和。

我这样写的时候系统提示错误说是没有给b矩阵赋值，不知道如何修改。
proc iml;
a={1 2 50 60,
        2 50 60 10};
n=0;
do while(n<10);
n=n+1;
k=ncol(a);
do i=1 to k;
b=loc(a[,i]>30);
c=sum(a[b,i]);
print a b c ;
end;
end;
quit;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：print while 各位大侠 quit End 如何 元素

我自己解决了：
由于在a矩阵的第一列所有的元素都是小于30的，这时候loc(a[,i]>30)无法返回数值。所以就会说b矩阵没有被赋值了
然后我用那个any函数来判断后再做计算便可以解决问题。
proc iml;
a={1 2 50 60,
   2 50 60 10};
n=0;
do while(n<10);
n=n+1;
k=ncol(a);
do i=1 to k;
b=any(a[,i]>30);
if b=1 then do;
d=loc(a[,i]>30);
c=sum(a[d,i]);
end;
else do;
d=loc(a[,i]<30);
c=sum(a[d,i]);
end;
print a b c ;
end;
end;
quit;

学习了

另外你这个问题如果不用IML
用SQL的话会很简单
因为你是求列的嘛

如果求行的话，用array也可以

playmore 发表于 2013-9-11 12:27
学习了

另外你这个问题如果不用IML

恩是的，用SQL的话会很简单。其实这个只是一个程序的开头，后面还要做其他的矩阵计算。
一开头就被卡着的感觉真不爽。主要还是IML用得少，对于它的函数运行结果不是很熟悉导致的。