schaefer.sim: schaefer.sim()

In Henning-Winker/FLRef: Reference point computation for advice rules

schaefer.sim()

Description

generates a Schafer surplus production model with process and observation error

Usage

schaefer.sim(
  k = 10000,
  r = 0.3,
  q = 0.5,
  pe = 0.1,
  oe = 0.2,
  bk = 0.9,
  years = 1980:2022,
  f0 = 0.2,
  fhi = 2.2,
  flo = 0.8,
  sigmaF = 0.15,
  iters = 1,
  blim = 0.3,
  bthr = 0.5,
  rel = FALSE
)

Arguments

k	carrying capacity
r	intrinsic rate of population increase
q	catchability coefficient
pe	process error
oe	process error
bk	initial fraction of b/k
years	time horizon
f0	factor for initial year as f0 = f/fmsy
fhi	factor for high F as fhi = f/fmsy
flo	factor for low F as flo = fbar/fmsy
sigmaF	variation on f trajectory
iters	number of iterations
rel	if TRUE metrics B/Bmsy and F/Fmsy are produced

Value

FLQuants

Examples

stk = schaefer.sim(iters=100,q=0.5) 
plotAdvice(stk)
plot(FLIndex(index=iter(stk@stock,1))) # index

Henning-Winker/FLRef documentation built on Dec. 16, 2024, 7:23 p.m.