Analyzing MCMC Output

Abstract

Markov chain Monte Carlo (MCMC) is a sampling-based method for estimating features of probability distributions. MCMC methods produce a serially correlated, yet representative, sample from the desired distribution. As such it can be difficult to know when the MCMC method is producing reliable results. We introduce some fundamental methods for ensuring a trustworthy simulation experiment. In particular, we present a workflow for output analysis in MCMC providing estimators, approximate sampling distributions, stopping rules, and visualization tools.

References in corpus (9)

• MCMC using Hamiltonian dynamics
• General state space Markov chains and MCMC algorithms
• Batch means and spectral variance estimators in Markov chain Monte Carlo
• Markov Chain Monte Carlo Estimation of Quantiles
• Unbiased Markov chain Monte Carlo with couplings
• Exact sampling for intractable probability distributions via a Bernoulli factory
• Kernel estimators of asymptotic variance for adaptive Markov chain Monte Carlo
• Quantifying Uncertainty in Transdimensional Markov Chain Monte Carlo Using Discrete Markov Models
• Estimating accuracy of the MCMC variance estimator: a central limit theorem for batch means estimators