choroplethr - Vignettes - j) 行政区画レベル1の地図を作る

Rのパッケージ`choroplethr`のビネットの翻訳です．
( 原文: http://cran.r-project.org/web/packages/choroplethr/index.html )

* a) [イントロダクション](http://qiita.com/nakamichi/items/23730349f31a377fb071)
* b) [アメリカの州単位コロプレス図](http://qiita.com/nakamichi/items/b8267674b610ef6e1cab)
* c) [アメリカの郡単位コロプレス図](http://qiita.com/nakamichi/items/2f63dcfb7b1cc407d87c)
* d) [アメリカの郵便番号に基づく地図](http://qiita.com/nakamichi/items/889d3883e9b2ff8e2cdf)
* e) [国家単位のコロプレス図](http://qiita.com/nakamichi/items/5090b4dadeed42c2e13c)
* f) [国勢調査の地図](http://qiita.com/nakamichi/items/de929da9335f96f6468c)
* g) [世界銀行データの地図](http://qiita.com/nakamichi/items/66fa782bb3104dd815bf)
* h) [コロプレス図の動画](http://qiita.com/nakamichi/items/02ca0f678933405cf18f)
* i) [自分の地図を作る](http://qiita.com/nakamichi/items/99ff93769f4f67b97d18)
* j) [行政区画レベル1の地図を作る](http://qiita.com/nakamichi/items/8fde40a12c6aca2abbd3)

----------

# j) 行政区画レベル1の地図を作る #


"行政区画レベル1" とは，国家の地方行政単位のうち最大のものを指す一般的な用語です．この単位は国によって名称が異なります．たとえば，アメリカでは行政区画レベル1の単位は「州（states）」と呼ばれますが，日本では「県（prefectures）」と呼ばれます．`choroplethr`には世界のすべての国の行政区画レベル1の地図が含まれており，これらの地図のコロプレス図を作成することができます．地図の種類について詳細を知るには次のようにタイプしてください．

```{r}
library(choroplethrAdmin1)
?admin1.map
```

このビネットでは，`choroplethr`を使って日本の行政区画レベル1のコロプレス図を作成してみます．使用するデータは，日本の総務省統計局が公開している2010年の人口推計値です．

```{r}
library(choroplethr)
?df_japan_census
data(df_japan_census)
head(df_japan_census)
```
```{r}
##    region pop_2010 percent_pop_change_2005_2010 pop_density_km2_2010
## 23  aichi  7411000                          2.2               1434.8
## 5   akita  1086000                         -5.2                 93.3
## 2  aomori  1373000                         -4.4                142.4
## 12  chiba  6216000                          2.6               1205.5
## 38  ehime  1431000                         -2.5                252.1
## 18  fukui   806000                         -1.9                192.4
```

## ステップ1: 県名が一致しているか確認する ##

このデータのコロプレス図を作成するためには，地図とデータの両方でどのように地域名が扱われているか理解しておく必要があります．地域名が地図とデータとで異なっている場合には，修正する必要があります．（例えば，地図では "Tōkyō" でデータでは "Tokyo" となっている場合には，地図が正しく描画されません．）

まずは国名のスペリングを見てみましょう．


```{r}
head(get_admin1_countries())
```
```{r}
## [1] "afghanistan"    "aland"          "albania"        "algeria"
## [5] "american samoa" "andorra"
```

この地図では日本は"japan"（すべて小文字）として記載されています．次に日本の県名がどうなっているか見てみましょう．


```{r}
head(get_admin1_regions("japan"))
```
```{r}
##      country region
## 1618   japan  aichi
## 1619   japan  akita
## 1620   japan aomori
## 1621   japan  chiba
## 1622   japan  ehime
## 1623   japan  fukui
```

全て小文字で，末尾に「-ken」はつかないようです．これは我々のデータと同じ記法であり，地図とデータで名前は一致していますから，次のステップに進むことができます．

日本の行政区画レベル1の地図を作成したい場合は次のようにタイプします．


```{r}
admin1_map("japan")
```

![admin-l1-jp.png](https://qiita-image-store.s3.amazonaws.com/0/58917/281be26b-8070-197b-8e0f-0c47c54fa2a9.png)

## ステップ2: データを用意する ##

`choroplethr`ではデータが指定の形式に従っている必要があります．データフレームは "region" と "value" という名前の列を持っていなければなりません．`df_japan_census`はregion列といくつかのデータを含む列を持っていますが，valueという名前の列はありません．例として使用するために， pop_2010 列を value 列にコピーしましょう．

```{r}
df_japan_census$value = df_japan_census$pop_2010
```

## ステップ3: データを探索する ##

データを`choroplethr`で利用する準備ができたので，`admin1_choropleth`関数でデータを見ることができます．


```{r}
library(choroplethr)
?admin1_choropleth
admin1_choropleth(country.name = "japan",
                  df           = df_japan_census,
                  title        = "2010 Japan Population Estimates",
                  legend       = "Population",
                  num_colors   = 1)
```
```{r}
## Warning in self$bind(): The following regions were missing and are being
## set to NA: NA
```

![jp-pop-est-2012.png](https://qiita-image-store.s3.amazonaws.com/0/58917/f47d9a23-1b20-41ad-cf95-59969bcad275.png)

どのコロプレス図の関数でも共通ですが，特定の地域にズームインしたり，値のスケールを変更したりすることもできます．例として，日本の関西地方を拡大して，連続的スケールで表示してみましょう．

```{r}
kansai = c("mie", "nara", "wakayama", "kyoto", "osaka", "hyogo", "shiga")
admin1_choropleth(country.name = "japan",
                  df           = df_japan_census,
                  title        = "2010 Japan Population Estimates - Kansai Region",
                  legend       = "Population",
                  num_colors   = 1,
                  zoom         = kansai)
```

![jp-kansai-pop-est-2012.png](https://qiita-image-store.s3.amazonaws.com/0/58917/84b157da-9fca-bf05-f04b-7ef52e99bb34.png)