統計学をRを使って学習するブログ

最初からつまづきました。
サンプルデータがどこにあるのかわからないというあれです。
見つけ出すのは時間の無駄でもあるので、あきらめてDATA_fileの中にあるファイルを使用する事にしました。

ts関数はRjpWikiも見たほうがいいかなと思います。

また、par(mfrow=c(4,2))のところで、下記のようなエラーがでます。

> ts.plot(GDP,type="l",main="GDP")
Error in plot.new() : figure margins too large

これはRStudioにありがちなplotsの領域（図が表示される領域）が小さいためにでるエラーになります。とりあえず、この領域をできるだけ大きくしてみるとエラーがでなくなりました。



練習問題のデータは、サンプルのDATA_fileの中の「IIP1998.csv」、「IIP1998_2.csv」が該当するようでした。

演習問題3-2はレポートの作成と書いてありますが、先進もうかと思います。

参考リンクをみて書くのかなとは思います。

練習問題は、もっとかっこいい書き方あるんじゃないかとは思いつつとりあえずやってみました。

##########練習問題2-1##########

ren21<-read.csv("DATA_file/IIP1998.csv",skip=2,header=TRUE)
head(ren21)
PROD<-ts(ren21$PROD,start=c(1998),frequency=4)
SHIP<-ts(ren21$SHIP,start=c(1998),frequency=4)
STOCK<-ts(ren21$STOCK,start=c(1998),frequency=4)
xlabel<-"年"
ylabel<-"鉱工業指数"
title<-"鉱工業指数傾向"
index<-c("生産指数","出荷指数","在庫指数(期末値)")
plot(PROD,type="o",xlab=xlabel,ylab=ylabel,main=title,col=1,lty=1)
lines(SHIP,type="o",col=2,lty=2)
lines(STOCK,type="o",col=3,lty=3)
legend("bottomleft","(x,y)",legend=index,pch=1,lty=c(1:3),col=c(1:3))
dev.off()

#y軸調整
plot(PROD,type="o",ylim=c(85,125),xlab=xlabel,ylab=ylabel,main=title,
col=1,lty=1)
lines(SHIP,type="o",col=2,lty=2)
lines(STOCK,type="o",col=3,lty=3)
legend("topright","(x,y)",legend=index,pch=1,lty=c(1:3),col=c(1:3))
dev.off()

#ts.plotの使用
ts.plot(PROD,SHIP,STOCK,type="o",xlab=xlabel,ylab=ylabel,main=title,
col=c(1:3),lty=c(1:3))
legend("top","(x,y)",legend=index,pch=1,lty=c(1:3),col=c(1:3))

##########練習問題2-2##########

#読み込みと確認
ren22<-read.csv("DATA_file/IIP1998_2.csv",skip=2,header=TRUE)
head(ren22)

#時系列オブジェクトの登録
CAPTL<-ts(ren22$CAPTL,start=c(1998),frequency=4)
CONST<-ts(ren22$CONST,start=c(1998),frequency=4)
COND<-ts(ren22$COND,start=c(1998),frequency=4)
CONNON<-ts(ren22$CONNON,start=c(1998),frequency=4)

#2002Q1から8期のデータに絞込み（20021Q時点を100とする）
ns1<-(2002-1998)*4+1;ne1<-ns1+8
x1<-CAPTL[ns1:ne1]/CAPTL[ns1]*100
x2<-CONST[ns1:ne1]/CONST[ns1]*100
x3<-COND[ns1:ne1]/COND[ns1]*100
x4<-CONNON[ns1:ne1]/CONNON[ns1]*100

#定数などなど
xlabel2<-"景気の谷(0)からの期数"
ylabel2<-"景気の谷=100の指数"
title2<-"鉱工業指数傾向の景気の谷との局所比較"
z<-0:8
index2<-c("資本財","建設財","耐久消費財","非耐久消費財")

#plot一回目
plot(z,x1,type="o",xlab=xlabel2,ylab=ylabel2,main=title2,
col=1,lty=1)
lines(z,x2,type="o",col=2,lty=2)
lines(z,x3,type="o",col=3,lty=3)
lines(z,x4,type="o",col=4,lty=4)
abline(h=100)

#軸調整
plot(z,x1,type="o",ylim=c(90,115),xlab=xlabel2,ylab=ylabel2,main=title2,
col=1,lty=1)
lines(z,x2,type="o",col=2,lty=2)
lines(z,x3,type="o",col=3,lty=3)
lines(z,x4,type="o",col=4,lty=4)
abline(h=100)
legend("topleft","(x,y)",legend=index2,pch=1,lty=c(1:4),col=c(1:4))

#Windowsメタファイル出力
win.metafile("DATA_file/Rplot32.emf")
plot(z,x1,type="o",ylim=c(90,115),xlab=xlabel2,ylab=ylabel2,main=title2,
col=1,lty=1)
lines(z,x2,type="o",col=2,lty=2)
lines(z,x3,type="o",col=3,lty=3)
lines(z,x4,type="o",col=4,lty=4)
abline(h=100)
legend("topleft","(x,y)",legend=index2,pch=1,lty=c(1:4),col=c(1:4))
dev.off()