【R】都道府県別　り災世帯数

日本は、地震や台風など自然災害による被害が何かと多い国です。 総務省統計局が発表している都道府県別のり災世帯数（2017年）を第六十九回日本統計年鑑　令和2年から取得し、地図上にプロットしてみます。

  library(leaflet)
  library(knitr)
  library(kableExtra)
  library(dplyr)
  library(tidyr)
  library(stringr)
  
  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
  
  table_df<-data.frame(都道府県=dat$都道府県, り災世帯数=dat$り災世帯数)
  
  datc_k <- cut(dat$り災世帯数, hist(dat$り災世帯数, plot=FALSE)$breaks, right=FALSE)
  datc_kcol <- rainbow(length(levels(datc_k)), start = col_start, end=col_end)[as.integer(datc_k)]

  windowsFonts(JP4=windowsFont("Biz Gothic"))
  windows(width=800, height=800)
  par(family="JP4")

  library(NipponMap)
  JapanPrefMap(datc_kcol, main="り災世帯数　(総務省統計局より) ")
  legend("bottomright", fill=rainbow(length(levels(datc_k)), start = col_start, end=col_end), legend=names(table(datc_k)))

  library(clipr)
  kable(table_df, align = "c") %>%
    kable_styling(full_width = F) %>%
    column_spec(1, bold = T) %>%
    collapse_rows(columns = 1, valign = "middle") %>%
    write_clip
  

2017年は台風や豪雨の災害が多くありました。り災世帯数もそれを反映しているようです。