ssa.otl {GillespieSSA}

Optimized tau-leap method (OTL)

Description

Optimized tau-leap method implementation of the SSA as described by Cao et al. (2006). Should be passed as method argument for ssa().

Usage

ssa.otl(epsilon = 0.03, nc = 10, hor = NA_real_, dtf = 10, nd = 100)

Arguments

epsilon

error control parameter.

nc

number of critical reactions threshold parameter.

hor

Highest order reaction vector. There must be one entry per species in x. Must be one of 1: first-order, 2: second-order or 22: homo-dimer. If hor is NA, defaults are all second-order.

dtf

Direct method threshold factor for temporarily suspending the OTL method.

nd

number of Direct method steps to perform during an OTL suspension.

Note

Third order-reactions (S_1+S_2+S_3 \rightarrow \ldots) 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

ssa.otl(
  hor = 1,
  nc = 10,
  epsilon = .03,
  dtf = 10,
  nd = 100
)