Calculates the Fbar value giving a maximum probability of ssb being below Blim of 5%

computeFp05(
  stock,
  sr,
  SBlim,
  range = c(0.01, 0.75),
  nyears = 3,
  sigmaR = 0.5,
  rho = 0.43,
  its = 500,
  verbose = TRUE
)

Arguments

stock: An FLStock over which the calculation is carried out.
sr: The stock-recruits relationship to use in fwd.
its: Number of iterations
verbose: Should progress be shown, TRUE.

Examples

data(ple4)
sr <- predictModel(model=bevholt, params=FLPar(a=1.4e6, b=1.5e5))
fp05 <- computeFp05(ple4, sr, SBlim=150000, its=300, range=c(0.40, 0.50))
#> [1] fbar: 0.4 - prob: 6e-04  - diff:  -0.0494 
#> [2] fbar: 0.5 - prob: 0.6111333  - diff:  0.5611333 
#> [3] fbar: 0.45 - prob: 0.1066667  - diff:  0.05666667 
#> [4] fbar: 0.425 - prob: 0.01606667  - diff:  -0.03393333 
#> [5] fbar: 0.4375 - prob: 0.04633333  - diff:  -0.003666667 
# RUN projection for obtained Fp.05 value
proj <- fwd(propagate(stf(ple4, nyears=100), 300), sr=sr,
  fbar=FLQuant(fp05, dimnames=list(year=2018:2117)),
  deviances=ar1rlnorm(rho=0.43, years=2018:2117, iters=300, meanlog=0,
  sdlog=0.5))
plot(ssb(proj), prob=c(0.01, 0.25, 0.50, 0.75, 0.99)) +
  geom_hline(yintercept=150000)