Skip to contents

Computes performance statistics of one or more projections (named list of ms_simulate() runs, one per management procedure or rule) and returns them as a long data.table in the format of the mseviz package: one row per run (mp), unit (biol: a stock code, a segment, or "All" for the whole system), statistic and iteration, with the value in data. The table can be passed to plot_bps(), plot_tradeoffs() and plot_quilt(), or to mseviz::plotBPs(), mseviz::plotTOs() and mseviz::perfTable().

Usage

ms_perf(runs, years = NULL, Blim = 0.25, Bpa = 0.5)

Arguments

runs

named list of ms_simulate() results (with keep = "all"), one per management procedure; all from the same system.

years

evaluation years (default the last 20).

Blim, Bpa

limit and precautionary biomass relative to Bmsy.

Value

A data.table with columns mp, biol, statistic, name

(plotmath label), desc, type ("value" or "prob"), year

(the evaluation period), iter and data.

Details

Statistics are computed per iteration over the evaluation years:

statisticunitdefinition
BBmsystockmedian B/Bmsy
FFmsystockmedian F/Fmsy
CMSYstockmean catch / MSY (MSY = rK/4)
PBlimstockshare of years with B < Blim
PBpastockshare of years with B < Bpa
PB80stockshare of years with B > 0.8 Bmsy
greenstockshare of years in the Kobe green zone (B > Bmsy, F < Fmsy)
Dayssegment, Allmean days fished
CAllmean total catch
RevAllmean revenue
ProfitAllmean profit
CMMSYAllmean total catch / sum of single-stock MSYs
AAVCAllmean annual absolute change of total catch (relative)
AAVEAllmean annual absolute change of total days (relative)
NBpaAllmean number of stocks below Bpa
PanyBlimAllshare of years with any stock below Blim
minBBmsyAlllowest B/Bmsy of any stock in any year

Probabilities (statistics starting with P and green) are best shown by their mean over iterations, which is the probability over years and iterations.