使用带n、p参数的SAS宏循环演示大数定律无输出问题求助

你的代码无输出的核心原因是仅定义了SAS宏，但从未调用执行它，同时代码还存在几处逻辑问题，以下是具体修正方案：

1. 核心问题：缺少宏调用 #

宏定义完成后，必须通过%LLN();（可指定参数，如%LLN(n=100, p=0.6);）触发执行，否则SAS不会运行宏内的data步代码。

2. 代码逻辑修正 #

原代码中的宏变量sum和mean完全多余，且data步内的sum未初始化会导致首次计算出现缺失值。修正后的完整代码如下：

%macro LLN(n=10, p=0.6);
   data LLN; 
      /* 初始化data步内的sum变量，避免缺失值 */
      sum = 0;
      %do i = 1 %to &n;
         x = rand("Bernoulli", &p); 
         sum = sum + x;
         mean = sum / &i; /* 用宏变量&i代替循环变量，避免data步变量混淆 */
         output; 
      %end;
   run;
   /* 打印结果到日志，直观查看大数定律效果 */
   proc print data=LLN;
      var x sum mean;
   run;
%mend LLN;

/* 调用宏，可自定义n和p参数 */
%LLN(n=100, p=0.6);

3. 修正说明 #

• 移除无用的宏变量sum和mean，在data步内直接初始化sum为0，解决缺失值问题。
• 循环内的mean计算改用宏变量&i，因为%do循环的i是宏变量，data步内默认不存在同名变量，不替换会报错。
• 添加proc print步骤，将结果直接输出到SAS日志，便于观察样本均值随观测数增加逐渐趋近于p的过程。

效果验证 #

调用宏后，SAS日志会输出对应n条观测，mean列数值会逐步逼近指定的p值（如示例中的0.6），清晰演示大数定律的核心结论。

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