【R】都道府県別　歩数

厚生労働省が、平成22年11月に実施した「国民健康・栄養調査」の結果で都道府県別の肥満及び主な生活習慣の状況から男性と女性（ともに20~歳）の一日当たりの平均歩数を表示してみました。

  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_f <- cut(dat$歩数_女性, hist(dat$歩数_女性, plot=FALSE)$breaks, right=FALSE)
  datc_fcol <- rainbow(length(levels(datc_f)), start = col_start, end=col_end)[as.integer(datc_f)]
  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_fcol, main="歩数_女性　(厚生労働省より) ")
  legend("bottomright", fill=rainbow(length(levels(datc_f)), start = col_start, end=col_end), legend=names(table(datc_f)))
  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
  

大都市を擁する都道府県の方が、歩数は少ないようですね。田舎は移動は自動車となりがちなのでしょう。