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

前回，前々回に続き，統計解析フリーソフトRで雨温図(Hythergraph)を書くというシリーズの3回目。

今回は4都市分の雨温図を並べて書いてみる。

前回，気象データとして，ASEAN各国の首都の月別最高気温・最低気温・雨量をまとめたcsvファイルを作ってみたのでそれを再び使用する：

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

前回も述べたが，データの出典は，カンボジアやラオスに関してはそれぞれの国の統計年鑑，その他の国については世界気象機関(WMO)のデータベースである。

スクリプトは以下の通りであるが，これまでとの大きな違いは，こんな感じ：

• ASEAN諸国のうち，カンボジア，ラオス，タイ，ベトナムの首都の雨温図を書くことにした
• 関数par()を使って画面を4分割し，4都市分がまとめてみられるようにした
• for文を使って描画作業を繰り返すことにした

#おまじない 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 <- c("Cambodia", "Lao PDR", "Thailand", "Vietnam") # # 画面を2×2に分割し，4つのグラフを表示できるようにする # また，軸表示用の余白を定義する： # ベクトルの数値は余白を表し，単位は"行"。底辺，左側，上側，右側の順で定義 # フォントのサイズはpsで指定 par(mfrow=c(2,2), mar=c(4, 4, 4, 4), ps = 9) # ##### # グラフの作成 # 対象国群の数だけ繰り返す ##### for (i in 1:length(countryname)) { # dtmp <- subset(dtmax, Country == countryname[i]) cityname <- data.matrix(dtmp$City) # # 雨量，月平均最高気温，月平均最低気温のデータを代入 tmax <- data.matrix(subset(dtmax, Country == countryname[i], xlabel)) tmin <- data.matrix(subset(dtmin, Country == countryname[i], xlabel)) rain <- data.matrix(subset(drain, Country == countryname[i], xlabel)) # # # 棒グラフを表示 # 戻り値は、各バーの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[i]) 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() # ##### # グラフの作成 # 対象国群の数だけ繰り返す ##### }

こんな図ができるはず。

今回のスクリプトを応用すれば，一画面に何都市分の雨温図でも表示できる。もちろん図が小さくなるので見えにくくなるが。