Extending the multi-level method for the simulation of stochastic biological systems

Abstract

The multi-level method for discrete state systems, first introduced by Anderson and Higham [Multiscale Model. Simul. 10:146--179, 2012], is a highly efficient simulation technique that can be used to elucidate statistical characteristics of biochemical reaction networks. A single point estimator is produced in a cost-effective manner by combining a number of estimators of differing accuracy in a telescoping sum, and, as such, the method has the potential to revolutionise the field of stochastic simulation. The first term in the sum is calculated using an approximate simulation algorithm, and can be calculated quickly but is of significant bias. Subsequent terms successively correct this bias by combining estimators from approximate stochastic simulations algorithms of increasing accuracy, until a desired level of accuracy is reached. In this paper we present several refinements of the multi-level method which render it easier to understand and implement, and also more efficient. Given the substantial and complex nature of the multi-level method, the first part of this work (Sections 2 - 5) is written as a tutorial, with the aim of providing a practical guide to its use. The second part (Sections 6 - 8) takes on a form akin to a research article, thereby providing the means for a deft implementation of the technique, and concludes with a discussion of a number of open problems.

38 pages

Cited by in corpus (12)

• Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art
• Multilevel rejection sampling for approximate Bayesian computation
• Multifidelity multilevel Monte Carlo to accelerate approximate Bayesian parameter inference for partially observed stochastic processes
• Tensor product approach to modelling epidemics on networks
• Uniformisation techniques for stochastic simulation of chemical reaction networks
• Robustly simulating biochemical reaction kinetics using multi-level Monte Carlo approaches
• Slow-scale split-step tau-leap method for stiff stochastic chemical systems
• Low variance couplings for stochastic models of intracellular processes with time-dependent rate functions
• Importance sampling for a robust and efficient multilevel Monte Carlo estimator for stochastic reaction networks
• Tensor product algorithms for inference of contact network from epidemiological data
• Coupling sample paths to the partial thermodynamic limit in stochastic chemical reaction networks
• Efficient simulation techniques for biochemical reaction networks