- ajout des indicateurs solagro dans les premiers traitements des réponses (au...

- ajout des indicateurs solagro dans les premiers traitements des réponses (au stade des redressement)
- ajout du graph de comparaison des gisements SRB/Solagro dans intrants valo

2_redressement_reponses.Rmd

@@ -30,7 +30,7 @@ load(paste0("1-exports_", params$campagne, ".RData"))
 load(paste0("0-donnees_ref", params$campagne, ".RData"))
 # menage dans l'environnement de travail : on supprime les palettes, les éléments pour les visualisation des intrants et les fonctions de formatage des tableaux
 rm(list = ls(pattern = "^pal_|^vec_col|_tot$|^palette_"))
-rm(famille_ordre, gismts_solagro, filigrane)
+rm(famille_ordre, filigrane)
 ```
 
 # Campagnes précédentes
@@ -369,6 +369,50 @@ indic_suivi_SRB
 
 ```
 
+### Calculs agrégats solagro
+
+```{r solagro}
+
+gsmt_sol_dep_mil <- expand_grid(Dep = ter, campagne  = mil_disp) %>% 
+  cross_join(gismts_solagro) 
+
+indic_suivi_solagro_base <- typo_dechets %>% 
+  select(cat_intrant = Categorie_intrants, Famille_intrants, type_gismnt_solagro) %>% 
+  distinct() %>% 
+  full_join(gsmt_sol_dep_mil, by = "type_gismnt_solagro", relationship = "many-to-many") %>% 
+  full_join(compil_pluriannuelle_intrants, .,  by = c("cat_intrant", "Famille_intrants", "Dep", "campagne")) %>% 
+  select(campagne, Dep, contains("solagro"), names(gismts_solagro), tonnage_total, gismt_reg = tMB_gismt_solagro) %>% 
+  group_by(campagne, ordre_solagro, type_gismnt_solagro, obj_srb_solagro, Dep, gismt_reg, gsmt_44, gsmt_49, gsmt_53, gsmt_72, gsmt_85) %>% 
+  summarise(tonnages_declares = sum(tonnage_total, na.rm = TRUE), .groups = "drop") %>% 
+  pivot_longer(cols = c(gsmt_44, gsmt_49, gsmt_53, gsmt_72, gsmt_85, gismt_reg), 
+               names_to = c(NA, "zone"), values_to = "gsmt_tmb", names_sep = "_") %>%
+  filter(Dep == zone | zone == "reg") %>% 
+  select(-Dep) %>% 
+  group_by(campagne, ordre_solagro, type_gismnt_solagro, obj_srb_solagro, zone, gsmt_tmb) %>% 
+  summarise(tonnages_declares = sum(tonnages_declares, na.rm = TRUE), .groups = "drop") %>% 
+  bind_rows(
+    # on écarte les cultures principales pour calculer le total
+    filter(.,  type_gismnt_solagro != "Hors gisements étude TEO Solagro") %>% 
+      group_by(campagne, zone) %>% 
+      summarise(ordre_solagro = 8.5, type_gismnt_solagro = "Total périmètre étude Solagro TEO", .groups = "drop", 
+                across(.cols = c(obj_srb_solagro, tonnages_declares, gsmt_tmb), .fns = ~sum(.x, na.rm = TRUE)))
+  ) %>% 
+  arrange(ordre_solagro) %>% 
+  mutate(type_gismnt_solagro = fct_inorder(type_gismnt_solagro), 
+         tx_mob_gsmt = round((tonnages_declares / gsmt_tmb * 100), 3))
+
+indic_suivi_solagro_xls_0 <- indic_suivi_solagro_base %>% 
+  pivot_wider(id_cols = c(zone, type_gismnt_solagro, obj_srb_solagro, gsmt_tmb), names_from = c(campagne), 
+              values_from = c(tonnages_declares, tx_mob_gsmt)) %>% 
+  arrange(zone, type_gismnt_solagro)
+ 
+indic_suivi_solagro_xls <-  bind_rows(
+  filter(indic_suivi_solagro_xls_0, zone == "reg"),
+  filter(indic_suivi_solagro_xls_0, zone != "reg") %>% mutate(obj_srb_solagro = NA)
+)
+rm(indic_suivi_solagro_xls_0)
+indic_suivi_solagro_xls
+```
 
 
 ## Tableau de suivi fin par catégorie d'intrants
@@ -441,13 +485,15 @@ writexl::write_xlsx(x = list(compil_reponses_ss_RGPD_redresse = compil_reponses_
                              reponse_complements_intrant = compil_reponses5,
                              comp_projet_annees_prcdtes = aide_redressement,
                              Suivi_intrants =  suivi_fin,
-                             Suivi_SRB = indic_suivi_SRB), 
+                             Suivi_SRB = indic_suivi_SRB,
+                             Suivi_SOLAGRO = indic_suivi_solagro_xls), 
                     col_names = TRUE,
                     path = paste0(getwd(), "/compilation/", lubridate::today(), "_indicateurs_intrants_", params$campagne, ".xlsx"))
                 
     
 save(compil_reponses5, compil_reponses_ssRGPD_redr, compil_pluriannuelle_intrants, mil_disp,
-     indic_suivi_SRB, suivi_fin, file = paste0("2-intrants_pour_valo_", params$campagne,".RData"))
+     indic_suivi_SRB, indic_suivi_solagro_base, suivi_fin, 
+     file = paste0("2-intrants_pour_valo_", params$campagne,".RData"))
 
 
 ```

4_intrants_valo.Rmd

@@ -321,41 +321,6 @@ indic_suivi_SRB %>%
 Taux de mobilisation du gisement SOLAGRO : tableaux
 
 ```{r indic group solagro}
