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.
Usage
getPredRate(
params,
n = initialN(params),
n_pp = initialNResource(params),
n_other = initialNOther(params),
t = 0,
...
)
Arguments
- 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
Value
A two dimensional array (predator species x prey size), where the prey size runs over fish community plus resource spectrum.
Your own predation rate function
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
<- setRateFunction(params, "PredRate", "myPredRate") params
Your function will then be called instead of mizerPredRate()
, with
the same arguments.
See also
Other rate functions:
getEGrowth()
,
getEReproAndGrowth()
,
getERepro()
,
getEncounter()
,
getFMortGear()
,
getFMort()
,
getFeedingLevel()
,
getMort()
,
getPredMort()
,
getRDD()
,
getRDI()
,
getRates()
,
getResourceMort()