GillespieSSA (version 0.6.2)

ssa.otl: 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.

References

Cao et al. (2006)

See Also

GillespieSSA-package, ssa()

Examples

# NOT RUN {
ssa.otl(
  hor = 1,
  nc = 10,
  epsilon = .03,
  dtf = 10,
  nd = 100
)
# }