我正在尝试使用R“ maps”包按县制作美国的热图/峰顶/专题图。不幸的是，数据库中似乎缺少一些县。
library(maps)
data(county.fips)
which(county.fips[,1] == 35006)
which(county.fips[,2] == 'new mexico,cibola')
结果未找到任何条目
> library(maps)
> data(county.fips)
> which(county.fips[,1] == 35006)
integer(0)
> which(county.fips[,2] == 'new mexico,cibola')
integer(0)
我知道这个fips代码存在！ （http://en.wikipedia.org/wiki/Cibola_County,_New_Mexico
http://en.wikipedia.org/wiki/List_of_counties_in_New_Mexico）
有人知道我可以如何充分利用自己的情节吗？我可以找到的所有信息都列在这里：http://grokbase.com/t/r/r-sig-geo/0964515wmd/missing-county-in-maps
谢谢你的帮助！！
最佳答案
啊哈，问题不在于地图，而在于county.fips。 Cibola县没有出现在该数据集中，因此当我将我的数据{变量至绘图，FIPS代码}与列表进行匹配时，Cibola始终会丢失。
有一种解决方法可显示西博拉县...为县写一条不涉及County.fips数据集的显式行。
countynames <- match(county.fips[,"fips"], mydatavector$fips)
cibola_value <- mydatavector[which(mydatavector$fips == 35006)]
m <- map("county")
map("county", countynames, col=rainbow(100)[mydatavector$var_to_plot], fill=TRUE, add=TRUE)
map("county", "new mexico,cibola", col=rainbow(100)[cibola_value], fill=TRUE, add=TRUE)
注意：我还发现map函数在绘制美国所有县区时遇到了麻烦，但是如果我将通话分为2组，它会起作用。
关于r - R map 包缺少县，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/13941586/