214

Counterexamples for optimal scaling of Metropolis-Hastings chains with rough target densities

Abstract

For sufficiently smooth targets of product form it is known that the variance of a single coordinate of the proposal in RWM (Random walk Metropolis) and MALA (Metropolis adjusted Langevin algorithm) should optimally scale as n1n^{-1} and as n13n^{-\frac{1}{3}} with dimension \(n\), and that the acceptance rates should be tuned to 0.2340.234 and 0.5740.574. We establish counterexamples to demonstrate that smoothness assumptions such as C1(R)\mathcal{C}^1(\mathbb{R}) for RWM and C3(R)\mathcal{C}^3(\mathbb{R}) for MALA are indeed required if these guidelines are to hold. The counterexamples identify classes of marginal targets (obtained by perturbing a standard Normal density at the level of the potential (or second derivative of the potential for MALA) by a path of fractional Brownian motion with Hurst exponent HH) for which these guidelines are violated. For such targets there is strong evidence that RWM and MALA proposal variances should optimally be scaled as n1Hn^{-\frac{1}{H}} and as n12+Hn^{-\frac{1}{2+H}} and will then obey anomalous acceptance rate guidelines. Useful heuristics resulting from this theory are discussed. The paper develops a framework capable of tackling optimal scaling results for quite general Metropolis-Hastings algorithms (possibly depending on a random environment).

View on arXiv
Comments on this paper