
arXiv: 1603.01136
This article considers the Sequential Monte Carlo (SMC) approximation of ratios of normalizing constants associated to posterior distributions which in principle rely on continuum models. Therefore, the Monte Carlo estimation error and the discrete approximation error must be balanced. A multilevel strategy is utilized to substantially reduce the cost to obtain a given error level in the approximation as compared to standard estimators. Two estimators are considered and relative variance bounds are given. The theoretical results are numerically illustrated for two Bayesian inverse problems arising from elliptic Partial Differential Equations (PDEs). The examples involve the inversion of observations of the solution of (i) a one-dimensional Poisson equation to infer the diffusion coefficient, and (ii) a two-dimensional Poisson equation to infer the external forcing.
FOS: Computer and information sciences, [MATH.MATH-PR] Mathematics [math]/Probability [math.PR], Multilevel Monte Carlo, FOS: Mathematics, Bayesian inverse problems, Mathematics - Numerical Analysis, Numerical Analysis (math.NA), Sequential Monte Carlo, Statistics - Computation, Computation (stat.CO)
FOS: Computer and information sciences, [MATH.MATH-PR] Mathematics [math]/Probability [math.PR], Multilevel Monte Carlo, FOS: Mathematics, Bayesian inverse problems, Mathematics - Numerical Analysis, Numerical Analysis (math.NA), Sequential Monte Carlo, Statistics - Computation, Computation (stat.CO)
| selected citations These citations are derived from selected sources. This is an alternative to the "Influence" indicator, which also reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | 16 | |
| popularity This indicator reflects the "current" impact/attention (the "hype") of an article in the research community at large, based on the underlying citation network. | Top 10% | |
| influence This indicator reflects the overall/total impact of an article in the research community at large, based on the underlying citation network (diachronically). | Top 10% | |
| impulse This indicator reflects the initial momentum of an article directly after its publication, based on the underlying citation network. | Top 10% |
