【R】都道府県別　農協数

2020年2月22日

総務省統計局が発表している都道府県別の農協数を第六十八回日本統計年鑑　平成３１年から取得し、地図上にプロットしてみます。

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_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=800, height=800)
par(family="JP4")

library(NipponMap)
JapanPrefMap(datc_mcol, main="農協数 （2016年）　　(総務省統計局より) ")
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

都道府県	農協数
北海道	233
青森	60
岩手	152
宮城	89
秋田	241
山形	23
福島	127
茨城	35
栃木	23
群馬	88
埼玉	26
千葉	63
東京	22
神奈川	16
新潟	27
富山	34
石川	20
福井	26
山梨	75
長野	88
岐阜	17
静岡	143
愛知	35
三重	23
滋賀	26
京都	12
大阪	37
兵庫	36
奈良	8
和歌山	15
鳥取	45
島根	19
岡山	23
広島	30
山口	26
徳島	31
香川	5
愛媛	19
高知	33
福岡	37
佐賀	5
長崎	24
熊本	39
大分	12
宮崎	69
鹿児島	27
沖縄	4