Regenerative Simulation for the Bayesian Lasso

06/06/2018
by   Y. -L. Chen, et al.
0

The Gibbs sampler of Park and Casella is one of the most popular MCMC methods for sampling from the posterior density of the Bayesian Lasso regression. As with many Markov chain samplers, their Gibbs sampler lacks a theoretically sound method of output analysis --- a method for estimating the variance of a given ergodic average and estimating how closely the chain is sampling from the stationary distribution, that is, the burn-in. In this paper, we address this shortcoming by identifying regenerative structure in the sampler of Park and Casella, thus providing a theoretically sound method of assessing its performance. The regenerative structure provides both a strongly consistent variance estimator, and an estimator of (an upper bound on) the total variation distance from the target posterior density. The result is a simple and theoretically sound way to assess the stationarity of the Park and Casella and, more generally, other MCMC samplers, for which regenerative simulation is possible. We perform a numerical study in which we validate the standard errors calculated by our regenerative method by comparing it with the standard errors calculated by an AR(1) heuristic approximation. Thus, we show that for the Bayesian Lasso model, the regenerative method is a viable and theoretically justified alternative to the existing ad-hoc MCMC convergence diagnostics.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/01/2022

Approximate Gibbs sampler for Bayesian Huberized lasso

The Bayesian lasso is well-known as a Bayesian alternative for lasso. Al...
research
04/08/2020

Posterior computation with the Gibbs zig-zag sampler

Markov chain Monte Carlo (MCMC) sampling algorithms have dominated the l...
research
04/06/2023

Variable-Complexity Weighted-Tempered Gibbs Samplers for Bayesian Variable Selection

Subset weighted-Tempered Gibbs Sampler (wTGS) has been recently introduc...
research
02/22/2018

Bayesian Lasso : Concentration and MCMC Diagnosis

Using posterior distribution of Bayesian LASSO we construct a semi-norm ...
research
01/25/2020

Particle-Gibbs Sampling For Bayesian Feature Allocation Models

Bayesian feature allocation models are a popular tool for modelling data...
research
06/28/2022

The split Gibbs sampler revisited: improvements to its algorithmic structure and augmented target distribution

This paper proposes a new accelerated proximal Markov chain Monte Carlo ...
research
10/06/2015

Bayesian Markov Blanket Estimation

This paper considers a Bayesian view for estimating a sub-network in a M...

Please sign up or login with your details

Forgot password? Click here to reset