library(maptools) gunma <- read.shape("gunma.shp") aged <- read.delim("agedprop.txt") gunma$att.data <- merge(aged,gunma$att.data,sort=F,by="JCODE") xy <- get.Pcent(gunma) x <- xy[,1] y <- xy[,2] DD <- gunma$att.data$AP2003 classes <- cut(DD,seq(min(DD),max(DD),length=5),include.lowest=T) # 等間隔に区分したい場合 # classes <- cut(DD,fivenum(DD),include.lowest=T) # 階級ごとの市町村数を等しくしたい場合 table(classes) cols <- c("azure","cyan","blue","navy") mappoly <- Map2poly(gunma,region.id=as.character(gunma$att.data$JCODE)) plot(mappoly,col=cols[ordered(classes)],xlab="",ylab="",axes=F) legend(min(x)-0.1*(max(x)-min(x)),min(y),legend=names(table(classes)),cex=0.4,fill=cols) title("群馬県市町村の65歳以上高齢者割合4区分(2003年)") text(x,y,gunma$att.data$CITY1,cex=0.3,pos=1,offset=0)