getRDD: Get density dependent recruitment

Description Usage Arguments Value Examples

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.

getRDD method with rdi argument.

getRDD method without rdi argument.

Usage

1
2
3
4
5
6
7
8
9
getRDD(object, n, n_pp, rdi, ...)

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

## S4 method for signature '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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
## 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,])

## End(Not run)

drfinlayscott/mizer documentation built on May 18, 2018, 9:20 a.m.