brpsson: Reference points for a seasonal model

In flr/FLBEIA: Bio-Economic Impact Assessment of Management Strategies using FLR

Reference points for a seasonal model

Description

Function to estimate reference points when having a seasonal model

Usage

brpsson(
  stk,
  B0,
  R0,
  rec.ss = 1,
  ssb.ss = 1,
  sr_model,
  sr_params,
  Fprop = rep(1/dim(stk)[4], dim(stk)[4]),
  Fscan = seq(0, 4, by = 0.01),
  oldest = 100,
  nrun = 200,
  tol = 0.01
)

Arguments

stk	An FLStock object.
B0	The value of the virgin biomass.
R0	The expected recrutiment in the virgin population.
rec.ss	The recruitment season (numeric). Default value = 1.
ssb.ss	The spawning season (numeric). Default value = 1.
sr_model	A character with the name of the model to simulate the recruitment process.
sr_params	A named vector with the SR parameter values.
Fprop	A vector with the same length as the number of seasons with the proportion of F by season. By default the same proportion in all the seasons is assumed.
Fscan	A vector with the F values to be simulated. By default, will scan values between 0 and 4 (by 0.01 increments).
oldest	The maximum age class to be considered. Default 100.
nrun	The maximum number of years to project forward until equilibrium.
tol	The desired accuracy.

See Also

plotBRPsson

Examples

library(FLBEIA)
data(multistk) # object with 2 seasons and 3 iterations

stk <- trim( multistk, year=1, iter = 1) 

# loop for different catch proportions by season
for (p in seq(0.1,0.9,0.1)) {
  fruns <- brpsson( stk, B0=1e+05, R0=27489766, rec.ss=2, ssb.ss=2, 
                    sr_model="bevholt", sr_params=c( a = 29988835.109, b = 9090.909), 
                    Fprop = c(p, 1-p))
} 

plotBRPsson( fruns, pdfnm="stk_Fbar_vs_SPR.pdf")
 

flr/FLBEIA documentation built on Feb. 22, 2024, 1:36 a.m.