rickerSurvModel: Generate recruit abundance with the Ricker Model with a...

In TESA-workshops/TESAsamSim: Closed-loop salmon simulator for 2021 TESA workshop

Description

This function calculates recruitment from Ricker curve with a marine survival covariate that is specific to brood year added in. Note that internal if statements prevent it from being vectorized so must be passed single values, i.e. all vectors for inputs and outputs are length 1.

Usage

rickerSurvModel(S, a, b, error, ppnAges, gamma, mSurvAtAge)

Arguments

S	A numeric of spawner abundance.
a	A numeric of alpha value, i.e. productivity at low spawner abundance.
b	A numeric vector of beta value, i.e. density dependence para- meter.
error	A numeric recruitment deviation, typically generated using rmvnorm() and relevant process variance estimates (sigma).
ppnAges	A numeric vector of proportion of spawner abundance at age
gamma	A numeric value represting the marine survival coefficient
mSurvAtAge	A numeric vector representing marine survival covariates at age

Value

A list of recruitment at each age for the modelled brood year (e.g., R2 = recruitment to age 2 from broodyear, R3 = recruitment to age 3, etc), as well as recruitment for all return ages combined from the modelled brood year (RecBY)

Examples

rickerSurvModel(S = 1000, a = 2.1, b = 0.00001, error = 0.8, ppnAges = c(0,0.83,0.17,0,0), gamma=0.4, mSurvAtAge=c(0,0.013, 0.015,0,0))

TESA-workshops/TESAsamSim documentation built on Feb. 6, 2021, 12:25 a.m.