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.

Intuitive Joint Priors for Bayesian Linear Multilevel Models: The R2D2M2 prior

Bayesian Regression Using a Prior on the Model Fit: The R2-D2 Shrinkage Prior

The R2D2 Prior for Generalized Linear Mixed Models

Bayesian Multilevel Structural Equation Modeling: An Investigation into Robust Prior Distributions for the Doubly Latent Categorical Model

Intuitive joint priors for variance parameters

Neuronized Priors for Bayesian Sparse Linear Regression

makemyprior: Intuitive Construction of Joint Priors for Variance Parameters in R

Optimal prior for Bayesian inference in a constrained parameter space

A flexible class of priors for orthonormal matrices with basis function-specific structure

Log-Scale Shrinkage Priors and Adaptive Bayesian Global-Local Shrinkage Estimation

Disjunct Support Spike and Slab Priors for Variable Selection in Regression under Quasi-sparseness

Empirical priors for prediction in sparse high-dimensional linear regression

Multi-Goal Prior Selection: A Way to Reconcile Bayesian and Classical Approaches for Random Effects Models

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

Bayesian Variable Selection for Mixed Effects Model with Shrinkage Prior.

Horseshoe Priors for Sparse Dirichlet-Multinomial Models

Bayesian High-dimensional Linear Regression with Sparse Projection-posterior

Structured prior distributions for the covariance matrix in latent factor models

Amortized Bayesian Multilevel Models

Structured Mixture of Continuation-ratio Logits Models for Ordinal Regression

Bayesian singular value regularization via a cumulative shrinkage process