#### 1_ Script création de la base de données pour faire un arbre html ##############################

####Chargement des packages ####
library(readxl)
library(writexl)
library(tidyverse)
library(here)


#### Chargement des bases ####

DONNEES_2020<- readRDS("~/Arborescence/Données/inosys_2020.rds")

DONNEES_2010<- readRDS("~/Arborescence/Données/inosys_2010.rds")


#####Chargements des fonctions 

source("~/Arborescence/Utilitaire/utilitaires.R")

#### création de la fonction pour calculer les effectifs par niveau avec le niveau précèdant
effectif<- function(data, group_var1, group_var2, group_var3) {
  data %>%
    mutate(nb_exploit = 1) %>%
    group_by({{ group_var1 }}, {{ group_var2 }}, {{ group_var3 }}) %>%
    summarise(NBEXP = secret_stat(nb_exploit, sum, na.rm = TRUE)) %>%
    rename(NIV = {{ group_var1 }}, NIV2 = {{ group_var2 }}) %>%
    filter(!grepl("^Exploitation_non_affectee_niv", NIV2) & NIV2 != 'Non_repechees')}


# Calcul des effectifs pour l'année 2010 ####
t<-DONNEES_2010%>%
  mutate(nb_exploit=1)%>%
  summarise(NBEXP=secret_stat(nb_exploit,sum,na.rm=TRUE))%>%
  mutate(NIV=NA)%>%
  mutate(NIV2="RA")

t0<-DONNEES_2010%>%
  mutate(nb_exploit=1)%>%
  group_by(TYPO_NATIONALE_APCA_NIV0)%>%
  summarise(NBEXP=secret_stat(nb_exploit,sum,na.rm=TRUE))%>%
  mutate(NIV="RA")%>%
  rename(NIV2=TYPO_NATIONALE_APCA_NIV0)

t1<-effectif(DONNEES_2010,TYPO_NATIONALE_APCA_NIV0,TYPO_NATIONALE_APCA_NIV1)

t2<-effectif(DONNEES_2010,TYPO_NATIONALE_APCA_NIV1,TYPO_NATIONALE_APCA_NIV2)

t3<-effectif(DONNEES_2010,TYPO_NATIONALE_APCA_NIV2,TYPO_NATIONALE_APCA_NIV3)


TOTAL_2010= rbind(t,t0,t1,t2,t3)%>%
  rename(NBEXP10=NBEXP)

# Suppression des objets temporaires
rm(t,t0,t1,t2,t3)

# Calcul des effectifs pour l'année 2020####

e<-DONNEES_2020%>%
  mutate(nb_exploit=1)%>%
  summarise(NBEXP=secret_stat(nb_exploit,sum,na.rm=TRUE))%>%
  mutate(NIV=NA)%>%
  mutate(NIV2="RA")

e0<-DONNEES_2020%>%
  mutate(nb_exploit=1)%>%
  group_by(TYPO_NATIONALE_APCA_NIV0)%>%
  summarise(NBEXP=secret_stat(nb_exploit,sum,na.rm=TRUE))%>%
  mutate(NIV="RA")%>%
  rename(NIV2=TYPO_NATIONALE_APCA_NIV0)

e1<-effectif(DONNEES_2020,TYPO_NATIONALE_APCA_NIV0,TYPO_NATIONALE_APCA_NIV1)

e2<-effectif(DONNEES_2020,TYPO_NATIONALE_APCA_NIV1,TYPO_NATIONALE_APCA_NIV2)

e3<-effectif(DONNEES_2020,TYPO_NATIONALE_APCA_NIV2,TYPO_NATIONALE_APCA_NIV3)


TOTAL_2020= rbind(e,e0,e1,e2,e3)%>%
  rename(NBEXP20=NBEXP)

# Suppression des objets temporaires
rm(e, e0, e1, e2, e3)

# Exportation vers Excel ####

write_xlsx(list("2010" = TOTAL_2010, "2020" = TOTAL_2020),
           here::here("Données", "arbre.xlsx"))

