Optimal acceptance rates for Metropolis algorithms: Moving beyond 0.234
Recent optimal scaling theory has produced a condition for the asymptotically optimal acceptance rate of Metropolis algorithms to be the well-known 0.234 when applied to certain multi-dimensional target distributions. These d-dimensional target distributions are formed of independent components, each of which is scaled according to its own function of d. We show that when the condition is not met the limiting process of the algorithm is altered, yielding an asymptotically optimal acceptance rate which might drastically differ from the usual 0.234. Specifically, we prove that as d-->[infinity] the sequence of stochastic processes formed by say the i*th component of each Markov chain usually converges to a Langevin diffusion process with a new speed measure [upsilon], except in particular cases where it converges to a one-dimensional Metropolis algorithm with acceptance rule [alpha]*. We also discuss the use of inhomogeneous proposals, which might prove to be essential in specific cases.
Year of publication: |
2008
|
---|---|
Authors: | Bédard, Mylène |
Published in: |
Stochastic Processes and their Applications. - Elsevier, ISSN 0304-4149. - Vol. 118.2008, 12, p. 2198-2222
|
Publisher: |
Elsevier |
Keywords: | Weak convergence Optimal scaling Langevin diffusion Generator Markov chain Monte Carlo |
Saved in:
Saved in favorites
Similar items by person
-
Scaling analysis of multiple-try MCMC methods
Bédard, Mylène, (2012)
- More ...