DTUAqua/spict

Plotting management scenarios when fractiles are used

Closed this issue · 1 comments

The plotting function plotspict.catch ignores catch fractile while plotting.

Should check if there are issues with other plotting functions and fractiles.

library(spict)

inp <- pol$lobster
inp$dteuler <- 1/4

inp$maninterval <- c(1992, 1993)
inp$maneval <- 1993

fit <- fit.spict(inp)

fit <- add.man.scenario(fit, "A", breakpointB = 0.5)
fit <- add.man.scenario(fit, "B", fractiles = list(catch = 0.25) , breakpointB = 0.5)

plotspict.catch(fit, stamp = "")

sumspict.manage(fit)
## C B/Bmsy F/Fmsy
## 1. A 1502.3   0.52   0.89
## 2. B 1258.0   0.52   0.87

Possible solution: retape with ffac corresponding to given catch percentile