Создание сюжета карты в R с использованием пакета plot_usmap - PullRequest
0 голосов
/ 16 ноября 2018

У меня есть набор данных, содержащий разные состояния, почтовые индексы и количество заявок, каждое в отдельных столбцах. Я пытаюсь создать график, показывающий общее количество претензий в соответствии с почтовыми индексами штата Массачусетс.

Dataset: enter image description here

Я использовал это, чтобы отфильтровать по МА:

MA_medicare <- medicare %>%
  filter(medicare$NPPES.Provider.State == "MA")

Затем я использовал это, чтобы установить код fips для plot_usmap:

MA_medicare$NPPES.Provider.State <- fips(MA_medicare$NPPES.Provider.State)
setnames(MA_medicare, old=c("NPPES.Provider.State"), new=c("fips"))

И последний пытался построить график (не уверен, почему это не работает):

plot_usmap(data = MA_medicare, values= c("Total.Claim.Count", "NPPES.Provider.Zip.Code"), include = c("MA")) + scale_fill_continuous(low= "white", high= "red") + theme(legend.position = "right") 

Ошибка: эстетика должна иметь длину 1 или совпадать с данными (4350838): заполнить

1 Ответ

0 голосов
/ 07 января 2019

Я разработчик usmap. plot_usmap принимает только один столбец значений для построения графика, поэтому вы, вероятно, ищете следующее:

plot_usmap(data = MA_medicare, values = "Total.Claim.Count", include = c("MA"))

Однако ваши данные представлены по почтовому индексу, и в настоящее время usmap не поддерживает карты почтовых индексов (только карты штатов и уровней графства). Он использует столбец FIPS для назначения цветов штатам / округам на карте. Поскольку вы определили коды FIPS по штатам, вы просто залите весь штат Массачусетс одним сплошным цветом.

...