Мне нужно создать набор из множества почтовых индексов США. Я в основном использую tigris и ggplot2, потому что я уже знаком с пакетом пакетов tidyverse. Моя проблема в том, что мне нужно ограничить свои результаты указанными c состояниями. Следующий код должен только изображать Орегон, но я получаю все США. Кто-нибудь может указать, что я делаю неправильно?
library(tigris)
library(ggmap)
library(tidyverse)
library(broom)
library(maptools)
library(rgeos)
library(sp)
options(tigris_use_cache = TRUE)
#download zctas from tigris package
zips <- zctas(cb = TRUE, state = 'oregon')
#turn shape file into dataframe
zip_df <- tidy(zips)
ggplot()+
geom_polygon(data = zip_df, aes(x=long, y=lat, group = group), color='red', fill = NA)