如何让SAS输出count=0时的所有age_groups？PROC FREQ字符型变量问题

问题背景 #

需要让SAS输出所有age_groups分类（包括count=0的情况）。由于部分age_groups对应的行缺失，手动在Excel文件间复制粘贴2000行数据极易导致统计错误。尝试过PROC FREQ的sparse选项但未生效，推测原因是age_group为字符型变量，缺失的分类完全未出现在频率统计结果中。

当前输出 #

仅包含有数据记录的age_group分类，每个mmwr_case（MMWR周）对应的年龄分组不完整，缺少count=0的行。

期望输出 #

每个mmwr_case下都完整列出所有age_group分类，即使该分组的count为0，保证每周的年龄分组全覆盖。

原SAS代码 #

PROC FREQ data=partial ORDER=INTERNAL ;
    tables mmwr_case*age_group / out=partial_2021 sparse list nocol nocum norow nopercent;
    by mmwr_case;
    where case_year = 2021; 
    title 'Cases, 2021, by MMWR Week';
RUN; 

解决方案 #

sparse选项仅能补全当前BY组中存在部分记录的分类，无法补全完全没有记录的字符型分类。需先构建包含所有age_group和mmwr_case的完整框架，再与原数据关联补全count值。

方法1：使用PROC SQL构建完整框架并统计

1. 创建包含所有age_group的完整列表

data age_groups_full;
    length age_group $50; /* 根据实际变量长度调整 */
    input age_group $;
    datalines;
0-4
5-17
18-29
30-49
50-64
65+
/* 替换为你实际的所有age_group分类 */
;
run;

2. 提取2021年所有的mmwr_case

proc sql;
    create table mmwr_weeks_2021 as
    select distinct mmwr_case
    from partial
    where case_year = 2021;
quit;

3. 生成笛卡尔积，得到每个周对应所有年龄组的空框架

proc sql;
    create table full_frame as
    select a.mmwr_case, b.age_group
    from mmwr_weeks_2021 a, age_groups_full b;
quit;

4. 左连接原数据，补全count值

proc sql;
    create table partial_2021 as
    select 
        f.mmwr_case, 
        f.age_group,
        coalesce(count(c.mmwr_case), 0) as count
    from full_frame f
    left join partial c
        on f.mmwr_case = c.mmwr_case 
        and f.age_group = c.age_group
        and c.case_year = 2021
    group by f.mmwr_case, f.age_group
    order by f.mmwr_case, f.age_group;
quit;

方法2：合并框架后用PROC FREQ统计

如果偏好使用PROC FREQ，可先合并完整框架与原数据，再通过权重变量统计：

/* 合并完整框架与原数据，标记是否有记录 */
data combined;
    merge full_frame partial(in=c);
    by mmwr_case age_group;
    count = ifn(c, 1, 0);
run;

/* 用PROC FREQ统计所有分组 */
proc freq data=combined order=internal;
    tables mmwr_case*age_group / out=partial_2021 list nocol nocum norow nopercent;
    by mmwr_case;
    weight count;
    title 'Cases, 2021, by MMWR Week';
run;

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