Mímisbrunnr知恵の泉

← 確率過程 一覧

🎓 レベル:標準 | 重要度:A(必須) 📎 前提:確率過程とは有限次元分布とコルモゴロフの拡張定理

要点(BLUF)

概念

過程を解析可能にする2つの「うまい構造」が定常性と独立増分です。定常性は「いつ観測を始めても統計的性質が同じ」。気温の日内変動のように平均が動くものは非定常、無限に続く安定した揺らぎは定常に近い。独立増分は「過去の変化と未来の変化が無関係」。ランダムウォークの1歩1歩が独立なのが典型です。

数式による定式化

強定常(strict stationarity):任意の時点列と任意のシフト hh について

(Xt1+h,,Xtn+h)=d(Xt1,,Xtn)(X_{t_1+h},\dots,X_{t_n+h}) \overset{d}{=} (X_{t_1},\dots,X_{t_n})

つまり有限次元分布が hh に依らない。弱定常(共分散定常):2次モーメントが有限で

E[Xt]=μ (一定),Cov(Xt,Xt+k)=γ(k) (ラグ k のみに依存)\mathbb{E}[X_t] = \mu \ (\text{一定}), \qquad \mathrm{Cov}(X_t, X_{t+k}) = \gamma(k) \ (\text{ラグ }k \text{ のみに依存})

独立増分:任意の t0<t1<<tnt_0 < t_1 < \dots < t_n について、増分

Xt1Xt0, Xt2Xt1, , XtnXtn1X_{t_1}-X_{t_0},\ X_{t_2}-X_{t_1},\ \dots,\ X_{t_n}-X_{t_{n-1}}

が互いに独立。さらに増分の分布がシフト hh で不変なら定常増分といいます。

直観

具体例

弱定常の例として AR(1) 過程 Xt=ϕXt1+εtX_t = \phi X_{t-1} + \varepsilon_tϕ<1|\phi|<1)。理論上、分散は σ2/(1ϕ2)\sigma^2/(1-\phi^2)、ラグ kk の自己相関は ϕk\phi^k です。独立増分の例として、ランダムウォークの重ならない増分の無相関を確認します。

import numpy as np
rng = np.random.default_rng(7)
phi, sigma, T = 0.7, 1.0, 120000
x = np.zeros(T)
for t in range(1, T):
    x[t] = phi*x[t-1] + sigma*rng.standard_normal()
x = x[2000:]                                   # バーンイン除去
var_th = sigma**2/(1-phi**2)
print(f"弱定常AR(1): 標本平均={x.mean():+.3f}(理論0) 標本分散={x.var():.3f}(理論{var_th:.3f})")
for k in [1, 2, 3]:
    ac = np.corrcoef(x[:-k], x[k:])[0, 1]
    print(f"  ラグ{k}自己相関={ac:+.3f} (理論 phi^k={phi**k:.3f})")
# 独立増分(ランダムウォーク)
W = np.cumsum(rng.standard_normal((40000, 40)), axis=1)
inc1 = W[:, 9] - W[:, 4]       # 区間(5,10]
inc2 = W[:, 29] - W[:, 19]     # 区間(20,30] 重ならない
print(f"重ならない増分の相関={np.corrcoef(inc1, inc2)[0,1]:+.4f} (独立なら0)")
# 弱定常AR(1): 標本平均=+0.000(理論0) 標本分散=1.967(理論1.961)
#   ラグ1自己相関=+0.702 (理論 phi^k=0.700)
#   ラグ2自己相関=+0.494 (理論 phi^k=0.490)
#   ラグ3自己相関=+0.348 (理論 phi^k=0.343)
# 重ならない増分の相関=-0.0022 (独立なら0)

AR(1) の分散と自己相関が理論値と一致(弱定常)。ランダムウォークの重ならない増分の相関はほぼ0(独立増分)。

他過程との関係

数式の直観的意味

独立増分があると、過程は「独立な小さな変化の和」になり、和の分散=分散の和(独立性より)から Var(Xt)t\mathrm{Var}(X_t)\propto t が出ます。さらに増分が定常なら中心極限定理が効いて、長い時間スケールでは増分の和が正規分布に近づく — これがブラウン運動が普遍的に現れる理由(自己相関と過程の特徴づけとあわせて)です。

⚠️ よくある誤解

対応シミュレーション

本文コードを再実行すれば AR(1) の弱定常性と増分の無相関が再現します。ϕ\phi を1に近づけると自己相関の減衰が遅くなり、定常性が崩れる境界(ϕ=1\phi=1 で単位根=ランダムウォーク)が観察できます。

関連