使用检验统计量(均值、标准偏差、计数)对R进行独立样本t检验



嗨,我正在尝试做一个独立的样本t检验,但没有完整的数据,只有汇总统计数据。有没有一个R包,我可以使用每个样本的平均值、标准偏差和计数?

DescTools包中的TTestA函数。

另外:你可以自己实现!来自维基百科关于t检验的文章,关于独立样本t检验、不等样本量的章节

示例数据:

means <- c(2.5, 3.7)
stddev <- c(1.5, 1.8)
n <- c(25,35)

计算汇总统计:

df <- sum(n)-2
stddev_pooled <- sqrt(sum((n-1)*stddev^2)/df)
tstat <- (means[1]-means[2])/(stddev_pooled*sqrt(sum(1/n)))

现在计算p值:

2*pt(abs(tstat),df=df,lower.tail=FALSE)

如果你想要不相等方差的情况,你可以在维基百科页面的下一节中实现这个公式(韦尔奇t检验(。

最新更新