ssa.otl.diag {GillespieSSA} | R Documentation |
Optimized tau-leap method for nu-diagonalized systems.
ssa.otl.diag(x,a,nu_tile,hor,nc,epsilon,dtf,nd)
x |
state vector. |
a |
vector of evaluated propensity functions. |
nu_tile |
state-change matrix. |
hor |
highest order reaction vector (one entry per species in x ) |
nc |
number of critical reactions threshold parameter. |
epsilon |
error control parameter. |
dtf |
Direct method threshold factor for temporarily suspending the OTL method. |
nd |
number of Direct method steps to perform during an OTL suspension. |
Performs one time step using the Explicit tau-leap method. It is usually called from within ssa
, but can be invoked directly, see ssa.otl
for Examples.
A list with three elements, 1) the time leap (tau
) and 2) the realized state change vector (nu_j
), and 3) a boolean value (suspendedTauLeapMethod
) indicating if the simulation should revert to the Direct method for nd
time steps.
Third order-reactions (S_1 + S_2 + S_3 —> ...) are not supported currently since they are approximations to sets of coupled first- and second-order reactions). See Cao et al. (2006) for more details.
## Not intended to be invoked stand alone