Generalised Gibbs sampler and multigrid Monte Carlo for Bayesian computation

Although Monte Carlo methods have frequently been applied with success, indiscriminate use of Markov chain Monte Carlo leads to unsatisfactory performances in numerous applications. We present a generalised version of the Gibbs sampler that is based on conditional moves along the traces of groups of transformations in the sample space. We explore its connection with the multigrid Monte Carlo method and its use in designing more efficient samplers. The generalised Gibbs sampler provides a framework encompassing a class of recently proposed tricks such as parameter expansion and reparameterisation. To illustrate, we apply this new method to Bayesian inference problems for nonlinear state-space models, ordinal data and stochastic differential equations with discrete observations.

[1]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[2]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[3]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[4]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[5]  Wang,et al.  Nonuniversal critical dynamics in Monte Carlo simulations. , 1987, Physical review letters.

[6]  Goodman,et al.  Multigrid Monte Carlo method. Conceptual foundations. , 1989, Physical review. D, Particles and fields.

[7]  Adrian F. M. Smith,et al.  Sampling-Based Approaches to Calculating Marginal Densities , 1990 .

[8]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[9]  J. Besag,et al.  Spatial Statistics and Bayesian Computation , 1993 .

[10]  Jun S. Liu,et al.  Detecting subtle sequence signals: a Gibbs sampling strategy for multiple alignment. , 1993, Science.

[11]  S. Chib,et al.  Bayesian analysis of binary and polychotomous response data , 1993 .

[12]  Walter R. Gilks,et al.  Adaptive Direction Sampling , 1994 .

[13]  Jun S. Liu,et al.  Covariance structure of the Gibbs sampler with applications to the comparisons of estimators and augmentation schemes , 1994 .

[14]  C. Geyer,et al.  Annealing Markov chain Monte Carlo with applications to ancestral inference , 1995 .

[15]  A. Pedersen A new approach to maximum likelihood estimation for stochastic differential equations based on discrete observations , 1995 .

[16]  G. Roberts,et al.  Updating Schemes, Correlation Structure, Blocking and Parameterization for the Gibbs Sampler , 1997 .

[17]  M. Pitt,et al.  Likelihood analysis of non-Gaussian measurement time series , 1997 .

[18]  Jun S. Liu,et al.  Simulated Sintering: Markov Chain Monte Carlo With Spaces of Varying Dimensions , 2007 .

[19]  Xiao-Li Meng,et al.  Seeking efficient data augmentation schemes via conditional and marginal augmentation , 1999 .

[20]  William H. Press,et al.  Numerical recipes in C , 2002 .