]>par(cex = 0.9)cutnum <- 16 jpn <- read.shape("C:/GISdata/japan_ver61.shp")jpndata <- jpn$att.data library(maptools)ライブラリ宣言y <- subset(jpnxy[, 2], jpndata$PREF == "群馬県")text(x, y, as.character(gunmadata$CITY1), cex = 0.7, pos = 1, offset = 0, col = "black")cols <- rev(heat.colors(cutnum))gunmadata <- merge(aged, gunmadata, sort = F, by = "JCODE") jpnpoly <- Map2poly(jpn, region.id = attr(jpn, "region.id"))gunmapoly <- subset(jpnpoly, jpndata$PREF == "群馬県")box() legend(max(x) - 1.1 * (max(x) - min(x)), min(y) + 1.37 * (max(y) - min(y)), legend = names(rev(table(classes))),cex = 0.7, fill = cols2, text.col = "black")gunmadata <- subset(jpndata, jpndata$PREF == "群馬県") plot(gunmapoly, col = cols[classes], xlab = "", ylab = "", axes = F)cols2 <- heat.colors(cutnum)text(x + 2e-04, y - 2e-04, as.character(gunmadata$CITY1), cex = 0.7, pos = 1, offset = 0, col = "white")x <- subset(jpnxy[, 1], jpndata$PREF == "群馬県")title("群馬県市町村65歳以上高齢者比率（2006年）", "出典/参考：中澤港「Rによる保健医療データ解析演習」, http://phi.ypu.jp/msb/")aged <- read.delim("C:/GISdata/agedprop.txt") DD <- gunmadata$AP2006jpnxy <- get.Pcent(jpn)classes <- cut(DD, seq(min(DD), max(DD), length = cutnum + 1), include.lowest = T)