# データの読み込み x <- read.delim("http://minato.sip21c.org/ldaR/japancensus2010tp.txt", encoding="CP932") # ふくらみ指数改を計算する関数 source("http://minato.sip21c.org/ldaR/PEI.R") # データから都道府県別の男性のふくらみ指数改を計算 malePEI <- as.vector(by(x$Males, x$Area, PEI, CLS=5))[2:48] # ふくらみ指数改から6段階の塗り分け色 classes <- cut(malePEI, c(9, 9.5, 10:14)*10, include.lowest=TRUE) cols <- cm.colors(7)[-1] layout(t(1:2)) # using mapdata library(mapdata) # ふくらみ指数改と参照都道府県名からなるデータフレーム定義 PN <- c("hokkaido","aomori","iwate","miyagi","akita","yamagata","fukushima", "ibaraki","tochigi","gunma","saitama","chiba","tokyo","kanagawa", "niigata","toyama","ishikawa","fukui","yamanashi","nagano","gifu", "shizuoka","aichi","mie","shiga","kyoto","osaka","hyogo","nara", "wakayama","tottori","shimane","okayama","hiroshima","yamaguchi", "tokushima","kagawa","ehime","kochi","fukuoka","saga","nagasaki", "kumamoto","oita","miyazaki","kagoshima","okinawa") map("japan", type="n") for (i in PN) { map("japan", region=i, fill=TRUE, add=TRUE, col=cols[classes[PN==i]]) } legend("bottomright", legend=names(table(classes)), cex=1, fill=cols) title("都道府県別ふくらみ指数(改)によるコロプレス図 (男性,2010年国勢調査) ") # using NipponMap (definition of PN is unnecessary) library(NipponMap) JapanPrefMap(cols[classes], inset=FALSE) legend("bottomright", legend=names(table(classes)), cex=1, fill=cols) title("都道府県別ふくらみ指数(改)によるコロプレス図 (男性,2010年国勢調査) ")