Abstract:Background: Hospital length of stay (LOS) is a key indicator of hospital care management efficiency, cost of care, and hospital planning. Hospital LOS is often used as a measure of a post-medical procedure outcome, as a guide to the benefit of a treatment of interest, or as an important risk factor for adverse events. Therefore, understanding hospital LOS variability is always an important healthcare focus. Hospital LOS data can be treated as count data, with discrete and non-negative values, typically right skewed, and often exhibiting excessive zeros. In this study, we compared the performance of the Poisson, negative binomial (NB), zero-inflated Poisson (ZIP), and zero-inflated negative binomial (ZINB) regression models using simulated and empirical data. Methods: Data were generated under different simulation scenarios with varying sample sizes, proportions of zeros, and levels of overdispersion. Analysis of hospital LOS was conducted using empirical data from the Medical Information Mart for Intensive Care database. Results: Results showed that Poisson and ZIP models performed poorly in overdispersed data. ZIP outperformed the rest of the regression models when the overdispersion is due to zero-inflation only. NB and ZINB regression models faced substantial convergence issues when incorrectly used to model equidispersed data. NB model provided the best fit in overdispersed data and outperformed the ZINB model in many simulation scenarios with combinations of zero-inflation and overdispersion, regardless of the sample size. In the empirical data analysis, we demonstrated that fitting incorrect models to overdispersed data leaded to incorrect regression coefficients estimates and overstated significance of some of the predictors. Conclusions: Based on this study, we recommend to the researchers that they consider the ZIP models for count data with zero-inflation only and NB models for overdispersed data or data with combinations of zero-inflation and overdispersion. If the researcher believes there are two different data generating mechanisms producing zeros, then the ZINB regression model may provide greater flexibility when modeling the zero-inflation and overdispersion.

Random Effects Modeling and the Zero-Inflated Poisson Distribution

Random effects compound Poisson model to represent data with extra zeros

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

A New Regression Model for the Analysis of Overdispersed and Zero-Modified Count Data

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

Zero truncated Poisson model: an alternative approach for analyzing count data with excess zeros

The consequences of checking for zero-inflation and overdispersion in the analysis of count data

Exponential Dispersion Models for Overdispersed Zero-Inflated Count Data

Semiparametric analysis of longitudinal zero-inflated count data

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

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

A simulation study comparing likelihood and non-likelihood approaches in analyzing overdispersed count data

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

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

Detecting overdispersion in count data: A zero-inflated Poisson regression analysis

Modeling data with zero inflation and overdispersion using GAMLSSs

Zero-inflated Beta distribution regression modeling

A Flexible Zero-Inflated Poisson-Gamma model with application to microbiome read counts

Semiparametric Bivariate Zero-Inflated Poisson Models with Application to Studies of Abundance for Multiple Species

Predicting the multivariate zero‐inflated counts: A novel model averaging method under Pearson loss

The Zero-Inflated Poisson - Probit regression model: a new model for count data