【R】都道府県別　労働災害

総務省統計局が発表している「統計でみる都道府県のすがた2020」から、都道府県別の労働災害発生の頻度と労働者災害補償保険給付率（％、対適用労働者数）を表示します。2017年のデータです。

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$労働災害発生頻度, 労働者災害補償保険給付率=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)]
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)]
windowsFonts(JP4=windowsFont("Biz Gothic"))
windows(width=1600, height=800)
par(family="JP4")
layout(matrix(1:2, 1, 2))

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)))
JapanPrefMap(datc_mcol, main="労働者災害補償保険給付率　％（対適用労働者数）")
legend("bottomright", fill=rainbow(length(levels(datc_m)), start = col_start, end=col_end), legend=names(table(datc_m)))

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