ssa.otl {GillespieSSA}R Documentation

Optimized tau-leap method (OTL)

Description

Optimized tau-leap method implementation of the SSA as described by Cao et al. (2006).

Usage

ssa.otl(x,a,nu,hor,nc,epsilon,dtf,nd)

Arguments

x state vector.
a vector of evaluated propensity functions.
nu 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.

Details

Performs one time step using the Explicit tau-leap method. Intended to be invoked by ssa.

Value

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.

Note

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.

References

Cao et al. (2006)

See Also

GillespieSSA-package, ssa

Examples

## Not intended to be invoked stand alone

[Package GillespieSSA version 0.5-0 Index]