如何将6×6的SAS数据集旋转90度？求Proc IML替代实现方法

需求说明 #

现有一个6行6列的SAS数据集have，需要将其对应的矩阵顺时针旋转90度，得到指定格式的输出（或数据集）。

原数据集代码 #

data have;
  input v1-v6;
  cards;
1 1 1 0 0 0
1 1 0 0 0 0
1 0 1 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1
;
run;

期望输出 #

0 0 0 1 1 1
0 0 0 0 1 1
0 0 0 1 0 1
0 0 1 0 0 0
0 1 0 0 0 0
1 0 0 0 0 0

PROC IML实现思路与代码 #

PROC IML是SAS的矩阵编程语言，天生适合处理这类矩阵变换需求，无需手动编写循环赋值，代码更简洁高效。

实现逻辑 #

顺时针旋转90度的矩阵运算等价于两个步骤：

1. 对原矩阵进行转置（行变列，列变行）
2. 将转置后的矩阵每行元素逆序排列（即逆序选取列）

完整代码 #

proc iml;
    /* 读取数据集`have`的所有数值变量到矩阵A */
    use have;
    read all var _NUM_ into A;
    close have;

    /* 执行顺时针90度旋转：转置后逆序取列 */
    rotated = T(A)[, 6:1];

    /* 在日志窗口输出旋转结果，匹配期望格式 */
    print rotated[format=best1. noheader rowname=""];

    /* 可选：将旋转后的矩阵保存为SAS数据集 */
    create rotated_data from rotated[colname=("v1":"v6")];
    append from rotated;
    close rotated_data;
quit;

代码说明 #

• use have; read all var _NUM_ into A;：将数据集have中的所有数值变量按行读取到矩阵A，矩阵维度与原数据集一致（6×6）
• T(A)：对矩阵A执行转置操作
• [,6:1]：保留所有行，按从第6列到第1列的顺序选取元素，实现每行逆序，完成旋转
• print语句的参数：format=best1.确保每个元素占1位，noheader和rowname=""去掉表头和行名，让输出格式与期望完全匹配
• 若需要将结果保存为数据集，使用create和append语句即可生成名为rotated_data的新数据集

内容的提问来源于stack exchange，提问作者whymath