make_results/supp_figure_5_6.R

########################################################################
#                                                                      #
#         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")
GaelMariani/NCSSDGproj documentation built on April 1, 2024, 9:18 a.m.