Блог

У меня есть следующие данные, которые я хочу нанести на карту штатов США:

 x = data.frame(state.x77[,"Income"])
 

Как я могу отобразить их в виде цветной заливки состояний на карте?

Я могу составить карту:

 map("state", boundary=TRUE, col=colorRampPalette(c("blue","green","yellow","red"))(50), fill=T)
 

Я могу сопоставить штаты:

 match(map("state",plot=F)$names, tolower(rownames(x)))
 

Я попытался изменить пример из ?map_data :

 if (require("maps")) {
states <- map_data("state")
arrests <- USArrests
names(arrests) <- tolower(names(arrests))
arrests$region <- tolower(rownames(USArrests))

choro <- merge(states, arrests, sort = FALSE, by = "region")
choro <- choro[order(choro$order), ]
qplot(long, lat, data = choro, group = group, fill = assault,
  geom = "polygon")
qplot(long, lat, data = choro, group = group, fill = assault / murder,
  geom = "polygon")
}
 

Но я не в состоянии это исправить. Спасибо за вашу помощь.

Вам не нужно объединять данные об арестах с данными карты. Вы можете передавать картографические данные отдельно в geom_map слой. Попробуй

 x = data.frame(region=tolower(rownames(state.x77)), 
    income=state.x77[,"Income"], 
    stringsAsFactors=F)

states_map <- map_data("state")
ggplot(x, aes(map_id = region))   
    geom_map(aes(fill = income), map = states_map)  
    scale_fill_gradientn(colours=c("blue","green","yellow","red"))   
    expand_limits(x = states_map$long, y = states_map$lat)