Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Try out changing R0 in the future_om_list #139

Open
k-doering-NOAA opened this issue Jul 21, 2022 · 1 comment
Open

Try out changing R0 in the future_om_list #139

k-doering-NOAA opened this issue Jul 21, 2022 · 1 comment
Labels
bug Something isn't working

Comments

@k-doering-NOAA
Copy link
Collaborator

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 k-doering-NOAA added the bug Something isn't working label Jul 21, 2022
@k-doering-NOAA
Copy link
Collaborator Author

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant