# encording with Shift-Jis

#
# world bank indicatorからのデータの取得
#

#
# 1. World bank open dataの利用
#

# World bank open data
# https://data.worldbank.org/

# World bank open data indicator
# https://data.worldbank.org/indicator


#
# 2. パッケージの読み込み
#

# World bank indicatorにアクセスするためのパッケージ

# Introduction to the wbstats R-package
# https://cran.r-project.org/web/packages/wbstats/vignettes/Using_the_wbstats_package.html


# install.packages("wbstats") #一度実行すればよい

#パッケージの読み込み
library(wbstats)

#
# 3. データのダウンロード
#

# データのIDは World bank indicatorsの中にある、変数を選択後、detailsを見て、調べることができる。

#NY.GDP.PCAP.KD : GDP per capita (constant 2010 US$)
GDPPC_2018_dat <- wb(indicator = "NY.GDP.PCAP.KD",country = "all",freq="Y", startdate=2018, enddate=2018)

# 国別データを取り出す
## ISO3c
GDPPC_JPN_2000_2018_dat <- wb(indicator = "NY.GDP.PCAP.KD",country = "JPN" ,freq="Y", startdate=2000, enddate=2018)

## ISO2c
GDPPC_JPN_2000_2018_dat <- wb(indicator = "NY.GDP.PCAP.KD",country = "JP",freq="Y", startdate=2000, enddate=2018)



#データフレームの要素は[行, 列]で指定
# 1から44行目を削除する。
GDPPC_2018_country_dat <-GDPPC_2018_dat[-1:-44,]

#記述統計
summary(GDPPC_2018_country_dat)
# ヒストグラムの作成
hist(GDPPC_2018_country_dat$value)

# 複数年にわたるデータの取得
# GDPPC_1993_2018_dat <- wb(indicator = 'NY.GDP.PCAP.PP.CD',country = 'all', startdate=1993, enddate=2018)

#
# 4. 複数のデータをダウンロード
#

#下記のデータを一括してダウンロード

#   'NY.GDP.MKTP.KD':'GDP (constant 2010 US$)',
#   'NY.GDP.PCAP.KD':'GDP per capita (constant 2010 US$)', 
#   'NE.CON.TOTL.KD':'Final consumption expenditure (constant 2010 US$)',  
#   'EN.ATM.CO2E.KD.GD':'CO2 emissions (kg per 2010 US$ of GDP)',
#   'SE.PRM.ENRR':'Gross enrollment ratio, primary, both sexes (%)', 
#   'SP.DYN.LE00.IN':'Life expectancy at birth, total (years)',
#   'SP.POP.TOTL':'Population,total',
#   'SL.UEM.TOTL.ZS':'Unemployment, total (% of total labor force) (modeled ILO estimate)',
#   'FR.INR.RINR':'Real interest rate (%)',
#   'IT.NET.USER.ZS':'Individuals using the Internet (% of population)',
#   'IT.CEL.SETS':'Mobile cellular subscriptions',
#   'IS.RRS.TOTL.KM':'Rail lines (total route-km)',
#   'DT.ODA.OATL.KD': 'Net official aid received (constant 2015 US$)'
#   'AG.LND.FRST.ZS': 'Forest area (% of land area)'


# ダウンロードするデータのリストの作成
ind = c("NY.GDP.MKTP.KD", "NY.GDP.PCAP.KD", "NE.CON.TOTL.KD", "EN.ATM.CO2E.KD.GD", "SE.PRM.ENRR", "SP.DYN.LE00.IN",
        "SP.POP.TOTL", "SL.UEM.TOTL.ZS", "FR.INR.RINR", "IT.NET.USER.ZS", "IT.CEL.SETS", "IS.RRS.TOTL.KM", "DT.ODA.OATL.KD","AG.LND.FRST.ZS")

# wide形式でのデータのダウンロード
wb_2000_2018_dat <- wb(indicator = ind ,country = 'all', startdate=2000, enddate=2018, return_wide = TRUE)

# C:\EconDat\data\wbに保存
write.csv(wb_200_2018_dat, "C:/EconDat/data/wb/wb_2000_2018_dat.csv")

# 発展的課題：各自が興味を持つデータをworld bank indicatorから、探してみてください。