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.

Asymptotic Bayes Optiamlity for Sparse Count Data

Posterior Contraction rate and Asymptotic Bayes Optimality for one-group shrinkage priors in sparse normal means problem

Asymptotic bayes optimality under sparsity for equicorrelated multivariate normal test statistics

Asymptotic Bayes' optimality under sparsity for exchangeable dependent multivariate normal test statistics

High-dimensional Grouped-regression using Bayesian Sparse Projection-posterior

Bayesian High-dimensional Linear Regression with Sparse Projection-posterior

Tail-adaptive Bayesian shrinkage

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

Bayesian Shrinkage Estimation of Credible Subgroups for Count Data with Excess Zeros

Testing the homogeneity of risk differences with sparse count data

Optimal Gaussian approximations to the posterior for log-linear models with Diaconis-Ylvisaker priors

The Bridged Posterior: Optimization, Profile Likelihood and a New Approach to Generalized Bayes

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

Optimal post-selection inference for sparse signals: a nonparametric empirical-Bayes approach

NAPA: Neighborhood-Assisted and Posterior-Adjusted Two-sample Inference

Consistent Group selection using Global-local prior in High dimensional setup

High-dimensional prediction for count response via sparse exponential weights

Empirical Bayes large-scale multiple testing for high-dimensional binary outcome data

Shrinkage with Robustness: Log-Adjusted Priors for Sparse Signals

On SURE-Type Shrinkage Estimation