Bayesian parameter inference for partially observed stopped processes
From MaRDI portal
Abstract: In this article we consider Bayesian parameter inference associated to partially-observed stochastic processes that start from a set B0 and are stopped or killed at the first hitting time of a known set A. Such processes occur naturally within the context of a wide variety of applications. The associated posterior distributions are highly complex and posterior parameter inference requires the use of advanced Markov chain Monte Carlo (MCMC) techniques. Our approach uses a recently introduced simulation methodology, particle Markov chain Monte Carlo (PMCMC) (Andrieu et. al. 2010 [1]), where sequential Monte Carlo (SMC) approximations (see Doucet et. al. 2001 [18] and Liu 2001 [27]) are embedded within MCMC. However, when the parameter of interest is fixed, standard SMC algorithms are not always appropriate for many stopped processes. In Chen et. al. [11] and Del Moral 2004 [15] the authors introduce SMC approximations of multi-level Feynman-Kac formulae, which can lead to more efficient algorithms. This is achieved by devising a sequence of nested sets from B0 to A and then perform the resampling step only when the samples of the process reach intermediate level sets in the sequence. Naturally, the choice of the intermediate level sets is critical to the performance of such a scheme. In this paper, we demonstrate that multi-level SMC algorithms can be used as a proposal in PMCMC. In addition, we propose a flexible strategy that adapts the level sets for different parameter proposals. Our methodology is illustrated on the coalescent model with migration.
Recommendations
- Bayesian static parameter estimation for partially observed diffusions via multilevel Monte Carlo
- Stopping-Time Resampling for Sequential Monte Carlo Methods
- Efficient \(\mathrm{SMC}^2\) schemes for stochastic kinetic models
- Bayesian parameter inference for partially observed stochastic differential equations driven by fractional Brownian motion
- Inference and rare event simulation for stopped Markov processes via reverse-time sequential Monte Carlo
Cites work
- scientific article; zbMATH DE number 3817476 (Why is no real title available?)
- scientific article; zbMATH DE number 2106098 (Why is no real title available?)
- A nonasymptotic theorem for unnormalized Feynman-Kac particle models
- Adaptive Multilevel Splitting for Rare Event Analysis
- Coalescence time for two genes from a subdivided population
- Coupling and Ergodicity of Adaptive Markov Chain Monte Carlo Algorithms
- Estimation in discretely observed diffusions killed at a threshold
- Exact Monte Carlo simulation of killed diffusions
- Genealogical particle analysis of rare events
- Importance sampling on coalescent histories. II: Subdivided population models
- Inference in Molecular Population Genetics
- Monte Carlo strategies in scientific computing
- Multilevel splitting for estimating rare event probabilities
- On large deviations theory and asymptotically efficient Monte Carlo estimation
- On particle methods for parameter estimation in state-space models
- Quantitative non-geometric convergence bounds for independence samplers
- Sampling per mode for rare event simulation in switching diffusions
- Sequential Monte Carlo Methods in Practice
- Sequential Monte Carlo Samplers
- Sequential Monte Carlo for rare event estimation
- Simulating probability distributions in the coalescent
- Splitting for rare event simulation: A large deviation approach to design and analysis
- Stopping-Time Resampling for Sequential Monte Carlo Methods
- The pseudo-marginal approach for efficient Monte Carlo computations
Cited in
(2)
This page was built for publication: Bayesian parameter inference for partially observed stopped processes
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q892438)
