Get density independent rate of egg production

Calculates the density independent rate of egg production \(R_{p.i}\) (units 1/year) before density dependence, by species. Used by getRDD to calculate the actual density dependent rate. See setReproduction for more details.

getRDI(
  params,
  n = params@initial_n,
  n_pp = params@initial_n_pp,
  n_other = params@initial_n_other,
  e_repro = getERepro(params, n = n, n_pp = n_pp, n_other = n_other),
  sex_ratio = 0.5
)

Arguments

params	A MizerParams object
n	A matrix of species abundances (species x size).
n_pp	A vector of the plankton abundance by size
n_other	A list of abundances for other dynamical components of the ecosystem
e_repro	The energy available for reproduction (optional). A matrix of size no. species x no. size bins. If not supplied, is calculated internally using getERepro.
sex_ratio	Proportion of the population that is female. Default value is 0.5.

Value

A numeric vector the length of the number of species

See also

getRDD

Other rate functions: getEGrowth(), getEReproAndGrowth(), getERepro(), getEncounter(), getFMortGear(), getFMort(), getFeedingLevel(), getMort(), getPlanktonMort(), getPredMort(), getPredRate(), getRDD(), getRates(), getStarvMort()

Examples

if (FALSE) {
data(NS_species_params_gears)
data(inter)
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 recruitment at a particular time step
getRDI(params,sim@n[21,,],sim@n_pp[21,])
}