使用metafor创建带有残差的漏斗图

Question

在R中使用metafor :对于不包含任何主持人的模型，如何在漏斗图的x轴上绘制残差？

Answer 1

虽然这不是故意的，但您可以覆盖rma对象内部的逻辑，该对象指示您已经拟合了一个'intercept-only‘模型。举例说明：

library(metafor)
dat <- escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=dat.bcg)
res <- rma(yi, vi, data=dat)
funnel(res)

这给了你一个标准的漏斗图。现在强制int.only为FALSE。然后你会得到残差的漏斗图：

res$int.only <- FALSE
funnel(res)

另一种方法是手动提取残差和相应的方差，并将它们直接传递给funnel()函数：

ei  <- resid(res)
vei <- diag(vcov(res, type="resid"))
funnel(ei, vei, xlab="Residual Value")