pith. sign in

arxiv: math/0407281 · v1 · submitted 2004-07-15 · 🧮 math.PR · math.ST· stat.TH

Improving Asymptotic Variance of MCMC Estimators: Non-reversible Chains are Better

classification 🧮 math.PR math.STstat.TH
keywords chainnon-reversibleasymptoticchainsreversiblevariancemcmcstate
0
0 comments X
read the original abstract

I show how any reversible Markov chain on a finite state space that is irreducible, and hence suitable for estimating expectations with respect to its invariant distribution, can be used to construct a non-reversible Markov chain on a related state space that can also be used to estimate these expectations, with asymptotic variance at least as small as that using the reversible chain (typically smaller). The non-reversible chain achieves this improvement by avoiding (to the extent possible) transitions that backtrack to the state from which the chain just came. The proof that this modification cannot increase the asymptotic variance of an MCMC estimator uses a new technique that can also be used to prove Peskun's (1973) theorem that modifying a reversible chain to reduce the probability of staying in the same state cannot increase asymptotic variance. A non-reversible chain that avoids backtracking will often take little or no more computation time per transition than the original reversible chain, and can sometime produce a large reduction in asymptotic variance, though for other chains the improvement is slight. In addition to being of some practical interest, this construction demonstrates that non-reversible chains have a fundamental advantage over reversible chains for MCMC estimation. Research into better MCMC methods may therefore best be focused on non-reversible chains.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Adaptive tuning of Hamiltonian Monte Carlo methods

    stat.CO 2025-06 conditional novelty 5.0

    ATune combines Gaussian theoretical analysis with burn-in simulation data to select system-specific splitting integrators and hyperparameter credible intervals for improved HMC stability and performance.