【R】都道府県別　使用電力量

2020年2月18日

総務省統計局が発表している都道府県別の使用電力量（電灯）を第六十四回日本統計年鑑（平成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$都道府県, 使用電力量100万kWh=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")
  layout(matrix(1:2, 1, 2))
  
  library(NipponMap)
  JapanPrefMap(datc_kcol, main="使用電力量（平成24年）100万kWh　(総務省統計局より) ")
  legend("bottomright", fill=rainbow(length(levels(datc_k)), start = col_start, end=col_end), legend=names(table(datc_k)))
  JapanPrefMap(datc_mcol, main="1人当たりの使用電力量 ")
  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

都道府県	使用電力量100万kWh
北海道	11818
青森	2911
岩手	2889
宮城	5006
秋田	2346
山形	2573
福島	4187
茨城	6285
栃木	4385
群馬	4499
埼玉	14675
千葉	12571
東京	29984
神奈川	18175
新潟	5242
富山	3076
石川	3322
福井	2354
山梨	1982
長野	5085
岐阜	4681
静岡	8189
愛知	16117
三重	4270
滋賀	3439
京都	6345
大阪	20039
兵庫	12957
奈良	3231
和歌山	2663
鳥取	1457
島根	1818
岡山	4907
広島	7071
山口	3514
徳島	1996
香川	2540
愛媛	3458
高知	1804
福岡	11401
佐賀	1938
長崎	3161
熊本	4026
大分	2791
宮崎	2491
鹿児島	3702
沖縄	2852