By explicitly representing well-stirred chemical reaction systems with independent, unit  Poisson processes we develop a new adaptive tau-leaping procedure. The procedure developed is novel in that we enforce any leap condition via a post-leap check as opposed to performing a pre-leap tau selection. Further, we perform the post-leap check in such a way that the statistics of the sample paths generated will not be skewed by the rejection of a leap. By performing a post-leap check to ensure accuracy, the method developed in this paper is guaranteed to never produce negative population values . The efficiency of the method developed here is demonstrated on a model of a decaying dimer .
By explicitly representing the reaction times of discrete chemical systems as the firing times of independent, unit rate Poisson processes we develop a new adaptive tau-leaping procedure. The procedure developed is novel in that accuracy is guaranteed by performing post-leap checks. Because the representation we use separates the randomness of the model from the state of the system, we are able to perform the post-leap checks in such a way that the statistics of the sample paths generated will not be skewed by the rejections of leaps. Further, since any leap condition is ensured with a probability of one, the simulation method naturally avoids negative population values .
R
well-stirred chemical reaction systems with
the reaction times of discrete chemical systems as the firing times of
A
rate
R
we enforce any leap condition via a
accuracy is guaranteed by performing
R
check as opposed to performing a pre-leap tau selection. Further, we
checks. Because the representation we use separates the randomness of the model from the state of the system, we are able to
R
check
checks
R
rejection of a leap. By performing a post-leap check to ensure accuracy, the method developed in this paper is guaranteed to never produce
rejections of leaps. Further, since any leap condition is ensured with a probability of one, the simulation method naturally avoids
D
. The efficiency of the method developed here is demonstrated on a model of a decaying dimer