Question

MERGED2014_15_PP <- read_csv("/MERGED2014_15_PP.csv")
a1 = subset(MERGED2014_15_PP, !is.na(SATMTMID), !is.na(SATVRMID))
a2 = subset(a1, !(SATMTMID=="NULL" | SATVRMID=="NULL"))
a3 = a2 %>% group_by(STABBR) %>% summarize(SATMTMIDavg = mean(as.numeric(SATMTMID)), SATVRMIDavg = mean(as.numeric(SATVRMID)))
a4 = subset(a3, !(STABBR=="PR" | STABBR=="VI" | STABBR=="DC"))
a5 = data.frame(state = tolower(abbr2state(a4$STABBR)), SATMTMIDavg=a4$SATMTMIDavg, SATVRMIDavg=a4$SATVRMIDavg)
# map_id creates the aesthetic mapping to the state name column in your data
a6 <- ggplot(a5, aes(group = 2, map_id = state)) + 
  # map points to the fifty_states shape data
  geom_map(aes(fill = SATMTMIDavg), map = fifty_states) + 
  expand_limits(x = fifty_states$long, y = fifty_states$lat) +
  coord_map() +
  scale_x_continuous(breaks = NULL) + 
  scale_y_continuous(breaks = NULL) +
  labs(x = "", y = "") +
  theme(legend.position = "right", 
        panel.background = element_blank())
a7 = a6 + aes(group = 2, fill2 = SATVRMIDavg) + scale_fill_colorplane() +
  theme(legend.position = "right")
a7

注意：部分代码来自http://www.et.bs.ehu.es/cran/web/packages/fiftystater/vignettes/fiftystater.html

我已经使用上面的代码生成了一个等值线，我现在正试图为每个机构添加MERGED2014_15_PP中的这一点。我能够生成具有经度，纬度和机构名称的数据框。但是，我无法使用geom_point添加它们，并且得到一个错误，即geom_point的aes与数据集的长度不同，即使我将它分开分组并且我有aes.inherit = FALSE。我还能尝试什么？

如何在这个等值的基础上添加点数？

0 个答案: