Я пытаюсь изобразить округа в континентальной части США по значению, которое округа с ними связали. Я использую пакет spmap
. Мои шейп-файлы взяты из: https://www.census.gov/geo/maps-data/data/cbf/cbf_counties.html. Я использую версию 500k.
capture ssc install spmap
capture ssc install shp2dta
Разархивирую файл. Затем я загружаю карту следующим образом:
shp2dta using "cb_2015_us_county_500k.shp", database(county2015db) coordinates(county2015coord) replace
use county2015db, clear
destring STATEFP, replace
destring GEOID, replace
Если я не хочу затенять округа, все работает отлично:
spmap using county2015coord if STATEFP<60 & STATEFP!=2 & STATEFP!=15, id(GEOID)
Однако заполняется только несколько округов (и карта намного меньше, и появляется Аляска), если я делаю следующее (просто используя _ID в качестве значений для этого игрушечного примера):
spmap _ID using county2015coord if STATEFP<60 & STATEFP!=2 & STATEFP!=15, id(GEOID) fcolor(Blues)
Что делать?