data data_set0;
input d_1 d_2 d_3;
datalines;
1 2 3
4 5 6
7 8 9
;
data data_set1;
input d_1 d_2 d_3;
datalines;
1 2 3
4 5 6
7 8 9
;
data data_set2;
input d_1 d_2 d_3;
datalines;
1 2 3
4 5 6
7 8 9
;
data data_set3;
input d_1 d_2 d_3;
datalines;
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_set0;
READ all var _ALL_ into data_set0_matrix[colname=varNames];
CLOSE data_set_0;
USE data_set1;
READ all var _ALL_ into data_set1_matrix[colname=varNames];
CLOSE data_set_1;
USE data_set2;
READ all var _ALL_ into data_set2_matrix[colname=varNames];
CLOSE data_set_2;
USE data_set3;
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; Just minor corrections will make your code work as above 🙂 EDIT: I see you made it work already. Thumbs up 🙂
... View more