Meta-analysis with the glmmTMB R package

TL;DR

This paper introduces a new covariance structure 'equalto' in the glmmTMB R package, enabling explicit modeling of sampling variances in meta-analyses, thus expanding the package's flexibility.

Contribution

It presents the 'equalto' covariance structure in glmmTMB, allowing users to incorporate known sampling variances directly into meta-analytical models.

Findings

The implementation produces estimates identical to the metafor package.

It enables modeling heteroscedasticity and dependence among sampling errors.

Demonstrated with meta-analyses from medicine, ecology, and social sciences.

Abstract

Meta-analytical models are typically formulated as a mixed-effects model where the sampling variances of the effect sizes are treated as known. In principle, such models could be fitted with standard mixed-modelling software such as the glmmTMB R package. This general-purpose package for generalized linear mixed models (GLMMs) provides flexibility in distributions and random effect covariance structures through the Template Model Builder (TMB). However, incorporating known sampling variances in the conventional inverse-variance formulation of meta-analysis was previously not easily accomplished in glmmTMB. Here, we introduce equalto, a new covariance structure in glmmTMB that allows users to supply a known sampling error variance-covariance matrix when fitting meta-analytic models. This enables explicit modelling of heteroscedasticity and dependence among sampling errors. The new implementation provides an alternative way to fit meta-analytic models, convenient for users already familiar with glmmTMB. Using simulations, we show that the new implementation produces model estimates identical to those from the established metafor package and illustrate its applicability with published meta-analyses in medicine, evolutionary ecology, and the social sciences. Further, this novel implementation in glmmTMB supports more flexible modelling of meta-analytical data, expanding the R toolkit available for evidence synthesis.