Calculates the potential rate (in units 1/year) at which a prey individual of
a given size \(w\) is killed by predators from species \(j\). In formulas
$${\tt pred\_rate}_j(w_p) = \int \phi_j(w,w_p) (1-f_j(w))
\gamma_j(w) N_j(w) \, dw.$$
This potential rate is used in getPredMort()
to
calculate the realised predation mortality rate on the prey individual.
getPredRate( params, n = initialN(params), n_pp = initialNResource(params), n_other = initialNOther(params), t = 0, ... )
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.) |
... | Unused |
A two dimensional array (predator species x prey size), where the prey size runs over fish community plus resource spectrum.
By default getPredRate()
calls mizerPredRate()
. However you can
replace this with your own alternative predation rate function. If
your function is called "myPredRate"
then you register it in a MizerParams
object params
with
params <- setRateFunction(params, "PredRate", "myPredRate")
Your function will then be called instead of mizerPredRate()
, with
the same arguments.
Other rate functions:
getEGrowth()
,
getEReproAndGrowth()
,
getERepro()
,
getEncounter()
,
getFMortGear()
,
getFMort()
,
getFeedingLevel()
,
getMort()
,
getPredMort()
,
getRDD()
,
getRDI()
,
getRates()
,
getResourceMort()