Commit 6a888af1 authored by Iago Mosqueira's avatar Iago Mosqueira
Browse files

Update AAP run, retros with rho

parent 7520c0bf
No preview for this file type
......@@ -60,14 +60,14 @@ retro <- FLStocks(parallel::mclapply(years, function(x) {
control=control, stdfile=btsgam@stdfile, wkdir=tempfile()) + window(stock, end=x)
}, mc.cores=length(years)))
# --- RUNS & STOCKS
# --- STOCKS & RESIDUALS
stocks <- FLStocks(lapply(runs, "+", stock))
residuals <- lapply(runs, residuals, stock)
residuals <- lapply(runs[!names(runs) %in% "McMC"], residuals, stock)
# --- SRR
srrs <- FLSRs(lapply(c("bevholt", "ricker", "segreg"), function(x)
fmle(as.FLSR(stocks[["BTS-GAM"]], model=x))))
srrs <- FLSRs(setNames(lapply(c("bevholt", "ricker", "segreg"), function(x)
fmle(as.FLSR(stocks[["BTS-GAM"]], model=x, name=x))), c("bevholt", "ricker", "segreg")))
# --- REFERENCE points
......
No preview for this file type
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment