【Ｒ】都道府県別　病院数

2020年2月16日

総務省統計局が発表している都道府県別の病院数を第六十四回日本統計年鑑（平成27年）から取得し、地図上にプロットしてみます。また、総人口100人当たりの病院数を計算し、併せてプロットします。

  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$単位病院数<-dat$病院数/dat$総人口*100
    
  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")
  par(family="JP4")
  layout(matrix(1:2, 1, 2))
  
  library(NipponMap)
  JapanPrefMap(datc_kcol, main="病院数（平成24年）　　(総務省統計局より) ")
  legend("bottomright", fill=rainbow(length(levels(datc_k)), start = col_start, end=col_end), legend=names(table(datc_k)))
  JapanPrefMap(datc_mcol, main="100人当たりの病院数 ")
  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

都道府県	病院数
北海道	574
青森	102
岩手	92
宮城	142
秋田	74
山形	68
福島	130
茨城	183
栃木	109
群馬	132
埼玉	347
千葉	278
東京	641
神奈川	345
新潟	131
富山	110
石川	100
福井	72
山梨	60
長野	131
岐阜	103
静岡	184
愛知	325
三重	102
滋賀	59
京都	173
大阪	535
兵庫	349
奈良	75
和歌山	89
鳥取	45
島根	54
岡山	171
広島	248
山口	148
徳島	114
香川	93
愛媛	143
高知	133
福岡	466
佐賀	108
長崎	159
熊本	214
大分	159
宮崎	140
鹿児島	261
沖縄	94