Abstract:Count data regression modeling has received much attention in several science fields in which the Poisson, Negative binomial, and Zero-Inflated models are some of the primary regression techniques. Negative binomial regression is applied to modeling count variables, usually when they are over-dispersed. A Poisson distribution is also utilized for counting data where the mean is equal to the variance. This situation is often unrealistic since the distribution of counts will usually have a variance that is not equal to its mean. Modeling it as Poisson distributed leads to ignoring under- or overdispersion, depending on if the variance is smaller or larger than the mean. Also, situations with outcomes having a larger number of zeros such as RNASeq data require Zero-inflated models. Variable selection through shrinkage priors has been a popular method to address the curse of dimensionality and achieve the identification of significant variables. We present a unified Bayesian hierarchical framework that implements and compares shrinkage priors in negative-binomial and zero-inflated negative binomial regression models. The key feature is the representation of the likelihood by a Polya-Gamma data augmentation, which admits a natural integration with a family of shrinkage priors. We specifically focus on the Horseshoe, Dirichlet Laplace, and Double Pareto priors. Extensive simulation studies address the efficiency of the model and mean square errors are reported. Further, the models are applied to data sets such as the Covid-19 vaccine, and Covid-19 RNA-Seq data among others. The models are robust enough to address variable selection, and MSE decreases as the sample size increases, having lower errors in p>n cases. The noteworthy results showed that the adverse events of Covid-19 vaccines were dependent on age, recovery, medical history, and prior vaccination with a remarkable reduction in MSE of the fitted values. No. of publications of Ph.D. students were dependent on the no. of children, and the no. of articles in the last three years. The models are robust enough to conduct both variable selections and produce effective fit because of their high shrinkage property and applicability to a broad range of biometric and public health high dimensional problems.

Empirical Analysis of Zipf's Law, Power Law, and Lognormal Distributions in Medical Discharge Reports

Power-Law Distributions in Empirical Data

Power-law distributions in binned empirical data

Large-Scale Analysis of Zipf’s Law in English Texts

A Novel Discrete Linear-Exponential Distribution for Modeling Physical and Medical Data

A new Bayesian regression model for counts in medicine

A heteroscedastic Bayesian generalized logistic regression model with application to scaling problems

Co-occurrence of the Benford-like and Zipf Laws Arising from the Texts Representing Human and Artificial Languages

Bayesian Analysis of Population Health Data

Bayesian Shrinkage Priors in Zero-Inflated and Negative Binomial Regression models with Real World Data Applications of COVID-19 Vaccine, and RNA-Seq

Degree distributions in networks: beyond the power law

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

Are there too many uncited articles? Zero inflated variants of the discretised lognormal and hooked power law distributions

A NEW flexible exponent power family of distributions with biomedical data analysis

Log-logistic Distribution as a Reliability Model: A Bayesian Analysis

POWER-LAW DISTRIBUTIONS BASED ON EXPONENTIAL DISTRIBUTIONS: LATENT SCALING, SPURIOUS ZIPF'S LAW, AND FRACTAL RABBITS

Mixed effect Poisson log-linear models for clinical and epidemiological sleep hypnogram data

Lognormals, Power Laws and Double Power Laws in the Distribution of Frequencies of Harmonic Codewords from Classical Music

Abstract 13965: Natural Language Processing of Hospitalization Discharge Summary to Predict 1-year Post-Discharge Mortality Among Patients With Acute Heart Failure

A Bayesian Analysis for the Parameters of the Exponential-Logarithmic Distribution

A Bayesian multivariate spatial approach for illness-death survival models