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

View source: R/SimRefs.R

FmmyR Documentation

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.