気象庁・過去データの一般的な落とし方とＲでの使い方＆超絶技法

昨日、データサイエンティスト養成講座: 初級編
いってきて、気象庁データ使うときに、おもわずハマったので、
もう一度手順の復習と、そこで教わった超絶技法のご紹介。

---

■気象庁・過去データの一般的な落とし方

＜＜課題＞＞
2015年1月1日から、1月31日までの大阪の最高気温をＣＳＶでダウンロード


＜＜手順＞＞

１．サイトアクセス

過去の気象データ・ダウンロード
http://www.data.jma.go.jp/gmd/risk/obsdl/index.php


２．地点を選ぶ
大阪をクリック

大阪をクリック

→チェックが付き、右側にエリアが選択されるはず


３．項目を選択

「項目を選ぶ」ボタンをクリック

日最高気温にチェック、右側に選択されているのを確認


４．期間を選ぶ
　「期間を選ぶ」ボタンをクリックし、
　希望期間をいれ、
　右側確認


５．ＣＳＶファイルをダウンロード
・もし、「表示オプションを選ぶ」で変えたいことがあれば、
クリックして、修正する。
その語、ＣＳＶファイルをダウンロードボタンをクリック

data.csvとかのファイル名でおちると思う。

---

■使い方

１．Ｒを立ち上げたら、

２．ファイル→ディレクトリの変更で、
　　今data.csvをダウンロードしてきたフォルダを指定

３．以下のコマンドでＶに読み込み


　　　　v<-read.csv("data.csv",skip=5)

　（頭６行がみだし。そこで、５行飛ばす。HEADERのデフォルトはTRUE)

４．２列目だけ取り出す
　　　　u<-v[[2]]

---

■超絶技法

Ｃｈｒｏｍｅで（IEではできない）

http://www.data.jma.go.jp/obd/stats/etrn/view/monthly_s3.php?prec_no=62&block_no=47772&year=&month=&day=&view=a2

のヘッダーから下を選択

これをctrl+Cで、クリップボードにコピー

ここで、Rから

x<-read.delim("clipboard",header=TRUE)

と打ち込み、Ｘをみると、

のように、年と温度がだんだんになっている。

ここで、１月の気温をとってくる方法


m<-nrow(x)
y<-x[seq(2,m,2),1]

おおお～SEQを使って２行目から、２行おきにとってくればよいのね

z<-mean(y)

で平均気温が求まり、

t.test(u,mu=z)

で上述１月の気温差のｔ検定が求まる。