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

Question

我一直在使用Metafor运行元分析。我已经能够从14篇文章中提取优势比、置信区间和样本量。

这是我一直在使用的代码：

data <- read.csv(file="OR.csv", check.names=F,sep=",")dat=datarandom = 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)