提问人：JoeBass 提问时间：8/13/2015 更新时间：8/13/2015 访问量：1091
R maps（） 中的县 FIPS 准确吗？
Are county FIPS accurate in R maps()?
问：
当我运行以下代码时...
require(maps)
colors <- data.frame(county=county.fips$polyname,color=rep("#FFFFFF",nrow(county.fips)), stringsAsFactors=FALSE)
colors[colors$county=="arizona,maricopa","color"] <- "#ABCABC"
map("county", col = colors$color, fill = TRUE)
我得到了一个不是马里科帕县的突出显示值......这是莫哈维县。
是我做错了什么，还是数据可疑？
我正在使用 maps_2.3-11
答：
2赞
Forrest R. Stevens
8/13/2015
#1
简短的回答是你做错了。您正在访问的名称与县数据库中的多边形的顺序不同。要执行所需的操作，应使用以下命令：
map("county")
map("county", "arizona,maricopa", col="#ABCABC", fill=T, add=T)
或者，如果确实需要按名称映射辅助值，则可以执行如下操作：state,county
## Get state,county names in the order they will be plotted:
c <- map("county", namesonly=T)
c1 <- rep("#FFFFFF", length(c))
c1[which(c=="arizona,maricopa")] <- "#FF0000"
map("county", col=c1, fill=T)
评论
1赞
hrbrmstr
8/13/2015
+1 一个更长的答案是该软件包有一个新的维护者，并且所有不准确之处（如“苏联”）都在即将发布的版本中得到处理。maps
0赞
JoeBass
8/13/2015
谢谢。不确定这是否是提出后续问题的正确方法，但是如果我有一个按 county.fips 顺序排列的 county.fips 值向量，我该如何将其应用于地图？
1赞
Forrest R. Stevens
8/13/2015
我已经修改了我的答案，希望能使用辅助数据为您提供一条前进的道路。不过，您需要将它们映射到提取的名称，如果您正在通过 FIPS 代码处理其他各种数据源，那么 chloroplethr 包最终可能会给您带来更少的麻烦和更漂亮的地图。
下一个：R 如何在组合子集上“拆分”？
评论