Displaying 1 result from an estimated 1 matches for "alphab".
Did you mean:
alpha
2012 Apr 25
1
Using apply() with a function involving ode()
...n for each
set of parameter values. I am having trouble getting the apply()
function to work.
Here is an example. Consider an epidemiology model of West Nile Virus:
wnv.model <- function(Time, State, Pars){
with(as.list(c(State, Pars)), {
dSB <- -(alphaB*betaB*SB*IM)/(SB + IB)
dIB <- (alphaB*betaB*SB*IM)/(SB + IB) - deltaB*IB
dSM <- bM*(SM + EM + IM) - (alphaM*betaB*IB*SM)/(SB
+ IB) - dM*SM
dEM <- (alphaM*betaB*SM*IB)/(SB + IB) - kappaM*EM -
dM*EM
dIM <...