東京に棲む日々

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

2014-11-15から1日間の記事一覧

状態空間モデル - 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]…