如何提取矩阵元素

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

corr.zip (9.84 KB) 本附件包括：

• corr.sas7bdat

2012-12-10 13:50:20 上传

数据文件

对于一个矩阵，如何把对角线1以下的下三角矩阵元素挑出来，并且按行排列？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：矩阵元素 如何提取 对角线 如何 元素 对角线

1. proc iml;                /* 重新排列相关系数矩阵 */
   
2.         use CorrMatrix;
   
3.         read all into CorrMatrix;
   
4.                         CorrPair=j(N*(N+1)/2,1,.);
   
5.         k=1;
   
6.         %do i=1 %to N;
   
7.                 %do j=1 %to &i;
   
8.                 CorrPair[k]=CorrMatrix[&i,&j];
   
9.                 k=k+1;
   
10.                 %end;
    
11.         %end;
    
12.         corr_name='Corr';
    
13.         create CorrList from CorrPair[colname=corr_name];
    
14.         append from CorrPair;
    
15.         close CorrList;
    
16. quit;

复制代码

proc iml;
a={1 2 3 ,1 2 3 , 1 2 3 };
temp1=1:ncol(a);
temp2=(1:nrow(a))`;
col=repeat(temp1,ncol(a),1);
row=repeat(temp2,1,nrow(a));
a[loc(row<col)]=j(1,ncol(loc(row<col)),0);
b=a[loc(row>=col)];
print a b;
quit;