4.6 Article

Manifold Markov chain Monte Carlo methods for Bayesian inference in diffusion models

出版社

WILEY
DOI: 10.1111/rssb.12497

关键词

constrained dynamics; Hamiltonian Monte Carlo; partially observed diffusions; stochastic differential equations

资金

  1. Lloyd's Register Foundation programme on data-centric engineering at the Alan Turing Institute
  2. Singapore Ministry of Education Tier 2 [MOE2016-T2-2-135]
  3. Leverhulme Trust Prize
  4. Young Investigator Award Grant (NUSYIA) [FY16 P16, R-155-000-180-133]

向作者/读者索取更多资源

This paper proposes a method for inferring the posterior distribution of nonlinear diffusion processes observed at discrete times, using ideas from statistical physics and computational chemistry. The method is highly automated and applicable in a wide range of settings, outperforming other approaches in the literature.
Bayesian inference for nonlinear diffusions, observed at discrete times, is a challenging task that has prompted the development of a number of algorithms, mainly within the computational statistics community. We propose a new direction, and accompanying methodology-borrowing ideas from statistical physics and computational chemistry-for inferring the posterior distribution of latent diffusion paths and model parameters, given observations of the process. Joint configurations of the underlying process noise and of parameters, mapping onto diffusion paths consistent with observations, form an implicitly defined manifold. Then, by making use of a constrained Hamiltonian Monte Carlo algorithm on the embedded manifold, we are able to perform computationally efficient inference for a class of discretely observed diffusion models. Critically, in contrast with other approaches proposed in the literature, our methodology is highly automated, requiring minimal user intervention and applying alike in a range of settings, including: elliptic or hypo-elliptic systems; observations with or without noise; linear or non-linear observation operators. Exploiting Markovianity, we propose a variant of the method with complexity that scales linearly in the resolution of path discretisation and the number of observation times. Python code reproducing the results is available at .

作者

我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。

评论

主要评分

4.6
评分不足

次要评分

新颖性
-
重要性
-
科学严谨性
-
评价这篇论文

推荐

暂无数据
暂无数据