Fusionner des attributs de données
Par défaut, merge() fusionne selon les colonnes qui portent le même nom dans les deux jeux de données. Dans votre cas, ce n'est pas approprié, puisque la colonne d'identifiants s'appelle tract dans un jeu de données et TRACTCE dans l'autre. Pour gérer cela, merge() offre les arguments optionnels by.x et by.y, où vous pouvez préciser le nom de la colonne sur laquelle faire la fusion, respectivement, dans les deux jeux de données.
merge() retourne un nouvel objet Spatial___DataFrame. Vous pouvez donc examiner le résultat en le représentant avec tmap.
Allons-y et fusionnons.
Cette activité fait partie du cours
Visualiser des données géospatiales avec R
Instructions de l’exercice
- Utilisez
merge(), en passant d'abord l'objet spatialnyc_tracts, puis la trame de donnéesnyc_income. Indiquezby.x = "TRACTCE"etby.y = "tract". Stockez le résultat dansnyc_tracts_merge. - Utilisez
summary()surnyc_tracts_mergepour vérifier que l'objet est unSpatialPolygonsDataFrameet que les données contiennent aussi la colonneestimaterequise denyc_income. - Utilisez
tm_shape()et ajoutez une couchetm_fill()pour créer une carte choroplèthe denyc_tracts_merge, en associant la couleur àestimate.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
library(sp)
library(tmap)
# Merge nyc_tracts and nyc_income: nyc_tracts_merge
# Call summary() on nyc_tracts_merge
# Choropleth with col mapped to estimate