期刊
STATISTICAL SCIENCE
卷 28, 期 3, 页码 424-446出版社
INST MATHEMATICAL STATISTICS
DOI: 10.1214/13-STS421
关键词
MCMC; Bayesian nonparametrics; algorithms; Gaussian random field; Bayesian inverse problems
资金
- EPSRC
- ERC [239870]
- St. Cross College
- ONR
- EPSRC [EP/F050798/1, EP/K014463/1] Funding Source: UKRI
- Engineering and Physical Sciences Research Council [EP/K014463/1, EP/F050798/1, EP/D002060/1] Funding Source: researchfish
- European Research Council (ERC) [239870] Funding Source: European Research Council (ERC)
Many problems arising in applications result in the need to probe a probability distribution for functions. Examples include Bayesian nonparametric statistics and conditioned diffusion processes. Standard MCMC algorithms typically become arbitrarily slow under the mesh refinement dictated by nonparametric description of the unknown function. We describe an approach to modifying a whole range of MCMC methods, applicable whenever the target measure has density with respect to a Gaussian process or Gaussian random field reference measure, which ensures that their speed of convergence is robust under mesh refinement. Gaussian processes or random fields are fields whose marginal distributions, when evaluated at any finite set of N points, are R-N-valued Gaussians. The algorithmic approach that we describe is applicable not only when the desired probability measure has density with respect to a Gaussian process or Gaussian random field reference measure, but also to some useful non-Gaussian reference measures constructed through random truncation. In the applications of interest the data is often sparse and the prior specification is an essential part of the overall modelling strategy. These Gaussian-based reference measures are a very flexible modelling tool, finding wide-ranging application. Examples are shown in density estimation, data assimilation in fluid mechanics, subsurface geophysics and image registration. The key design principle is to formulate the MCMC method so that it is, in principle, applicable for functions; this may be achieved by use of proposals based on carefully chosen time-discretizations of stochastic dynamical systems which exactly preserve the Gaussian reference measure. Taking this approach leads to many new algorithms which can be implemented via minor modification of existing algorithms, yet which show enormous speed-up on a wide range of applied problems.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据