Fmmy: Fmmy() Uses opt.bisect to derive the F at Maximum Median...

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

Fmmy() Uses opt.bisect to derive the F at Maximum Median Yield from stochastic simulations

Description

Fmmy() Uses opt.bisect to derive the F at Maximum Median Yield from stochastic simulations

Usage

Fmmy(
  brp,
  sigmaR = 0.5,
  rho = 0,
  nyears = 100,
  iters = 250,
  yrs.eval = NULL,
  range = "missing",
  tol = 0.001,
  maxit = 15,
  verbose = TRUE
)

Arguments

brp	output object from computeFbrp() of class FLBRP
sigmaR	lognormal recruitment standard deviation
rho	AR1 recruitment autocorrelation coefficient
nyears	number of simulation years
iters	number simulation iterations
yrs.eval	last years to be used evaluation period, default nyears/2
range	range of Fbar value to be evaluated
tol	tolerance
maxit	number of steps
verbose	cat comments

Value

list of FLPar, FLStock and FLBRP objects

Examples

data(ple4)
bh = srrTMB(as.FLSR(ple4,model=bevholtSV),spr0=spr0y(ple4))
brp = computeFbrp(ple4,bh,proxy=c("bx","msy"),x=35,blim=0.1)
fmmy = Fmmy(brp,sigmaR=0.7,rho=0.3)
getF(fmmy) # FMMY value 
plotFsim(fmmy)
brpfmmy = computeFbrp(ple4,bh,proxy=getF(fmmy),blim=0.1)
fsim = Fsim(brpfmmy,sigmaR=0.7,rho=0.3)
plotFsim(fsim)

Henning-Winker/FLRef documentation built on July 14, 2024, 7:28 a.m.