A practical guide to pseudo-marginal methods for computational inference in systems biology

by   David J. Warne, et al.

For many stochastic models of interest in systems biology, such as those describing biochemical reaction networks, exact quantification of parameter uncertainty through statistical inference is intractable. Likelihood-free computational inference techniques enable parameter inference when the likelihood function for the model is intractable but the generation of many sample paths is feasible through stochastic simulation of the forward problem. The most common likelihood-free method in systems biology is approximate Bayesian computation that accepts parameters that result in low discrepancy between stochastic simulations and measured data. However, it can be difficult to assess how the accuracy of the resulting inferences are affected by the choice of acceptance threshold and discrepancy function. The pseudo-marginal approach is an alternative likelihood-free inference method that utilises a Monte Carlo estimate of the likelihood function. This approach has several advantages, particularly in the context of noisy, partially observed, time-course data typical in biochemical reaction network studies. Specifically, the pseudo-marginal approach facilitates exact inference and uncertainty quantification, and may be efficiently combined with particle filters for low variance, high-accuracy likelihood estimation. In this review, we provide a practical introduction to the pseudo-marginal approach using inference for biochemical reaction networks as a series of case studies. Implementations of key algorithms and examples are provided using the Julia programming language; a high performance, open source programming language for scientific computing.


Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes

Models of stochastic processes are widely used in almost all fields of s...

Particle Methods for Stochastic Differential Equation Mixed Effects Models

Parameter inference for stochastic differential equation mixed effects m...

Bayesian Computation with Intractable Likelihoods

This article surveys computational methods for posterior inference with ...

Correlated pseudo-marginal schemes for time-discretised stochastic kinetic models

Performing fully Bayesian inference for the reaction rate constants gove...

Augmented pseudo-marginal Metropolis-Hastings for partially observed diffusion processes

We consider the problem of inference for nonlinear, multivariate diffusi...

Generalised likelihood profiles for models with intractable likelihoods

Likelihood profiling is an efficient and powerful frequentist approach f...

Rapid Bayesian inference for expensive stochastic models

Almost all fields of science rely upon statistical inference to estimate...

Please sign up or login with your details

Forgot password? Click here to reset