0

我一直在绘制等值线图，当我将绘制的颜色与它们分配给它的数字进行比较时，它不合适。

所以这是我的数据

    zip latitude    longitude   count2.x  count2.freq   reg   colorBuckets  colma
99501   61.21680    -149.87828    AK          67      Alaska     1         #EDF8FB
35010  32.90343      -85.92669    AL        1582      Alabama    3         #99D8C9
90001   33.97291    -118.24878    CA       20970    California   6         #006D2C
20001   38.90771    -77.01732     DC         952         NA      2         #CCECE6

所以我从一开始就使用的代码如下

    library("zipcode")
library("maps")
library("plyr")
library("RColorBrewer")
colors=brewer.pal(6, "BuGn")
data(zipcode)
merge1<-merge(zipcode, tel2, by.x='zip', by.y='zip_code', all.y=TRUE)
result<- ddply(merge1, .(state), transform, count2 = count(state))
#remove NA's#
final<- result[complete.cases(result),]
#remove duplicates#
nodupl <- final[!duplicated(final$state),]
#add state to abbreviations#
nodupl$reg<-state.name[match(nodupl$count2.x, state.abb)]
#intervalle bestimmen#
nodupl$colorBuckets<- as.numeric(cut(nodupl$count2.freq, c(0,500,1000,5000,10000,15000,22000)))
#intervall legend#
text<- c("< 500", "500 - 999", "1000 - 4999","5000 - 9999", "10000 - 14999", "15000 - 22000")
#see what color is assign to where#
nodupl$colma<- colors[nodupl$colorBuckets]
map("state",regions=nodupl$reg, exact=FALSE, col = colors[nodupl$colorBuckets],  fill = TRUE,resolution = 0,lty = 0)
map("state",col = "black",fill=FALSE,add=TRUE,lty=1,lwd=1)
#legend plotten#
par(xpd=TRUE)
legend("bottomleft", text, horiz=FALSE, fill=colors)

因此，如果我看到图例，问题又是正确分配的颜色，但是如果我仔细检查列（行 count2.freq）中的数字与地图上不适合的颜色。例如，加利福尼亚的颜色很浅，但应该是深色的。有没有人看到做错了什么。另外，我在地图上定位图例时遇到了一些麻烦，所以我再也看不到地图了。我能做些什么呢？感谢您的帮助，即使是星期六

4

1 回答 1

0

阿拉斯加地区在该map功能中不可用。nodupl$reg因此，您的地图显示了（即阿拉巴马州和加利福尼亚州）的第二个和第三个条目。但是使用了您的第一种和第二种颜色。

要以所需颜色打印这些状态，请使用命令

map("state", regions=nodupl$reg[2:3], exact=FALSE,
    col = colors[nodupl$colorBuckets][2:3],  fill = TRUE,resolution = 0,lty = 0)

但我也建议搜索阿拉斯加地图。

于 2012-10-20T14:41:09.530 回答