SSMSE icon indicating copy to clipboard operation
SSMSE copied to clipboard

Try out changing R0 in the future_om_list

Open k-doering-NOAA opened this issue 2 years ago • 1 comments

Thanks @rwildermuth for trying out some different specifications for changing recruitment regimes. Trying to make R0 time varying did not work, and I think this is a bug.

Works:

template_mod_change <- create_future_om_list(example_type = "model_change")
rec_dev_specify <- template_mod_change[[1]]
rec_dev_specify$pars <- "rec_devs"
rec_dev_specify$scen <- c("replicate", "all")
rec_dev_specify$input <- data.frame(first_yr_averaging = NA, #rep(datfile$styr, 1),
                                    last_yr_averaging = NA, #rep(2019, 1),
                                    last_yr_orig_val = c(2019, 2024),
                                    first_yr_final_val = c(2020, 2025),
                                    ts_param = c("sd", "mean"),
                                    method = c("absolute", "additive"),
                                    value = c(1.25, -1))

Does not work:

template_mod_change <- create_future_om_list(example_type = "model_change")
rec_dev_specify <- template_mod_change[[1]]
rec_dev_specify$pars <- c("SR_LN(R0)")#, "SR_regime_BLK1repl_2000")
rec_dev_specify$scen <- c("replicate", "all") # note: could change this to c("random", "all") if did not want to replicate the same recdevs sequences across scenarios
rec_dev_specify$input <- data.frame(first_yr_averaging = NA, #rep(datfile$styr, 1),
                                    last_yr_averaging = NA, #rep(2019, 1),
                                    last_yr_orig_val = rep(2024, 1),
                                    first_yr_final_val = rep(2025, 1),
                                    ts_param = rep("mean", 1),
                                    method = rep("absolute", 1),
                                    value = c(15.3516))#, -2.25458))

k-doering-NOAA avatar Jul 21 '22 20:07 k-doering-NOAA