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

【问题标题】：How to select the 5 minimum values with SAS Proc IML?如何使用 SAS Proc IML 选择 5 个最小值？
【发布时间】：2015-02-04 00:11:38
【问题描述】：

我想知道是否可以使用 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 ;

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

如果你有任何想法！ :) 如果不是 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

【问题讨论】：

我认为您被否决了，因为您的问题有一些非常通用的标签。我已经为您删除了这些标签，以便将这个问题保留在那些遵循“sas”标签的人中......
我不完全确定您在寻找什么。一个包含数据和预期输出的示例会很有帮助。
我用一个例子编辑了这个问题，如果不清楚，对不起

标签： sas sas-iml

【解决方案1】：

阅读文章"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 模块。

【讨论】：

【解决方案2】：

您可以在 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;

【讨论】：