在SAS中计算均值所用观测数：无需转宽格式生成新列n

开发者阿泽

2026-7-29

解决方案：无需转宽格式即可计算均值涉及的层级数量

完全不需要将长格式数据转换为宽格式，用SAS的BY组DATA步或者PROC SQL分组聚合就能直接生成目标列n，以下是两种可行的实现方式：

方法1：DATA步BY组处理（适合保留原数据行顺序的场景）

首先确保数据按entity排序（BY组处理的前提），然后通过retain变量跨观测计数非缺失的value：

/* 先按entity排序，保证BY组分组正确 */
proc sort data=have;
    by entity;
run;

/* 计算每个entity的有效观测数n，并赋值给该entity的所有行 */
data want;
    set have;
    by entity;
    retain n 0; /* 保留n的数值，避免每次观测都重置 */
    if first.entity then n = 0; /* 每个entity的第一行重置计数 */
    if not missing(value) then n + 1; /* 非缺失value则计数+1 */
run;

代码说明：

first.entity是SAS BY组自动生成的变量，遍历到每个entity的第一行时触发重置操作
retain n 0让n的值在不同观测间传递，不会默认重置为缺失
最终每个entity的所有行都会带上该组的有效观测数（即计算均值用到的层级数量）

方法2：PROC SQL分组聚合（代码更简洁，无需提前排序）

通过子查询统计每个entity的非缺失value数量，再左连接回原数据集：

proc sql;
    create table want as
        select a.*, b.n
        from have as a
        left join (
            /* 子查询：按entity分组，统计非缺失value的数量 */
            select entity, count(value) as n
            from have
            group by entity
        ) as b
        on a.entity = b.entity;
quit;