R/project_methods.R
mizerPredMort.Rd
Calculates the total predation mortality rate \(\mu_{p,i}(w_p)\) (in units
of 1/year) on each prey species by prey size:
$$\mu_{p.i}(w_p) = \sum_j {\tt pred\_rate}_j(w_p)\, \theta_{ji}.$$
You would not usually call this
function directly but instead use getPredMort()
, which then calls this
function unless an alternative function has been registered, see below.
mizerPredMort(params, n, n_pp, n_other, t, pred_rate, ...)
params | A MizerParams object |
---|---|
n | A matrix of species abundances (species x size). |
n_pp | A vector of the resource abundance by size |
n_other | A list of abundances for other dynamical components of the ecosystem |
t | The time for which to do the calculation (Not used by standard mizer rate functions but useful for extensions with time-dependent parameters.) |
pred_rate | A two dimensional array (predator species x predator size) with the feeding level. |
... | Unused |
A two dimensional array (prey species x prey size) with the predation mortality
By default getPredMort()
calls mizerPredMort()
. However you can
replace this with your own alternative predation mortality function. If
your function is called "myPredMort"
then you register it in a MizerParams
object params
with
params <- setRateFunction(params, "PredMort", "myPredMort")
Your function will then be called instead of mizerPredMort()
, with the
same arguments.
Other mizer rate functions:
mizerEGrowth()
,
mizerEReproAndGrowth()
,
mizerERepro()
,
mizerEncounter()
,
mizerFMortGear()
,
mizerFMort()
,
mizerFeedingLevel()
,
mizerMort()
,
mizerPredRate()
,
mizerRDI()
,
mizerRates()
,
mizerResourceMort()