使用二分水平与PRT的SAS Proc Means：双水平CLASS的P值疑问

当你用Proc Means做双水平CLASS变量的均值比较时，双侧检验的P值理论上不会随组顺序变化，仅检验统计量符号反转。你遇到P值差异，大概率是以下几个问题导致：

1. CLASS变量并非严格双水平 #

如果CLUSTER变量存在缺失值，或者有隐藏的第三水平，Proc Means会自动把缺失值/额外水平作为单独一组，此时均值比较会从t检验变成单因素ANOVA，P值自然会变化。

• 排查代码：

  Proc freq data=prelim_notsnp2;
    tables CLUSTER / missing;  /* 显示缺失值对应的组 */
  Run;
  

2. 数据或变量范围有误 #

• 检查Var q1–q24;中的连字符：如果是全角的–而非半角-，SAS会把它当成变量名的一部分，而非变量范围，导致实际分析的变量不对，结果异常。直接把全角符号改成半角-即可。
• 确认两次运行的数据集是否一致：如果prelim_notsnp2数据集在两次运行间有更新（比如新增/删除观测、修改变量值），P值必然会变。

3. 缺失值处理影响样本量 #

Proc Means默认对每个变量单独计算有效样本量（即只排除该变量缺失的观测），如果某个变量在两组中的有效样本量发生变化，t检验的自由度和标准误会改变，最终影响P值。

• 排查代码：

  Proc means data=prelim_notsnp2 n nmiss;
    Var q1-q24;
    Class CLUSTER;
  Run;
  

  对比两次运行中各变量的有效样本量（N）和缺失数（NMISS）是否一致。

验证方法 #

用Proc Ttest做双水平检验，对比Proc Means的结果：

Proc ttest data=prelim_notsnp2;
  Class CLUSTER;
  Var q1-q24;
Run;

Proc Ttest会明确输出双侧P值，正常情况下和Proc Means的prt结果完全一致。如果不一致，说明你的CLUSTER变量实际不是双水平（比如包含缺失组）。

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