廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

R {Nippon} パッケージのコロプレス図(塗り分け地図)に凡例をつける

全国最低賃金

地域別最低賃金の全国一覧 |厚生労働省 の表をプロットしてみる.

カラーパレットには RColorBrewer パッケージを使います.

f:id:abrahamcow:20161110012343p:plain

library(rvest)
library(dplyr)
url1 <-"http://www.mhlw.go.jp/stf/seisakunitsuite/bunya/koyou_roudou/roudoukijun/minimumichiran/"
tab1 <-read_html(url1) %>%
  html_table(header = TRUE)

nuriwake_chizu <-function(d,pal.name="Oranges"){
  n.FD <-nclass.FD(d)
  n.bins <-min(n.FD,9)
  colclass=cut(d,seq(min(d),max(d),length.out = n.bins+1),include.lowest = TRUE,dig.lab = 4)
  col_ <-RColorBrewer::brewer.pal(n.bins,name=pal.name)
  col_tab <-data.frame(class=factor(levels(colclass)),col=col_,stringsAsFactors = FALSE)
  outcol <-merge(data.frame(id=1:47,class=colclass),col_tab,by="class")
  outcol <- outcol[order(outcol$id),]
  old_par =par(no.readonly = TRUE)
  par(mai=c(0,0,0,0))
  p=Nippon::JapanPrefecturesMap(col=outcol$col)
  legend("bottomright",legend=col_tab$class,fill=col_tab$col,bty = "n")
  par(old_par)
  invisible(p)
}
nuriwake_chizu(d =tab1[[1]][-48,2])

Google アナリティクスデータ

最低賃金のデータは JIS コード順(北海道, 青森, ... , 沖縄の順番)に並んでいたからそのままプロットできた.

そうでないデータは JIS コード順にソートしてやる必要があります.

それをやるには例えば ISO 3166-2:JP - Wikipedia から表をとってきてマージすればよい.

f:id:abrahamcow:20161110012839p:plain

library(RGA)
authorize()
prof <-list_profiles()
dat_ga <-get_ga(profileId = prof$id[1],
              start.date = "2016-10-01",
              end.date = "2016-10-31",
              dimensions = "ga:region",
              filters = "ga:country==Japan",
              metrics = "ga:users")
dat_ga <- dat_ga %>% 
  mutate(name=sub(" Prefecture","",region))

url2 <- "https://en.wikipedia.org/wiki/ISO_3166-2:JP"
tab2 <-read_html(url2)%>%
  html_table(header = TRUE,fill=TRUE)

pref_jis <-tab2[[1]]
colnames(pref_jis) <- c("Code","namae","name")
dat_ga2 <-left_join(pref_jis,dat_ga,by="name") %>% arrange(Code)
nuriwake_chizu(dat_ga2$users)