getRDD: Get density dependent reproduction rate

In mizer: Dynamic Multi-Species Size Spectrum Modelling

Get density dependent reproduction rate

Description

Calculates the density dependent rate of egg 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 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 getRDI().
...	Unused

Value

A numeric vector the length of the number of species.

See Also

getRDI()

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)

mizer documentation built on Oct. 17, 2024, 9:07 a.m.