# The metafor Package

A Meta-Analysis Package for R

### Site Tools

plots:contour_enhanced_funnel_plot_2

# Differences

This shows you the differences between two versions of the page.

 — plots:contour_enhanced_funnel_plot_2 [2019/06/03 18:43] (current) Line 1: Line 1: + ===== Contour-Enhanced Funnel Plot 2 ===== + + ==== Description ==== + + Here is another example of a contour-enhanced funnel plot (Peters et al., 2008) (see [[plots:​contour_enhanced_funnel_plot|here]] for a brief discussion of the idea behind such plots). This one is a very close re-creation of Figure 2(b) from Peters et al. (2008). Note that the values for ''​ylim''​ technically have to be greater than 0 when using the inverse standard error on the y-axis, but we can get essentially the same limits by using a very small value for the lower limit. Also, the dataset used by Peters et al. (2008) contains an extra study, which isn't included in the Cochrane review by Graves et al. (2010). + + ==== Plot ==== + + {{plots:​contour_enhanced_funnel_plot_2.png?​nolink}} + + ==== Code ==== + + + library(metafor) + + ### decrease margins so the full space is used + par(mar=c(5,​4,​1,​2)) + + ### calculate log odds ratios and corresponding sampling variances + dat <- escalc(measure="​OR",​ ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat.graves2010) + + ### create contour enhanced funnel plot + funnel(dat$yi,​ dat$vi, yaxis="​seinv",​ + ​xlim=c(-3,​2),​ ylim=c(.00001,​8),​ xaxs="​i",​ yaxs="​i",​ las=1, + ​level=c(.10,​ .05, .01), shade=c("​white",​ "​gray55",​ "​gray75"​),​ + ​legend=TRUE,​ back="​gray90",​ ylab="​Precision (1/​se)"​) + ​ + + ==== References ==== + + Graves, P. M., Deeks, J. J., Demicheli, V., & Jefferson, T. (2010). Vaccines for preventing cholera: Killed whole cell or other subunit vaccines (injected). //Cochrane Database of Systematic Reviews, 8//, CD000974. + + Peters, J. L., Sutton, A. J., Jones, D. R., Abrams, K. R., & Rushton, L. (2008). Contour-enhanced meta-analysis funnel plots help distinguish publication bias from other causes of asymmetry. //Journal of Clinical Epidemiology,​ 61//(10), 991--996.