English

An Adaptive Interacting Wang-Landau Algorithm for Automatic Density Exploration

Computation 2012-06-15 v3 Applications Methodology

Abstract

While statisticians are well-accustomed to performing exploratory analysis in the modeling stage of an analysis, the notion of conducting preliminary general-purpose exploratory analysis in the Monte Carlo stage (or more generally, the model-fitting stage) of an analysis is an area which we feel deserves much further attention. Towards this aim, this paper proposes a general-purpose algorithm for automatic density exploration. The proposed exploration algorithm combines and expands upon components from various adaptive Markov chain Monte Carlo methods, with the Wang-Landau algorithm at its heart. Additionally, the algorithm is run on interacting parallel chains -- a feature which both decreases computational cost as well as stabilizes the algorithm, improving its ability to explore the density. Performance is studied in several applications. Through a Bayesian variable selection example, the authors demonstrate the convergence gains obtained with interacting chains. The ability of the algorithm's adaptive proposal to induce mode-jumping is illustrated through a trimodal density and a Bayesian mixture modeling application. Lastly, through a 2D Ising model, the authors demonstrate the ability of the algorithm to overcome the high correlations encountered in spatial models.

Keywords

Cite

@article{arxiv.1109.3829,
  title  = {An Adaptive Interacting Wang-Landau Algorithm for Automatic Density Exploration},
  author = {Luke Bornn and Pierre Jacob and Pierre Del Moral and Arnaud Doucet},
  journal= {arXiv preprint arXiv:1109.3829},
  year   = {2012}
}

Comments

33 pages, 20 figures (the supplementary materials are included as appendices)

R2 v1 2026-06-21T19:06:32.208Z