ssa.otl.diag {GillespieSSA}R Documentation

Optimized tau-leap method (OTL) for nu-diagonalized systems

Description

Optimized tau-leap method for nu-diagonalized systems.

Usage

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

Arguments

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.

Details

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.

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.

See Also

ssa.otl,

Examples

## Not intended to be invoked stand alone

[Package GillespieSSA version 0.5-3 Index]