########################################################################
# #
# produce supp. Fig 5 and 6 - sensitivity analyses 10p #
# #
########################################################################
rm(list = ls(), envir = .GlobalEnv)
### ----- LOAD DATA
pos <- get(load(here::here("results", "sensitivity_analysis_res_pos0.1.RData")))
neg <- get(load(here::here("results", "sensitivity_analysis_res_neg0.1.RData")))
### ----- FORMAT DATA
## ---- Positive values
sensit_ana_res_pos <- pos[["res_analyses"]]
null_values_pos <- pos[["null_values"]]
## ---- Negative values
sensit_ana_res_neg <- neg[["res_analyses"]]
null_values_neg <- neg[["null_values"]]
### ----- PLOT DATA
## ---- Positive values
NCSSDGproj::supp_fig_null_hist(null_vals = null_values_pos,
res_null_mod = sensit_ana_res_pos,
name = "Supp_fig5_sensit_pos_10p")
## ---- Negative values
NCSSDGproj::supp_fig_null_hist(null_vals = null_values_neg,
res_null_mod = sensit_ana_res_neg,
name = "Supp_fig6_sensit_neg_10p")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.