addPolygons() 関数
所得データが欠損している郵便番号を確認し、なぜデータがないのか仮説を立ててみましょう。
今回マッピングするのは実際のZIPコードではなくZCTAなので、NC(ノースカロライナ州)のすべての地域に境界があるわけではありません。ファイルサイズを小さくするために簡略化しているため、境界が重なる場合もあります。ポリゴンをマッピングする際には、こうしたトレードオフを考慮する必要があります。
前置きはこのくらいにして、地図を作りましょう。shp データは、追加の引数を指定しなくても、そのまま leaflet()、addTiles()、addPolygons() にパイプして、ノースカロライナ州のZIPコードをマップできます。
スタートしやすいように、IRSの所得変数を含む shp の SpatialPolygonsDataFrame はすでに読み込まれています。
この演習はコースの一部です
R で作る leaflet のインタラクティブマップ
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# map the polygons in shp
___ %>%
leaflet() %>%
addTiles() %>%
___()