Conditional Monte Carlo for reaction networks
From MaRDI portal
Publication:5075688
Abstract: Reaction networks are often used to model interacting species in fields such as biochemistry and ecology. When the counts of the species are sufficiently large, the dynamics of their concentrations are typically modeled via a system of differential equations. However, when the counts of some species are small, the dynamics of the counts are typically modeled stochastically via a discrete state, continuous time Markov chain. A key quantity of interest for such models is the probability mass function of the process at some fixed time. Since paths of such models are relatively straightforward to simulate, we can estimate the probabilities by constructing an empirical distribution. However, the support of the distribution is often diffuse across a high-dimensional state space, where the dimension is equal to the number of species. Therefore generating an accurate empirical distribution can come with a large computational cost. We present a new Monte Carlo estimator that fundamentally improves on the "classical" Monte Carlo estimator described above. It also preserves much of classical Monte Carlo's simplicity. The idea is basically one of conditional Monte Carlo. Our conditional Monte Carlo estimator has two parameters, and their choice critically affects the performance of the algorithm. Hence, a key contribution of the present work is that we demonstrate how to approximate optimal values for these parameters in an efficient manner. Moreover, we provide a central limit theorem for our estimator, which leads to approximate confidence intervals for its error.
Recommendations
- Variance-reduced simulation of lattice discrete-time Markov chains with applications in reaction networks
- Robustly simulating biochemical reaction kinetics using multi-level Monte Carlo approaches
- Multilevel Monte Carlo for continuous time Markov chains, with applications in biochemical kinetics
- Conditional expectations in network traffic estimation
- Decay rates of higher-order norms of solutions to the Navier-Stokes-Landau-Lifshitz system
Cites work
- A multiple time interval finite state projection algorithm for the solution to the chemical master equation
- A solver for the stochastic master equation applied to gene regulatory networks
- Accurate chemical master equation solution using multi-finite buffers
- Complexity of Multilevel Monte Carlo Tau-Leaping
- Computational Complexity Analysis for Monte Carlo Approximations of Classically Scaled Population Processes
- Error analysis of tau-leap simulation methods
- scientific article; zbMATH DE number 1999206 (Why is no real title available?)
- Integer Programming
- Introduction to probability
- Markov Chains
- Multilevel Monte Carlo for continuous time Markov chains, with applications in biochemical kinetics
- Multiscale Modeling of Chemical Kinetics via the Master Equation
- Non-explosivity of stochastically modeled reaction networks that are complex balanced
- Product-form stationary distributions for deficiency zero chemical reaction networks
- Representations of Markov processes as multiparameter time changes
- Sliding Window Abstraction for Infinite Markov Chains
- Solving chemical master equations by adaptive wavelet compression
- Solving the chemical master equation by a fast adaptive finite state projection based on the stochastic simulation algorithm
- Solving the chemical master equation for monomolecular reaction systems analytically
- Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems, I. Effect of Inequality of Variance in the One-Way Classification
- Spectral approximation of solutions to the chemical master equation
- Stochastic analysis of biochemical systems
- Stochastic modelling for systems biology.
- Synthesis of variance
- The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior
- Time-dependent product-form Poisson distributions for reaction networks with higher order complexes
- Weak error analysis of numerical methods for stochastic models of population processes
Cited in
(8)- Automated generation of conditional moment equations for stochastic reaction networks
- A reaction network scheme which implements inference and learning for hidden Markov models
- Bootstrapping least-squares estimates in biochemical reaction networks
- Method of conditional moments (MCM) for the chemical master equation
- Synthetic likelihood method for reaction network inference
- Searching for Rare Growth Factors Using Multicanonical Monte Carlo Methods
- Variance reduction in stochastic reaction networks using control variates
- Stochastic filtering of reaction networks partially observed in time snapshots
This page was built for publication: Conditional Monte Carlo for reaction networks
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5075688)
