漏斗图的变化

说明

漏斗图显示了x轴上观察到的效应大小或结果，与y轴上观察到的效应大小或结果的精度进行了比较。基于Sterne和Egger(2001)，y轴的推荐选择是标准误差(按降序)，这也是metafor包中漏斗()函数的默认设置。在没有发表偏差和异质性的情况下，人们可能会看到这些点形成漏斗形状，大多数点落入伪置信度区域内，界限为^θ±1.96SEθ^±1.96SE，其中^θθ^是基于固定效应模型估计的效果或结果，SESE是来自y轴的标准误差值。对于y轴的其他精度度量，漏斗的预期形状可能会有很大不同。下图显示了y轴的各种选择，以及这如何影响漏斗图的形状(以及伪置信度区域的形式)。

图例

代码

library(metafor)
 
### fit fixed-effects model
res <- rma(yi, vi, data=dat.hackshaw1998, measure="OR", method="FE")
 
### set up 2x2 array for plotting
par(mfrow=c(2,2))
 
### draw funnel plots
funnel(res, main="Standard Error")
funnel(res, yaxis="vi", main="Sampling Variance")
funnel(res, yaxis="seinv", main="Inverse Standard Error")
funnel(res, yaxis="vinv", main="Inverse Sampling Variance")

参考文献

Sterne, J. A. C., & Egger, M. (2001). Funnel plots for detecting bias in meta-analysis: Guidelines on choice of axis. Journal of Clinical Epidemiology, 54(10), 1046–1055.