Displaying 1 result from an estimated 1 matches for "dpredat".
Did you mean:
predate
2016 Apr 26
0
vectors of equations in ode / desolve
...<- rIng * Prey * Predator
GrowthPrey <- rGrow * Prey * (1 - Prey/K)
MortPredator <- rMort * Predator
dPrey <- GrowthPrey - Ingestion
# the growth of predator would depend on the ingestion of all prays
regardless of species
sumIngestion<-sum(Ingestion)
dPredator <- sumIngestion * assEff - MortPredator
# returning an array of preys does produce an error
return(list(c(dPrey, dPredator)))
})
}
# Example is made for three different species of prey
pars <- c(rIng = 0.2, # /day, rate of ingestion
# I do not know how to pass on these pa...