SAS中实现R语言循环分层建模及rbind式结果合并的方法咨询

你熟悉用R实现「分层创建子集→子集建模→提取结果→循环合并」的逻辑，现在需要转成SAS代码，核心有两个疑问：

1. 如何批量选择分组变量并循环创建子集（不用逐个列变量）
2. 如何循环合并模型结果（类似R的rbind）

一、批量选择分组变量+循环创建子集 #

SAS可以通过**系统字典表DICTIONARY.COLUMNS**批量获取变量名，再用宏循环完成遍历：

• 先从字典表提取指定范围的变量（比如原数据第4到第6列），存为宏变量列表
• 用%DO循环遍历每个分组变量，再通过PROC SQL获取该变量的所有类别，生成类别宏变量列表
• 嵌套循环每个类别，用WHERE语句筛选子集

二、循环合并结果（类似R的rbind） #

SAS里可以先定义好结果数据集的结构，再在每次循环中用PROC APPEND将当前模型结果追加进去，这是最贴合rbind逻辑的实现方式；也可以用DATA步的SET语句合并临时数据集。

完整SAS代码示例（对应你提供的R示例数据） #

/* 1. 创建对应R示例的数据集dd */
data dd;
    do i = 1 to 10;
        y = 0.1 * i;
        base = 3;
        treat = (i > 5);
        group1 = mod(i, 2); /* 对应R的rep(1:0,5) */
        group2 = 1 - mod(i, 2); /* 对应R的rep(0:1,5) */
        group3 = (i in (3,4,7,8,10)); /* 对应R的group3 */
        output;
    end;
    drop i;
run;

/* 2. 批量获取分组变量（第4到第6列，对应R的colnames(dd)[4:6]） */
proc sql noprint;
    select name into :group_vars separated by ' '
    from dictionary.columns
    where libname = 'WORK' and memname = 'DD' and varnum between 4 and 6;
quit;

/* 3. 初始化结果数据集，提前定义列结构 */
data results;
    length group_var $30 group_val $20 intercept base treat 8;
    stop; /* 只生成表结构，不产生空行 */
run;

/* 4. 嵌套宏循环实现分层建模+结果合并 */
%macro loop_groups;
    %let n_groups = %sysfunc(countw(&group_vars));
    %do i = 1 %to &n_groups;
        %let current_group = %scan(&group_vars, &i);
        
        /* 获取当前分组变量的所有类别 */
        proc sql noprint;
            select distinct strip(&current_group) into :group_vals separated by ' '
            from dd;
        quit;
        %let n_vals = %sysfunc(countw(&group_vals));
        
        /* 内层循环：遍历每个类别建模 */
        %do j = 1 %to &n_vals;
            %let current_val = %scan(&group_vals, &j);
            
            /* 对筛选后的子集运行线性回归，输出系数到临时数据集 */
            proc reg data=dd(where=(&current_group = &current_val)) noprint outest=model_result;
                model y = base treat;
            run;
            
            /* 整理结果格式 */
            data temp_result;
                set model_result;
                group_var = "&current_group";
                group_val = "&current_val";
                keep group_var group_val intercept base treat;
            run;
            
            /* 追加到总结果数据集（对应R的rbind） */
            proc append base=results data=temp_result force;
            run;
            
            /* 清理临时数据集，节省内存 */
            proc datasets lib=work nolist; delete temp_result model_result; run; quit;
        %end;
    %end;
%mend;
%loop_groups;

/* 查看最终合并后的结果 */
proc print data=results;
run;

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