Coupled Markov Chain Monte Carlo for High-Dimensional Regression with Half-T Priors

• Activity Number: 376 - SBSS Paper Competition Winners (Part 2)
• Type: Topic-Contributed
• Date/Time: Thursday, August 12, 2021 : 12:00 PM to 1:50 PM
• Sponsor: Section on Bayesian Statistical Science
• Author(s): Niloy Biswas* and Anirban Bhattacharya and Pierre Etienne Jacob and James Johndrow
• Companies: Harvard University and Texas A&M University and Harvard University and University of Pennsylvania
• Keywords: Bayesian inference; Horseshoe prior; Couplings; Parallel computation; High dimensions; MCMC

Abstract

Continuous shrinkage priors are commonly used in Bayesian analysis of high-dimensional data, due to both their computational advantages and favorable statistical properties. We develop coupled Markov chain Monte Carlo (MCMC) algorithms for Bayesian shrinkage regression in high dimensions. Following Glynn & Rhee (2014), these couplings can be used in parallel computation strategies and practical diagnostics of convergence. Focusing on a class of shrinkage priors which includes the Horseshoe, we demonstrate the scalability of the proposed couplings with data from a genome-wide association study with 2000 rows and 100,000 covariates. The results highlight the impact of the shrinkage prior on the computational efficiency of the coupling procedure, and motivate priors where the local precisions are Half-t distributed with degree of freedom larger than one, which are statistically justifiable in terms of posterior concentration, and lead to practical computations.