ggplot2

日銀買入額(2010-12-15～2016-06-10)

ggplot2,xts,Quandl,gdata パッケージ

(過去の記事)
日銀買入額
日銀買入額(2010-12-15～2015-07-09)
日銀買入額(2010-12-15～2016-02-19)

(使用したデータ)

(作成したデータファイル)
データファイル:2010_2015 (コードの中で読み込む)
d2010_2015.RData

(コードの中で読み込むデータ:2016年分)
2016年3月以前の買入結果
2016年4月以降の買入結果

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
39
40
41
42
43
44
45
library(xts)
library(gdata)
library(ggplot2)
library(grid) 
library(Quandl)
#
#作成したデータファイル読み込み
load(url("http://statrstart.github.io/data/d2010_2015.RData"))
#
#日銀のサイトからデータを取り込む
#2016年3月以前の買入結果
temp <- tempfile()
download.file("http://www3.boj.or.jp/market/jp/pastetfreit.zip",temp)
con <- unzip(temp, "2016.xls")
d201601_03 <- read.xls(con,header=F,skip=8)
unlink(temp)
#2016年4月以降の買入結果
temp <- tempfile()
download.file("http://www3.boj.or.jp/market/jp/etfreit.zip",temp)
con <- unzip(temp, "2016.xls")
d201604 <- read.xls(con,header=F,skip=8)
unlink(temp)
#
d2016<-rbind(d201601_03[,2:4],d201604[,2:4]) 
names(d2016)<-c("date","ETF","REIT")
d2016$ETF[is.na(d2016$ETF)]<-0
d2016$REIT[is.na(d2016$REIT)]<-0
#rbind
d2010_2016<-rbind(d2010_2015,d2016)
#
x <- read.zoo(d2010_2016)
#ETF累計
a<-ggplot(data = fortify(cumsum(x[,1]), melt = TRUE),aes(x = Index, y = Value) ) +
geom_bar(stat="identity",position=position_dodge(),fill=rgb(1,0,0,alpha=0.7)) +
labs(title=paste("日銀買い入れ:累計(ETF)  ",d2010_2016[1,"date"],"_",d2010_2016[nrow(d2010_2016),"date"]), x="", y="")
#REIT累計
b<-ggplot(data = fortify(cumsum(x[,2]), melt = TRUE),aes(x = Index, y = Value) ) +
geom_bar(stat="identity",position=position_dodge(),fill=rgb(0,0,1,alpha=0.7)) +
labs(title=paste("日銀買い入れ:累計(REIT)  ",d2010_2016[1,"date"],"_",d2010_2016[nrow(d2010_2016),"date"]), x="", y="")
#png("boj20160610.png")
grid.newpage() #空の画面を作る
pushViewport(viewport(layout=grid.layout(2, 1))) #画面を区切る
print(a, vp=viewport(layout.pos.row=1)) 
print(b, vp=viewport(layout.pos.row=2)) 
#dev.off()

日経平均株価データの取得にQuandlパッケージを用います。

買入日の株価指数の終値でETF(日経平均株価)の買入を行ったと仮定する

1
2
3
4
5
#日経平均株価
N225=Quandl("YAHOO/INDEX_N225", collapse="daily", start_date="2010-12-15", type="xts")
#買入日の株価指数の終値でETF(日経平均株価)の買入を行ったと仮定する
dat<-merge(d2010_2016,data.frame(date=index(N225),coredata(N225[,"Close"])))
sum(dat$ETF)/sum(dat$ETF/dat$Close)

14621.61

日経平均株価と日銀買い入れ:累計(ETF)

1
2
3
4
5
6
7
8
9
10
#png("boj20160610_02.png",width=800,height=640)
par(mar=c(5,7, 3,7))  
plot.zoo(merge(N225[,"Close"],cumsum(x[,1]))[,2],type="h",lend=1,col="pink",xlab="",ylab="",yaxt = "n")
mtext("日銀買い入れ:累計(ETF)",side = 4, line =4)  # 右の縦軸のラベル
axis(4,las=1)                        # 右の縦軸を表示
par(new=T)
plot.zoo(merge(N225[,"Close"],cumsum(x[,1]))[,1],xlab="",ylab="",xaxt = "n",las=1)
mtext("N225(終値)",side = 2, line = 4) 
title("N225(終値) & 日銀買い入れ:累計(ETF)")
#dev.off()

(参考)
yahoo ファイナンス : 日本銀行

日本銀行が保有する国債の銘柄別残高
日本銀行による国庫短期証券の銘柄別買入額

2016年5月31日現在 （単位、億円）

日本銀行が保有する国債の残高
3,119,139
日本銀行による国庫短期証券買入額
403,023
合計
3,522,162億円(約３５２兆円)