Rによる空間データの統計分析(1)

空間データの分類と可視化

library(maptools)
library(classInt)

#shpファイル読み込み(IDvarには、データのID的なものを指定すると後のテーブル合成などの時に役立つ)
jpn_pref<-readShapePoly("jpn_pref.shp",IDvar="PREF_CODE")
jpn_COD<-read.csv("COD.csv")
#属性テーブルと空間データテーブルとの結合(IDvarで指定した列で自動マッチング?)
jpn_pref_code<-spCbind(jpn_pref,jpn_COD)

pa10<-c("white","grey","grey2") #カラーパレット
#等量分類
q_pref<-classIntervals(round(jpn_pref_code$malignant,2),n=5,style="quantile")
plot(q_pref,pal=pa10,cex.axis=1.3,cex.lab=1.2,lwd=2,main="") #累積密度分布図
q_pref_COL<-findColours(q_pref,pa10)
plot(jpn_pref_code,col=q_pref_COL)

参考資料