computeFp05 | R Documentation |
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
)
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. |
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))
# 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.