イタリアの州の相対的な生産性に関するデータがあります。生産性レベルに応じて 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)
ここで、簡単にするために、5 つのグループは 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) と同じスタイルで塗りつぶしたいかもしれませんが、線は逆方向または垂直方向に進んでいます。この場合でも、私はそれに対処することができませんでした。
ありがとうございました