commerce.R

library(tidyverse)

# Source https://www.douane.gouv.fr/la-douane/opendata?f%5B0%5D=categorie_opendata_facet%3A459&recherche_opendata=&op=
# https://www.douane.gouv.fr/sites/default/files/2022-03/14/202201-stat-national-ce-import.zip
"data/stat-national-ce-import/Libelle_PAYS.txt" %>%
  read_csv2()

"data/stat-national-ce-import/NATIONAL_NC8PAYSI.txt" %>%
  read_csv2(
    ., 
    col_names = c("flux", "mois", "annee", "code_cpf6", 
                  "code_a129", "code_nc8", "code_pays", 
                  "valeur", "masse", "usup")
  ) %>%
  glimpse()

"data/stat-national-ce-import/NATIONAL_NC8PAYSI.txt" %>%
  read_csv2(
    ., 
    col_names = c("flux", "mois", "annee", "code_cpf6", 
                  "code_a129", "code_nc8", "code_pays", 
                  "valeur", "masse", "usup")
  ) %>%
  filter(code_pays == "RU")

df_nc8 <-  read_csv2(
  file = "data/stat-national-ce-import/Libelle_NC8_2021.txt", 
  locale = locale(encoding = "latin1")
  ) 

"data/stat-national-ce-import/NATIONAL_NC8PAYSI.txt" %>%
  read_csv2(
    ., 
    col_names = c("flux", "mois", "annee", "code_cpf6", 
                  "code_a129", "code_nc8", "code_pays", 
                  "valeur", "masse", "usup")
  ) %>%
  filter(code_pays == "RU") %>%
  group_by(mois, annee, code_cpf6, code_nc8) %>%
  summarise(valeur = sum(valeur)) %>%
  filter(mois == "01", annee == 2022) %>%
  arrange(desc(valeur)) %>%
  left_join(y = df_nc8, by = c("code_nc8" = "NC8")) %>%
  write_csv(x = ., file = "imports_russie_janvier2022.csv")

"data/stat-national-ce-import/indicateur_prevalence_NR012022_COMEXT012021.txt" %>%
  read_csv2() %>%
  group_by(flux) %>%
  summarise(valeur = sum(valeur)) 

"data/stat-national-ce-import/indicateur_prevalence_NR012022_COMEXT012021.txt" %>%
  read_csv2() %>%
  group_by(flux, ) %>%
  summarise(valeur = sum(valeur))