Contexte

On va s’intéresser aux données sur les accidents corporels de circulation mises à disposition par le ministère de l’Intérieur, ici en version 2018 : https://www.data.gouv.fr/fr/datasets/base-de-donnees-accidents-corporels-de-la-circulation/

Les objectifs vont être d’importer ses données, de les jointer avec des données sur la population, par département et commune, et de calculer à différents niveaux des indicateurs d’accidentalité.

On représentera ensuite ces indicateurs sur une ou plusieurs cartes qui seront exportées.

Pour cela, il faut télécharger le dossier https://github.com/tvroylandt/r-cartostat qui contient aussi les slides du cours.

Packages

Import des données

Les données sont disponibles dans le dossier data/dataset/exo.

Elles se composent :

Données accidents

  • Importer le fichier sur les accidents dans un objet appelé df_accidents18.
  • Combien de colonnes comporte-t-il ?
  • Combien de lignes ?
  • Assurez vous que le type des colonnes est cohérent avec leur contenu.
  • Sinon, recodez les colonnes qui sont mal importées.
  • Visualisez le jeu de données.

Données population

  • Importer les données départementales de population dans un objet df_pop_dep17.
  • Idem pour les communales dans df_pop_com17.

Bonus - offres

On a construit un jeu de données sur les offres par département : offres_stmt19.sas7bdat.

  • Importer ce jeu de données.

Calculs départementaux

On va regarder ce que ça donne au niveau des départements.

Regroupement

  • Regrouper les données par département.
  • Pour cela on somme les indicateurs.
  • Calculer aussi un nombre d’accidents.

Jointure

  • Jointer sur les données de population.

Calculer des indicateurs

  • Calculer le nombre d’accidents pour 10 000 habitants pour chaque département.

  • Calculer le nombre de tués pour 10 000 habitants pour chaque département.

  • Découper ces deux parts en classes (à choisir manuellement).

  • Ne conserver que les variables d’intérêt.

  • Enlever les départements à blanc (20 et 976).

  • Refaire les opérations en conservant la variable mois.

Luminosité

  • Calculer une part d’accident qui ont lieu le jour, par département.
  • Découper cette part en classes selon les quantiles.
  • Découper cette part en classes définies manuellement.
  • Garder code_dep et les variables de pourcentages.

Bonus - Agglomération

  • Construire un score de gravité pour chaque accident :

    • Indemne = 0.1 ;
    • Blessés légers = 1 ;
    • Blessés hospitalisés = 2 ;
    • Tués = 5.

  • Calculer le score moyen dans et hors agglomération, par département.

  • Calculer la différence de score entre dans et hors agglomération, par département (utiliser pivot_wider puis pivot_longer).

Recoupement avec les données spatiales

Chargement

  • Charger le fond de carte départemental qui se trouve dans data/shp.

Jointure

  • On jointe les différentes données avec les fichiers spatiaux :

    • par département ;
    • par département/mois ;
    • pour les accidents de jour ;
    • pour le score agglomération/hors agglomération.

Représentation cartographique

Aplats de couleur

  • Représenter la part des accidents qui arrivent le jour en continu.

  • Mettre un thème blanc et changer la légende.

  • Faire de même en zoomant sur les Hauts-de-France

  • Qu’observe-t-on sur la légende ?

  • Faire de même avec la variable découpée en quantile.

Ajouter le texte

  • Ajouter la valeur, arrondie à l’unité, sur les départements.
  • Si possibles, pour l’Île-de-France, ne le faire que sur le zoom.

Région

  • Charger le fond de carte régions.
  • Faire apparaître en plus gros les contours des régions.

Cercle proportionnels

  • Représenter le nombre d’accidents par des cercles proportionnels.

  • En colorer l’intérieur par la part d’accidents par habitant, en continu.

  • En colorer l’intérieur par la part de tués par habitant, en classes.

Effet de la saisonnalité

  • Faire une carte du nombre d’accidents par habitatn par mois et département.
  • On peut facetter avec facet_wrap ou facet_grid.

Bonus - Agglomération

  • Réprésenter le score moyen en et hors agglomération, avec une facette.

Niveau bassin d’emploi

On se penche maintenant sur le niveau bassin d’emploi.

Regroupement en bassins d’emploi

  • Importer la table de passage commune -> bassin d’emploi.

  • Jointer sur les accidents pour construire le niveau bassin.

  • Construire le nombre d’accidents par bassin.

  • A-t-on des trous ?

  • Sur quelles communes ? Regarder d’où ça vient sur le site de l’Insee.

  • Importer le fond de carte au niveau bassin d’emploi.

Cartes

  • Jointer sur le fond de carte.
  • Réaliser une carte du nombre d’accidents par bassin d’emploi, en cercles proportionnels.
  • Avec en surimpression les régions.

Export