make_results/supp_figure_3.R

In GaelMariani/NCSSDGproj: What the Package Does (One Line, Title Case)

####################################################################
#                                                                  #
# produce supp. Fig 3 - Correspondance Analysis for negative links #
#                                                                  #
####################################################################
rm(list = ls(), envir = .GlobalEnv)


### ----- LOAD DATA

  ## ---- Data of links between NCS and SDG
  sheets  <- NCSSDGproj::read_all_sheets()
  

### ----- FORMAT DATA

  ## ---- From sheets to df
  matrix_all <- NCSSDGproj::sheets_to_df(sheets_list = sheets, binary = TRUE)
  
  ## ---- From dataframes to contingency matrices 
  
    # --- binary data
    matrix_conting_bin <- lapply(matrix_all, NCSSDGproj::contingency_mat_targets, binary = TRUE)
    
    # --- 0-1-2 score
    matrix_conting <- lapply(matrix_all, NCSSDGproj::contingency_mat_targets, binary = FALSE)
  
  ## ---- Informations on the NCS
  info_NCS <- NCSSDGproj::NCS_info(matrix_cont = matrix_conting_bin[[1]])
  

### ----- ANALYSIS  

  ## ---- Correspondance analysis for POSITIVE links with the contribution of each target and NCS to the variance of each axis
  contri_CA_neg <- NCSSDGproj::CA_contri_vars(matrix_cont  = matrix_conting_bin[["score_neg"]],
                                              colNCS_ter   = "#228B22",
                                              colNCS_coast = "#5EA9A2",
                                              colNCS_mar   = "#1134A6")

  
### PLOT DATA 

  ## ---- Figure with all panels for POSITIVE links
  NCSSDGproj::supp_fig2_3(data           = contri_CA_neg[["CorresAna"]],
                          arrow          = FALSE,
                          data_arrow     = contri_CA_neg[["data_arrow"]],
                          colNCS_ter     = "#228B22", 
                          colNCS_coast   = "#5EA9A2",
                          colNCS_mar     = "#1134A6",
                          save           = TRUE,
                          name           = "Supp_fig3")