考虑下面两个独立的虚构正常数据集。
set.seed(910)
x1 = rnorm(100, 50, 7)
x2 = rnorm(90, 60, 8)
mean(x1)
[1] 49.86995
mean(x2)
[1] 60.04863
样本均值显示X¯1<X¯2建议μ1可能“显着”小于μ2.R 中的 Welch 两样本 t 检验表明,原假设H0:μ1=μ2被拒绝支持Ha:μ1<μ2.P值接近0.
更具体地说,单侧 CI:显示 95% 的置信度
μ1−μ2≤−8.412.
t.test(x1, x2, alt="less")
Welch Two Sample t-test
data: x1 and x2
t = -9.5241, df = 187.95, p-value < 2.2e-16
alternative hypothesis:
true difference in means is less than 0
95 percent confidence interval:
-Inf -8.412074
sample estimates:
mean of x mean of y
49.86995 60.04863
这里没有什么是“证明”μ1<μ2.当然,因为我模拟了这些数据,我们知道 μ1−μ2=−10.然而,在真正的统计应用程序中,我们不会有如此准确的信息。我们能做的最好的就是报告样本意味着显示X¯1=49.87<X¯2=60.05
并说我们“有理由确定”μ1−μ2≤−8.412.