SAS中嵌套Do循环生成多层模拟数据运行报错问题求助

报错核心原因 #

1. 混淆了宏编译和数据步运行的时序：将误差生成函数用%let赋值给宏变量&e，宏变量是纯文本替换，仅在宏编译阶段解析，无法在数据步运行时动态生成随机数，直接触发语法或逻辑错误。
2. 存在多余循环层：原代码中do I=1 to (10000/&N)的循环会生成远超需求的冗余观测，与1000×4×5的目标样本结构不符。
3. 宏循环使用冗余：遍历样本量分组、误差分布分组不需要用宏循环实现，用普通数据步循环即可，避免不必要的宏解析问题。

修正后完整代码 #

%macro MLM;
data sim_data;
    call streaminit(525600);
    * 外层：1000次模拟循环;
    do rep=1 to 1000;
        * 第二层：4组样本量分组;
        do a=1 to 4;
            select(a);
                when(1) N=100;
                when(2) N=250;
                when(3) N=500;
                when(4) N=1000;
            end;
            * 第三层：生成对应数量的组水平样本;
            do J=1 to N;
                u0 = RAND('NORMAL',0,1);
                w1 = RAND('BERNOULLI',.5);
                * 第四层：5种误差分布遍历;
                do b=1 to 5;
                    select(b);
                        when(1) e = RAND('NORMAL',0,1);
                        when(2) e = RAND('UNIFORM');
                        when(3) e = RAND('CHISQUARE', N-1);
                        when(4) e = RAND('LOGNORMAL');
                        when(5) e = RAND('BETA',1,1);
                    end;
                    * 生成个体水平变量与结果;
                    x1 = RAND('BERNOULLI',.5);
                    b0 = 2.5 + 1*w1 + u0;
                    b1 = 0.15;
                    y = b0 + b1*x1 + e;
                    output;
                end;
            end;
        end;
    end;
run;
%mend MLM;
%MLM;

调整说明 #

• 所有宏循环、宏判断均替换为数据步原生循环与select条件判断，规避宏解析与数据步运行的时序冲突，更适配新手使用场景
• 删除冗余的I层循环，最终生成的数据集完全匹配1000次模拟×4组样本量×5种误差分布的结构要求
• 误差e改为数据步普通变量，每次循环动态生成对应分布的随机数，保证模拟逻辑正确
• 保留原设定的随机种子与变量生成规则，确保模拟结果可复现

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