IML 如何在矩阵中新生成n列

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

初学IML，想在矩阵中产生新的一列，用来计算第一列与第二列对应元素的差，如果原矩阵中不存在第三列直接计算，
会报错 ERROR: (execution) Invalid subscript or subscript out of range. 我用x3=0新建了一列可以解决，以下是程序：

1. data a;
   
2. input x1 x2;
   
3. cards;
   
4. 1 2
   
5. 9 5
   
6. 5 6
   
7. 6 8
   
8. ;
   
9. /*不先生成x3会报错*/
   
10. data a;
    
11. set a;
    
12. x3=0;
    
13. run;
    
14. /*****************/
    
15. proc iml;
    
16. use work.a;
    
17. read all into b;
    
18. reset print;
    
19. b[,3]=b[,1]-b[,2];
    
20. quit;

复制代码

但是由于需要计算得出的列比较多，每列都用类似于x3=0的先声明太麻烦，有没有简单的办法能让矩阵产生新的一列，或者生成新的一个矩阵用于计算也行（比如生成一个m*n的零矩阵）？
先谢过了！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：EXECUTION invalid script scrip RESET 矩阵 execution ERROR 如何

b=b||(b[,1]-b[,2]);

Do it simply as,

    data a;
    input x1 x2;
    cards;
    1 2
    9 5
    5 6
    6 8
    ;
    /*不先生成x3会报错*/
    data a;
    set a;
    run;
    /*****************/
    proc iml;
    use work.a;
    read all into b;
    *reset print;
    *b[,3]=b[,1]-b[,2];

        x3=b[ ,+];
        print b x3;
    quit;

bobguy 发表于 2012-7-6 09:51
Do it simply as,

    data a;

  原来x3=b[ ,+];可以，谢谢LS