Metafor:使用提供的置信区间

问题描述 投票:0回答:1

我一直在使用 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

r analysis medical metafor
1个回答
0
投票

rma()
函数的第二个参数用于采样方差。但是,您将其设置为
SE_X
,这似乎是 标准错误。所以,你应该使用
rma(OR, sei=SE_X, data=dat)

但是,如果

OR
真的代表优势比,那么这无论如何都是不对的。人们可以为对数优势比构建置信区间,然后对区间边界进行反向变换(即求幂)。假设您的标准误差也适用于对数优势比,那么您会想要这样做:

random <- rma(log(OR), sei=SE_X, data=dat)
forest(random, atransf=exp)
© www.soinside.com 2019 - 2024. All rights reserved.