In pacific-hake/hake-assessment: Stock Assessment for Pacific Hake

Model estimates

Recruitment (numbers of fish)

plot_recruitment(list(base_model),
                 list(base_model_name),
                 inc_means = TRUE,
                 leg_pos = "none") +
  theme(axis.title.x = element_text(size = 12),
        axis.title.y = element_text(size = 12),
        legend.key.size = unit(0.5, units = "cm"),
        legend.text = element_text(size = 10))

Spawning Biomass

plot_biomass(list(base_model),
             list(base_model_name),
             ylim = c(0, 5),
             point_shape = ts_single_model_pointshape,
             leg_pos = "none") +
  theme(axis.title.x = element_text(size = 12),
        axis.title.y = element_text(size = 12),
        legend.key.size = unit(0.5, units = "cm"),
        legend.text = element_text(size = 10))

Relative Spawning Biomass

plot_rel_biomass(model_lst = list(base_model),
                 model_names = list(base_model_name),
                 ylim = c(0, 2.5),
                 alpha = 0.2,
                 leg_pos = "none") +
  theme(axis.title.x = element_text(size = 12),
        axis.title.y = element_text(size = 12),
        legend.key.size = unit(0.5, units = "cm"),
        legend.text = element_text(size = 10))

Fishing intensity and exploitation rate

:::::: {.columns} ::: {.column width="50%"}

plot_fishing_intensity(base_model,
                       line_width = 0.5,
                       point_size = 1.5,
                       point_stroke = 1,
                       fspr40_text_x = 1970) +
  theme(axis.title.x = element_text(size = 10),
        axis.title.y = element_text(size = 10),
        axis.text.x = element_text(size = 8),
        legend.key.size = unit(0.5, units = "cm"),
        legend.text = element_text(size = 10))

:::

::: {.column width="50%"}

plot_exploitation_fraction(base_model,
                           line_width = 0.5,
                           point_size = 1.5,
                           point_stroke = 1) +
  theme(axis.title.x = element_text(size = 10),
        axis.title.y = element_text(size = 10),
        axis.text.x = element_text(size = 8),
        legend.key.size = unit(0.5, units = "cm"),
        legend.text = element_text(size = 10))