我有来自三个不同来源(医院)的配对t检验的数据摘要,这意味着我没有全部样本数据。我要进行配对t检验的是平均值(每个配对样本数据的差异)。
数据看起来像下面的表格(每行都是一个病人),我有mean(diff)和sd(diff)。
有什么功能,我可以很容易地进行配对t检验,并得到95%的CI?由于病人的隐私问题,没有办法获得整个数据集。
3 months 6 months diff
1 3 -2
2 1 1
5 9 -4发布于 2018-12-04 02:08:24
配对t检验与t检验完全等价,t检验用零的零假设检验差异的平均值。
diffs <- c(-2,1,-4)
m <- mean(diffs)
s <- sd(diffs)
n <- length(diffs)现在假设您不了解diffs,只知道m、s和n (这三者都是必需的)
pval <- 2*pt(abs(m/s),lower.tail=FALSE,df=n)
## 0.555
qq <- qt(0.975,df=n)
ci <- m + c(-1,1)*qq*s
## [1] -9.675648 6.342314https://stackoverflow.com/questions/53604064
复制相似问题