从SAS / IML中的四个单独矩阵中创建一个平均矩阵

Question

我首次在SAS Enterprise指南中使用IML/SAS，并希望执行以下操作：

1. 将一些数据集读取到IML矩阵
2. 平均矩阵
3. 将结果IML矩阵返回SAS数据集

我的输入数据集看起来像以下内容（这是虚拟数据 - 实际集合更大）。输入数据集的格式也是我从输出数据集中获得的格式。

data_set0:     d_1    d_2    d_3
               1      2      3
               4      5      6 
               7      8      9 

我进行以下操作：

proc iml;
    /* set the names of the migration matrix columns */
    varNames = {"d_1","d_2","d_3"};
    /* 1. transform input data set into matrix 
    USE data_set_0;
    READ all var _ALL_ into data_set0_matrix[colname=varNames];
    CLOSE data_set_0;
    USE data_set_1;
    READ all var _ALL_ into data_set1_matrix[colname=varNames];
    CLOSE data_set_1;
    USE data_set_2;
    READ all var _ALL_ into data_set2_matrix[colname=varNames];
    CLOSE data_set_2;
    USE data_set_3;
    READ all var _ALL_ into data_set3_matrix[colname=varNames];
    CLOSE data_set_3;
    /* 2. find the average matrix */
    matrix_sum = (data_set0_matrix + data_set1_matrix + 
                  data_set2_matrix + data_set3_matrix)/4;
    /* 3. turn the resulting IML matrix back into a SAS data set */ 
    create output_data from matrix_sum[colname=varNames];
    append from matrix_sum; 
    close output_data; 
 quit; 

我一直在尝试大量的东西，但似乎对我无济于事。我目前会读取的错误：

ERROR: Matrix matrix_sum has not been set to a value

我在做什么错？感谢您的帮助。

Answer 1

上述代码有效。在此代码的完整版本中（这是可读性的简化），我误为我的一个变量之一。

如果其他人想使用SAS/IML找到平均矩阵，我将把问题放在上面。