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

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 ;

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

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

阅读文章"计算 SAS 中的第 k 个最小数据值"定义模块，如本文所示。然后使用以下命令：

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;

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

Answer 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;