Accelerating Asymptotically Exact MCMC for Computationally Intensive Models via Local Approximations
Author(s)
Conrad, Patrick R.; Marzouk, Youssef M.; Pillai, Natesh S.; Smith, Aaron
Download1402.1694v4.pdf (2.613Mb)
OPEN_ACCESS_POLICY
Open Access Policy
Creative Commons Attribution-Noncommercial-Share Alike
Terms of use
Metadata
Show full item recordAbstract
We construct a new framework for accelerating Markov chain Monte Carlo in posterior sampling problems where standard methods are limited by the computational cost of the likelihood, or of numerical models embedded therein. Our approach introduces local approximations of these models into the Metropolis-Hastings kernel, borrowing ideas from deterministic approximation theory, optimization, and experimental design. Previous efforts at integrating approximate models into inference typically sacrifice either the sampler’s exactness or efficiency; our work seeks to address these limitations by exploiting useful convergence characteristics of local approximations. We prove the ergodicity of our approximate Markov chain, showing that it samples asymptotically from the exact posterior distribution of interest. We describe variations of the algorithm that employ either local polynomial approximations or local Gaussian process regressors. Our theoretical results reinforce the key observation underlying this paper: when the likelihood has some local regularity, the number of model evaluations per MCMC step can be greatly reduced without biasing the Monte Carlo average. Numerical experiments demonstrate multiple order-of-magnitude reductions in the number of forward model evaluations used in representative ODE and PDE inference problems, with both synthetic and real data.
Date issued
2015-10Department
Massachusetts Institute of Technology. Department of Aeronautics and AstronauticsJournal
Journal of the American Statistical Association
Publisher
American Statistical Association
Citation
Conrad, Patrick R., Youssef M. Marzouk, Natesh S. Pillai, and Aaron Smith. “Accelerating Asymptotically Exact MCMC for Computationally Intensive Models via Local Approximations.” Journal of the American Statistical Association (October 21, 2015): 00–00.
Version: Author's final manuscript
ISSN
0162-1459
1537-274X