Get density dependent rate of larvae production

Calculates the density dependent rate of larvae production \(R_i\) (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 dependent "stock-recruitment" relationship function. See setReproduction for more details.

getRDD(
  params,
  n = params@initial_n,
  n_pp = params@initial_n_pp,
  n_other = params@initial_n_other,
  sex_ratio = 0.5,
  rdi = getRDI(params, n = n, n_pp = n_pp, n_other = n_other, sex_ratio = sex_ratio)
)

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
sex_ratio	Proportion of the population that is female. Default value is 0.5.
rdi	A vector of density independent recruitment for each species. If not specified rdi is calculated internally using getRDI.

Value

A numeric vector the length of the number of species.

See also

getRDI

Other rate functions: getEGrowth(), getEReproAndGrowth(), getERepro(), getEncounter(), getFMortGear(), getFMort(), getFeedingLevel(), getMort(), getPlanktonMort(), getPredMort(), getPredRate(), getRDI(), 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 rate at a particular time step
getRDD(params,sim@n[21,,],sim@n_pp[21,])
}