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, ...)A two dimensional array (prey species x prey size) with the predation mortality
A MizerParams object
A matrix of species abundances (species x size).
A vector of the resource abundance by size
A list of abundances for other dynamical components of the ecosystem
The time for which to do the calculation (Not used by standard mizer rate functions but useful for extensions with time-dependent parameters.)
A two dimensional array (predator species x predator size) with the feeding level.
Unused
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(),
mizerERepro(),
mizerEReproAndGrowth(),
mizerEncounter(),
mizerFMort(),
mizerFMortGear(),
mizerFeedingLevel(),
mizerMort(),
mizerPredRate(),
mizerRDI(),
mizerRates(),
mizerResourceMort()