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生成的置信区间。

有没有人能帮我一下？

谢谢

S

Answer 1

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)