31
6
v1v2 (latest)

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 nn, and that the acceptance rates should be tuned to 0.2340.234 and 0.5740.574. We establish counterexamples to demonstrate that smoothness assumptions of the order of C1(R)\mathcal{C}^1(\mathbb{R}) for RWM and C3(R)\mathcal{C}^3(\mathbb{R}) for MALA are indeed required if these scaling rates are to hold. The counterexamples identify classes of marginal targets for which these guidelines are violated, obtained by perturbing a standard Normal density (at the level of the potential for RWM and the second derivative of the potential for MALA) using roughness generated by a path of fractional Brownian motion with Hurst exponent HH. 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