getRDI {mizer} | R Documentation |
Get density independent rate of egg production
Description
Calculates the density-independent rate of total egg production
(units 1/year) before density dependence, by species.
Usage
getRDI(
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 |
Details
This rate is obtained by taking the per capita rate at
which energy is invested in reproduction, as calculated by
getERepro()
,
multiplying it by the number of individuals and integrating over
all sizes
and then multiplying by the reproductive efficiency
and dividing by the egg size
w_min
, and by a factor of two
to account for the two sexes:
Used by getRDD()
to calculate the actual, density dependent rate.
See setReproduction()
for more details.
Value
A numeric vector the length of the number of species.
Your own reproduction function
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. For an example of an alternative reproduction function
see constantEggRDI()
.
See Also
Other rate functions:
getEGrowth()
,
getERepro()
,
getEReproAndGrowth()
,
getEncounter()
,
getFMort()
,
getFMortGear()
,
getFeedingLevel()
,
getMort()
,
getPredMort()
,
getPredRate()
,
getRDD()
,
getRates()
,
getResourceMort()
Examples
params <- NS_params
# 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)