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

View source: R/stockRecruitModels.R

rickerSurvModelR Documentation

Generate recruit abundance with the Ricker Model with a marine survival co-variate

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,
  sig = NULL,
  biasCor = NULL
)

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

sig

A numeric vector of Ricker sigma values

biasCor

A logical TRUE/FALSE vector indicating if log-normal bias correction should be applied. If NULL, then default is FALSE

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


CamFreshwater/samSim documentation built on Sept. 25, 2023, 10:22 a.m.