東京に棲む日々

データ分析、統計、ITを勉強中。未だ世に出ず。

2014-11-01から1ヶ月間の記事一覧

状態空間モデル - R dlm (4)

ARIMA(p,d,q)モデルを動的線形モデルで表現してみる。 ARIMAモデルの状態空間表現は一意に決定されないらしい、ここでは文献に従った形式を見て行く。 I(d)は差分を表し、通常のモデリングアプローチでは、データがトレンドを持ち非定常な場合は、その差分を…

状態空間モデル - R dlm (3)

今回は季節要素モデル。 例えば、データが四半期(周期=4)だった場合を考える。 Y[t-1] = α1 + v[t-1] Y[t] = α2 + v[t] Y[t+1] = α3 + v[t+1] Y[t+2] = α4 + v[t+2] α1+α2+α3+α4 = w[t], w[t] ~ N(0, w_seas) α1,2,3,4が各四半期成分となり、和が平均0で分…

状態空間モデル - R dlm (2)

前回より続く。 線形成長モデルは以下に定義される。 Y[t] = μ[t] + v[t], v[t]~N(0, V) μ[t] = μ[t-1] + β[t-1] + w1[t], w1[t]~N(0, W1) β[t] = β[t-1] + w2[t], w2[t]~N(0, W2) 動的線形モデルで表現する場合、各要素は以下になる。 t=0, θ[0]~N(m[0]…