-code_reg <- "52"
-ter <- list_dep_in_reg(code_reg) %>%
-  c(., "reg") 
-noms_ter <- nom_zone(type_zone = "Départements", code_zone = list_dep_in_reg(code_reg)) %>%
-  c(., nom_zone(type_zone = "Régions", code_zone = code_reg)) %>% 
-  tibble(ter, nom_ter = .)
-  
-gsmt_sol_dep_mil <- expand_grid(Dep = ter, campagne  = mil_disp) %>% 
-  cross_join(gismts_solagro)
-
-indic_suivi_solagro_base <- typo_dechets %>% 
-  select(cat_intrant = Categorie_intrants, Famille_intrants, type_gismnt_solagro) %>% 
-  distinct() %>% 
-  full_join(gsmt_sol_dep_mil, by = "type_gismnt_solagro") %>% 
-  full_join(compil_pluriannuelle_intrants, .,  by = c("cat_intrant", "Famille_intrants", "Dep", "campagne")) %>% 
-  arrange(ordre_solagro) %>% 
-  mutate(type_gismnt_solagro = fct_inorder(type_gismnt_solagro)) %>% 
-  select(campagne, Dep, contains("solagro"), names(gismts_solagro), tonnage_total, gismt_reg = tMB_gismt_solagro) %>% 
-  group_by(campagne, type_gismnt_solagro, obj_srb_solagro, Dep, gismt_reg, gsmt_44, gsmt_49, gsmt_53, gsmt_72, gsmt_85) %>% 
-  summarise(tonnages_declares = sum(tonnage_total, na.rm = TRUE), .groups = "drop") %>% 
-  pivot_longer(cols = c(gsmt_44, gsmt_49, gsmt_53, gsmt_72, gsmt_85, gismt_reg), 
-               names_to = c(NA, "zone"), values_to = "gsmt_tmb", names_sep = "_") %>%
-  filter(Dep == zone | zone == "reg") %>% 
-  select(-Dep) %>% 
-  group_by(campagne, type_gismnt_solagro, obj_srb_solagro, zone, gsmt_tmb) %>% 
-  summarise(tonnages_declares = sum(tonnages_declares, na.rm = TRUE), .groups = "drop") %>% 
-  bind_rows(
-    # on écarte les cultures principales pour calculer le total
-    filter(.,  type_gismnt_solagro != "Hors gisements étude TEO Solagro") %>% 
-      group_by(campagne, zone) %>% 
-      summarise(type_gismnt_solagro = "Total périmètre étude Solagro TEO", .groups = "drop", 
-                across(.cols = c(obj_srb_solagro, tonnages_declares, gsmt_tmb), .fns = ~sum(.x, na.rm = TRUE)))
-  ) %>% 
-  mutate(tx_mob_gsmt = round((tonnages_declares / gsmt_tmb * 100), 3)) 
-
 
 
 creer_tab_solagro_zone_evol <- function(terr = "reg") {
@@ -459,6 +424,35 @@ map(ter, ~creer_taux_mob_solagro_ter(terr = .x, typ = "tmb"))
 ```
 
 
+```{r passage SRB solagro}
+
+
+# TO DO : graph de comparaison des gisements
+gr_sol_1 <- gismts_solagro %>% 
+  arrange(ordre_solagro) %>% 
+  filter(type_gismnt_solagro != "Hors gisements étude TEO Solagro") %>% 
+  mutate(lib_fam = forcats::fct_inorder(str_wrap(type_gismnt_solagro, 30))) %>% 
+  select(lib_fam, obj_2030_tonnes = obj_srb_solagro, SRB = tMB_gismt_SRB, Solagro = tMB_gismt_solagro) %>% 
+  pivot_longer(cols = c(Solagro, SRB), names_to = "type_gisement", values_to = "tonnes") %>% 
+  ggplot(aes(x = lib_fam)) + 
+  scale_x_discrete(limits = rev, labels = ~str_wrap(.x, 30)) +
+  coord_flip() + 
+  scale_y_continuous(labels = ~french_formatting(.x/1000), limits = c(0, 25000000)) 
+
+
+gr_sol_1 +
+  geom_col(aes(y = tonnes, fill = type_gisement), position = "dodge") +
+  scale_fill_manual(values = c("SRB" = gouv_colors("b1")[[1]], "Solagro" = gouv_colors("f2")[[1]])) + 
+  theme(legend.position = "right" ) +
+  geom_errorbar(aes(ymin = obj_2030_tonnes-0.00001, ymax = obj_2030_tonnes+ 0.00001 , color = gouv_colors("p1")), # gouv_colors("m2")
+                position = position_dodge(width = 3)) +
+  scale_color_identity(labels = "", guide = "legend") +
+  guides(fill = guide_legend("Estimation du\ngisement", order = 1, reverse = TRUE), 
+         error_bar = guide_legend("Objectifs 2030", order = 2)) + 
+  labs(title = "Comparaison des gisements estimés aux objectifs 2030", x ="", y = "milliers de tonnes", 
+       caption = filigrane, color = "Objectifs 2030 SRB") 
+
+```
 
 
 # 1- Tonnages `r params$campagne`