Abstract: This work focuses on sampling from hidden Markov models (Cappe et al, 2005) whose observations have intractable density functions. We develop a new sequential Monte Carlo (Doucet et al, 2000 and Gordon et al, 1993) algorithm and a new particle marginal Metropolis-Hastings (Andrieu et al, 2010) algorithm for these purposes. We build from Jasra, et al (2013) and Whiteley, et al (2013) to construct the sequential Monte Carlo (SMC) algorithm (which we call the alive twisted particle filter). Like the alive particle filter of Jasra, et al (2013), our new SMC algorithm adopts an approximate Bayesian computation (Tavare et al, 1997) estimate of the HMM. Our alive twisted particle filter also uses a twisted proposal as in Whiteley, et al (2013) to obtain a low-variance estimate of the HMM normalising constant. We demonstrate via numerical examples that, in some scenarios, this estimate has a much lower variance than that of the estimate obtained via the alive particle filter. The low variance of this normalising constant estimate encourages the implementation of our SMC algorithm within a particle marginal Metropolis-Hastings (PMMH) scheme, and we call the resulting methodology ``alive twisted PMMH. We numerically demonstrate on a stochastic volatility model how our alive twisted PMMH can converge faster than the standard alive PMMH of Jasra, et al (2013).
Recommendations
- Twisted particle filters
- The alive particle filter and its use in particle Markov chain Monte Carlo
- Adaptive particle allocation in iterated sequential Monte Carlo via approximating meta-models
- \(\mathrm{SMC}^2\): an efficient algorithm for sequential analysis of state space models
- Particle filtering for continuous-time hidden Markov models
Cites work
- scientific article; zbMATH DE number 5919872 (Why is no real title available?)
- scientific article; zbMATH DE number 44406 (Why is no real title available?)
- scientific article; zbMATH DE number 2106098 (Why is no real title available?)
- A sequential particle algorithm that keeps the particle system alive
- A variant of importance splitting for rare event estimation: fixed number of successes
- An adaptive sequential Monte Carlo method for approximate Bayesian computation
- Approximate Bayesian Computation for Smoothing
- Bayesian non-parametric hidden Markov models with applications in genomics
- Establishing some order amongst exact approximations of MCMCs
- Filtering via approximate Bayesian computation
- Inference in hidden Markov models.
- Likelihood analysis of non-Gaussian measurement time series
- Parameter estimation for hidden Markov models with intractable likelihoods
- Particle Markov Chain Monte Carlo Methods
- Sequential Monte Carlo for rare event estimation
- Subadditive ergodic theory
- The alive particle filter and its use in particle Markov chain Monte Carlo
- Twisted particle filters
Cited in
(3)
This page was built for publication: Twisting the alive particle filter
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q292346)