Abstract:Background: Two characteristics of commonly used outcomes in medical research are zero inflation and non-negative integers; examples include the number of hospital admissions or emergency department visits, where the majority of patients will have zero counts. Zero-inflated regression models were devised to analyze this type of data. However, the performance of zero-inflated regression models or the properties of data best suited for these analyses have not been thoroughly investigated. Methods: We conducted a simulation study to evaluate the performance of two generalized linear models, negative binomial and zero-inflated negative binomial, for analyzing zero-inflated count data. Simulation scenarios assumed a randomized controlled trial design and varied the true underlying distribution, sample size, and rate of zero inflation. We compared the models in terms of bias, mean squared error, and coverage. Additionally, we used logistic regression to determine which data properties are most important for predicting the best-fitting model. Results: We first found that, regardless of the rate of zero inflation, there was little difference between the conventional negative binomial and its zero-inflated counterpart in terms of bias of the marginal treatment group coefficient. Second, even when the outcome was simulated from a zero-inflated distribution, a negative binomial model was favored above its ZI counterpart in terms of the Akaike Information Criterion. Third, the mean and skewness of the non-zero part of the data were stronger predictors of model preference than the percentage of zero counts. These results were not affected by the sample size, which ranged from 60 to 800. Conclusions: We recommend that the rate of zero inflation and overdispersion in the outcome should not be the sole and main justification for choosing zero-inflated regression models. Investigators should also consider other data characteristics when choosing a model for count data. In addition, if the performance of the NB and ZINB regression models is reasonably comparable even with ZI outcomes, we advocate the use of the NB regression model due to its clear and straightforward interpretation of the results.

Semiparametric analysis of longitudinal zero-inflated count data

The Quadratic Inference Functions in Measurement Error Model for Longitudinal Data

Variable selection approach for zero-inflated count data via adaptive lasso

On performance of parametric and distribution‐free models for zero‐inflated and over‐dispersed count responses

Random Effects Modeling and the Zero-Inflated Poisson Distribution

A simulation study of the performance of statistical models for count outcomes with excessive zeros

Joint modeling of zero-inflated longitudinal measurements and time-to-event outcomes with applications to dynamic prediction

Parameter Estimations on Zero-Inflated Negative Binomial Model with Incomplete Data

Evaluation of negative binomial and zero-inflated negative binomial models for the analysis of zero-inflated count data: application to the telemedicine for children with medical complexity trial

Modeling Zero-inflated Count Data Using Generalized Poisson and Ordinal Logistic Regression Models in Medical Research

Multilevel modeling in single-case studies with zero-inflated and overdispersed count data

Zero-inflated Poisson Models with Censored Data and Its Application

Statistical model for overdispersed count outcome with many zeros: an approach for direct marginal inference

A bivariate zero-inflated negative binomial model and its applications to biomedical settings

A Semiparametric Change-Point Regression Model for Longitudinal Observations

Modeling Time-Dependent Overdispersion in Longitudinal Count Data

Semiparametric Bayesian Analysis for Longitudinal Mixed Effects Models with Non-Normal AR(1) Errors

A comparison of statistical methods for modeling count data with an application to hospital length of stay

Semiparametric and Nonparametric Regression Analysis of Longitudinal Data - Rejoinder

Semiparametric Regression Analysis of Longitudinal Data with Informative Drop-Outs

A nonparametric maximum likelihood estimation for biased-sampling data with zero-inflated truncation