Calculates the density-independent rate of total egg production \(R_{p.i}\)
(units 1/year) before density dependence, by species. This is obtained by
taking the total rate at which energy is invested in reproduction,
multiplying by the reproductive efficiency erepro and dividing by the egg
size w_min, and by a factor of two to account for the two sexes.
You would not usually call this
function directly but instead use getRDI(), which then calls this
function unless an alternative function has been registered, see below.
mizerRDI(params, n, n_pp, n_other, t, e_growth, mort, e_repro, ...)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.)
An array (species x size) with the energy available for
growth as calculated by getEGrowth(). Unused.
An array (species x size) with the mortality rate as calculated
by getMort(). Unused.
An array (species x size) with the energy available for
reproduction as calculated by getERepro().
Unused
A numeric vector with the rate of egg production for each species.
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.
Used by getRDD() to calculate the actual, density dependent rate.
See setReproduction() for more details.
Other mizer rate functions:
mizerEGrowth(),
mizerEReproAndGrowth(),
mizerERepro(),
mizerEncounter(),
mizerFMortGear(),
mizerFMort(),
mizerFeedingLevel(),
mizerMort(),
mizerPredMort(),
mizerPredRate(),
mizerRates(),
mizerResourceMort()