Conjugate Bayes for probit regression via unified skew-normal distributions
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel pith:R5KRYF4Srecord.jsonopen to challenge →
read the original abstract
Regression models for dichotomous data are ubiquitous in statistics. Besides being useful for inference on binary responses, these methods serve also as building blocks in more complex formulations, such as density regression, nonparametric classification and graphical models. Within the Bayesian framework, inference proceeds by updating the priors for the coefficients, typically set to be Gaussians, with the likelihood induced by probit or logit regressions for the responses. In this updating, the apparent absence of a tractable posterior has motivated a variety of computational methods, including Markov Chain Monte Carlo routines and algorithms which approximate the posterior. Despite being routinely implemented, Markov Chain Monte Carlo strategies face mixing or time-inefficiency issues in large p and small n studies, whereas approximate routines fail to capture the skewness typically observed in the posterior. This article proves that the posterior distribution for the probit coefficients has a unified skew-normal kernel, under Gaussian priors. Such a novel result allows efficient Bayesian inference for a wide class of applications, especially in large p and small-to-moderate n studies where state-of-the-art computational methods face notable issues. These advances are outlined in a genetic study, and further motivate the development of a wider class of conjugate priors for probit models along with methods to obtain independent and identically distributed samples from the unified skew-normal posterior.
This paper has not been read by Pith yet.
Forward citations
Cited by 2 Pith papers
-
End-to-end probabilistic hierarchical forecasting of large hierarchies via probabilistic top-down
e2eTD forecasts a small subset of aggregate series (~0.3% of hierarchy) and propagates joint samples via probabilistic top-down disaggregation to produce coherent probabilistic forecasts, achieving lowest weighted sca...
-
Model selection with proper scoring rules on data sets of time series: prefer the mean scaled score
Mean scaled score is recommended over rank-based aggregation for model selection on time series datasets because skewness causes non-mean criteria to select misspecified models with short tests.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.