事实上,你可以使用宏函数来轻松地做到这一点:
1)设置一个空数据集:
data result;
run;
2)创建一个宏函数来处理您的数据集。VAR 宏变量将是列表 B1、B2、B3、B4 中的数据集数量。
%macro processDatasets(var);
proc iml ;
use B&var;
read all var _ALL_ into B&var;
close B&var;
g= B&var.[1,1];
ExG = B&var.[3,1];
Res = B&var.[5,1];
E =2;
R =3;
您在这里使用调用 symputx("g", B&var.[1,1], 'G'); 存储值。然后将全局宏变量添加到数据集 to_store 中。
call symputx("g", B&var.[1,1], 'G');.
....
test = g/(g+(ExG/E)+(Res/(R*E)));
print g ExG res test [format = 10.6];
create try var {g ExG res test};
append;
close try;
run; quit;
这是如何将宏变量存储到当前数据集中的方法。
data to_store;
g=call symget("g");
ExG=call symget("ExG");
...
run;
3) 这是将多个数据集连接成一个结果的技巧。
data result;
set result to_store;
run;
%mend processDataset;
4)您可以对函数进行多次调用:
%processDataset(1);
%processDataset(2);
%processDataset(3);
...
5) 或者你可以遍历一个数字列表:
%macro processVarList;
%let var_list=1|2|3|4;
%let k=1;
%do %while (%qscan(&var_list, &k,|) ne );
%let VAR = %scan(&var_list, &k,|);
%processDataset(&var);
%let k = %eval(&k + 1);
%end;
%mend processVarList;
%processVarList;
6)您可以根据要存储在结果数据集中的内容来调整此代码。我没有测试代码,但想法就在那里。我想您必须将所有变量 G、ExG、Res、E、R 存储到数据集中才能合并到最终结果数据集中。
最后,您将所有价值都放入结果数据集中。