代码不适用于其他 shp 文件

Question

here 中提出的问题已解决，但不知何故不适用于其他 shp 文件。知道为什么吗？ 我一直在尝试使用来自Coloring the states according to a given variable using ggplot（第一个链接）的 shp 文件

代码：

library(maptools)
library(ggplot2)
download.file("http://geocommons.com/overlays/204743.zip", destfile="204743.zip")
dir.create("ind")
unzip("204743.zip", exdir="ind")
india <- readShapeSpatial("ind/india_state.shp")
f_india <- fortify(india)
i <- sapply(india@data$NAME, function(x) agrep(x, data$Row.Labels, max.distance=.3, ignore.case=T)[1]) 
india@data$maj <- data$MAJORITY[i]
f_india <- merge(x=f_india, y=india@data, by.x="id", by.y="ID")
ggplot(f_india, aes(x=long, y=lat, group=group, fill=maj)) + 
  geom_polygon(colour="black") 

错误的输出：

任何帮助将不胜感激。

Answer 1

您的合并不止一次地匹配每一行，从而在数据中创建重复的行。此外，还有一个排序问题，当没有匹配的行时，您会丢失多边形。

替换

f_india <- merge(x=f_india, y=india@data, by.x="id", by.y="ID")

与

f_india <- merge(x=f_india, y=unique(india@data), by.x="id", by.y="ID",all.x=T)
f_india <- f_india[with(f_india, order(id, order)), ]

编辑：

这是我得到的输出。灰色多边形是由不匹配的名称引起的，但这是我能看到的唯一问题。