Кажется, я только что сделал самую странную хороплет в мире.
Код прост:
df$cut <- cut_interval(df$Count, 5)
ggplot(df, aes(lon, lat)) +
geom_polygon(aes(fill = df$cut), colour = alpha("white", 1/2), size = 0.2) +
geom_polygon(data = df , colour = "white", fill = NA) +
scale_fill_brewer(palette = "PuRd")
(см. Pastebin для dput кадра данных)
Но карта сумасшедшая:
Где я ошибся?
Глядя на другие фрагменты code, я не могу понять, как это сделать.
Вы можете построить хороплет следующим образом:
library(ggplot2)
library(XML)
tab <- readHTMLTable("http://www.50states.com/abbreviations.htm", which=1)
df$region <- tolower(tab[match(df$State, tab[, 2]), 1])
states <- map_data("state")
choro <- merge(states, df, sort = FALSE, by = "region")
choro <- choro[order(choro$order), ]
ggplot(choro, aes(long, lat.x)) +
geom_polygon(aes(group = group, fill = Count), colour="grey55") +
coord_quickmap() +
facet_wrap(~Year)
Во-первых, вы должны сопоставить названия регионов в df$State с названиями в states$region. Затем данные карты объединяются с вашим фреймом данных и упорядочиваются. Наконец, сюжет — это кусок пирога с использованием ggplot.
mapsили пакетеcoroplethr. Другое дело, что у вас есть несколько лет для некоторых штатов. Интересно, нужно ли вам агрегировать данные перед созданием карты. В противном случае, возможно, вам нужно было бы выбрать конкретные годы, которые вам нужны? - person jazzurro schedule 06.02.2016