Unire gli attributi dei dati
Per impostazione predefinita, merge() unisce in base alle colonne con lo stesso nome in entrambi gli insiemi di dati. Nel tuo caso non va bene, perché la colonna degli ID si chiama tract in un dataset e TRACTCE nell’altro. Per gestire questo caso, merge() ha gli argomenti facoltativi by.x e by.y, con cui puoi specificare rispettivamente i nomi delle colonne su cui effettuare l’unione nei due dataset.
merge() restituisce un nuovo oggetto Spatial___DataFrame, quindi puoi dare un’occhiata al risultato tracciandolo con tmap.
Procediamo e facciamo l’unione.
Questo esercizio fa parte del corso
Visualizzare dati geospaziali in R
Istruzioni dell'esercizio
- Usa
merge(), passando prima l’oggetto spazialenyc_tractse poi il data framenyc_income. Specificaby.x = "TRACTCE"eby.y = "tract". Salva il risultato innyc_tracts_merge. - Usa
summary()sunyc_tracts_mergeper verificare che l’oggetto sia unSpatialPolygonsDataFramee che i dati contengano anche la colonnaestimatenecessaria danyc_income. - Usa
tm_shape()e aggiungi un layertm_fill()per creare una mappa coropletica dinyc_tracts_merge, mappando il colore suestimate.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
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