View source: R/stockRecruitModels.R
rickerModel | R Documentation |
This function calculates recruitment from Ricker curve with AR(1) process
(according to Peterman et al. 2003; modified to take more recent parameter-
ization). Uses parameters from arima.mle (a, -b, sig, rho in log space) with
multivariate normally distributed errors. 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. Note that by default
prevErr and rho are NULL, resulting in a standard Ricker model.
rickerModel(
S,
a,
b,
error,
rho = NULL,
prevErr = NULL,
sig = NULL,
biasCor = NULL
)
S |
A numeric vector of spawner abundances. |
a |
A numeric vector of alpha values, i.e. productivity at low spawner abundance. |
b |
A numeric vector of beta values, i.e. density dependence para- meter. |
error |
A numeric vector of recruitment deviations, typically generated
using |
rho |
A numeric vector of rho values, i.e. AR1 coefficient. outside of model using multivariate normal (or equivalent) distribution. |
prevErr |
A numeric vector representing recruitment deviations from previous brood year. |
sig |
A numeric vector of Ricker sigma values |
biasCor |
A logical TRUE/FALSE indicating if log-normal bias correction should be applied. If NULL, then default is FALSE |
Returns a list of R, a numeric representing recruit abundance, and
errNext
which is used to generate subsequent process error (i.e. next
year's prevErr.
#Spawner and recruit values represent millions of fish, stock-recruit
parameters approximate those of Fraser River sockeye salmon Chilko CU.
#without autoregressive error
rickerModel(S = 1.1, a = 1.8, b = 1.2, error = 0.3)
#with autoregressive error
rickerModel(S = 1.1, a = 1.8, b = 1.2, error = 0.3, rho = 0.2,
prevErr = 0.7)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.