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.

Estimation of mean using under-reported and overdispersed count data

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

Accurate inference in negative binomial regression

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

Growth Estimators and Confidence Intervals for the Mean of Negative Binomial Random Variables with Unknown Dispersion

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

Two-Sample Inference in Highly Dispersed Negative Binomial Models

A Penalized Empirical Likelihood Approach for Estimating Population Sizes under the Negative Binomial Regression Model

Estimating the Population Mean Under Contaminated Normal Distribution Based on Deciles

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

Data Dispersion: Now You See It… Now You Don't

Exponential Dispersion Models for Overdispersed Zero-Inflated Count Data

A Bayesian nonparametric approach to correct for underreporting in count data

The Tilted Beta-Binomial Distribution in Overdispersed Data: Maximum Likelihood and Bayesian Estimation

A robust regression model for bounded count health data

Analysis of Frequency Count Data Using the Negative Binomial Distribution

A marginal moment matching approach for fitting endemic-epidemic models to underreported disease surveillance counts

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

Modeling Time-Dependent Overdispersion in Longitudinal Count Data

Small area prediction of counts under machine learning-type mixed models

A Hierarchical Framework for Correcting Under-Reporting in Count Data