Minimax Mixing Time of the Metropolis-Adjusted Langevin Algorithm for Log-Concave Sampling

09/27/2021
by   Keru Wu, et al.
0

We study the mixing time of the Metropolis-adjusted Langevin algorithm (MALA) for sampling from a log-smooth and strongly log-concave distribution. We establish its optimal minimax mixing time under a warm start. Our main contribution is two-fold. First, for a d-dimensional log-concave density with condition number κ, we show that MALA with a warm start mixes in Õ(κ√(d)) iterations up to logarithmic factors. This improves upon the previous work on the dependency of either the condition number κ or the dimension d. Our proof relies on comparing the leapfrog integrator with the continuous Hamiltonian dynamics, where we establish a new concentration bound for the acceptance rate. Second, we prove a spectral gap based mixing time lower bound for reversible MCMC algorithms on general state spaces. We apply this lower bound result to construct a hard distribution for which MALA requires at least Ω̃(κ√(d)) steps to mix. The lower bound for MALA matches our upper bound in terms of condition number and dimension. Finally, numerical experiments are included to validate our theoretical results.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/08/2023

A Simple Proof of the Mixing of Metropolis-Adjusted Langevin Algorithm under Smoothness and Isoperimetry

We study the mixing time of Metropolis-Adjusted Langevin algorithm (MALA...
research
12/23/2020

Optimal dimension dependence of the Metropolis-Adjusted Langevin Algorithm

Conventional wisdom in the sampling literature, backed by a popular diff...
research
04/10/2023

When does Metropolized Hamiltonian Monte Carlo provably outperform Metropolis-adjusted Langevin algorithm?

We analyze the mixing time of Metropolized Hamiltonian Monte Carlo (HMC)...
research
04/03/2023

Improved Bound for Mixing Time of Parallel Tempering

In the field of sampling algorithms, MCMC (Markov Chain Monte Carlo) met...
research
06/10/2021

Lower Bounds on Metropolized Sampling Methods for Well-Conditioned Distributions

We give lower bounds on the performance of two of the most popular sampl...
research
02/14/2020

Local continuity of log-concave projection, with applications to estimation under model misspecification

The log-concave projection is an operator that maps a d-dimensional dist...
research
06/13/2022

On the Computational Complexity of Metropolis-Adjusted Langevin Algorithms for Bayesian Posterior Sampling

In this paper, we study the computational complexity of sampling from a ...

Please sign up or login with your details

Forgot password? Click here to reset