0

我正在使用 R 包 GWmodel。在绘制 GWR 的输出时，我发现自己无法编辑 ssplot 函数中的侧栏。

我正在使用的几行代码（数据已经在包中）就足够了。您只需要安装 GWmodel 库：

library(GWmodel)
data("DubVoter")

bw.gwr.1 <- bw.gwr(GenEl2004 ~ DiffAdd + LARent + SC1 + Unempl + LowEduc + Age18_24 + Age25_44 + Age45_64, data = Dub.voter, approach = "CV", kernel = "gaussian", adaptive = FALSE)

gwr.res <- gwr.basic(GenEl2004 ~ DiffAdd + LARent + SC1 + Unempl + LowEduc + Age18_24 + Age25_44 + Age45_64, data = Dub.voter, bw = bw.gwr.1,  kernel = "gaussian", adaptive = FALSE, F123.test = TRUE)



Greens <- c("#238B45","#74C476","#BAE4B3","#EDF8E9")

spplot(gwr.res$SDF, "Unempl", cuts = 4, at = c(fivenum(gwr.res$SDF$Unempl)), col.regions = Greens)

结果如下：

在此处输入图像描述

正如您在 spplot 函数中看到的，参数 at 的值是 c(fivenum(gwr.res$SDF$Unempl))，它对应于 gwr 中的最小值、第 1 个四分位数、中位数、第 3 个四分位数、最大值。 res$SDF$Unnempl。查看图表，我们可以看到颜色将此划分为 4 个类（使用参数 cut = 4 at = c(fivenum(gwr.res$SDF$Unempl))）。

我想做的是用调色板编辑侧栏，并用 c(fivenum(gwr.res$SDF$Unempl)) 中的值添加替换 6,4,2,0,-2，它们是：

-2.427，-0.916，-0。6.706, -0.494, 7.566

4

1 回答 1

1

这是有效的代码：

labelat = fivenum(gwr.res$SDF$Unempl)
labeltext = labelat 


spplot(gwr.res$SDF, "Unempl", cuts = 4, at = c(fivenum(gwr.res$SDF$Unempl)), col.regions = Greens,
colorkey=list(width=0.3,      
              space="right", 
              tick.number=5, 
              labels=list(  
                at=labelat, 
                labels=labeltext )))

在此处输入图像描述

请注意，侧栏上的值与 -2.427、-0.916、-0 不同。6.706、-0.494、7.566，因为数据不同。

于 2015-09-02T16:46:57.230 回答