2009
DOI: 10.1063/1.3078490
|View full text |Cite
|
Sign up to set email alerts
|

An exact accelerated stochastic simulation algorithm

Abstract: An exact method for stochastic simulation of chemical reaction networks, which accelerates the stochastic simulation algorithm ͑SSA͒, is proposed. The present "ER-leap" algorithm is derived from analytic upper and lower bounds on the multireaction probabilities sampled by SSA, together with rejection sampling and an adaptive multiplicity for reactions. The algorithm is tested on a number of well-quantified reaction networks and is found experimentally to be very accurate on test problems including a chaotic re… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1
1

Citation Types

0
36
0

Year Published

2009
2009
2019
2019

Publication Types

Select...
7

Relationship

3
4

Authors

Journals

citations
Cited by 16 publications
(36 citation statements)
references
References 18 publications
0
36
0
Order By: Relevance
“…For example, the bounds that Cao et al [7] have proposed may provide a promising avenue for this particular endeavour. Ongoing work aims to extend the present algorithm to spatially developing systems [21] with multiresolution capabilities [4] with delays and in extending the recently proposed exact R-leaping algorithm [18] to systems with delays. …”
Section: Discussionmentioning
confidence: 98%
“…For example, the bounds that Cao et al [7] have proposed may provide a promising avenue for this particular endeavour. Ongoing work aims to extend the present algorithm to spatially developing systems [21] with multiresolution capabilities [4] with delays and in extending the recently proposed exact R-leaping algorithm [18] to systems with delays. …”
Section: Discussionmentioning
confidence: 98%
“…Equations (19) and (20) imply that smaller values of p result in much better approximation. However, a value for p that is too small can make the geometric sampling described by (14) inefficient, as the cumulative probabilities of the form P r(Y s k ≤ n) will involve too many terms.…”
Section: Quality Of Approximationmentioning
confidence: 96%
“…Several accelerated methods have been proposed that are either exact or approximate. Exact methods typically involve optimisations over the standard algorithm, such as the next reaction method [8], the optimised DM [5], the logarithmic DM [17] and ER-leap [19]. Most of these approaches involve the use of appropriate data structures in order to generate the simulation events efficiently.…”
Section: Introductionmentioning
confidence: 99%
“…However, this method has the following drawbacks, such as, it is easily influenced by the random sample size, low efficient and inaccurate calculation results (Peer and Sharma 2007;Mjolsness et al 2009;Sakalauskas 2002). In this article, to solve stochastic models accurately and precisely, a novel two-phase approach is proposed, that is, initially, according to the feature of different removal probability density functions, disassembly probability density functions of feasible disassembly paths are determined by a time-domain method or frequency method, and additionally, after disassembly probability density functions have been obtained, typical stochastic models in the different disassembly decision-making are solved by a numerical solution method.…”
Section: The Minimum Expected Value Model Of Disassembly Timementioning
confidence: 99%