ResearchTrend.AI
  • Papers
  • Communities
  • Events
  • Blog
  • Pricing
Papers
Communities
Social Events
Terms and Conditions
Pricing
Parameter LabParameter LabTwitterGitHubLinkedInBlueskyYoutube

© 2025 ResearchTrend.AI, All rights reserved.

  1. Home
  2. Papers
  3. 2203.05351
54
4
v1v2 (latest)

Multi-index Sequential Monte Carlo ratio estimators for Bayesian Inverse problems

10 March 2022
K. Law
N. Walton
Shan Yang
Ajay Jasra
ArXiv (abs)PDFHTML
Abstract

We consider the problem of estimating expectations with respect to a target distribution with an unknown normalizing constant, and where even the unnormalized target needs to be approximated at finite resolution. This setting is ubiquitous across science and engineering applications, for example in the context of Bayesian inference where a physics-based model governed by an intractable partial differential equation (PDE) appears in the likelihood. A multi-index Sequential Monte Carlo (MISMC) method is used to construct ratio estimators which provably enjoy the complexity improvements of multi-index Monte Carlo (MIMC) as well as the efficiency of Sequential Monte Carlo (SMC) for inference. In particular, the proposed method provably achieves the canonical complexity of MSE−1^{-1}−1, while single level methods require MSE−ξ^{-\xi}−ξ for ξ>1\xi>1ξ>1. This is illustrated on examples of Bayesian inverse problems with an elliptic PDE forward model in 111 and 222 spatial dimensions, where ξ=5/4\xi=5/4ξ=5/4 and ξ=3/2\xi=3/2ξ=3/2, respectively. It is also illustrated on a more challenging log Gaussian process models, where single level complexity is approximately ξ=9/4\xi=9/4ξ=9/4 and multilevel Monte Carlo (or MIMC with an inappropriate index set) gives ξ=5/4+ω\xi = 5/4 + \omegaξ=5/4+ω, for any ω>0\omega > 0ω>0, whereas our method is again canonical.

View on arXiv
Comments on this paper