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

全国最低賃金

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

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

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から表をとってきてマージすればよい.

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)