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

在R中使用metafor：

---

对于不包含任何调节器的模型，如何在漏斗图的x轴上绘制残差

尽管这不是有意的，但您可以覆盖指示您已安装“仅拦截”模型的逻辑内部

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)

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

漏斗（）

函数：

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

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