如何使用不同的数据子集迭代运行 SAS 过程？

Question

我想使用现有 SAS 数据集的不同子集重复运行 PROC REG。这是一个简单的示例数据集：

DATA data_main;
  input trt depth year response;
  cards;
1 1 2014 1.1
1 2 2014 1.2
2 1 2014 1.3
2 2 2014 1.4
1 1 2013 2.2
1 2 2013 2.4
2 1 2013 2.6
2 2 2013 2.8
;
run;

对于 trt 和 depth 的每个组合，我想运行这个过程，其中current_data 是trt 和depth 的当前组合：

PROC REG data = current_data;
    model response = year;
run;

我想在一个数据集或文本文件中捕获所有迭代的回归系数和 p 值。

在我的实际数据集中，输入和 trt 的级别数要多得多，所以我试图避免手动编码每个组合。有人可以向我解释如何做到这一点吗？

Answer 1

考虑运行一个宏迭代 trt 和 depth 的组合。下面的嵌套循环迭代地重新创建current_data 数据集，并在输出相应结果表的回归过程中使用它。根据需要为所有组合调整循环限制中的值范围：

%macro loopregression;
%do j = 1 %to 2;                 * TRT VALUES;
    %do i = 2013 %to 2014;       * DEPTH VALUES; 
        DATA current_data; 
            SET data_main;  
            if trt = &j;    
            if depth = &i;
        run;

        PROC REG data = current_data noprint outest=results&i&j;;
            model response = year;
        run;
    %end;
%end;
%mend loopregression;

%loopregression;