Add missing Canary assessment

[seananderson]

https://waves-vagues.dfo-mpo.gc.ca/Library/338602.pdf

[seananderson]

Digitize this:

And find B0 in doc.

0.2 B0 was equivalent of LRP... update shortly after as 0.4 Bmsy, but values not shown

From: https://waves-vagues.dfo-mpo.gc.ca/Library/338492.pdf

[seananderson]

RAM legacy code... but I think we're best working with the docs above so we know exactly what we're working with, plus we want uncertainty.

library(ggplot2)
library(dplyr)
library(ramlegacy)
download_ramlegacy()
d <- load_ramlegacy()
names(d)
d$stock[grep("Canary", d$stock$commonname),]

x <- d$timeseries %>% dplyr::filter(stockid == "CROCKWCVANISOGQCI")
unique(x$tsid)

names(d)

d$timeseries_ids_views %>% head
d$timeseries_ids_views %>% filter(stockid == "CROCKWCVANISOGQCI")

filter(x, tsid == "BdivBmsytouse-dimensionless") %>% 
  ggplot(aes(tsyear, tsvalue)) + geom_line()

filter(x, tsid == "BdivBmsypref-dimensionless") %>% 
  ggplot(aes(tsyear, tsvalue)) + geom_line()

[seananderson]

https://cran.r-project.org/web/packages/digitize/index.html

[seananderson]

The initial plot in this thread is for run 11:

Res Doc p189 Fig. J26.