RIPS()中的县FIPS是否准确?

当我运行以下代码时……

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)

我得到了一个非马里科帕县的突出价值……这是莫哈维县.

《RIPS()中的县FIPS是否准确?》

我做错了什么,或者数据是否可疑?

我正在使用maps_2.3-11

最佳答案 简短的回答是你做错了.您访问的名称与县数据库中的多边形的顺序不同.要做你想做的事,你应该使用以下内容:

map("county")
map("county", "arizona,maricopa", col="#ABCABC", fill=T, add=T)

《RIPS()中的县FIPS是否准确?》

作为替代方案,如果您确实需要按州,县名映射辅助值,则可以执行以下操作:

##  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)

《RIPS()中的县FIPS是否准确?》

点赞