Skip to contents

Get density dependent reproduction rate

Source: R/rate_functions.R
getRDD.Rd

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(object, ...)

Arguments

object

A MizerParams or MizerSim object.

...

Additional arguments that depend on the class of object.

For a MizerParams object:

n

A matrix of species abundances (species x size). Defaults to the initial abundances stored in object.

n_pp

A vector of the resource abundance by size. Defaults to the initial resource abundance stored in object.

n_other

A named list of the abundances of other dynamical components. Defaults to the initial values stored in object.

t

The time for which to do the calculation. Defaults to 0.

rdi

A vector of density-independent reproduction rates for each species. If not specified, it is calculated internally using getRDI().

For a MizerSim object:

time_range

The time range over which to return the rates. Either a vector of values, a vector of min and max time, or a single value. Defaults to the whole time range of the simulation.

Value

  • MizerParams: A numeric vector the length of the number of species.

  • MizerSim: An ArrayTimeBySpecies object (time x species).

See also

getRDI()

Other rate functions: getDiffusion(), getEGrowth(), getERepro(), getEReproAndGrowth(), getEncounter(), getFMort(), getFMortGear(), getFeedingLevel(), getFlux(), getMort(), getPredMort(), getPredRate(), getRDI(), getRates(), getResourceMort()

Examples

# \donttest{
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)
#>        Sprat      Sandeel       N.pout      Herring          Dab      Whiting 
#> 7.273519e+11 4.098406e+11 9.405199e+12 1.103948e+12 1.115039e+10 5.394960e+11 
#>         Sole      Gurnard       Plaice      Haddock          Cod       Saithe 
#> 3.855063e+10 7.536194e+11 2.654321e+13 1.733106e+12 8.259403e+09 1.117118e+11 
# }