如何用SAS批量填充60列缺失值？基于列运算均值规则

问题背景 #

需按特定规则填充B_01至B_60共60列的缺失值：

• B_01缺失值：(所有非缺失B_01/A的平均值) × 当前行A值
• B_n（n≥2）缺失值：(所有非缺失B_n/A的平均值) × 当前行A值 - 该行生成的new_B_01值

手动逐列编写Proc SQL效率极低，可通过SAS宏实现批量生成代码。

批量解决方案：SAS宏循环 #

利用SAS宏的循环功能，自动生成60列的计算语句，无需手动重复编写：

%macro fill_missing_cols;
Proc sql;
  create table new_table as
  Select *
  /* 处理第一列B_01 */
  , sum(B_01/A)/count(B_01)*A as new_B_01
  /* 循环生成B_02到B_60的计算逻辑 */
  %do col_num=2 %to 60;
  , sum(B_&col_num/A)/count(B_&col_num)*A - new_B_01 as new_B_&col_num
  %end;
  from table_one;
quit;
%mend fill_missing_cols;

/* 调用宏执行批量填充 */
%fill_missing_cols;

代码说明 #

1. 宏循环逻辑：%do col_num=2 %to 60自动遍历从2到60的列编号，生成每一列的计算语句
2. 统计简化：count(B_n)等价于手动编写的sum(case when B_n is missing then . else 1 end)，用于统计B_n的非缺失值数量，代码更简洁
3. 规则匹配：sum(B_n/A)/count(B_n)直接对应所有非缺失B_n/A的平均值，后续按需求完成乘A、减new_B_01的计算

示例验证 #

针对你提供的示例数据，宏会自动生成：

• new_B_01：(2/5 + 3/7 + 3/9)/3 × A，完全匹配需求公式
• new_B_02：(3/5 + 1.2/7 + 0.3/9)/3 × A - new_B_01，与需求逻辑一致

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