【R】都道府県別　世帯数を日本地図に表示

総務省統計局が公表したデータを基に、平成27年度の世帯数と世帯数当たりの人員数を日本地図上に表示します。

library(leaflet)


dat <- read.csv("http://www.dinov.tokyo/Data/JP_Pref/Pref_population.csv", header = TRUE,
                 fileEncoding="CP932")
dat$人員_世帯<-dat$総人口/dat$世帯数
datc_m <- cut(dat$世帯数/1000000, hist(dat$世帯数/1000000, plot=FALSE)$breaks, right=FALSE)
datc_mcol <- rainbow(length(levels(datc_m)), start=0, end=0.9, alpha=0.6)[as.integer(datc_m)]
datc_h <- cut(dat$人員_世帯, hist(dat$人員_世帯, plot=FALSE)$breaks, right=FALSE)
datc_hcol <- rainbow(length(levels(datc_h)), start=0, end=0.9, alpha=0.6)[as.integer(datc_h)]
windowsFonts(JP1=windowsFont("MS Gothic"),
             JP2=windowsFont("MS Mincho"),
             JP3=windowsFont("Meiryo"),
             JP4=windowsFont("Biz Gothic"))
windows(width=1600, height=800)
par(family="JP4")
layout(matrix(1:2, 1, 2))

library(NipponMap)
JapanPrefMap(datc_mcol, main="平成27年都道府県別世帯数(百万世帯) ")
legend("bottomright", fill=rainbow(length(levels(datc_m)), alpha=0.6), legend=names(table(datc_m)))
JapanPrefMap(datc_hcol, main="平成27年都道府県別世帯当たりの人員数(人／世帯)")
legend("bottomright", fill=rainbow(length(levels(datc_h)), alpha=0.6), legend=names(table(datc_h)))

やはり東京を含めた大都市が世帯当たりの人員数が少なく、核家族化、少子化が進んでいるようです。