我有意大利各省相对生产力的数据。我根据生产力水平建立了5个组。最后,我想使用意大利各省的地图来表示这种情况,因为我不能使用颜色,因为地图的打印将是灰色的。我使用了以下内容:
library(spdep)
library(spectralGP)
library(maptools)
library(gstat)
library (RColorBrewer)
plot(submap.a,col="beige",border="Black",axes=T, xlim=c(7,20), ylim=c(36,48))
plot(submap.b,col="darkgrey",border="Black",add=T)
plot(submap.c,col="black",border="Black",add=T,density=24)
plot(submap.d, col="grey30",border="Black",add=T)
plot(submap.e,col="grey13",border="Black", add=T)
其中,为简单起见,五个组是 a,b,c,d,e
结果在这里http://s9.postimg.org/m5co5g3nz/Rplot.png
我想做两件我无法弄清楚的事情。
首先我需要添加图例,我尝试过
legend_map <- c("a", "b","c", "d", "e")
legend(17,47,legend=legend_map,fill=c("white","darkgrey", ??,"grey30","grey13"),bty="n",cex=0.6)
我无法找到一种方法来插入组 (c) 的代码来代替 ??。我试过了
"black"[density=24]
"black"$density=24
但他们不工作
其次,为了让事情看起来更清晰,我可能想用与 (c) 组相同的样式填充 (e) 组省份,但线条走向相反方向甚至垂直。即使在这种情况下,我也无法处理它。
谢谢