Displaying 1 result from an estimated 1 matches for "lvmod0d".
2010 Dec 11
2
Predator Prey Models
...00). The disturbance can be the introduction of 40 prey
(N=40) and 10 predators (Pred = 10). I would like to see my model go from a
state of equilibrium (up to t = 99), show this disturbance (at t = 100) and
then slowly work its way back to equilibrium. Does anybody know if this
could be done?
LVmod0D <- function(Time, State, Pars) {
with(as.list(c(State, Pars)), {
IngestPred <- rI * N * Pred
GrowthN <- rG * N * (1 - N/K)
MortPred <- rM * Pred
dN <- GrowthN - IngestPred
dPred <- IngestPred * AE - MortPred
return(list(c(dN, dPred)))
})
}
pars <- c(rI = 0.1, rG...