Displaying 1 result from an estimated 1 matches for "minmod".
Did you mean:
mikmod
2010 Apr 06
1
estimating the starting value within a ODE using nls and lsoda
All-
I am interested in estimating a parameter that is the starting value for an ODE model.
That is, in the typical combined fitting procedure using nls and lsoda (alternatively rk4), I first defined the ODE model:
minmod <- function(t, y, parms) {
G <- y[1]
X <- y[2]
with(as.list(parms),{
I_t <- approx(time, I.input, t)$y
dG <- -1*(p1 + X)*G +p1*G_b
dX <- -1*p2*X + p3*(I_t-I_b)
list(c(dG, dX))
})
}
Then I estimated the parameters of the model using nls:
fit.rk4 <- nls(no...