Rで雨温図(Hythergraph)を書く（その2）

前回からの続き。

今回は気象データを読み込んで雨温図(Hythergraph)を書くことにする。

まず，気象データであるが，ASEAN諸国の首都の月別最高気温，最低気温，雨量をまとめたcsvファイルを作ってみた：

「aseanclimate.csv」をダウンロード

データの出典として，カンボジアやラオスに関してはそれぞれの国の統計年鑑，その他の国については世界気象機関(WMO)のデータベースを用いた。

スクリプトは次の通りで，前回との違いは，csvファイルの読み込みを行うところである。

あと，ブルネイではなく，カンボジアのプノンペンの雨温図を書くことにする。

ちなみに，小生の場合，csvファイルはcドライブ直下に作った"rtemp"というフォルダーに置いている："c:\\rtemp\\aseanclimate.csv"

rm(list=ls(all=TRUE)) # # # csvデータの読み込み d <- read.csv("c:\\rtemp\\aseanclimate.csv") attach(d) names(d) # # データフレームを雨量，月平均最高気温，月平均最低気温の3つに分割 dtmax <- subset(d, Item == "tempmax") dtmin <- subset(d, Item == "tempmin") drain <- subset(d, Item == "rainfall") # # 棒グラフの横軸タイトル＝ラベル（月名）の定義 xlabel <- c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec") # # 対象国の指定 countryname <- "Cambodia" dtmp <- subset(dtmax, Country == countryname) cityname <- data.matrix(dtmp$City) # # 雨量，月平均最高気温，月平均最低気温のデータを代入 tmax <- data.matrix(subset(dtmax, Country == countryname, xlabel)) tmin <- data.matrix(subset(dtmin, Country == countryname, xlabel)) rain <- data.matrix(subset(drain, Country == countryname, xlabel)) # # # 軸表示用の余白の定義 # ベクトルの数値は余白を表し，単位は"行"。底辺，左側，上側，右側の順で定義 # フォントのサイズはpsで指定 par(mar=c(4, 4, 4, 4), ps = 9) # # 棒グラフを表示 # 戻り値は、各バーのx座標位置 barx <- barplot(rain, xlim=c(0, 14), ylim=c(0, 800), yaxp = c(0, 800, 8), names.arg=xlabel, las=3) # # グラフの上書きを指定し， # ラベルや座標軸を非表示のまま月平均最高気温をプロット par(new=T) plot(barx, tmax, xlim=c(0, 14), ylim=c(0, 40), type="o", pch=1, axes=F, xlab="", ylab="") # # グラフの上書きを指定し， # ラベルや座標軸を非表示のまま月平均最低気温をプロット # 最後に折れ線グラフのy軸を表示 par(new=T) plot(barx, tmin, xlim=c(0, 14), ylim=c(0, 40), type="o", pch=19, axes=F, xlab="", ylab="") gtitle <- paste(cityname, ", ", countryname) mtext(gtitle, side=1, line=2) mtext("Mean Total Rainfall [mm]",side=2, line=2) mtext("Mean Daily Max and Min Temp. [C]",side=4, line=2) axis(side=4, yaxp = c(0, 40, 8)) box()

こんな図ができるはず。

今回，csvファイルの読み込みを行ったのは，複数都市の雨温図を並べて書くという作業を行うための伏線である。

次回（その3）では4都市分の雨温図を書くことにする。