inverse a square matrix in data step

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

Several call routine in SAS/IML are now can be defined in FCMP and used in data step and other proccedures. Here is an example for inverse a matrix(array) in data step.

proc fcmp outlib=sasuser.funcs.temp;;
     subroutine invert (mat[*,*],invmat[*,*] );
      outargs invmat;
      call inv(mat, invmat);
      endsub;
run;
quit;

options cmplib=sasuser.funcs;

data _null_;
  array a(3,3) _temporary_   ( 1 2 3,
                               2 2 4,
                               3 4 5
                             );
  array b(3,3) _temporary_ ;
  call invert (a, b);

  do i=1 to dim1(b);
    do j=1 to dim2(b);
      put b[i,j] @;
    end;
    put;
  end;
run;

***check with iml;
proc iml;
a={ 1 2 3,
     2 2 4,
     3 4 5
    };
b=inv(a);
print b;
quit;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Data step inverse Square matrix Data mat null