Calculates the density-independent rate of total egg production \(R_{di}\) (units 1/year) before density dependence, by species.

getRDI( params, n = initialN(params), n_pp = initialNResource(params), n_other = initialNOther(params), t = 0, ... )

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.) |

... | Unused |

A numeric vector the length of the number of species.

This rate is obtained by taking the per capita rate \(E_r(w)\psi(w)\) at
which energy is invested in reproduction, as calculated by `getERepro()`

,
multiplying it by the number of individuals\(N(w)\) and integrating over all sizes
\(w\) and then multiplying by the reproductive efficiency \(\epsilon\)
and dividing by the egg size `w_min`

, and by a factor of two to account for
the two sexes:
$$R_{di} = \frac{\epsilon}{2 w_{min}} \int N(w) E_r(w) \psi(w) \, dw$$

Used by `getRDD()`

to calculate the actual, density dependent rate.
See `setReproduction()`

for more details.

By default `getRDI()`

calls `mizerRDI()`

. However you can
replace this with your own alternative reproduction function. If
your function is called `"myRDI"`

then you register it in a MizerParams
object `params`

with

params <- setRateFunction(params, "RDI", "myRDI")

Your function will then be called instead of `mizerRDI()`

, with the
same arguments.

Other rate functions:
`getEGrowth()`

,
`getEReproAndGrowth()`

,
`getERepro()`

,
`getEncounter()`

,
`getFMortGear()`

,
`getFMort()`

,
`getFeedingLevel()`

,
`getMort()`

,
`getPredMort()`

,
`getPredRate()`

,
`getRDD()`

,
`getRates()`

,
`getResourceMort()`

if (FALSE) { 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 density-independent reproduction rate at a particular time step getRDI(params, n = N(sim)[15, , ], n_pp = NResource(sim)[15, ], t = 15) }