DescTools包中的TTestA函数。
嗨,我正在尝试做一个独立的样本t检验,但没有完整的数据,只有汇总统计数据。有没有一个R包,我可以使用每个样本的平均值、标准偏差和计数?
另外:你可以自己实现!来自维基百科关于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检验(。