mizer (version 1.0.1)

getRDD: Get density dependent recruitment

Description

Calculates the density dependent recruitment (total egg production) \(R_i\) for each species. This is the flux entering the smallest size class of each species. The density dependent recruitment is the density independent recruitment after it has been put through the density dependent stock-recruitment relationship function. This method is used by the project method for performing simulations.

Usage

getRDD(object, n, n_pp, rdi, ...)
# S4 method for MizerParams,matrix,numeric,matrix
getRDD(object, n, n_pp, rdi,
  sex_ratio = 0.5)
# S4 method for MizerParams,matrix,numeric,missing
getRDD(object, n, n_pp,
  sex_ratio = 0.5)

Arguments

object

An MizerParams object

n

A matrix of species abundance (species x size)

n_pp

A vector of the background abundance by size

rdi

A matrix of density independent recruitment (optional) with dimensions no. sp x 1. If not specified rdi is calculated internally using the getRDI method.

...

Other arguments (currently unused).

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.

Examples

# NOT RUN {
data(NS_species_params_gears)
data(inter)
params <- MizerParams(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 energy at a particular time step
getRDD(params,sim@n[21,,],sim@n_pp[21,])
# }