R

日銀買入額

zoo,xts,quantmod,gdata パッケージ

日銀ETF,J-REIT買入
2014年4月以降の買入結果,2014年3月以前の買入結果
(エクセル形式 + 圧縮)
2014年3月以前の買入結果を解凍し、すべてのxlsファイルを作業ディレクトリへ
Rで読み込み、結合、d2010_201403.dat　保存。
xlsファイルを作業ディレクトリから移動。 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
#library(gdata)
#d2010<-read.xls("2010.xls",sheet=1,header=F,skip=8)
#d2010<-d2010[,2:4]
#d2011<-read.xls("2011.xls",sheet=1,header=F,skip=8)
#d2011<-d2011[,2:4]
#d2012<-read.xls("2012.xls",sheet=1,header=F,skip=8)
#d2012<-d2012[,2:4]
#d2013<-read.xls("2013.xls",sheet=1,header=F,skip=8)
#d2013<-d2013[,2:4]
#d2014<-read.xls("2014.xls",sheet=1,header=F,skip=8)
#d2014<-d2014[,2:4]
#d2010_201403<-rbind(d2010,d2011,d2012,d2013,d2014)
#names(d2010_201403)<-c("date","ETF","REIT")
#save(d2010_201403, file="d2010_201403.dat")

作成したデータファイル　
d2010_201403.dat

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
library(gdata)
library(xts)
library(quantmod)
load("d2010_201403.dat")
#(Rのコマンドで行う)2014年4月以降の買入結果を解凍して2014.xlsを作業ディレクトリへ
#d2014<-read.xls("2014.xls",sheet=1,header=F,skip=8)
#Rのコマンドでzipファイルをダウンロードし、解凍する。
temp <- tempfile()
download.file("http://www3.boj.or.jp/market/jp/etfreit.zip",temp)
con <- unzip(temp, "2014.xls")
d2014 <- read.xls(con,header=F,skip=8)
unlink(temp)
d2014<-d2014[,2:4]
names(d2014)<-c("date","ETF","REIT")
#tail(d2010_201403);head(d2014)
d<-rbind(d2010_201403,d2014)
d.xts <- as.xts(zoo(d[,-1]),as.POSIXct(d[,1]))
#head(d.xts)
#png("nichi01.png", width=1000, height=800)
par(mfrow=c(2,1),mar = c(4, 5, 3, 3))
plot(d.xts[,1],type="h",main="日本銀行 ETF 買入額　単位：億円")
plot(d.xts[,2],type="h",main="日本銀行 J-REIT 買入額　単位：億円")
par(mfrow=c(1,1))
#dev.off()

日経平均株価 & 日本銀行 ETF 買入額 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
getSymbols("^N225",from="2010-10-15")
#head(N225);tail(N225)
N225["2014-12-01",c(1:6)]<-c(17475.10,17649.02,17474.27,17590.10,0,17590.10)
#chartSeries(N225)
#tail(N225$N225.Close)
#グラフにする期間を指定
start<-"2010-10-15" ; end<-"2015-01-09"
r<-as.POSIXct(c(start,end))
t<-paste(start,"::",end,sep="")
#png("nichi02.png", width=1000, height=800)
par(oma = c(0,3, 0, 3))
plot(N225$N225.Close[t],main=paste("日経平均株価 & 日本銀行 ETF 買入額(",start,"～",end,")"),xlim=c(r[1],r[2]),ylim=c(5000,19000),las=1)
mtext("日経平均株価",side = 2, line = 4) 
par(new=T)
plot(d.xts[t,1],type="h",xlim=c(r[1],r[2]),ylim=c(0,1000),xaxt="n",xlab="",yaxt="n",ylab="",main="",col="red")
axis(4,las=1)  
mtext("ETF 買入額(億円)",side = 4, line = 3) 
#dev.off()

1991_2014までのtopixのデータは
株価時系列データ取得
を利用してyahooファイナンスから入手しました。

topix1991_2014.dat

（2015年のデータは最新のグラフを作るときに入手）

【TOPIX】株価指数 時系列データ 日足
2015 を選択 -> csvダウンロード

ダウンロードしたcsvファイルを作業フォルダに置いて

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
library(gdata)
library(xts)
library(quantmod)
topix2015 <- read.csv("/home/user/R/work/I102_2015.csv", header=F,skip=2)
names(topix2015)<-c("Date","Open","High","Low","Close")
#topix1991_2014.datを読み込んで結合
load("topix1991_2014.dat")
topix<-rbind(topix2015,topix1991_2014)
topix.xts <- as.xts(zoo(topix[,-1]), as.POSIXct(topix[,1]))
#chartSeries(topix.xts)
#plot(topix.xts[,4])
#topix & 日本銀行 ETF 買入額
#日本銀行 ETF 買入額の読み込みは上と同じ
d2014<-read.xls("2014.xls",sheet=1,header=F,skip=8)
d2014<-d2014[,2:4]
names(d2014)<-c("date","ETF","REIT")
load("d2010_201403.dat")
d<-rbind(d2010_201403,d2014)
d.xts <- as.xts(zoo(d[,-1]),as.POSIXct(d[,1]))
#グラフにする期間を指定
start<-"2010-10-15" ; end<-"2015-01-09"
r<-as.POSIXct(c(start,end))
t<-paste(start,"::",end,sep="")
#png("nichi03.png", width=1000, height=800)
par(oma = c(0,3, 0, 3))
plot(topix.xts$Close[t],main=paste("TOPIX & 日本銀行 ETF 買入額(",start,"～",end,")"),xlim=c(r[1],r[2]),ylim=c(500,1500),las=1)
mtext("TOPIX",side = 2, line = 4) 
par(new=T)
plot(d.xts[t,1],type="h",xlim=c(r[1],r[2]),ylim=c(0,1000),xaxt="n",xlab="",yaxt="n",ylab="",main="",col="red")
axis(4,las=1)  
mtext("ETF 買入額(億円)",side = 4, line = 3) 
#dev.off()