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

Description Usage Arguments Note References See Also Examples

View source: R/ssa.otl.R

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

1
2
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

1
2
3
4
5
6
7
ssa.otl(
  hor = 1,
  nc = 10,
  epsilon = .03,
  dtf = 10,
  nd = 100
)

GillespieSSA documentation built on July 27, 2019, 1:02 a.m.