A HCR to set total catch based on SSB depletion level

A HCR to set total catch based on SSB depletion level

catchSSB.hcr(
  stk,
  dtarget = 0.4,
  dlimit = 0.1,
  lambda = 1,
  MSY,
  dtaclow = 0.85,
  dtacupp = 1.15,
  yrs = 1,
  metric = "ssb",
  args,
  tracking
)

Arguments

stk

The perceived FLStock.

MSY

Assumed or estimated MSY.

args

MSE arguments, class list.

tracking

Structure for tracking modules outputs.

dtarget=0.40

Depletion level from which catch is decreased.

dlimit=0.10

Depletion level at which fishing is stopped.

lambda=1

Multiplier for MSY level.

dtaclow=0.85

Maximum proportional decrease in allowable catch.

dtacupp=1.15

Maximum proportional increase in allowable catch.

Value

A list with elements ctrl, of class fwdControl, and tracking.

Examples

data(sol274)
#> Warning: namespace ‘DBI’ is not available and has been replaced
#> by .GlobalEnv when processing object ‘om’
catchSSB.hcr(stock(om), MSY=140000, tracking=FLQuant(),
args=list(ay=2018, data_lag=1, management_lag=1, frq=1))
#> $ctrl
#> An object of class "fwdControl"
#>  (step) year quant relYear   min             value   max
#>       1 2019 catch      NA    NA 140000.000(0.000)    NA
#>       2 2019 catch    2018 0.850                NA 1.150
#>    iters:  100 
#> 
#> 
#> $tracking
#> An object of class "FLQuant"
#> , , unit = unique, season = all, area = unique
#> 
#>      year
#> quant 1 
#>   all NA
#> 
#> units:  NA 
#> 
# APPLY hcr over a range of dtarget values
lapply(seq(0.30, 0.80, by=0.1), function(x) {
  catchSSB.hcr(stock(om), MSY=140000, dtarget=x,
  args=list(ay=2018, data_lag=1, management_lag=1, frq=1),
  tracking=FLQuant())$ctrl } )
#> [[1]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min      value   max
#>       1 2019 catch      NA    NA 140000.000    NA
#>       2 2019 catch    2018 0.850         NA 1.150
#>    iters:  100 
#> 
#> 
#> [[2]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min             value   max
#>       1 2019 catch      NA    NA 140000.000(0.000)    NA
#>       2 2019 catch    2018 0.850                NA 1.150
#>    iters:  100 
#> 
#> 
#> [[3]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min                 value   max
#>       1 2019 catch      NA    NA 121577.715(20663.032)    NA
#>       2 2019 catch    2018 0.850                    NA 1.150
#>    iters:  100 
#> 
#> 
#> [[4]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min                value   max
#>       1 2019 catch      NA    NA 97262.172(16530.426)    NA
#>       2 2019 catch    2018 0.850                   NA 1.150
#>    iters:  100 
#> 
#> 
#> [[5]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min                value   max
#>       1 2019 catch      NA    NA 81051.810(13775.355)    NA
#>       2 2019 catch    2018 0.850                   NA 1.150
#>    iters:  100 
#> 
#> 
#> [[6]]
#> An object of class "fwdControl"
#>  (step) year quant relYear   min                value   max
#>       1 2019 catch      NA    NA 69472.980(11807.447)    NA
#>       2 2019 catch    2018 0.850                   NA 1.150
#>    iters:  100 
#> 
#>