R: Штат и графство Очертания для choroplethrZip

Question

У меня есть набор данных с почтовыми индексами, которые я использовал choroplethrZip для построения графика данных. Я рассматриваю данные о состоянии штата и уездах. Однако почтовые индексы не обязательно соответствуют состояниям и линиям округа. Я попытался использовать reference_map = TRUE, но у него нет линий графства, и он начинает выглядеть немного занятым. Есть ли способ изменить карту ссылок по умолчанию на ту, у которой есть состояния и линии округа, без остальной части детали карты? Другими словами, мне не нужны улицы и топография.

Вот мой код с образцом данных, похожий на данные, которые я использую. Вы можете видеть проблему с состоянием границ Техаса.

#zip.regions metadata file for choroplethrZip 
data(zip.regions) 
head(zip.regions) 

#Test data file:A data.frame containing population estimates 
# for US Zip Code Tabulated Areas (ZCTAs) in 2012. 
data(df_pop_zip) 

#Create a choropleth of US Zip Codes 
zip_choropleth(df_pop_zip, 
       state_zoom="texas", 
       title="2012 Texas State ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) 

#Zoom County 
dd_fips = c(48113, 48121) 
zip_choropleth(df_pop_zip, 
       county_zoom=dd_fips, 
       title="2012 Denton & Dallas ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) 

TexasPlot

DentonDallasPlot

Answer 1

Я использовал этот блог и узнал, как сделать эту работу: http://www.arilamstein.com/blog/2015/07/02/exploring-the-demographics-of-ferguson-missouri/

Вот мой окончательный код:

library(choroplethrZip) 
library(ggplot2) 
library(choroplethr) 


#Pull in zip.regions metadata file for choroplethrZip 
data(zip.regions) 
head(zip.regions) 

#Test data file:A data.frame containing population estimates 
# for US Zip Code Tabulated Areas (ZCTAs) in 2012. 
data(df_pop_zip) 

# highlight a county 
highlight_county = function(county_fips) 
{ 
    library(choroplethrMaps) 
    data(county.map, package="choroplethrMaps", envir=environment()) 
    df = county.map[county.map$region %in% county_fips, ] 
    geom_polygon(data=df, aes(long, lat, group = group), color = "yellow", fill = NA, size = 1) 
} 

#Zoom County 
dd_fips = c(48113, 48121) 
zip_choropleth(df_pop_zip, 
       county_zoom=dd_fips, 
       title="2012 Denton & Dallas ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) + 
       highlight_county(dd_fips) 

Rplot Counties

Мне просто нужно добавить функцию highlight_county, и она отлично поработала. Когда я тестировал свои данные (вместо общих данных о населении), он также работал.