sprat-xsa-experimental

1
library(ggplot2)
2

3
# build fishing mortality uncertainty data
4
pics <- FLQuants(
5
  Rec = rnorm(1000, rec(result.xsa.stock), result.xsa.sigma * rec(result.xsa.stock)),
6
  SSB = rnorm(1000, ssb(result.xsa.stock), result.xsa.sigma * ssb(result.xsa.stock)),
7
  F = rnorm(1000, fbar(result.xsa.stock), result.xsa.sigma * fbar(result.xsa.stock)),
8
  Catch = catch(result.xsa.stock)
9
)
10

11
# get ssb
12
xsa.ssb <- pics$SSB
13
# read cmsy output from table
14
cmsy.b <- read_clip_tbl()
15

16
plot(xsa.ssb) + 
17
  geom_line(data = cmsy.b, aes(y = b), linetype = "dashed") + 
18
  geom_ribbon(data = cmsy.b, aes(y = b, ymin = b_lo, ymax = b_hi), alpha = 0.15, fill = "blue") + 
19
  theme_minimal()