Description Usage Arguments Details Value See Also Examples
This function creates an object of class rlk
(subclass of ode
), which holds the basic information of the Rational Law Kinetics system in question.
1 2 |
A |
The matrix of powers (pxd). Here d the number of species. |
B |
The matrix of powers (pxd). Here d the number of species. |
s |
|
r |
An object of class |
rx0 |
An object of class |
Rational Law Kinetics is a class of ODE systems, having the following vector field:
\frac{dx}{dt} = θ (x^A / (1 + x^B))
with x^A = (∏_{i=1}^dx_i^{A_{ji}})_{j=1}^p and θ an estimatable parameter matrix of dimension dxp. By convention theta will only be reported as a vector (concatinated column-wise).
An object with S3 class "rlk" and "ode".
ode, numsolve, field
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | # Rational law kinetics
A <- matrix(
c(1, 0, 0, 0,
0, 1, 2, 0,
0, 0, 0, 1), ncol = 4, byrow = TRUE)
theta <- matrix(
c(0, 2, -0.5, 0,
1, 0, 0, 1,
-1, -1, -1, -1), ncol = 3, byrow = TRUE)
x0 <- c(X = 1, Y = 4, Z = 0.1, W = 0.1)
time <- seq(0, 1, by = .1)
r <- rlk(A, A[c(2, 1, 3), ])
# Solve system
numsolve(o = r, time = time, x0 = x0, param = theta)
# Evaluate field
field(r, x0, theta)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.