addPolygons() Function
Let's look at those zip codes with missing data to hypothesize why they do not have income data.
We are mapping ZCTAs (not actual zip codes) so not every part of NC will have a boundary. Our boundaries may overlap because the file was simplified to reduce size. These are trade offs to consider when mapping polygons.
Enough nuance, let's make a map. You can pipe the shp
data directly into our calls to leaflet()
, addTiles()
, and addPolygons()
without supplying any additional arguments to map North Carolina's zip codes.
To get you started, the shp
SpatialPolygonsDataFrame including the IRS income variables has been loaded for you.
This exercise is part of the course
Interactive Maps with leaflet in R
Hands-on interactive exercise
Have a go at this exercise by completing this sample code.
# map the polygons in shp
___ %>%
leaflet() %>%
addTiles() %>%
___()