图例标签多次显示相同的间隔

Question

我正在 R 中制作阴影图。图例为零和缺失创建单独的类别，但将它们都标记为“[0,0)”

知道如何解决这个问题吗？

library(gpclib)
library(maptools)     # loads sp library too
library(RColorBrewer) # creates nice color schemes
library(classInt)     # finds class intervals for continuous variables


CAcounty.shp <- readShapePoly(file.choose(),proj4string=CRS("+proj=longlat"))

 plotvar <- CAcounty.shp@data$mediresp
 nclr <- 8
 plotclr <- brewer.pal(nclr,"BuPu")
 class <- classIntervals(plotvar, nclr, style="quantile")
 colcode <- findColours(class, plotclr)

 plot(CAcounty.shp, xlim=c(-122.5, -117), ylim=c(32,42))
 plot(CAcounty.shp, col=colcode, add=T) 

XXXXX <- read.csv("XXXXXXXX.csv")
plotvar <- XXXXX$XXXXXX 

title(main=“XXXXXXXXXXXXXXXX",
sub=“Internal Data")
legend(-119, 42, legend=names(attr(colcode, "table")),
    fill=attr(colcode, "palette"), cex=0.6, bty="n")
 plotclr <- brewer.pal(nclr,"RdYlGn")
 colcode <- findColours(class, plotclr)
points(KPhosp$lon, KPhosp$lat, pch=16, col=colcode, cex=1)

Answer 1

这对我来说是这样的，但你的假设是错误的：

plotvar=c(rep(NA,20),rep(0,20),1:10)
nclr <- 8
plotclr <- brewer.pal(nclr,"BuPu")
class <- classIntervals(plotvar, nclr, style="quantile")
colcode <- findColours(class, plotclr)

plot(1:10)

legend(5,5, legend=names(attr(colcode, "table")),
    fill=attr(colcode, "palette"), cex=0.6, bty="n")

生成具有大量 [0,0) 值的图例，这些值实际上都是零（不是缺失值 - 通过从向量中删除 NA 进行检查，没有任何变化）。

这是因为我有 8 个分位数中断，但是太多的零填充了不止一个分位数。

如果这是您的问题，那么您需要考虑另一种将值映射到颜色的方法，这与缺失值无关！

这里有一些类似的代码，展示了如何向图例添加一个额外的项目，以指示颜色为 NA 的项目会发生什么：

set.seed(310366)
x=1:30
y=30*runif(30)
plotvar=sample(5+10*c(rep(NA,20),runif(10)))
nclr <- 8
plotclr <- brewer.pal(nclr,"Spectral")
class <- classIntervals(plotvar, nclr, style="quantile")
colcode <- findColours(class, plotclr)

plot(x,y,col="black",pch=21,bg=colcode)

legend(5,10, legend=c("Missing",names(attr(colcode, "table"))),
    fill=c("white",attr(colcode, "palette")), cex=0.6, bty="n")

我认为您必须手动执行此操作。问题不大。