【R】都道府県別　そば収穫量

農林水産省が公表している 作物統計調査 から平成29年（2017年）都道府県別のそば収穫量を地図に描画します。収穫量が不明の都道府県は0にしています。

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$都道府県, 出荷量t=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="平成２９年 都道府県別 そば収穫量　t　　(農林水産省より) ")
legend("bottomright", fill=rainbow(length(levels(datc_m)), start = col_start, end=col_end), legend=names(table(datc_m)))

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

北海道が圧倒的ですが、全国的に広く作られているようです。

出典：農林水産省　 2017年　作物統計調査　 平成29年産野菜生産出荷統計