SAS proc report嵌套表头按分组值合并行消除空单元格问题咨询

现有代码出现行错位、大量空单元格的核心原因有两个：

• proc report中cylinders被设置为display类型，该类型会为数据集内每一条观测单独生成一行。当前频数数据集里每个气缸数对应3条不同产地的观测，最终每个气缸数会占3行，非对应产地的列自动留空，形成错位。
• 用固定值var=1做SQL内连接的写法冗余，一旦两边数据集行数变化很容易生成笛卡尔积脏数据，进一步放大错位问题。

不需要调整proc freq的核心逻辑，只需要优化合并逻辑、修正proc report的变量类型设置即可：

1. 去掉冗余的关联键字段，直接在数据步用_n_=1的判断只读一次卡方结果，将总样本量、卡方P值retain到所有频数行，避免笛卡尔积风险
2. 将cylinders设置为group类型，proc report会自动合并相同气缸数的行，每个气缸数仅占1行，各产地对应的频数、行百分比会自动对齐到嵌套表头的对应列下
3. 给统计值加上适配格式，提升报表可读性

修正后的可运行代码如下：

proc freq data=sashelp.cars;    
    tables origin*cylinders / chisq outpct out=freqpct(keep=cylinders origin count pct_row); 
    output out=chisq(keep=N P_PCHI) chisq;
run;

data report_data;
    set freqpct;
    if _n_ = 1 then set chisq;
    retain N P_PCHI;
run;

proc report data=report_data nowd headline headskip;
    column cylinders origin,(count pct_row) N P_PCHI;
    /* 关键设置：cylinders为group类型自动合并同值行 */
    define cylinders / group '气缸数' order=internal;
    define origin / across '产地' order=internal;
    define count / display '频数' format=comma6.;
    define pct_row / display '行百分比' format=percent8.1;
    define N / group '总样本量' format=comma6.;
    define P_PCHI / group '卡方检验P值' format=pvalue6.4;
run;

运行后生成的报表无空单元格，每个气缸数仅显示一次，各产地的频数、行百分比完全对齐嵌套表头，最右侧的总样本量、卡方P值也不会出现错位。

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