Description
Sensitivity Analysis for Publication Bias in Meta-Analyses.
Description
Performs sensitivity analysis for publication bias in meta-analyses (per Mathur & VanderWeele, 2020 [<doi:10.31219/osf.io/s9dp6>]). These analyses enable statements such as: "For publication bias to shift the observed point estimate to the null, 'significant' results would need to be at least 30-fold more likely to be published than negative or 'nonsignificant' results." Comparable statements can be made regarding shifting to a chosen non-null value or shifting the confidence interval. Provides a worst-case meta-analytic point estimate under maximal publication bias obtained simply by conducting a standard meta-analysis of only the negative and "nonsignificant" studies.
README.md
PublicationBias
PublicationBias is an R package that provides bias correction and sensitivity analysis for publication bias in meta-analyses (per Mathur & VanderWeele, 2020).
Installation
You can install PublicationBias from CRAN with:
install.packages("PublicationBias")
You can install the development version of PublicationBias from GitHub with:
# install.packages("devtools")
devtools::install_github("mathurlabstanford/PublicationBias")
Example
Start by generating some example data from the metafor package.
library(PublicationBias)
dat <- metafor::escalc(measure = "RR", ai = tpos, bi = tneg, ci = cpos,
di = cneg, data = dat.bcg)
Calculate the meta-analytic effect size estimate, correcting for an assumed selection ratio of 5 (i.e., affirmative results are 5x more likely to be published than nonaffirmative ones).
pubbias_meta(yi = dat$yi, vi = dat$vi, selection_ratio = 5,
model_type = "fixed", favor_positive = FALSE)
Calculate how high the selection ratio would need to be to attenuate the effect size estimate to the null.
pubbias_svalue(yi = dat$yi, vi = dat$vi, q = 0,
model_type = "fixed", favor_positive = FALSE)