powered by
This is a helper function to tranform GillesieSSA-style paramters to GillespieSSA2.
port_reactions(x0, a, nu)
A set of reaction()s to be used by ssa().
reaction()
ssa()
The x0 parameter of GillespieSSA::ssa().
x0
GillespieSSA::ssa()
The a parameter of GillespieSSA::ssa().
a
The nu parameter of GillespieSSA::ssa().
nu
x0 <- c(Y1 = 1000, Y2 = 1000) a <- c("c1*Y1","c2*Y1*Y2","c3*Y2") nu <- matrix(c(+1,-1,0,0,+1,-1),nrow=2,byrow=TRUE) port_reactions(x0, a, nu)
Run the code above in your browser using DataLab