ComtradeFromR.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Comtradeからのデータ取得法（API&R）\n",
    "> UNによる解説ページ  \n",
    "> * [APIの解説](http://comtrade.un.org/data/doc/api/)  \n",
    "> * [Rから利用する方法の解説](http://comtrade.un.org/data/Doc/api/ex/r)  \n",
    "\n",
    "## 1. Rの関数（作成）"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "get.Comtrade <- function(url=\"http://comtrade.un.org/api/get?\"\n",
    "                         ,maxrec=250000\n",
    "                         ,type=\"C\"\n",
    "                         ,freq=\"A\"\n",
    "                         ,px=\"HS\"\n",
    "                         ,ps=\"now\"\n",
    "                         ,r\n",
    "                         ,p=\"all\"\n",
    "                         ,rg=\"all\"\n",
    "                         ,cc=\"TOTAL\"\n",
    "                         ,fmt=\"json\"\n",
    ")\n",
    "{\n",
    "  string<- paste(url\n",
    "                 ,\"max=\",maxrec,\"&\" #maximum no. of records returned\n",
    "                 ,\"type=\",type,\"&\" #type of trade (c=commodities)\n",
    "                 ,\"freq=\",freq,\"&\" #frequency\n",
    "                 ,\"px=\",px,\"&\" #classification\n",
    "                 ,\"ps=\",ps,\"&\" #time period\n",
    "                 ,\"r=\",r,\"&\" #reporting area\n",
    "                 ,\"p=\",p,\"&\" #partner country\n",
    "                 ,\"rg=\",rg,\"&\" #trade flow\n",
    "                 ,\"cc=\",cc,\"&\" #classification code\n",
    "                 ,\"fmt=\",fmt        #Format\n",
    "                 ,sep = \"\"\n",
    "  )\n",
    "  \n",
    "  if(fmt == \"csv\") {\n",
    "    raw.data<- read.csv(string,header=TRUE)\n",
    "    return(list(validation=NULL, data=raw.data))\n",
    "  } else {\n",
    "    if(fmt == \"json\" ) {\n",
    "      raw.data<- fromJSON(file=string)\n",
    "      data<- raw.data$dataset\n",
    "      validation<- unlist(raw.data$validation, recursive=TRUE)\n",
    "      ndata<- NULL\n",
    "      if(length(data)> 0) {\n",
    "        var.names<- names(data[[1]])\n",
    "        data<- as.data.frame(t( sapply(data,rbind)))\n",
    "        ndata<- NULL\n",
    "        for(i in 1:ncol(data)){\n",
    "          data[sapply(data[,i],is.null),i]<- NA\n",
    "          ndata<- cbind(ndata, unlist(data[,i]))\n",
    "        }\n",
    "        ndata<- as.data.frame(ndata)\n",
    "        colnames(ndata)<- var.names\n",
    "      }\n",
    "      return(list(validation=validation,data =ndata))\n",
    "    }\n",
    "  }\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1.1. 関数の引数について  \n",
    "___\n",
    "#### maxrec\n",
    "最大データ数.デフォルトは250,000に設定（認証ユーザーの場合の最大数）.ゲストユーザーは50,000まで.  \n",
    "___\n",
    "#### type\n",
    "貿易のタイプ  \n",
    "* \"C\":財貿易（デフォルト）\n",
    "* \"S\":サービス貿易\n",
    "___\n",
    "#### freq\n",
    "データの周期  \n",
    "* \"A\":年次（デフォルト）\n",
    "* \"M\":月次\n",
    "___\n",
    "#### px\n",
    "品目分類の種別.デフォルトは\"HS\"  \n",
    "* \"HS\":Harmonized System (HS), as reported (e.g. if data was originally submitted to UN Comtrade in HS1996 then HS1996 is displayed) (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"H0\":HS 1992 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"H1\":HS 1996 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"H2\":HS 2002 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"H3\":HS 2007 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"H4\":HS 2012 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"ST\":Standard International Trade Classification (SITC), as reported (e.g. if data was originally submitted to UN Comtrade in SITC Rev. 1 then SITC Rev. 1 is displayed) (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"S1\":SITC Revision 1 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"S2\":SITC Revision 2 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"S3\":SITC Revision 3 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"S4\":SITC Revision 4 (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "* \"BEC\":Broad Economic Categories (BEC) (full classification code list in JSON, formatted for select2.js, pass \"id\" to cc parameter) \n",
    "___\n",
    "#### ps\n",
    "データの取得期間（5期間まで、コンマで区切ってもOK?)\n",
    "* \"now\":直近1期間\n",
    "* \"recent\":直近5期間\n",
    "* YYYY:西暦を指定（2013など）\n",
    "* YYYYMM:年月を指定（月次データの場合）\n",
    "___\n",
    "#### r\n",
    "報告側の地域・国（入力必須、5地域まで、コンマで区切ってもOK?）\n",
    "* \"all\":すべての国\n",
    "* [その他の場合](http://comtrade.un.org/data/cache/reporterAreas.json)\n",
    "___\n",
    "#### p\n",
    "相手側の地域・国（5地域まで、コンマで区切ってもOK?）\n",
    "* \"all\":すべての国（デフォルト）\n",
    "* [その他の場合](http://comtrade.un.org/data/cache/partnerAreas.json)\n",
    "___\n",
    "#### rg\n",
    "貿易のタイプ（コンマで区切ってもOK?)\n",
    "* \"all\":すべてのタイプ（デフォルト）\n",
    "* \"1\":輸入\n",
    "* \"2\":輸出\n",
    "* \"4\":再輸入\n",
    "* \"3\":再輸出\n",
    "___\n",
    "#### cc\n",
    "品目分類コード（20までなら、コンマで区切ってもOK?)\n",
    "* \"TOTAL\":合計値（デフォルト）\n",
    "* \"AG1\",\"AG2\",\"AG3\",\"AG4\",\"AG5\",\"AG6\":分類の桁（\"AG2\"なら2桁分類）\n",
    "* \"ALL\":全コード\n",
    "___\n",
    "#### fmt\n",
    "データ・フォーマット\n",
    "* \"json\":json形式（デフォルト）\n",
    "* \"csv\":csv形式\n",
    "___\n",
    "### 1.2. データの使用制限\n",
    "___\n",
    "#### リクエスト制限\n",
    "* 認証ユーザーは1時間10,000リクエストまで\n",
    "* ゲストユーザーは1時間100リクエストまで\n",
    "___\n",
    "#### パラメータの組み合わせの制限\n",
    "* ps（データの取得期間）, r（報告側の地域・国） and p（相手側の地域・国） are limited to 5 codes each\n",
    "* Only one of the above codes may use the special \"ALL\" value.\n",
    "* Classification codes (cc) are limited to 20 items. \"ALL\" is always a valid classification code. \n",
    "___"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. データ取得\n",
    "### 2.1. JSON形式\n",
    "> `rjson`ライブラリを予めインストールすること"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<dl>\n",
       "\t<dt>$validation</dt>\n",
       "\t\t<dd><dl class=dl-horizontal>\n",
       "\t<dt>status.name</dt>\n",
       "\t\t<dd>\"Ok\"</dd>\n",
       "\t<dt>status.value</dt>\n",
       "\t\t<dd>\"0\"</dd>\n",
       "\t<dt>status.category</dt>\n",
       "\t\t<dd>\"0\"</dd>\n",
       "\t<dt>status.description</dt>\n",
       "\t\t<dd>\"\"</dd>\n",
       "\t<dt>status.helpUrl</dt>\n",
       "\t\t<dd>\"\"</dd>\n",
       "\t<dt>count.value</dt>\n",
       "\t\t<dd>\"6\"</dd>\n",
       "\t<dt>count.started</dt>\n",
       "\t\t<dd>\"2016-07-25T11:14:49.0335334-04:00\"</dd>\n",
       "\t<dt>count.finished</dt>\n",
       "\t\t<dd>\"2016-07-25T11:14:49.392338-04:00\"</dd>\n",
       "\t<dt>count.durationSeconds</dt>\n",
       "\t\t<dd>\"0.3588046\"</dd>\n",
       "\t<dt>datasetTimer.started</dt>\n",
       "\t\t<dd>\"2016-07-25T11:14:49.0335334-04:00\"</dd>\n",
       "\t<dt>datasetTimer.finished</dt>\n",
       "\t\t<dd>\"2016-07-25T11:14:50.3751506-04:00\"</dd>\n",
       "\t<dt>datasetTimer.durationSeconds</dt>\n",
       "\t\t<dd>\"1.3416172\"</dd>\n",
       "</dl>\n",
       "</dd>\n",
       "\t<dt>$data</dt>\n",
       "\t\t<dd><table>\n",
       "<thead><tr><th></th><th scope=col>pfCode</th><th scope=col>yr</th><th scope=col>period</th><th scope=col>periodDesc</th><th scope=col>aggrLevel</th><th scope=col>IsLeaf</th><th scope=col>rgCode</th><th scope=col>rgDesc</th><th scope=col>rtCode</th><th scope=col>rtTitle</th><th scope=col>...</th><th scope=col>qtAltCode</th><th scope=col>qtAltDesc</th><th scope=col>TradeQuantity</th><th scope=col>AltQuantity</th><th scope=col>NetWeight</th><th scope=col>GrossWeight</th><th scope=col>TradeValue</th><th scope=col>CIFValue</th><th scope=col>FOBValue</th><th scope=col>estCode</th></tr></thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>H4          </td><td>2015        </td><td>2015        </td><td>2015        </td><td>0           </td><td>0           </td><td>1           </td><td>Import      </td><td>842         </td><td>USA         </td><td>...         </td><td>NA          </td><td>            </td><td>NA          </td><td>NA          </td><td>NA          </td><td>NA          </td><td>301035596398</td><td>NA          </td><td>NA          </td><td>0           </td></tr>\n",
       "\t<tr><th scope=row>2</th><td>H4          </td><td>2015        </td><td>2015        </td><td>2015        </td><td>0           </td><td>0           </td><td>2           </td><td>Export      </td><td>842         </td><td>USA         </td><td>...         </td><td>NA          </td><td>            </td><td>NA          </td><td>NA          </td><td>NA          </td><td>NA          </td><td>279990083749</td><td>NA          </td><td>NA          </td><td>0           </td></tr>\n",
       "\t<tr><th scope=row>3</th><td>H4         </td><td>2015       </td><td>2015       </td><td>2015       </td><td>0          </td><td>0          </td><td>3          </td><td>Re-Export  </td><td>842        </td><td>USA        </td><td>...        </td><td>NA         </td><td>           </td><td>NA         </td><td>NA         </td><td>NA         </td><td>NA         </td><td>46422145927</td><td>NA         </td><td>NA         </td><td>0          </td></tr>\n",
       "\t<tr><th scope=row>4</th><td>H4          </td><td>2015        </td><td>2015        </td><td>2015        </td><td>0           </td><td>0           </td><td>1           </td><td>Import      </td><td>842         </td><td>USA         </td><td>...         </td><td>NA          </td><td>            </td><td>NA          </td><td>NA          </td><td>NA          </td><td>NA          </td><td>297499549233</td><td>NA          </td><td>NA          </td><td>0           </td></tr>\n",
       "\t<tr><th scope=row>5</th><td>H4          </td><td>2015        </td><td>2015        </td><td>2015        </td><td>0           </td><td>0           </td><td>2           </td><td>Export      </td><td>842         </td><td>USA         </td><td>...         </td><td>NA          </td><td>            </td><td>NA          </td><td>NA          </td><td>NA          </td><td>NA          </td><td>236377320352</td><td>NA          </td><td>NA          </td><td>0           </td></tr>\n",
       "\t<tr><th scope=row>6</th><td>H4         </td><td>2015       </td><td>2015       </td><td>2015       </td><td>0          </td><td>0          </td><td>3          </td><td>Re-Export  </td><td>842        </td><td>USA        </td><td>...        </td><td>NA         </td><td>           </td><td>NA         </td><td>NA         </td><td>NA         </td><td>NA         </td><td>50441857206</td><td>NA         </td><td>NA         </td><td>0          </td></tr>\n",
       "</tbody>\n",
       "</table>\n",
       "</dd>\n",
       "</dl>\n"
      ],
      "text/latex": [
       "\\begin{description}\n",
       "\\item[\\$validation] \\begin{description*}\n",
       "\\item[status.name] \"Ok\"\n",
       "\\item[status.value] \"0\"\n",
       "\\item[status.category] \"0\"\n",
       "\\item[status.description] \"\"\n",
       "\\item[status.helpUrl] \"\"\n",
       "\\item[count.value] \"6\"\n",
       "\\item[count.started] \"2016-07-25T11:14:49.0335334-04:00\"\n",
       "\\item[count.finished] \"2016-07-25T11:14:49.392338-04:00\"\n",
       "\\item[count.durationSeconds] \"0.3588046\"\n",
       "\\item[datasetTimer.started] \"2016-07-25T11:14:49.0335334-04:00\"\n",
       "\\item[datasetTimer.finished] \"2016-07-25T11:14:50.3751506-04:00\"\n",
       "\\item[datasetTimer.durationSeconds] \"1.3416172\"\n",
       "\\end{description*}\n",
       "\n",
       "\\item[\\$data] \\begin{tabular}{r|lllllllllllllllllllllllllllllllllll}\n",
       "  & pfCode & yr & period & periodDesc & aggrLevel & IsLeaf & rgCode & rgDesc & rtCode & rtTitle & ... & qtAltCode & qtAltDesc & TradeQuantity & AltQuantity & NetWeight & GrossWeight & TradeValue & CIFValue & FOBValue & estCode\\\\\n",
       "\\hline\n",
       "\t1 & H4           & 2015         & 2015         & 2015         & 0            & 0            & 1            & Import       & 842          & USA          & ...          & NA           &              & NA           & NA           & NA           & NA           & 301035596398 & NA           & NA           & 0           \\\\\n",
       "\t2 & H4           & 2015         & 2015         & 2015         & 0            & 0            & 2            & Export       & 842          & USA          & ...          & NA           &              & NA           & NA           & NA           & NA           & 279990083749 & NA           & NA           & 0           \\\\\n",
       "\t3 & H4          & 2015        & 2015        & 2015        & 0           & 0           & 3           & Re-Export   & 842         & USA         & ...         & NA          &             & NA          & NA          & NA          & NA          & 46422145927 & NA          & NA          & 0          \\\\\n",
       "\t4 & H4           & 2015         & 2015         & 2015         & 0            & 0            & 1            & Import       & 842          & USA          & ...          & NA           &              & NA           & NA           & NA           & NA           & 297499549233 & NA           & NA           & 0           \\\\\n",
       "\t5 & H4           & 2015         & 2015         & 2015         & 0            & 0            & 2            & Export       & 842          & USA          & ...          & NA           &              & NA           & NA           & NA           & NA           & 236377320352 & NA           & NA           & 0           \\\\\n",
       "\t6 & H4          & 2015        & 2015        & 2015        & 0           & 0           & 3           & Re-Export   & 842         & USA         & ...         & NA          &             & NA          & NA          & NA          & NA          & 50441857206 & NA          & NA          & 0          \\\\\n",
       "\\end{tabular}\n",
       "\n",
       "\\end{description}\n"
      ],
      "text/plain": [
       "$validation\n",
       "                        status.name                        status.value \n",
       "                               \"Ok\"                                 \"0\" \n",
       "                    status.category                  status.description \n",
       "                                \"0\"                                  \"\" \n",
       "                     status.helpUrl                         count.value \n",
       "                                 \"\"                                 \"6\" \n",
       "                      count.started                      count.finished \n",
       "\"2016-07-25T11:14:49.0335334-04:00\"  \"2016-07-25T11:14:49.392338-04:00\" \n",
       "              count.durationSeconds                datasetTimer.started \n",
       "                        \"0.3588046\" \"2016-07-25T11:14:49.0335334-04:00\" \n",
       "              datasetTimer.finished        datasetTimer.durationSeconds \n",
       "\"2016-07-25T11:14:50.3751506-04:00\"                         \"1.3416172\" \n",
       "\n",
       "$data\n",
       "  pfCode   yr period periodDesc aggrLevel IsLeaf rgCode    rgDesc rtCode\n",
       "1     H4 2015   2015       2015         0      0      1    Import    842\n",
       "2     H4 2015   2015       2015         0      0      2    Export    842\n",
       "3     H4 2015   2015       2015         0      0      3 Re-Export    842\n",
       "4     H4 2015   2015       2015         0      0      1    Import    842\n",
       "5     H4 2015   2015       2015         0      0      2    Export    842\n",
       "6     H4 2015   2015       2015         0      0      3 Re-Export    842\n",
       "  rtTitle rt3ISO ptCode ptTitle pt3ISO ptCode2 ptTitle2 pt3ISO2 cstCode cstDesc\n",
       "1     USA    USA    124  Canada    CAN    <NA>                                 \n",
       "2     USA    USA    124  Canada    CAN    <NA>                                 \n",
       "3     USA    USA    124  Canada    CAN    <NA>                                 \n",
       "4     USA    USA    484  Mexico    MEX    <NA>                                 \n",
       "5     USA    USA    484  Mexico    MEX    <NA>                                 \n",
       "6     USA    USA    484  Mexico    MEX    <NA>                                 \n",
       "  motCode motDesc cmdCode        cmdDescE qtCode      qtDesc qtAltCode\n",
       "1                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "2                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "3                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "4                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "5                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "6                   TOTAL All Commodities      1 No Quantity      <NA>\n",
       "  qtAltDesc TradeQuantity AltQuantity NetWeight GrossWeight   TradeValue\n",
       "1                    <NA>        <NA>      <NA>        <NA> 301035596398\n",
       "2                    <NA>        <NA>      <NA>        <NA> 279990083749\n",
       "3                    <NA>        <NA>      <NA>        <NA>  46422145927\n",
       "4                    <NA>        <NA>      <NA>        <NA> 297499549233\n",
       "5                    <NA>        <NA>      <NA>        <NA> 236377320352\n",
       "6                    <NA>        <NA>      <NA>        <NA>  50441857206\n",
       "  CIFValue FOBValue estCode\n",
       "1     <NA>     <NA>       0\n",
       "2     <NA>     <NA>       0\n",
       "3     <NA>     <NA>       0\n",
       "4     <NA>     <NA>       0\n",
       "5     <NA>     <NA>       0\n",
       "6     <NA>     <NA>       0\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "library(\"rjson\")\n",
    "s1 <- get.Comtrade(r=\"842\", p=\"124,484\")\n",
    "s1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "### 2.2. CSV形式\n",
    "> JSONと異なり、validationが出力されない"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "$validation\n",
       "NULL\n",
       "\n",
       "$data\n",
       "  Classification Year Period Period.Desc. Aggregate.Level Is.Leaf.Code\n",
       "1             H4 2015   2015         2015               0            0\n",
       "2             H4 2015   2015         2015               0            0\n",
       "3             H4 2015   2015         2015               0            0\n",
       "4             H4 2015   2015         2015               0            0\n",
       "5             H4 2015   2015         2015               0            0\n",
       "6             H4 2015   2015         2015               0            0\n",
       "  Trade.Flow.Code Trade.Flow Reporter.Code Reporter Reporter.ISO Partner.Code\n",
       "1               1     Import           842      USA          USA          124\n",
       "2               2     Export           842      USA          USA          124\n",
       "3               3  Re-Export           842      USA          USA          124\n",
       "4               1     Import           842      USA          USA          484\n",
       "5               2     Export           842      USA          USA          484\n",
       "6               3  Re-Export           842      USA          USA          484\n",
       "  Partner Partner.ISO X2nd.Partner.Code X2nd.Partner X2nd.Partner.ISO\n",
       "1  Canada         CAN                NA           NA               NA\n",
       "2  Canada         CAN                NA           NA               NA\n",
       "3  Canada         CAN                NA           NA               NA\n",
       "4  Mexico         MEX                NA           NA               NA\n",
       "5  Mexico         MEX                NA           NA               NA\n",
       "6  Mexico         MEX                NA           NA               NA\n",
       "  Customs.Proc..Code Customs Mode.of.Transport.Code Mode.of.Transport\n",
       "1                 NA      NA                     NA                NA\n",
       "2                 NA      NA                     NA                NA\n",
       "3                 NA      NA                     NA                NA\n",
       "4                 NA      NA                     NA                NA\n",
       "5                 NA      NA                     NA                NA\n",
       "6                 NA      NA                     NA                NA\n",
       "  Commodity.Code       Commodity Qty.Unit.Code    Qty.Unit Qty\n",
       "1          TOTAL All Commodities             1 No Quantity  NA\n",
       "2          TOTAL All Commodities             1 No Quantity  NA\n",
       "3          TOTAL All Commodities             1 No Quantity  NA\n",
       "4          TOTAL All Commodities             1 No Quantity  NA\n",
       "5          TOTAL All Commodities             1 No Quantity  NA\n",
       "6          TOTAL All Commodities             1 No Quantity  NA\n",
       "  Alt.Qty.Unit.Code Alt.Qty.Unit Alt.Qty Netweight..kg. Gross.weight..kg.\n",
       "1                NA           NA      NA             NA                NA\n",
       "2                NA           NA      NA             NA                NA\n",
       "3                NA           NA      NA             NA                NA\n",
       "4                NA           NA      NA             NA                NA\n",
       "5                NA           NA      NA             NA                NA\n",
       "6                NA           NA      NA             NA                NA\n",
       "  Trade.Value..US.. CIF.Trade.Value..US.. FOB.Trade.Value..US.. Flag\n",
       "1      301035596398                    NA                    NA    0\n",
       "2      279990083749                    NA                    NA    0\n",
       "3       46422145927                    NA                    NA    0\n",
       "4      297499549233                    NA                    NA    0\n",
       "5      236377320352                    NA                    NA    0\n",
       "6       50441857206                    NA                    NA    0\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "library(\"rjson\")\n",
    "s2 <- get.Comtrade(r=\"842\", p=\"124,484\", fmt=\"csv\")\n",
    "s2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2.3. データの繰り返し取得\n",
    "上記（2.1.か2.2.）を`for`ループなどで繰り返す."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.3.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}