【R】都道府県別　乗用車保有台数

自動車検査登録情報協会が発表している都道府県別の乗用車保有台数（刑事土砂を含む）を描画してみました。また、人口（15～64歳）一人当たりの保有台数も計算して描画しました。

library(leaflet)

dat <- read.csv("http://www.dinov.tokyo/Data/JP_Pref/Pref_data.csv", header = TRUE,
                fileEncoding="UTF-8")
col_start <- 0.2
col_end <- 0.0

dat$乗用車率<-dat$乗用車 / dat$人口15_64

datc_m <- cut(dat$乗用車, hist(dat$乗用車, plot=FALSE)$breaks, right=FALSE)
datc_mcol <- rainbow(length(levels(datc_m)), start = col_start, end=col_end)[as.integer(datc_m)]
datc2_m <- cut(dat$乗用車率, hist(dat$乗用車率, plot=FALSE)$breaks, right=FALSE)
datc2_mcol <- rainbow(length(levels(datc2_m)), start = col_start, end=col_end)[as.integer(datc2_m)]
windowsFonts(JP1=windowsFont("MS Gothic"),
             JP2=windowsFont("MS Mincho"),
             JP3=windowsFont("Meiryo"),
             JP4=windowsFont("Biz Gothic"))
windows(width=1700, height=600)
par(family="JP4")
layout(matrix(1:2, 1, 2))

library(NipponMap)
JapanPrefMap(datc_mcol, main="令和1年 都道府県別 乗用車保有台数　台　(自動車検査登録情報協会より) ")
legend("bottomright", fill=rainbow(length(levels(datc_m)), start = col_start, end=col_end), legend=names(table(datc_m)))
JapanPrefMap(datc2_mcol, main="令和1年 都道府県別 15～64歳人口あたりの乗用車保有率　台/人 ")
legend("bottomright", fill=rainbow(length(levels(datc2_m)), start = col_start, end=col_end), legend=names(table(datc2_m)))

単純な保有台数で比べると人口の多い大都市が多い傾向にあるように見えますが、人口一人当たりの保有台数を比べるとむしろ交通の便があまりよくない地方が多いようです。

出典：一般財団法人　自動車検査登録情報協会　都道府県別・車種別保有台数表（令和1年8月分）