posted on 2020-08-24, 10:46authored byAxel Finke, Adam M. Johansen, Dario Spanò
We develop particle Gibbs samplers for static-parameter estimation in discretely observed piecewise deterministic process (PDPs). PDPs are stochastic processes that jump randomly at a countable number of stopping times but otherwise evolve deterministically in continuous time. A sequential Monte Carlo (SMC) sampler for filtering in PDPs has recently been proposed. We first provide new insight into the consequences of an approximation inherent within that algorithm. We then derive a new representation of the algorithm. It simplifies ensuring that the importance weights exist and also allows the use of variance-reduction techniques known as backward and ancestor sampling. Finally, we propose a novel Gibbs step that improves mixing in particle Gibbs samplers whose SMC algorithms make use of large collections of auxiliary variables, such as many instances of SMC samplers. We provide a comparison between the two particle Gibbs samplers for PDPs developed in this paper. Simulation results indicate that they can outperform reversible-jump MCMC approaches.
Funding
Engineering and Physical Sciences Research Council (EPSRC) Doctoral Training Grant EP/J500586/1.
EPSRC Grant EP/I017984/1.
CRiSM, an EPSRC/HEFCE-funded grant.
History
School
Science
Department
Mathematical Sciences
Published in
Annals of the Institute of Statistical Mathematics
This is a post-peer-review, pre-copyedit version of an article published in Annals of the Institute of Statistical Mathematics. The final authenticated version is available online at: https://doi.org/10.1007/s10463-014-0455-z.