景気ウオッチャー調査のプロットスクリプト

上記の記事に用いたRのプログラムは以下の通り

watch <- read.csv("keikiWatch.csv")

sumTs <- ts(watch$sum, start=c(2000,1), freq=12)

houseTs <- ts(watch$house, start=c(2000,1), freq=12)

firmTs <- ts(watch$firm, start=c(2000,1), freq=12)

empTs <- ts(watch$emp, start=c(2000,1), freq=12)
plot(sumTs, lwd=3, col="black", ylim=c(25, 65),

main="景気の現状判断DI：景気ウオッチャー（内閣府）",

 ylab="D I", xlab="2000年1月-2007年12月")

lines(houseTs, lwd=3, lty=2, col="blue")

lines(firmTs, lwd=3, lty=3, col="red")

lines(empTs, lwd=3, lty=4, col="green")

abline(v=2006); abline(v=2007); abline(h=50)

legend(locator(1), legend=c("合計", "家計", "企業","雇用"),

lty=1:4, lwd=3, col=c("black", "blue", "red", "green"))