Chômage par département

Nous reprenons la table chomage_par_dep pour afficher le pourcentage de chômages par département.

library(maps)
library(rgdal)
library(sp)
library(plotrix)   # Créer des échelles de couleurs
library(classInt)

# Nous reprenons le fichier que nous avons créé précédemment
chomage_par_dep <- read.csv(file = "chomage_par_dep.csv", sep=";", header = TRUE,dec = ',')

# Nous enlevons la première colonne qui est inutile.
chomage_par_dep <- chomage_par_dep[,-1]

chomage_par_dep[,1] <- as.character(chomage_par_dep[,1])

# Nous ajoutons des 0 pour les codes départements comportant un seul chiffre
for (i in 1:length(chomage_par_dep[,1])){
  if(nchar(chomage_par_dep[i,1])==1){
    chomage_par_dep[i,1]=paste(c("0"),chomage_par_dep[i,1],sep="")
  }
}


# Lecture des communes
departement <- readOGR(dsn="C:/Users/Pauline/Desktop/GEOFLA_2-1_DEPARTEMENT_SHP_LAMB93_FXX_2015-12-01/GEOFLA/1_DONNEES_LIVRAISON_2015/GEOFLA_2-1_SHP_LAMB93_FR-ED152/DEPARTEMENT", layer="DEPARTEMENT")

# Jointure entre departement et chomage_par_dep pour récupérer les coordonnées géographiques
names(chomage_par_dep)[1] <- c("CODE_DEPT")
departement <- merge(departement,chomage_par_dep[,c(1,5)], by.x='CODE_DEPT',by.y='CODE_DEPT')

# Lecture des limites des communes
frontiere_departement <- readOGR(dsn="C:/Users/Pauline/Desktop/GEOFLA_2-1_DEPARTEMENT_SHP_LAMB93_FXX_2015-12-01/GEOFLA/1_DONNEES_LIVRAISON_2015/GEOFLA_2-1_SHP_LAMB93_FR-ED152/DEPARTEMENT", layer="LIMITE_DEPARTEMENT")
View(frontiere_departement)

frontiere <- frontiere_departement[frontiere_departement$NATURE %in% 'Limite de département',]

europe <- readOGR(dsn="C:/Users/Pauline/Desktop/ne_110m_admin_0_countries", layer="ne_110m_admin_0_countries")
europe <- europe[europe$region_un=="Europe",]


col <- findColours(classIntervals(
  departement@data[,12], 100, style="quantile"),
  smoothColors("white",98,"purple"))

# Légende
leg <- findColours(classIntervals(
  round(departement@data[,12]),6,style="quantile"),
  smoothColors("white",3,"purple"),
  under="moins de", over="plus de", between="–",
  cutlabels=FALSE)

# Projection en Lambert 93
europe <- spTransform(europe, CRS("+init=epsg:2154"))

# Traçage de la carte

plot(frontiere,  col="#FFFFFF")
plot(europe,  col="#E6E6E6", border="#AAAAAA",lwd=1, add=TRUE)
plot(frontiere,  col="#D8D6D4", lwd=6, add=TRUE)
plot(departement,col=col, border=col,lwd=.1, add=TRUE)
title(main = "Taux de chômage par département pour le premier tour")

legend("topright",fill=attr(leg, "palette"),
       legend=names(attr(leg,"table")),
       title = "Taux de chômage :", cex=0.75)