Visualizations for Assessing Convergence and Mixing of Markov Chain Monte Carlo Simulations

Jaakko Peltonen, Jarkko Venna, Samuel Kaski

Research output: Contribution to journalArticleScientificpeer-review

9 Citations (Scopus)

Abstract

Bayesian inference often requires approximating the posterior distribution by Markov chain Monte Carlo sampling. The samples come from the true distribution only after the simulation has converged, which makes detecting convergence a central problem. Commonly, several simulation chains are started from different points, and their overlap is used as a measure of convergence. Convergence measures cannot tell the cause of convergence problems; it is suggested that complementing them with proper visualization will help. A novel connection is pointed out: linear discriminant analysis (LDA) minimizes the overlap of the simulation chains measured by a common multivariate convergence measure. LDA is thus justified for visualizing convergence. However, LDA makes restrictive assumptions about the chains, which can be relaxed by a recent extension called discriminative component analysis (DCA). Lastly, methods are introduced for unidentifiable models and model families with variable number of parameters, where straightforward visualization in the parameter space is not feasible.
Original languageEnglish
Pages (from-to)4453-4470
JournalCOMPUTATIONAL STATISTICS AND DATA ANALYSIS
Volume53
Issue number12
DOIs
Publication statusPublished - 7 Jul 2009
Externally publishedYes
Publication typeA1 Journal article-refereed

Fingerprint

Dive into the research topics of 'Visualizations for Assessing Convergence and Mixing of Markov Chain Monte Carlo Simulations'. Together they form a unique fingerprint.

Cite this