Calculates the density-independent rate of total egg production \(R_{di}\) (units 1/year) before density dependence, by species.
getRDI( 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 numeric vector the length of the number of species.
This rate is obtained by taking the per capita rate \(E_r(w)\psi(w)\) at
which energy is invested in reproduction, as calculated by getERepro(),
multiplying it by the number of individuals\(N(w)\) and integrating over all sizes
\(w\) and then multiplying by the reproductive efficiency \(\epsilon\)
and dividing by the egg size w_min, and by a factor of two to account for
the two sexes:
$$R_{di} = \frac{\epsilon}{2 w_{min}} \int N(w) E_r(w) \psi(w) \, dw$$
Used by getRDD() to calculate the actual, density dependent rate.
See setReproduction() for more details.
By default getRDI() calls mizerRDI(). However you can
replace this with your own alternative reproduction function. If
your function is called "myRDI" then you register it in a MizerParams
object params with
params <- setRateFunction(params, "RDI", "myRDI")
Your function will then be called instead of mizerRDI(), with the
same arguments.
Other rate functions:
getEGrowth(),
getEReproAndGrowth(),
getERepro(),
getEncounter(),
getFMortGear(),
getFMort(),
getFeedingLevel(),
getMort(),
getPredMort(),
getPredRate(),
getRDD(),
getRates(),
getResourceMort()
if (FALSE) { params <- newMultispeciesParams(NS_species_params_gears, inter) # Project with constant fishing effort for all gears for 20 time steps sim <- project(params, t_max = 20, effort = 0.5) # Get the density-independent reproduction rate at a particular time step getRDI(params, n = N(sim)[15, , ], n_pp = NResource(sim)[15, ], t = 15) }