R

FRED

plotrix、quantmod、ggplot2、gridExtra パッケージ

(参考)
ASSET MANAGEMENT CONSULTING CORP.

FRED(Federal Reserve Economic Data) には日本に関するデータも多い。
FRED
「japan」で検索
Search Results
2,441 Series(2015/3/1 現在)

ドル円、日経２２５、日銀：総資産残高

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
library(quantmod)
library(ggplot2)
library(reshape2)
library(grid)
library(gridExtra)
#mergeするため、as.yearmon日付を統一する
#ドル円と日経２２５は月平均をとる。na.rm=Tを忘れずに
#ドル円
getSymbols("DEXJPUS", src = "FRED")
DY <- DEXJPUS["1999::",]
DYm<-apply.monthly(DY,mean,na.rm=T)
DY<-as.xts(coredata(DYm),as.yearmon(index(DYm)))
#日経２２５
getSymbols("NIKKEI225",src="FRED")
nikkei<- NIKKEI225["1999::",]
nikkeim<-apply.monthly(nikkei,mean,na.rm=T)
nikkei<-as.xts(coredata(nikkeim),as.yearmon(index(nikkeim)))
#日銀：総資産残高
getSymbols("JPNASSETS",src="FRED")
jpa<- JPNASSETS["1999::",]
JB<-as.xts(coredata(jpa),as.yearmon(index(jpa)))
#mergeする。 all=Tを忘れずに
dat3<-merge(DY,nikkei,JB,all=T)
#関数を定義
xtsplot2<-function(x) { 
library(ggplot2)
library(reshape2)
data<-data.frame(date=index(x),coredata(x))
data.melt <- melt(data,id.vars=c("date"),measure.vars=c(colnames(data[2:ncol(data)])))
ggplot(data.melt, aes(x=as.Date(date),y=value,colour=variable,group=variable)) + geom_line(colour = "blue",size =0.5,linetype=1,alpha = 1)+labs(x="", y="")
}
#グラフ
#png("FRED01.png",width=1000,height=800)
p1<-xtsplot2(dat3[,1])+ labs(title="USD/Yen")
p2<-xtsplot2(dat3[,2])+ labs(title="NIKKEI225")
p3<-xtsplot2(dat3[,3])+ labs(title="Bank of Japan: Total Assets for Japan")
grid.arrange(p1,p2,p3,ncol=1,nrow=3)
#dev.off()

plotrix::twoord.plot 関数を使う

1
2
3
4
5
6
#2軸のグラフ
library(plotrix)
#png("FRED02.png",width=1000,height=800)
twoord.plot(index(dat3),coredata(dat3[,1]),index(dat3),coredata(dat3[,2]),type="l",xlab="",
ylab="USD/Yen",rylab="Nikkei225",main="USD/Yen & Nikkei225",do.first="plot_bg();grid(col=\"white\",lty=1)")
#dev.off()