scholarly journals Asymptotic variance for random walk Metropolis chains in high dimensions: logarithmic growth via the Poisson equation

2019 ◽  
Vol 51 (4) ◽  
pp. 994-1026
Author(s):  
Aleksandar Mijatović ◽  
Jure Vogrinc

AbstractThere are two ways of speeding up Markov chain Monte Carlo algorithms: (a) construct more complex samplers that use gradient and higher-order information about the target and (b) design a control variate to reduce the asymptotic variance. While the efficiency of (a) as a function of dimension has been studied extensively, this paper provides the first results linking the efficiency of (b) with dimension. Specifically, we construct a control variate for a d-dimensional random walk Metropolis chain with an independent, identically distributed target using the solution of the Poisson equation for the scaling limit in [30]. We prove that the asymptotic variance of the corresponding estimator is bounded above by a multiple of $\log(d)/d$ over the spectral gap of the chain. The proof hinges on large deviations theory, optimal Young’s inequality and Berry–Esseen-type bounds. Extensions of the result to non-product targets are discussed.

2017 ◽  
Vol 54 (4) ◽  
pp. 1233-1260 ◽  
Author(s):  
Alain Durmus ◽  
Sylvain Le Corff ◽  
Eric Moulines ◽  
Gareth O. Roberts

Abstract In this paper we consider the optimal scaling of high-dimensional random walk Metropolis algorithms for densities differentiable in the Lp mean but which may be irregular at some points (such as the Laplace density, for example) and/or supported on an interval. Our main result is the weak convergence of the Markov chain (appropriately rescaled in time and space) to a Langevin diffusion process as the dimension d goes to ∞. As the log-density might be nondifferentiable, the limiting diffusion could be singular. The scaling limit is established under assumptions which are much weaker than the one used in the original derivation of Roberts et al. (1997). This result has important practical implications for the use of random walk Metropolis algorithms in Bayesian frameworks based on sparsity inducing priors.


Sign in / Sign up

Export Citation Format

Share Document