Ознакомьтесь с документацией Census API . ZCTA вложен в состояния в иерархии географии. Кажется немного странным, так как обычно ZCTA должны быть независимыми от штатов или округов. Более «чистый» ZCTA имеет уровень 860; тот, который вложен в состояния, имеет уровень 871. Я понятия не имею, насколько они могут отличаться.
Таким образом, для ZCTA вам нужно присвоить ему состояние в regionin
, например:
censusapi::getCensus("sf1",
vars = c("P0010001"),
region = "zip code tabulation area:*",
vintage = 2010,
key = Sys.getenv("CENSUS_KEY"),
regionin = "state:09")
Если вам нужны почтовые индексы для каждого штата, вы можете застрять apply
ing или map
пинговать коды FIPS всех штатов и подключиться к regionin
.