retrospective_AM: retrospective_AM (retrospective of Assessment model in MSE)

Description

Plots the true retrospective of an assessment model during the MSE. A series of time series estimates of SSB, F, and VB are plotted over the course of the MSE are plotted against the operating model (true) values (in black).

Usage

retrospective_AM(MSE, sim = 1, MP, MSE_Hist = NULL, plot_legend = FALSE)

Arguments

MSE

An object of class MSE created by runMSE with PPD = TRUE.

sim

Integer between 1 and MSE@nsim. The simulation number for which the retrospectives will be plotted.

Character. The name of the management procedure created by make_MP containing the asssessment model.

MSE_Hist

Optional. The list containing historical data for the MSE, created by runMSE with argument Hist = TRUE. Currently only used to plot operating model vulnerable biomass in historical period.

plot_legend

Logical. Whether to plot legend to reference year of assessment in the MSE.

Value

A series of figures for spawning stock biomass (SSB, including absolute magnitude and relative to MSY and virgin), fishing mortality (F, including absolute magnitude and relative to MSY), and vulnerable biomass (VB) estimates over the course of the MSE are plotted against the operating model (true) values (in black).

Details

For assessment models that utilize annual harvest rates (u), the instantaneous fishing mortality rates are obtained as F = -log(1 - u).

Examples

Run this code

# NOT RUN {
DD_MSY <- makeMP(DD_TMB, HCR_MSY, diagnostic = "full")
myMSE_hist <- DLMtool::runMSE(DLMtool::testOM, Hist = TRUE)
myMSE <- DLMtool::runMSE(DLMtool::testOM, MPs = "DD_MSY", PPD = TRUE)
retrospective_AM(myMSE, sim = 1, MP = "DD_MSY")
retrospective_AM(myMSE, sim = 1, MP = "DD_MSY", Hist = myMSE_hist)
# }