recruitment: Recruitment

In vquennessen/densityratio: Using Density Ratios to Help Manage Fisheries In And Around Marine Reserves

Recruitment

Description

recruitment calculates the number of new recruits entering the population

Usage

recruitment(
  t,
  cr,
  fdr,
  SSB,
  A = 5,
  R0 = 1e+05,
  H,
  B0,
  Eps,
  Sigma_R,
  Rec_age,
  Recruitment_mode,
  LDP = 0.1
)

Arguments

t	temporary numeric value, the current time step .
cr	temporary numeric value, the current control rule .
fdr	temporary numeric value, the current final target density ratio.
SSB	numeric array, the spawning stock biomass of the whole stock for each area, at each timestep, under each control rule, and for each estimate of natural mortality, in kg.
A	numeric value, the number of total areas in the model. Default value is 5.
R0	numeric value, set arbitrarily, the unfished recruitment. Default value is 1e+5.
H	numeric value, the steepness of the stock-recruitment curve.
B0	numeric value, the unfished biomass, in kg.
Eps	numeric matrix, the recruitment error terms.
Sigma_R	numeric value, the recruitment standard deviation.
Rec_age	numeric value, the age at recruitment, in years.
Recruitment_mode	character value, values can be: 'closed' - the recruits in each area originate from adults in that area. 'pool' - the recruits in each area come from a pool of larvae produced by adults in all areas. 'regional_DD' - larvae experience regional density dependence before settling evenly across all areas 'local_DD' - larvae experience local density dependence before settling evely across all areas Default value is 'pool'.
LDP	numeric value, the larval drift proportion, the proportion of larvae that drift from one area to an adjacent area before settling. Default value is 0.1.

Value

a numeric value representing the number of new recruits coming into the population in area a, at timestep t, under control rule cr.

Examples

A = 5; TimeT = 70; CR = 6; FDR = 4
SSB <- array(rep(10, A*TimeT*CR*FDR), c(A, TimeT, CR, FDR))
NuR <- array(rnorm(A*TimeT*CR*FDR, 0, 0.5), c(A, TimeT, CR, FDR))
Eps <- epsilon(A, TimeT, CR, FDR, NuR, Rho_R = 0)
recruitment(t = 3, cr = 1, fdr = 1, SSB, A, R0 = 1e+5, H = 0.65,
   B0 = 1e+5/1.1, Eps, Sigma_R = 0.5, Rec_age = 2, Recruitment_mode = 'pool',
   LDP = 0.1)

vquennessen/densityratio documentation built on Aug. 28, 2022, 5:36 p.m.