Sampling latent states for high-dimensional non-linear state space models with the embedded HMM method

Alexander Y. Shestopaloff , Dept. of Statistical Sciences, University of Toronto
Radford M. Neal, Dept. of Statistical Sciences and Dept. of Computer Science, University of Toronto

We propose a new scheme for selecting pool states for the embedded Hidden Markov Model (HMM) Markov Chain Monte Carlo (MCMC) method. This new scheme allows the embedded HMM method to be used for efficient sampling in state space models where the state can be high-dimensional. Previously, embedded HMM methods were only applicable to low-dimensional state-space models. We demonstrate that using our proposed pool state selection scheme, an embedded HMM sampler can have similar performance to a well-tuned sampler that uses a combination of Particle Gibbs with Backward Sampling (PGBS) and Metropolis updates. The scaling to higher dimensions is made possible by selecting pool states locally near the current value of the state sequence. The proposed pool state selection scheme also allows each iteration of the embedded HMM sampler to take time linear in the number of the pool states, as opposed to quadratic as in the original embedded HMM sampler.

Bayesian Analysis, 2017, 26 pages. Available here.


Associated references:

This is a revised version of the following technical report:

Shestopaloff, A. Y. and Neal, R. M. (2016) ``Sampling latent states for high-dimensional non-linear state space models with the embedded HMM method'', Technical Report, 21 pages: abstract, pdf.

This technical report is also available from arXiv.org.