内閣府の消費者動向調査

URLは、http://www.esri.cao.go.jp/jp/stat/shouhi/shouhi.html

中の長期推移のデータをみたけど、元号表記での全角混じりのデータ、、、、書きこむ人も半角を全角に直してそうだ、、、

面倒なので、データ部分だけ持ってきて、dataframeにする

  1: str(con.bhv)

  2: 'data.frame': 86 obs. of 5 variables:

   3: $ index : num 42.7 45.4 48.3 44.9 48.7 49.2 46.1 47.7 48 44 ...

   4: $ life : num 42.4 44.4 47.1 43.6 47.1 47.5 45.1 46.6 46.4 43 ...

   5: $ income : num 40.1 42.1 45.9 41.9 46 46.6 43.3 46.1 46.2 41.8 ...

   6: $ employemnt: num 41.3 46.4 49.6 45.6 50.2 50.7 46.9 48.2 49.1 44 ...

   7: $ spending : num 47.1 48.5 50.7 48.3 51.6 51.9 49.2 50 50.1 47.3 ...

xtsオブジェクトにする

日付を作る。

   1: dates <- paste(rep(2004:2011, each=12), 1:12, “01”, sep=”-“)

   2: head(dates)

   3: [1] "2004-1-01" "2004-2-01" "2004-3-01" "2004-4-01" "2004-5-01" "2004-6-01"

上下を引く

1: dates <- dates[c(-1,-2)]

2: dates <-　dates[1:(length(dates)-8)]

3: ret <- as.xts(con.bhv,　order.by=as.Date(dates))

4: head(ret)

5: index life income employemnt spending

6: 2004-03-01 42.7 42.4　40.1 41.3　　47.1

7: 2004-04-01 45.4 44.4 42.1 46.4　48.5

8: 2004-05-01 48.3 47.1 45.9 49.6　50.7

9: 2004-06-01 44.9 43.6 41.9 45.6　48.3

10: 2004-07-01 48.7 47.1 46.0 50.2　51.6

11: 2004-08-01 49.2 47.5 46.6 50.7　51.9

実は、xtsオブジェクトを作ってみたかっただけ、、、

data.frameに戻して、グラフ化。

   1: ret <- data.frame(date=index(ret), ret)

   2:  R> rownames(ret) <- NULL

   3: ret1 <- melt(ret, id="date")

   4: ggplot(ret1, aes(date, value, colour=variable))+geom_line()

雇用の振れ幅が大きいのが分かる。



２００９年以降は、インデックス（態度指数を抜く)

   1: ggplot(ret1[ret1$date>"2008-12-30",], aes(date, value, colour=variable))+geom_line()

よく分からない。
２０１１年３月以降に落ち込みがあるくらい。急激と言えるとは思うけど。。。

4つの指標を比較しないといけないけど、よく分からないので、また今度。