Section: New Results
Advanced Monte Carlo methods.
Metropolis Hastings algorithm in large dimension. With T. Lelièvre and B. Miasojedow, B. Jourdain considers the Random Walk Metropolis algorithm on with Gaussian proposals, and when the target probability measure is the -fold product of a one dimensional law. It is well-known that, in the limit tends to infinity, starting at equilibrium and for an appropriate scaling of the variance and of the timescale as a function of the dimension , a diffusive limit is obtained for each component of the Markov chain. They generalize this result when the initial distribution is not the target probability measure (  ). The obtained diffusive limit is the solution to a stochastic differential equation nonlinear in the sense of McKean. In  , they prove convergence to equilibrium for this equation. They also discuss practical counterparts in order to optimize the variance of the proposal distribution to accelerate convergence to equilibrium. The analysis confirms the interest of the constant acceptance rate strategy (with acceptance rate between 1/4 and 1/3).