Efficient Low-Order Approximation of First-Passage Time Distributions

Abstract

We consider the problem of computing first-passage time distributions for reaction processes modelled by master equations. We show that this generally intractable class of problems is equivalent to a sequential Bayesian inference problem for an auxiliary observation process. The solution can be approximated efficiently by solving a closed set of coupled ordinary differential equations (for the low-order moments of the process) whose size scales with the number of species. We apply it to an epidemic model and a trimerisation process, and show good agreement with stochastic simulations.

5 pages, 3 figures

References in corpus (6)

• First-passage times in complex scale-invariant media
• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Disease extinction in the presence of non-Gaussian noise
• Validity conditions for moment closure approximations in stochastic chemical kinetics
• First passage time statistics for two-channel diffusion
• Rare events in stochastic populations under bursty reproduction