Delayed acceptance particle MCMC for exact inference in stochastic kinetic models

Abstract

Recently-proposed particle MCMC methods provide a flexible way of performing Bayesian inference for parameters governing stochastic kinetic models defined as Markov (jump) processes (MJPs). Each iteration of the scheme requires an estimate of the marginal likelihood calculated from the output of a sequential Monte Carlo scheme (also known as a particle filter). Consequently, the method can be extremely computationally intensive. We therefore aim to avoid most instances of the expensive likelihood calculation through use of a fast approximation. We consider two approximations: the chemical Langevin equation diffusion approximation (CLE) and the linear noise approximation (LNA). Either an estimate of the marginal likelihood under the CLE, or the tractable marginal likelihood under the LNA can be used to calculate a first step acceptance probability. Only if a proposal is accepted under the approximation do we then run a sequential Monte Carlo scheme to compute an estimate of the marginal likelihood under the true MJP and construct a second stage acceptance probability that permits exact (simulation based) inference for the MJP. We therefore avoid expensive calculations for proposals that are likely to be rejected. We illustrate the method by considering inference for parameters governing a Lotka-Volterra system, a model of gene expression and a simple epidemic process.

Statistics and Computing (to appear)

References in corpus (1)

• The pseudo-marginal approach for efficient Monte Carlo computations

Cited by in corpus (11)

• Speeding Up MCMC by Delayed Acceptance and Data Subsampling
• Importance sampling type estimators based on approximate marginal MCMC
• Rapid Bayesian inference for expensive stochastic models
• Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes
• Importance sampling correction versus standard averages of reversible MCMCs in terms of the asymptotic variance
• Computing Bayes: From Then 'Til Now'
• Expectation propagation for continuous time stochastic processes
• A survey of Monte Carlo methods for noisy and costly densities with application to reinforcement learning and ABC
• A linear noise approximation for stochastic epidemic models fit to partially observed incidence counts
• Accelerating sequential Monte Carlo with surrogate likelihoods
• A persistent-homology-based Bayesian prior for potential coefficient reconstruction in an elliptic PDE