metafor-package | R Documentation |

The metafor package provides a comprehensive collection of functions for conducting meta-analyses in **R**. The package can be used to calculate various effect size or outcome measures and then allows the user to fit equal-, fixed-, and random-effects models to these data. By including study-level variables (‘moderators’) as predictors in these models, (mixed-effects) meta-regression models can also be fitted. For meta-analyses of \mjeqn2 \times 22x2 tables, proportions, incidence rates, and incidence rate ratios, the package also provides functions that implement specialized methods, including the Mantel-Haenszel method, Peto's method, and a variety of suitable generalized linear mixed-effects models (i.e., mixed-effects logistic and Poisson regression models). For non-independent effects/outcomes (e.g., due to correlated sampling errors, correlated true effects or outcomes, or other forms of clustering), the package also provides a function for fitting multilevel and multivariate models to meta-analytic data.

Various methods are available to assess model fit, to identify outliers and/or influential studies, and for conducting sensitivity analyses (e.g., standardized residuals, Cook's distances, leave-one-out analyses). Advanced techniques for hypothesis testing and obtaining confidence intervals (e.g., for the average effect or outcome or for the model coefficients in a meta-regression model) have also been implemented (e.g., the Knapp and Hartung method, permutation tests, cluster-robust inference methods / robust variance estimation).

The package also provides functions for creating forest, funnel, radial (Galbraith), normal quantile-quantile, L'Abbé, Baujat, bubble, and GOSH plots. The presence of publication bias (or more precisely, funnel plot asymmetry or ‘small-study effects’) and its potential impact on the results can be examined via the rank correlation and Egger's regression test, the trim and fill method, the test of excess significance, and by applying a variety of selection models.

[`escalc`

] Before a meta-analysis can be conducted, the relevant results from each study must be quantified in such a way that the resulting values can be further aggregated and compared. The `escalc`

function can be used to compute a wide variety of effect size or outcome measures (and the corresponding sampling variances) that are often used in meta-analyses (e.g., risk ratios, odds ratios, risk differences, mean differences, standardized mean differences, response ratios / ratios of means, raw or r-to-z transformed correlation coefficients). Measures for quantifying some characteristic of individual groups (e.g., in terms of means, proportions, or incidence rates and transformations thereof), measures of change (e.g., raw and standardized mean changes), and measures of variability (e.g., variability ratios and coefficient of variation ratios) are also available.

[`rma.uni`

] The various meta-analytic models that are typically used in practice are special cases of the general linear (mixed-effects) model. The `rma.uni`

function (with alias `rma`

) provides a general framework for fitting such models. The function can be used in combination with any of the effect size or outcome measures computed with the `escalc`

function or, more generally, any set of estimates (with corresponding sampling variances or standard errors) one would like to analyze. The notation and models underlying the `rma.uni`

function are explained below.

For a set of \mjseqni = 1, ..., k independent studies, let \mjseqny_i denote the observed value of the effect size or outcome measure in the \mjeqni\textrmthith study. Let \mjseqn\theta_i denote the corresponding (unknown) true effect/outcome, such that \mjdeqny_i | \theta_i \sim N(\theta_i, v_i).y_i | \theta_i ~ N(\theta_i, v_i). In other words, the observed effect sizes or outcomes are assumed to be unbiased and normally distributed estimates of the corresponding true effects/outcomes with sampling variances equal to \mjseqnv_i (where \mjseqnv_i is just the square of the standard errors of the estimates). The \mjseqnv_i values are assumed to be known. Depending on the outcome measure used, a bias correction, normalizing, and/or variance stabilizing transformation may be necessary to ensure that these assumptions are (at least approximately) true (e.g., the log transformation for odds/risk ratios, the bias correction for standardized mean differences, Fisher's r-to-z transformation for correlations; see `escalc`

for more details).

According to the **random-effects model**, we further assume that \mjeqn\theta_i \sim N(\mu, \tau^2)\theta_i ~ N(\mu, \tau^2), that is, the true effects/outcomes are normally distributed with \mjseqn\mu denoting the average true effect/outcome and \mjseqn\tau^2 the variance in the true effects/outcomes (\mjseqn\tau^2 is therefore often referred to as the amount of ‘heterogeneity’ in the true effects/outcomes). The random-effects model can also be written as \mjsdeqny_i = \mu + u_i + \epsilon_i, where \mjeqnu_i \sim N(0, \tau^2)u_i ~ N(0, \tau^2) and \mjeqn\epsilon_i \sim N(0, v_i)\epsilon_i ~ N(0, v_i). The fitted model provides estimates of \mjseqn\mu and \mjseqn\tau^2, that is, \mjdeqn\hat\mu = \frac\sum_i=1^k w_i y_i\sum_i=1^k w_i,\mu-hat = \sum w_i y_i / \sum w_i, where \mjeqnw_i = 1/(\hat\tau^2 + v_i)w_i = 1/(\tau-hat^2 + v_i) and \mjeqn\hat\tau^2\tau-hat^2 denotes an estimate of \mjseqn\tau^2 obtained with one of the many estimators that have described in the literature for this purpose (this is the standard ‘inverse-variance’ method for random-effects models).

A special case of the model above is the **equal-effects model** (also sometimes called the common-effects model) which arises when \mjseqn\tau^2 = 0. In this case, the true effects/outcomes are homogeneous (i.e., \mjeqn\theta_1 = \theta_2 = ... = \theta_k \equiv \theta\theta_1 = \theta_2 = ... = \theta_k = \theta) and hence we can write the model as \mjsdeqny_i = \theta + \epsilon_i, where \mjseqn\theta denotes *the* true effect/outcome in the studies, which is estimated with \mjdeqn\hat\theta = \frac\sum_i=1^k w_i y_i\sum_i=1^k w_i,\theta-hat = \sum w_i y_i / \sum w_i, where \mjeqnw_i = 1/v_iw_i = 1/v_i (again, this is the standard ‘inverse-variance’ method as described in the meta-analytic literature). Note that the commonly-used term ‘fixed-effects model’ is not used here - for an explanation, see here.

Study-level variables (often referred to as ‘moderators’) can also be included as predictors in meta-analytic models, leading to so-called ‘meta-regression’ analyses (to examine whether the effects/outcomes tend to be larger/smaller under certain conditions or circumstances). When including moderator variables in a random-effects model, we obtain a **mixed-effects meta-regression model**. This model can be written as \mjdeqny_i = \beta_0 + \beta_1 x_i1 + \beta_2 x_i2 + ... + \beta_p' x_ip' + u_i + \epsilon_i,y_i = \beta_0 + \beta_1 x_i1 + \beta_2 x_i2 + ... + \beta_p' x_ip' + u_i + \epsilon_i, where \mjeqnu_i \sim N(0, \tau^2)u_i ~ N(0, \tau^2) and \mjeqn\epsilon_i \sim N(0, v_i)\epsilon_i ~ N(0, v_i) as before and \mjeqnx_ijx_ij denotes the value of the \mjeqnj\textrmthjth moderator variable for the \mjeqni\textrmthith study (letting \mjseqnp = p' + 1 denote the total number of coefficients in the model including the model intercept). Therefore, \mjseqn\beta_j denotes how the average true effect/outcome changes for a one-unit increase in \mjeqnx_ijx_ij and the model intercept \mjseqn\beta_0 denotes the average true effect/outcome when the values of all moderator variables are equal to zero. The value of \mjseqn\tau^2 in the mixed-effects model denotes the amount of ‘residual heterogeneity’ in the true effects/outcomes (i.e., the amount of variability in the true effects/outcomes that is not accounted for by the moderators included in the model).

[`rma.mh`

] The Mantel-Haenszel method provides an alternative approach for fitting equal-effects models when dealing with studies providing data in the form of \mjeqn2 \times 22x2 tables or in the form of event counts (i.e., person-time data) for two groups (Mantel & Haenszel, 1959). The method is particularly advantageous when aggregating a large number of studies with small sample sizes (the so-called sparse data or increasing strata case). The Mantel-Haenszel method is implemented in the `rma.mh`

function. It can be used in combination with risk ratios, odds ratios, risk differences, incidence rate ratios, and incidence rate differences.

[`rma.peto`

] Yet another method that can be used in the context of a meta-analysis of \mjeqn2 \times 22x2 table data is Peto's method (see Yusuf et al., 1985), implemented in the `rma.peto`

function. The method provides an estimate of the (log) odds ratio under an equal-effects model. The method is particularly advantageous when the event of interest is rare, but see the documentation of the function for some caveats.

[`rma.glmm`

] Dichotomous outcomes and event counts (based on which one can calculate outcome measures such as odds ratios, incidence rate ratios, proportions, and incidence rates) are often assumed to arise from binomial and Poisson distributed data. Meta-analytic models that are directly based on such distributions are implemented in the `rma.glmm`

function. These models are essentially special cases of generalized linear mixed-effects models (i.e., mixed-effects logistic and Poisson regression models). For \mjeqn2 \times 22x2 table data, a mixed-effects conditional logistic model (based on the non-central hypergeometric distribution) is also available. Random/mixed-effects models with dichotomous data are often referred to as ‘binomial-normal’ models in the meta-analytic literature. Analogously, for event count data, such models could be referred to as ‘Poisson-normal’ models.

[`rma.mv`

] Standard meta-analytic models assume independence between the observed effect sizes or outcomes obtained from a set of studies. This assumption is often violated in practice. Dependencies can arise for a variety of reasons. For example, the sampling errors and/or true effects/outcomes may be correlated in multiple treatment studies (e.g., when multiple treatment groups are compared with a common control/reference group, such that the data from the control/reference group is used multiple times to compute the observed effect sizes or outcomes) or in multiple endpoint studies (e.g., when more than one effect size estimate or outcome is calculated based on the same sample of subjects due to the use of multiple endpoints or response variables). Correlations in the true effects/outcomes can also arise due to other forms of clustering (e.g., when multiple effects/outcomes derived from the same author, lab, or research group may be more similar to each other than effects/outcomes derived from different authors, labs, or research groups). In ecology and related fields, the shared phylogenetic history among the organisms studied (e.g., plants, fungi, animals) can also induce correlations among the effects/outcomes. The `rma.mv`

function can be used to fit suitable meta-analytic multivariate/multilevel models to such data, so that the non-independence in the effects/outcomes is accounted for. Network meta-analyses (also called multiple/mixed treatment comparisons) can also be carried out with this function.

The metafor package is a work in progress and is updated on a regular basis with new functions and options. With `metafor.news()`

, you can read the ‘NEWS’ file of the package after installation. Comments, feedback, and suggestions for improvements are always welcome.

To cite the package, please use the following reference:

Viechtbauer, W. (2010). Conducting meta-analyses in R with the metafor package. *Journal of Statistical Software*, **36**(3), 1-48. doi: 10.18637/jss.v036.i03

The paper mentioned above is a good starting place for those interested in using the package. The purpose of the article is to provide a general overview of the package and its capabilities (as of version 1.4-0). Not all of the functions and options are described in the paper, but it should provide a useful introduction to the package. The paper can be freely downloaded from the URL given above or can be directly loaded with the command `vignette("metafor")`

.

In addition to reading the paper, carefully read this page and then the help pages for the `escalc`

and the `rma.uni`

functions (or the `rma.mh`

, `rma.peto`

, `rma.glmm`

, and/or `rma.mv`

functions if you intend to use these methods). The help pages for these functions provide links to many additional functions, which can be used after fitting a model. You can also read the entire documentation online at https://wviechtb.github.io/metafor/ (where it is nicely formatted and the output from all examples is provided).

A (pdf) diagram showing the various functions in the metafor package (and how they are related to each other) can be opened with the command `vignette("diagram")`

.

Finally, additional information about the package, several detailed analysis examples, examples of plots and figures provided by the package (with the corresponding code), some additional tips and notes, and a FAQ can be found on the package website at https://www.metafor-project.org.

Wolfgang Viechtbauer wvb@metafor-project.org

package website: https://www.metafor-project.org

author homepage: `https://www.wvbauer.com`

Suggestions on how to obtain help with using the package can found on the package website at: https://www.metafor-project.org/doku.php/help

Cooper, H., Hedges, L. V., & Valentine, J. C. (Eds.) (2009). *The handbook of research synthesis and meta-analysis* (2nd ed.). New York: Russell Sage Foundation.

Hedges, L. V., & Olkin, I. (1985). *Statistical methods for meta-analysis*. San Diego, CA: Academic Press.

Mantel, N., & Haenszel, W. (1959). Statistical aspects of the analysis of data from retrospective studies of disease. *Journal of the National Cancer Institute*, **22**(4), 719–748. `https://doi.org/10.1093/jnci/22.4.719`

Viechtbauer, W. (2010). Conducting meta-analyses in R with the metafor package. *Journal of Statistical Software*, **36**(3), 1–48. `https://doi.org/10.18637/jss.v036.i03`

Yusuf, S., Peto, R., Lewis, J., Collins, R., & Sleight, P. (1985). Beta blockade during and after myocardial infarction: An overview of the randomized trials. *Progress in Cardiovascular Disease*, **27**(5), 335–371. `https://doi.org/10.1016/s0033-0620(85)80003-7`

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.