我一直在使用 Metafor 进行元分析。 我已经能够从 14 篇文章中提取比值比、置信区间、样本量。
这是我一直在使用的代码:
data <- read.csv(file="OR.csv", check.names=F,sep=",")
dat=data
random = rma(OR,SE_X,data=dat)
forest (random)
funnel (random)
我的问题:rma函数生成的置信区间与作者提供的有很大不同。 这显然是一个问题,我想使用文章提供的置信区间,而不是 RMA 生成的置信区间。
有人可以帮我吗?
谢谢你
S
rma()
函数的第二个参数用于采样方差。但是,您将其设置为 SE_X
,这似乎是 标准错误。所以,你应该使用rma(OR, sei=SE_X, data=dat)
。
但是,如果
OR
真的代表优势比,那么这无论如何都是不对的。人们可以为对数优势比构建置信区间,然后对区间边界进行反向变换(即求幂)。假设您的标准误差也适用于对数优势比,那么您会想要这样做:
random <- rma(log(OR), sei=SE_X, data=dat)
forest(random, atransf=exp)