Abstract:Abstract Background In meta-analyses of a binary outcome, double zero events in some studies cause a critical methodology problem. The generalized linear mixed model (GLMM) has been proposed as a valid statistical tool for pooling such data. Three parameter estimation methods, including the Laplace approximation (LA), penalized quasi-likelihood (PQL) and adaptive Gauss–Hermite quadrature (AGHQ) were frequently used in the GLMM. However, the performance of GLMM via these estimation methods is unclear in meta-analysis with zero events. Methods A simulation study was conducted to compare the performance. We fitted five random-effects GLMMs and estimated the results through the LA, PQL and AGHQ methods, respectively. Each scenario conducted 20,000 simulation iterations. The data from Cochrane Database of Systematic Reviews were collected to form the simulation settings. The estimation methods were compared in terms of the convergence rate, bias, mean square error, and coverage probability. Results Our results suggested that when the total events were insufficient in either of the arms, the GLMMs did not show good point estimation to pool studies of rare events. The AGHQ method did not show better properties than the LA estimation in terms of convergence rate, bias, coverage, and possibility to produce very large odds ratios. In addition, although the PQL had some advantages, it was not the preferred option due to its low convergence rate in some situations, and the suboptimal point and variance estimation compared to the LA. Conclusion The GLMM is an alternative for meta-analysis of rare events and is especially useful in the presence of zero-events studies, while at least 10 total events in both arms is recommended when employing GLMM for meta-analysis. The penalized quasi-likelihood and adaptive Gauss–Hermite quadrature are not superior to the Laplace approximation for rare events and thus they are not recommended.

Approximations to the Log-Likelihood Function in the Nonlinear Mixed-Effects Model

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

Approximate Inferences for Nonlinear Mixed Effects Models with Scale Mixtures of Skew-Normal Distributions

Maximum likelihood estimation of nonlinear mixed-effects models with crossed random effects by combining first-order conditional linearization and sequential quadratic programming

Essential formulae for restricted maximum likelihood and its derivatives associated with the linear mixed models

On the maximum likelihood estimation in general log-linear models

Restricted maximum likelihood estimation in generalized linear mixed models

Maximum Approximated Likelihood Estimation

Asymptotic Properties of Estimation in Generalized Partial Linear Mixed Models

Linear mixed models with marginally symmetric nonparametric random effects

An MCMC method for the evaluation of the Fisher information matrix for non-linear mixed effect models

Flexible Bayesian quantile regression for nonlinear mixed effects models based on the generalized asymmetric Laplace distribution

Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations

Approximate Marginal Likelihood Inference in Mixed Models for Grouped Data

Linear Mixed-Effects Models for Non-Gaussian Repeated Measurement Data

An Approximate Restricted Likelihood Ratio Test for Variance Components in Generalized Linear Mixed Models

Laplace approximation, penalized quasi-likelihood, and adaptive Gauss–Hermite quadrature for generalized linear mixed models: towards meta-analysis of binary outcome with sparse data

Computing log-likelihood and its derivatives for restricted maximum likelihood methods

Inference for High-Dimensional Linear Mixed-Effects Models: A Quasi-Likelihood Approach

Approximate Laplace approximations for scalable model selection

Fast selection of nonlinear mixed effect models using penalized likelihood