ssa.otl: Optimized tau-leap method (OTL)

In GillespieSSA: Gillespie's Stochastic Simulation Algorithm (SSA)

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 ---> ...) 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
)

GillespieSSA documentation built on March 18, 2022, 7:55 p.m.