SP_gf: PBS groundfish surplus production wrapper function

In pbs-assess/ggmse: Tools for Working with openMSE (Including DLMtool, MSEtool, and SAMtool)

PBS groundfish surplus production wrapper function

Description

PBS groundfish surplus production wrapper function

Usage

SP_gf(
  x,
  Data,
  reps = 1,
  LRP = 0.4,
  TRP = 0.6,
  RP_type = "SSB_SSBMSY",
  start = list(r_prior = c(0.3, 0.1)),
  use_r_prior = TRUE,
  tac_max_increase = 1.2,
  tac_max_decrease = 0.5,
  tac_floor = 0.1,
  tac_increase_buffer = 1.05,
  initial_tac = NULL
)

SP6040_gf(x, Data, reps = 1, ...)

SP8040_gf(x, Data, reps = 1, ...)

SP4010_gf(x, Data, reps = 1, ...)

add_SP_prior(
  mp,
  r_prior,
  tac_max_increase = 1.2,
  other_start = NULL,
  tac_max_decrease = 0.5,
  tac_floor = 0.1,
  tac_increase_buffer = 1.05,
  initial_tac = NULL,
  ...
)

Arguments

x	Index
Data	Data
reps	Reps
LRP	Lower reference point
TRP	Target reference point
RP_type	Reference point type
start	Starting list for the model
use_r_prior	Logical but whether to use a prior
tac_max_increase	Maximum proportional increase in TAC
tac_max_decrease	Maximum proportional decrease in TAC
tac_floor	Floor for TAC
tac_increase_buffer	Proportional buffer below which TAC won't change
initial_tac	An initial-year TAC to compare the HCR "meta" rules to instead of using the last year's catch.
...	Other args
mp	MP to wrap
r_prior	Mean and SD of r prior
other_start	A named list of other elements to pass to SAMtool::SP().

Examples

library(SAMtool)
my_mp <- add_SP_prior(SP4010_gf, c(0.3, 0.05))
om <- MSEtool::testOM
om@nsim <- 5
om@proyears <- 10
## Not run: 
mse <- runMSE(om, MPs = "my_mp")

## End(Not run)

pbs-assess/ggmse documentation built on Nov. 21, 2023, 8:06 p.m.