The online home for the publications of the American Statistical Association
139
Views
7
CrossRef citations
0
Altmetric

Theory and Methods

Auxiliary Variable Methods for Markov Chain Monte Carlo with Applications

Pages 585-595 | Received 01 Jul 1996, Published online: 17 Feb 2012
 
Translator disclaimer

Suppose that one wishes to sample from the density π(x) using Markov chain Monte Carlo (MCMC). An auxiliary variable u and its conditional distribution π(u|x) can be defined, giving the joint distribution π(x, u) = π(x)π(u|x). A MCMC scheme that samples over this joint distribution can lead to substantial gains in efficiency compared to standard approaches. The revolutionary algorithm of Swendsen and Wang is one such example. Besides reviewing the Swendsen-Wang algorithm and its generalizations, this article introduces a new auxiliary variable method called partial decoupling. Two applications in Bayesian image analysis are considered: a binary classification problem in which partial decoupling out performs Swendsen-Wang and single-site Metropolis methods, and a positron emission tomography (PET) reconstruction that uses the gray level prior of Geman and McClure. A generalized Swendsen–Wang algorithm is developed for this problem, which reduces the computing time to the point where MCMC is a viable method of posterior exploration.

Key Words: Please supply.
  • PrivacyBadger has replaced this AddThis button.

People also read

Original text