getRDD {mizer} | R Documentation |
Get density dependent reproduction rate
Description
Calculates the density dependent rate of egg production (units
1/year) for each species. This is the flux entering the smallest size class
of each species. The density dependent rate is the density independent
rate obtained with
getRDI()
after it has been put through the
density dependence function. This is the Beverton-Holt function
BevertonHoltRDD()
by default, but this can be changed. See
setReproduction()
for more details.
Usage
getRDD(
params,
n = initialN(params),
n_pp = initialNResource(params),
n_other = initialNOther(params),
t = 0,
rdi = getRDI(params, n = n, n_pp = n_pp, n_other = n_other, t = t),
...
)
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.) |
rdi |
A vector of density-independent reproduction rates for each
species. If not specified, rdi is calculated internally using
|
... |
Unused |
Value
A numeric vector the length of the number of species.
See Also
Other rate functions:
getEGrowth()
,
getERepro()
,
getEReproAndGrowth()
,
getEncounter()
,
getFMort()
,
getFMortGear()
,
getFeedingLevel()
,
getMort()
,
getPredMort()
,
getPredRate()
,
getRDI()
,
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 rate at a particular time step
getRDD(params, n = N(sim)[15, , ], n_pp = NResource(sim)[15, ], t = 15)