SAS中Do Loop迭代异常：r1与r未达预期最大值

首先明确：你看到的r1、r最大值偏小，可能是循环实际跑了全范围，但只有符合统计条件的记录被输出，也可能是循环本身的上限写错了，按以下步骤排查：

• 先验证循环是否真的未迭代到预期范围
  在各循环块内添加put语句，输出当前循环变量值，比如：

  do n = 2 to &usern;
    put "当前n=" n;
    do n1 = 1 to (n - 1);
      put "  当前n1=" n1;
      do r1 = 0 to n1;
        put "    当前r1=" r1;
        do r = (r1+1) to &usern;
          put "      当前r=" r;
          /* 原有代码 */
        end;
      end;
    end;
  end;
  

  运行后查看日志，确认n、n1、r1、r是否真的没跑到你预期的范围（比如n到46时，n1到45，r1到n1，r到n）。

• 修正r的循环上限（明显错误）
  你预期r的范围是r1+1到n，但代码里写的是to &usern，这完全不符合统计逻辑（r是n个样本中的响应数，不可能超过n）。把r的循环改成：

  do r = (r1+1) to n;
  

  这个错误会导致r跑到比n大的无意义值，且这些值对应的P0_s、P1_s大概率不满足输出条件，所以不会出现在结果里。

• 检查内层条件对输出的过滤
  代码里有if r - l > 0 and n - r > 0 then do;，这个条件会过滤掉：

  • r=l的情况（即第二阶段没有响应）
  • r=n的情况（即所有样本都响应）
    如果你的统计逻辑需要保留这些情况，把条件改成if r - l >= 0 and n - r >= 0 then do;，或者根据实际试验设计的需求调整。

• 验证统计条件的过滤作用
  输出的记录必须满足P0_s <= &alpha and P1_s >= 1 - &beta，当r1或r过大时：

  • r1越大，PET_P0（早期终止概率）越高，但P0_s=1 - cdf('BINOMIAL', r, &p0, n)会越小，可能满足alpha，但P1_s=1 - cdf('BINOMIAL', r, &p1, n)会急剧下降，无法达到1-beta（0.9），所以不会被输出。
    可以在if判断前添加put语句，输出P0_s和P1_s的值，看哪些循环变量对应的统计量不满足条件，从而理解为什么某些值没出现在结果里。

• 检查变量类型和数值精度
  确认所有变量都是数值型，没有隐性类型转换导致的截断。SAS默认的数值精度足够处理这类计算，但可以通过format语句强制输出更多小数位，比如P0_s P1_s 10.6;，避免因精度问题导致条件判断出错。

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