【R】都道府県別　犬の登録頭数

厚生労働省が「都道府県別の犬の登録頭数と予防注射頭数等（平成25年度～平成30年度）」を公表しています。このデータから、平成30年度の都道府県別の犬の登録頭数を表示してみます。

library(leaflet)

dat <- read.csv("http://www.dinov.tokyo/Data/JP_Pref/Pref_Dog2018.csv", header = FALSE,
                 fileEncoding="CP932")

datc_m <- cut(dat$V2, hist(dat$V2, plot=FALSE)$breaks, right=FALSE)
datc_mcol <- rainbow(length(levels(datc_m)))[as.integer(datc_m)]
windowsFonts(JP1=windowsFont("MS Gothic"),
             JP2=windowsFont("MS Mincho"),
             JP3=windowsFont("Meiryo"),
             JP4=windowsFont("Biz Gothic"))
windows(width=800, height=800)
par(family="JP4")

library(NipponMap)
JapanPrefMap(datc_mcol, main="平成30年度都道府県別 犬の登録頭数 　千頭　(厚生労働省より) ")
legend("bottomright", fill=rainbow(length(levels(datc_m))), legend=names(table(datc_m)))

当然といえば当然ですが、人口の多い地域では犬の登録頭数が多そうです。意外にも北関東も多い傾向がありそうです。