Question

我想知道是否可以使用IML按行选择5个最小值或最大值？

这是我的代码：

Proc iml ;
  use table;
  read all var {&varlist} into matrix ; 

n=nrow(matrix) ; /* n=369 here*/
p=ncol(matrix);  /* p=38 here*/

test=J(n,5,.) ; 

Do i=1 to n ;
    test[i,1]=MIN(taux[i,]);
End;
Quit ;

所以我想获得一个矩阵测试，其中包含第一列的最大最小值，然后第二列包含我的行的最小值，除了第一个值等等...

如果您有任何想法！ :) 事件，如果不是IML（但使用SAS：base，sql ..）

例如：

    Data test;                                                                                                                input x1-x10 ;                                                                                                              cards;     
1 9 8 7 3 4 2 6
9 3 2 1 4 7 12 -2
;run;

我想获得按行排序的结果：

1 2 3 4 6 7 8 9
-2  1 2 3 4 7 12

为了在另一个表中选择我的5个最小值：

y1 y2 y3 y4 y5
1  2 3 4  6
-2 1 2 3 4

Answer 1

您可以在PROC IML中使用call sort()对列进行排序。因为要分离列而不对整个矩阵进行排序，请提取列，对其进行排序，然后更新原始列。

您想对行进行排序，因此转置矩阵，进行排序，然后转置回来。

proc iml;

have = {1 9 8 7 3 4 2 6,
9 3 2 1 4 7 12 -2};

print have;

n = nrow(have);

have = have`; /*Transpose because sort works on columns*/

do i=1 to n;
    tmp = have[,i];
    call sort(tmp,1);
    have[,i]=tmp;
end;

have = have`;

want = have[,1:5];
print want;
quit;

Answer 2

阅读文章"Compute the kth smallest data value in SAS" 像文章中一样定义模块。然后使用以下内容：

have = {1 9 8 7 3 4 2 6,
        9 3 2 1 4 7 12 -2};
x = have`;    /* transpose */

ord = j(5,ncol(x));
do j = 1 to ncol(x);
   ord[,j] = ordinal(1:5, x[,j]);
end;
print ord;

如果数据中缺少值并想要排除它们，请使用SMALLEST模块而不是ORDINAL模块。

如何使用SAS Proc IML选择5个最小值？

2 个答案: