r - Metafor：使用提供的置信区间

Question

我一直在使用 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 生成的置信区间。

任何人都可以帮我解决这个问题吗？

谢谢

小号

score 0 · Accepted Answer

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

但是，如果OR真的代表优势比，那么这无论如何都是不对的。可以为对数优势比构建 CI，然后对区间边界进行反向变换（即取幂）。假设您的标准误差也适用于对数优势比，那么您需要执行以下操作：

random <- rma(log(OR), sei=SE_X, data=dat)
forest(random, atransf=exp)

1 回答 1