R语言双样本t检验：比较方向如何判定？

这个问题确实容易搞懵，毕竟t.test的帮助文档没把这个点写得特别直白——我来给你拆解清楚：

核心逻辑：因子水平顺序决定对比方向 #

当你用公式格式y ~ group调用t.test()时，它的对比方向完全由分组变量（这里是smoking）的因子水平顺序决定：

• R会把分组变量的第一个水平当作「参考组」，第二个水平当作「对比组」
• 检验计算的差值是：对比组均值 - 参考组均值
• alt="greater"对应的备择假设是：这个差值 > 0，也就是「对比组均值 > 参考组均值」

针对你的具体场景 #

你的smoking变量取值是yes/no，分两种情况：

1. 如果smoking是字符向量（没转成因子）：
   R会自动把它转成因子，并且按字母顺序排序水平——"no"在"yes"前面（n的字母顺序比y早），所以因子水平是c("no", "yes")。
   此时t.test(lungCap~smoking, alt="greater")检验的是：mean(yes组) > mean(no组)，这和你想要验证的「非吸烟者（no）肺活量更大」正好相反！

2. 如果你手动指定了因子水平顺序：
   比如你执行smoking <- factor(smoking, levels=c("yes", "no"))，把yes设为第一个水平（参考组），no设为第二个水平（对比组）。
   此时检验的差值是mean(no组) - mean(yes组)，alt="greater"就对应你要的备择假设：mean(no组) > mean(yes组)。

怎么快速验证？ #

你可以用一小段测试代码确认：

# 构造模拟数据：no组均值8，yes组均值6（符合你预期的非吸烟者肺活量大）
set.seed(123)
lungCap <- c(rnorm(50, 8), rnorm(50, 6))
smoking <- rep(c("no", "yes"), each=50)

# 运行你原来的命令
t.test(lungCap~smoking, alt="greater")

看输出里的difference in means：会显示yes - no，结果是负数，而且p值接近1——这说明检验在尝试证明yes组>no组，但实际数据是反过来的，所以不显著。

为什么帮助文档没写？ #

因为这个逻辑是R中因子对比的通用规则，不是t.test()特有的——所有基于公式的统计函数（比如lm()）都遵循「后水平减前水平」的对比逻辑，所以帮助文档默认你已经了解这个基础规则。

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