编辑 7:
经过相当多的帮助,我已经能够得到一张接近我需要的结果的地图。但是我仍然需要在地图上显示州界,但我无法弄清楚。为了制作一个合适的可重现示例,我需要链接到数据集,因为 dput 太大了。
为了简单起见,我只对三个状态进行了子集化,但边界线没有出现。我希望能够让边界线像白线一样穿过情节,就像它们在地图的其余部分一样。谢谢你的帮助。
数据集:
https://www.dropbox.com/s/0evuvrlm49ab9up/PRISM_1895_db.csv?dl=0
代表代码:
PRISM_1895_db <- read.csv("PRISM_1895_db.csv")
regions<- c("north dakota","south dakota","nebraska","kansas","oklahoma","texas","minnesota","iowa","missouri","arkansas", "illinois", "indiana", "wisconsin")
ggplot() +
geom_polygon(data=subset(map_data("state"), region %in% regions), aes(x=long, y=lat, group=group), col="white") +
geom_point(data = PRISM_1895_db2, aes(x = longitude, y = latitude, color = APPT), alpha = .5, size = 3.5)
图: