hcr,biodyn,FLPar-method | R Documentation |
Harvest Control Rule, calculates F, or Total Allowable Catch (TAC) based on a hockey stock harvest control rule.
Harvest Control Rule, calculates F, or Total Allowable Catch (TAC) based on a hockey stock harvest control rule.
## S4 method for signature 'biodyn,FLPar'
hcr(
object,
refs = hcrParam(ftar = 0.7 * mpb:::fmsy(refs), btrig = 0.8 * mpb:::bmsy(refs), fmin =
0.01 * mpb:::fmsy(refs), blim = 0.4 * mpb:::bmsy(refs)),
params = refs,
stkYrs = max(as.numeric(dimnames(stock(object))$year)),
refYrs = max(as.numeric(dimnames(catch(object))$year)),
hcrYrs = max(as.numeric(dimnames(stock(object))$year)),
tac = TRUE,
bndF = NULL,
bndTac = NULL,
maxF = 2,
...
)
## S4 method for signature 'FLStock,FLBRP'
hcr(
object,
refs,
params = hcrParam(ftar = 0.7 * fmsy(refs), btrig = 0.8 * bmsy(refs), fmin = 0.01 *
fmsy(refs), blim = 0.4 * bmsy(refs)),
stkYrs = max(as.numeric(dimnames(stock(object))$year)),
refYrs = max(as.numeric(dimnames(catch(object))$year)),
hcrYrs = max(as.numeric(dimnames(stock(object))$year)),
tac = TRUE,
tacMn = TRUE,
bndF = NULL,
bndTac = NULL,
maxF = 2,
...
)
object |
an object of class |
... |
other parameters, i.e.
refs |
FLPar
object with value(s) for F or TAC if tac==TRUE
FLPar
object with value(s) for F or TAC if tac==TRUE
## Not run:
bd =sim()
bd=window(bd,end=29)
for (i in seq(29,49,1))
bd=fwd(bd,harvest=hcr(bd,yr=i,yr=i+1)$hvt)
## End(Not run)
## Not run:
bd =sim()
bd=window(bd,end=29)
for (i in seq(29,49,1))
bd=fwd(bd,harvest=hcr(bd,yr=i,yr=i+1)$hvt)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.