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

Question

在 R 中使用 metafor：对于不包含任何调节器的模型，如何在漏斗图的 x 轴上绘制残差？

Answer 1

尽管这不是有意的，但您可以覆盖rma表明您已安装“仅拦截”模型的逻辑内部对象。为了显示：

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")