У меня есть данные об относительной продуктивности итальянских провинций. Я построил 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
Я хотел бы сделать еще 2 вещи, которые я не смог понять.
Сначала мне нужно добавить легенду, я пробовал с
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
но они не работают
Во-вторых, чтобы все выглядело яснее, я могу захотеть заполнить провинции группы (e) тем же стилем, что и группа (c), но с линиями, идущими в обратном направлении или даже вертикально. Даже в этом случае я не смог с этим справиться.
Спасибо