2014年6月11日、消費税増税後としては2回目に当たる2014年5月分の国内企業物価指数が公表されました。

2014年6月11日、消費税増税後としては2回目に当たる2014年5月分の国内企業物価指数が公表されました。

・総平均の指数値、106.1は民主党政権以降の最高値。
・2014年5月の前月比、増税後である2014年4月との前月比はプラス0.28%。
・工業製品において2014年5月の前月比の最も高い類別は非鉄金属のプラス0.83%。最も低いそれは情報通信機器のマイナス0.78%。


なお同指数に関する基本統計量等は以下の通り。グラフは2004年5月から1直近まで、10年間の時系列推移。折線の緑部分が民主党政権期間、青部分は自民党への政権交代以降。
140611_01

> library(tseries)
> library(ggplot2)

> dataset$date < - as.Date(dataset$date)

> jimin < - dataset[which(dataset$date>=as.Date("2013-1-1") & dataset$date< = as.Date("2014-5-1")),] 

> minshu < - dataset[which(dataset$date>=as.Date("2009-9-1") & dataset$date< = as.Date("2012-12-1")),] 

> p < - ggplot(NULL)

> p < - p+geom_line(data=dataset,aes(date,dgcpi))

> p < - p+geom_line(data=minshu,aes(date,dgcpi),colour="3")

> p < - p+geom_line(data=jimin,aes(date,dgcpi),colour="4")

> p < - p+labs(x="日付",y="国内企業物価指数")

> p

> quantile(dataset$dgcpi)
   0%   25%   50%   75%  100% 
 95.3  99.7 100.3 101.9 109.3 

> quantile(minshu$dgcpi)
   0%   25%   50%   75%  100% 
 99.3 100.0 100.3 101.1 102.2 

> quantile(jimin$dgcpi)
   0%   25%   50%   75%  100% 
100.5 101.6 102.5 102.8 106.1 

> mean(dataset$dgcpi)
[1] 100.6355

> mean(minshu$dgcpi)
[1] 100.575

> mean(jimin$dgcpi)
[1] 102.5118

> adf.test(dataset$dgcpi,k=0)

        Augmented Dickey-Fuller Test

data:  dataset$dgcpi
Dickey-Fuller = -1.387, Lag order = 0, p-value = 0.8312
alternative hypothesis: stationary

> adf.test(diff(dataset$dgcpi),k=0)

        Augmented Dickey-Fuller Test

data:  diff(dataset$dgcpi)
Dickey-Fuller = -6.3999, Lag order = 0, p-value = 0.01
alternative hypothesis: stationary

 警告メッセージ: 
In adf.test(diff(dataset$dgcpi), k = 0) :
  p-value smaller than printed p-value

> adf.test(minshu$dgcpi,k=0)

        Augmented Dickey-Fuller Test

data:  minshu$dgcpi
Dickey-Fuller = -1.1971, Lag order = 0, p-value = 0.8886
alternative hypothesis: stationary

> adf.test(diff(minshu$dgcpi),k=0)

        Augmented Dickey-Fuller Test

data:  diff(minshu$dgcpi)
Dickey-Fuller = -4.8405, Lag order = 0, p-value = 0.01
alternative hypothesis: stationary

 警告メッセージ: 
In adf.test(diff(minshu$dgcpi), k = 0) :
  p-value smaller than printed p-value

> adf.test(jimin$dgcpi,k=0)

        Augmented Dickey-Fuller Test

data:  jimin$dgcpi
Dickey-Fuller = -1.3829, Lag order = 0, p-value = 0.8075
alternative hypothesis: stationary

> adf.test(diff(jimin$dgcpi),k=0)

        Augmented Dickey-Fuller Test

data:  diff(jimin$dgcpi)
Dickey-Fuller = -4.1813, Lag order = 0, p-value = 0.01693
alternative hypothesis: stationary

> time1 < - seq(1,length(minshu$date))

> length(time1)
[1] 40

> time2 < - seq(1,length(jimin$date))

> length(time2)
[1] 17

> kaiki.minshu < - lm(minshu$dgcpi~time1)

> kaiki.minshu

Call:
lm(formula = minshu$dgcpi ~ time1)

Coefficients:
(Intercept)        time1  
   99.99423      0.02833  

> kaiki.jimin < - lm(jimin$dgcpi~time2)

> kaiki.jimin

Call:
lm(formula = jimin$dgcpi ~ time2)

Coefficients:
(Intercept)        time2  
   100.2044       0.2564  

アプリケーション
R Core Team (2013). R: A language and environment for statistical computing.
R Foundation for Statistical Computing, Vienna, Austria.
URL http://www.R-project.org/.