On the utility of Metropolis-Hastings with asymmetric acceptance ratio

03/26/2018
by   Christophe Andrieu, et al.
0

The Metropolis-Hastings algorithm allows one to sample asymptotically from any probability distribution π. There has been recently much work devoted to the development of variants of the MH update which can handle scenarios where such an evaluation is impossible, and yet are guaranteed to sample from π asymptotically. The most popular approach to have emerged is arguably the pseudo-marginal MH algorithm which substitutes an unbiased estimate of an unnormalised version of π for π. Alternative pseudo-marginal algorithms relying instead on unbiased estimates of the MH acceptance ratio have also been proposed. These algorithms can have better properties than standard PM algorithms. Convergence properties of both classes of algorithms are known to depend on the variability of the estimators involved and reduced variability is guaranteed to decrease the asymptotic variance of ergodic averages and will shorten the burn-in period, or convergence to equilibrium, in most scenarios of interest. A simple approach to reduce variability, amenable to parallel computations, consists of averaging independent estimators. However, while averaging estimators of π in a pseudo-marginal algorithm retains the guarantee of sampling from π asymptotically, naive averaging of acceptance ratio estimates breaks detailed balance, leading to incorrect results. We propose an original methodology which allows for a correct implementation of this idea. We establish theoretical properties which parallel those available for standard PM algorithms and discussed above. We demonstrate the interest of the approach on various inference problems. In particular we show that convergence to equilibrium can be significantly shortened, therefore offering the possibility to reduce a user's waiting time in a generic fashion when a parallel computing architecture is available.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
12/29/2020

Metropolis-Hastings with Averaged Acceptance Ratios

Markov chain Monte Carlo (MCMC) methods to sample from a probability dis...
research
06/26/2018

Large Sample Asymptotics of the Pseudo-Marginal Method

The pseudo-marginal algorithm is a variant of the Metropolis-Hastings al...
research
12/07/2017

Exact active subspace Metropolis-Hastings, with applications to the Lorenz-96 system

We consider the application of active subspaces to inform a Metropolis-H...
research
12/10/2021

Comparison of Markov chains via weak Poincaré inequalities with application to pseudo-marginal MCMC

We investigate the use of a certain class of functional inequalities kno...
research
05/27/2021

Unbiased estimation of equilibrium, rates, and committors from Markov state model analysis

Markov state models (MSMs) have been broadly adopted for analyzing molec...
research
10/29/2019

Asymptotically Unbiased Generative Neural Sampling

We propose a general framework for the estimation of observables with ge...

Please sign up or login with your details

Forgot password? Click here to reset