2軸グラフの作成

AllData02

#各種条件設定
FirstDate<-as.Date("2007/9/1")
LastDate<-as.Date("2014/8/1")
period<-"months"
path01<-paste("C:/Users/",username,"/Desktop/R_Folder/",sep="")
path02<-paste("C:/Users/",username,"/Desktop/R_Graph/",sep="")
setwd(path01)
#パッケージ読込
library(urca)
library(tseries)
library(forecast)
#分析開始
AllData01<-data.frame(date=seq(FirstDate,LastDate,by=period),dummy=1)
AllData02<-AllData01
TSeriesDataFileList<-dir(path01)
for(iii in 1:2)#本スクリプトサンプルは系列が2つの場合
{
dataset<-read.table(file=TSeriesDataFileList[iii],header=TRUE,sep=",")
#日付列書式の統一
dataset[,1]<-as.Date(dataset[,1])
colnames(dataset)[1]<-"date"
dataset<-subset(dataset,FirstDate<=dataset[,1])
dataset<-subset(dataset,LastDate>=dataset[,1])
for(ccc in 2:ncol(dataset))
{
#以下は非数値データ行の削除及びデータの実数化。オリジナルデータ書式、分析目的に合わせて適宜追加、変更。
dataset<-subset(dataset,dataset[,ccc]!=".")
dataset<-subset(dataset,is.na(dataset[,ccc])==F)
dataset[,ccc]<-as.double(as.character(dataset[,ccc]))
}
AllData01<-merge(AllData01,dataset,by="date",all=T,sort=T) #欠損はNAとして処理
AllData02<-merge(AllData02,dataset,by="date",sort=T) #全系列に数値データが存在する日付のみ
#グラフ作成
if(iii==2){
setwd(path02)
png(file="AllData02.png",width=900,height=900)
par(mar=c(5,5,5,5),ps=22,cex.axis=1,cex.lab=1,cex.main=1,cex.sub=1)
date<-AllData02[,1]
plot(date,AllData02[,3],type="l",xaxt="n",ylab=colnames(AllData02)[3],col=2,main=paste("対象期間:",format(FirstDate,"%Y/%b/%d"),"-",format(LastDate,"%Y/%b/%d"),sep=""))
par(new=T)
plot(date,AllData02[,4],type="l",xaxt="n",yaxt="n",ylab="",col=4,xlab="date")
axis(side=4)
mtext(side=4,line=3,text=colnames(AllData02)[4])
axis.Date(side=1,at=seq(FirstDate,LastDate,period),format="%Y%m")
legend(
"topleft",
legend=c(colnames(AllData02)[3],colnames(AllData02)[4]),
lty=c(1,1),
col=c(2,4),
bty="n"
)
dev.off()
}
}

参考文献
福地純一郎、伊藤有希(2011).『Rによる計量経済分析』.朝倉書店.186pp.

アプリケーション
R Core Team (2013). R: A language and environment for statistical computing.
R Foundation for Statistical Computing, Vienna, Austria.
URL http://www.R-project.org/.