Exact active subspace Metropolis–Hastings, with applications to the Lorenz-96 system

  • We consider the application of active subspaces to inform a Metropolis-Hastings algorithm, thereby aggressively reducing the computational dimension of the sampling problem. We show that the original formulation, as proposed by Constantine, Kent, and Bui-Thanh (SIAM J. Sci. Comput., 38(5):A2779-A2805, 2016), possesses asymptotic bias. Using pseudo-marginal arguments, we develop an asymptotically unbiased variant. Our algorithm is applied to a synthetic multimodal target distribution as well as a Bayesian formulation of a parameter inference problem for a Lorenz-96 system.

Export metadata

Additional Services

Share in Twitter Search Google Scholar
Author:Ingmar Schuster, Paul Constantine, T. J. Sullivan
Document Type:Article
Year of first publication:2017
ArXiv Id:http://arxiv.org/abs/1712.02749