acpass

モデリング公式

ひとことで言うと

AR(1)の予測は『平均回帰』が本質です。kk 期先の点予測は Y^t+k=ϕk(Ytμ)+μ\hat Y_{t+k}=\phi^k(Y_t-\mu)+\mu で、ϕk\phi^k00 に向かうため遠い将来ほど無条件平均 μ\mu に吸い寄せられます。同時に予測誤差の分散は増えていき、最終的に無条件分散 γ0\gamma_0 に収束する——遠い未来ほど『平均に近いが不確実』になるのです。

AR(1)でφ=0.8,μ=0,σ²=0.36,現在Y_t=2のk期先予測。点予測がφ^k・Y_tで0へ平均回帰する(緑線)。予測誤差帯(薄緑)がk=1で±0.6からk→∞で±1(無条件分散γ_0=1の平方根)へ広がる点予測は平均μ=0へ回帰, 帯は無条件分散γ_0へ拡大Y_t=2μ=0先k期

AR(1) ϕ=0.8,μ=0,σ2=0.36\phi=0.8,\mu=0,\sigma^2=0.36、現在 Yt=2Y_t=2 の予測。点予測 Y^t+k=ϕkYt=1.6,1.28,\hat Y_{t+k}=\phi^k Y_t=1.6,1.28,\dots は平均 00 へ回帰(緑線)。予測誤差の帯(薄緑)は ±0.6±1\pm0.6\to\pm1 と無条件分散 γ0=1\gamma_0=1 へ広がる。

数式で表すと

Y^t+k=ϕkYt, Var=σ21ϕ2k1ϕ2\hat Y_{t+k}=\phi^k Y_t,\ \mathrm{Var}=\sigma^2\tfrac{1-\phi^{2k}}{1-\phi^2}

AR(1)の kk 期先予測は ϕk\phi^k 倍に減衰。予測誤差分散は徐々に増える。

AR(1) モデル Ytμ=ϕ(Yt1μ)+εtY_t-\mu=\phi(Y_{t-1}-\mu)+\varepsilon_tϕ<1|\phi|<1 で定常、εt\varepsilon_t は分散 σ2\sigma^2 のホワイトノイズ)に基づく将来予測を考えます。最小平均二乗誤差の意味で最適な点予測は条件付き期待値で、AR(1) では Y^t+kt=ϕk(Ytμ)+μ\hat Y_{t+k\mid t}=\phi^k(Y_t-\mu)+\mu となります(μ=0\mu=0 なら Y^t+k=ϕkYt\hat Y_{t+k}=\phi^k Y_t)。漸化式を繰り返し代入すれば導けます:Y^t+1=μ+ϕ(Ytμ)\hat Y_{t+1}=\mu+\phi(Y_t-\mu)Y^t+2=μ+ϕ2(Ytμ)\hat Y_{t+2}=\mu+\phi^2(Y_t-\mu)、…と ϕ\phi が1回ずつ掛かっていくためです。ϕ<1|\phi|<1 なので ϕk0\phi^k\to0、したがって kk\to\infty で点予測は無条件平均 μ\mu に収束します。これが『平均回帰』——遠い将来ほど予測は現在値の情報を忘れ、長期平均に吸い寄せられるという AR(1) 予測の本質です。 次に予測の不確実性です。kk 期先の予測誤差分散は Var(予測誤差)=σ21ϕ2k1ϕ2\mathrm{Var}(\text{予測誤差})=\sigma^2\dfrac{1-\phi^{2k}}{1-\phi^2} です。k=1k=1 では σ2\sigma^2 ですが、kk が増えるにつれて単調に増え、kk\to\inftyσ2/(1ϕ2)=γ0\sigma^2/(1-\phi^2)=\gamma_0 つまり無条件分散に収束します。直感的には、遠い未来の値については現在値の手がかりがほぼ効かなくなり、予測の不確実性が『何も知らないときの分散(無条件分散)』にまで膨らむということです。点予測は平均 μ\mu へ、予測誤差分散は無条件分散 γ0\gamma_0 へ——どちらも『無条件の世界』に向かって収束するのが対になっています。 具体例:ϕ=0.8,μ=0,σ2=0.36\phi=0.8,\mu=0,\sigma^2=0.36γ0=0.36/(10.64)=1\gamma_0=0.36/(1-0.64)=1。現在 Yt=2Y_t=2。 点予測: k=1k=11.61.6k=2k=21.281.28k=5k=50.655\approx0.655kk\to\infty00(平均回帰)。 予測誤差分散: σk2=10.64k\sigma_k^2=1-0.64^kk=1k=10.360.36k=2k=20.59040.5904kk\to\infty1=γ01=\gamma_0

試験に出る性質

k期先の点予測

Y^t+kt=ϕk(Ytμ)+μ\hat Y_{t+k\mid t}=\phi^k(Y_t-\mu)+\muμ=0\mu=0 なら ϕkYt\phi^k Y_t)。漸化式の繰り返し代入で ϕ\phikk 回掛かる。

平均回帰

ϕ<1|\phi|<1 より ϕk0\phi^k\to0kk\to\infty で点予測は無条件平均 μ\mu に収束。遠い将来ほど現在値を忘れる。

予測誤差分散

σ2(1ϕ2k)/(1ϕ2)\sigma^2(1-\phi^{2k})/(1-\phi^2)k=1k=1 では σ2\sigma^2kk とともに単調増加。

無条件分散へ収束

kk\to\inftyσ2/(1ϕ2)=γ0\sigma^2/(1-\phi^2)=\gamma_0。不確実性が無条件分散に収束する。

二重の収束

点予測は平均 μ\mu へ、予測誤差分散は γ0\gamma_0 へ。遠い未来は『平均に近いが最も不確実』。

例で見る

ϕ=0.8,μ=0,σ2=0.36\phi=0.8,\mu=0,\sigma^2=0.36γ0=1\gamma_0=1。現在 Yt=2Y_t=2。 点予測: k=1:1.6k=1{:}1.6k=2:1.28k=2{:}1.28k:0k\to\infty{:}0(平均回帰)。 予測誤差分散: k=1:0.36k=1{:}0.36k=2:0.5904k=2{:}0.5904k:1=γ0k\to\infty{:}1=\gamma_0

つまずきポイント

  • 点予測が現在値 YtY_t をいつまでも保つと思う(ϕk0\phi^k\to0 で無条件平均 μ\mu に回帰。長期予測は現在値の情報を忘れる)
  • 予測誤差分散が無限に発散すると思う(定常 AR では γ0\gamma_0 に収束する。発散するのは単位根 ϕ=1\phi=1 の非定常な場合)
  • 1期先分散 σ2\sigma^2 をk期先でも使う(k期先は σ2(1ϕ2k)/(1ϕ2)\sigma^2(1-\phi^{2k})/(1-\phi^2)kk とともに増える)

定着クイズ

AR(1) μ=0\mu=0kk 期先点予測は?

ϕ=0.8,μ=0,Yt=2\phi=0.8,\mu=0,Y_t=2 のとき k=2k=2 期先の点予測は?

AR(1) で kk\to\infty のとき予測誤差分散は?

この用語を扱う問題(2