Regime Beverton-Holt stock recruitment relationship — bevholt

Defines a two-regime recruitment function where each regime is defined by a Beverton-Holt stock recruitment relationship with differeing levels of R0 or sigR.

Usage

bevholt_regimes(
  h,
  sbpr,
  R0s,
  sigRs,
  nyears,
  regime_length,
  starting_regime,
  seed
)

Arguments

h: steepness (0.2 <= h <= 1.0)
sbpr: spawning biomass per recruit at unfished conditions
R0s: vector of unfished recruitments (one for each regime, max 2)
sigRs: vector of recruitment variability (one for each regime, max 2)
nyears: total number of years of recruitments to generate
regime_length: vector of length of each regime
starting_regime: regime to start with (0 = first regime, 1 = second regime)
seed: random seed for reproducability