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.

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

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

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

Random Effects Modeling and the Zero-Inflated Poisson Distribution

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

Exponential Dispersion Models for Overdispersed Zero-Inflated Count Data

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

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

Time Series Regression for Zero-Inflated and Overdispersed Count Data: A Functional Response Model Approach

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

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

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

A non-parametric model to address overdispersed count response in a longitudinal data setting with missingness.

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

Too many zeros and/or highly skewed? A tutorial on modelling health behaviour as count data with Poisson and negative binomial regression

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

Modeling Time-Dependent Overdispersion in Longitudinal Count Data

Semiparametric analysis of longitudinal zero-inflated count data

A general averaging method for count data with overdispersion and/or excess zeros in biomedicine

Estimation of mean using under-reported and overdispersed count data