plot_scenario: Plot a single time series for various scenarios and MPs

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

Plot a single time series for various scenarios and MPs

Description

Plot historical and projected time series for various scenarios (B, B/BMSY, catch, etc.) either (1) separate panels by MP that overlap by OM scenario, or (2) as a full grid by OM scenario and MP. See also plot_scenario_index().

Usage

plot_scenario_BMSY(
  object_list,
  MP = NULL,
  quantiles = c(0.025, 0.975),
  facet = c("wrap", "grid"),
  palette = "Dark2",
  n_samples = 0,
  seed = 42,
  french = isTRUE(getOption("french"))
)

plot_scenario_SSB(
  object_list,
  MP = NULL,
  quantiles = c(0.025, 0.975),
  facet = c("wrap", "grid"),
  palette = "Dark2",
  n_samples = 2,
  seed = 42,
  french = isTRUE(getOption("french"))
)

plot_scenario_F(
  object_list,
  MP = NULL,
  quantiles = c(0.025, 0.975),
  facet = c("wrap", "grid"),
  palette = "Dark2",
  n_samples = 2,
  seed = 42,
  french = isTRUE(getOption("french"))
)

plot_scenario_FMSY(
  object_list,
  MP = NULL,
  quantiles = c(0.025, 0.975),
  facet = c("wrap", "grid"),
  palette = "Dark2",
  n_samples = 2,
  seed = 42,
  french = isTRUE(getOption("french"))
)

plot_scenario_catch(
  object_list,
  MP = NULL,
  quantiles = c(0.025, 0.975),
  facet = c("wrap", "grid"),
  palette = "Dark2",
  n_samples = 2,
  seed = 42,
  french = isTRUE(getOption("french"))
)

Arguments

object_list	A named list of MSE objects from MSEtool. Names become scenario names.
MP	Optional subset of management procedures
quantiles	Length 2 vector for the quantiles to plot as ribbons
facet	Character to indicate a wrapped facet plot by MP or grid by OM and MP. See description.
palette	The name of an RColorBrewer colour palette
n_samples	The number of timeseries samples to illustrate. Only used if facet = "grid".
seed	The seed to set before drawing samples
french	Logical to use French translation.

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