Description Usage Arguments Details Value See Also Examples
This function creates an object of class plk
(subclass of ode
), which holds the basic information of the Power Law Kinetics system in question.
1 2 
A 
The matrix of powers (pxd). Here d the number of species. 
s 

r 
An object of class 
rx0 
An object of class 
Power Law Kinetics is a class of ODE systems, having the following vector field:
\frac{dx}{dt} = θ x^A
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 columnwise).
An object with S3 class "plk" and "ode".
ode, numsolve, field
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19  # Power law system
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)
p < plk(A)
# Solve system
numsolve(p, Time, x0, theta)
# Evaluate field
field(p, x0, theta)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.