如何用SAS按分数阈值增量统计符合条件的成员观测数？

针对需求，这里提供两种SAS实现方式，解决你用数组和DO循环遇到的问题：

方法1：DATA步+临时数组+DO循环 #

这种方式用临时数组存储每个阈值的计数，遍历所有观测完成累加后统一输出结果：

/* 假设你的原始数据集名为member_data，包含ID和score变量 */
data threshold_counts;
  set member_data end=eof;
  /* 临时数组：对应21个阈值（0.20到0.40，步长0.01）的计数器 */
  array counts[21] _temporary_;
  
  /* 遍历每个阈值，对当前观测的score判断后累加计数 */
  do i = 1 to 21;
    current_threshold = 0.20 + (i-1)*0.01;
    if score > current_threshold then counts[i] + 1;
  end;
  
  /* 处理完所有观测后，输出每个阈值及对应的达标人数 */
  if eof then do i = 1 to 21;
    threshold = 0.20 + (i-1)*0.01;
    count = counts[i];
    output;
  end;
  
  /* 只保留需要的变量 */
  keep threshold count;
run;

/* 查看结果 */
proc print data=threshold_counts noobs;
  format threshold 5.2;
run;

关键说明： #

• 临时数组_temporary_仅用于内存中存储计数，不会写入输出数据集，避免冗余。
• 不需要用put nobs=;，这里通过counts[i] + 1直接累加符合条件的观测数，最后统一输出即可。
• end=eof标记最后一条观测，确保所有数据处理完成后再输出统计结果。

方法2：PROC SQL简洁实现 #

如果更习惯SQL语法，可直接生成阈值序列并分组统计：

proc sql;
  create table threshold_counts as
    select 
      0.20 + (i-1)*0.01 as threshold format=5.2,
      sum(case when score > calculated threshold then 1 else 0 end) as count
    from member_data,
         /* 生成1到21的序列，对应21个阈值 */
         (select i from generate(1 to 21)) as thresholds
    group by threshold;
quit;

这种方式无需手动处理循环和数组，SAS会自动完成分组统计。

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