Displaying 1 result from an estimated 1 matches for "pkmod".
Did you mean:
kmod
2013 Apr 21
1
lsoda question from deSolve package
...following data and function:
require(deSolve)
times <- c(0:24)
tin <- 0.5
D <- 400
V <- 26.3
k <-0.056
k12 <- 0.197118
k21 <- 0.022665
yini <- c(dy1 = 0,dy2 = 0)
events <- data.frame(var = "dy1",time = c(10,15),value = c(200,100),method = "add")
pkmod <- function(t, y, p) {
if (t < tin) R <- (D/tin) else R <- 0
dy1 <- R - (p["k12"]+p["k"])* y[1] + p["k21"]*y[2]
dy2 <- p["k12"]* y[1] - p["k21"] *y[2]
list(c(dy1, dy2))
}
p <- c(k=k, k12=k12, k21=k21, V=V)
result <-...