Kalman Filterの予測・innovation・更新
Kalman Filterの予測更新サイクル、gainが不確実性を調整する仕組み、最適性の条件、金融応用の失敗点を解説します。
Mark編集 · 一次情報は記事末尾に掲載
要点
Kalman Filterは線形 Gaussian 状態空間モデルで潜在状態の条件付き平均と不確実性を再帰推定します。まず遷移式で状態と共分散を予測し、新観測と予測を比べます。この予測誤差が innovation です。Kalman gain は予測状態の不確実性と測定不確実性を比べ、どれだけ状態を修正するか決めます。モデル内では厳密でも、誤った動態やノイズ共分散、非線形性、外れ値、未来データでの調整は精密な誤答を作ります。
フィルターは確率的な記述から始まります
各時点で持つのは滑らかな線上の一点だけでなく、状態推定値と誤差共分散です。
線形性は状態と観測を過去状態と攪乱のアフィン関数にし、Gaussian 仮定は条件分布を Gaussian に保ちます。
そのモデルでは条件付き平均が平均二乗誤差を最小にします。モデル外でも有用ですが同じ最適性は主張できません。
予測は平均と不確実性を動かします
遷移行列は以前の状態推定を前方へ写し、既知の制御や入力が予測平均を移動できます。
共分散も伝播し、プロセスノイズ分だけ増えます。速く変わる状態は観測間で不確実になります。
この事前予測は現在観測を使いません。更新と分けることでリアルタイム情報を監査できます。
Innovationが新しい情報を測ります
Innovation は観測値から観測式の予測値を引いた、事前状態が予想しなかった部分です。
その分散は投影した状態不確実性と測定ノイズを合わせます。分散が大きければ大きな生の驚きも普通かもしれません。
仕様が正しければ標準化 innovation は予測不能なノイズに近く、系列パターンは未吸収情報を示します。
Kalman gainが驚きを配分します
Gain は innovation を状態修正へ写し、測定に比べて状態予測が不確実なほど大きくなります。
観測ノイズが大きければ gain は小さく、精密な測定なら驚きの多くが状態修正になります。
一般に固定平滑パラメータではなく、共分散、欠測、時変システム行列に応じて変わります。
共分散更新も答えの一部です
観測後は不確実性が解消され状態共分散が縮み、次の予測でプロセスノイズにより再び広がります。
丸め誤差が非対称や非半正定値行列を作るため、数値的に安定した共分散式が重要です。
フィルター共分散だけの信頼帯はパラメータ不確実性を除く場合があります。遷移やノイズの推定にも誤差があります。
初期化と尤度が経路全体に効きます
初期状態は既知、定常、推定、diffuse にできます。初期フィルター値は選択に強く依存し得ます。
予測誤差と分散から Gaussian 尤度を作り未知パラメータを推定します。見た目の滑らかさだけでノイズを調整しません。
初期化、保留予測、境界推定を比較します。ゼロ近くの分散や一近くの持続性は弱い識別の兆候です。
拡張法は別々の問題を解きます
Extended Filter は非線形関数を局所線形化し、Unscented Filter は選択点を伝播し、Particle Filter は状態分布を模擬します。
自動的に精度が上がるわけではありません。強い非線形、厚い裾、多峰状態には各手法の診断が要ります。
全市場状態に一つの線形 Gaussian 再帰を強制するより、外れ値頑健モデルやレジーム転換が適切な場合があります。
金融フィルターは予言装置ではありません
ノイズのある気配値から効率価格、時変 beta、イールドカーブ因子、変動ヘッジ比率を推定できます。
推定値はモデルとパラメータを条件にします。「予測」という語が明日のリターンを教えるわけではありません。
各日付に実際に使えた情報で一段階 innovation、ローリング予測、状態修正、回転率、取引費用を評価します。
よくある質問
Kalman Filterは移動平均ですか?
いいえ。重みは固定窓でなく状態遷移モデルと変化する不確実性から決まります。
Kalman gainが高いとは?
予測状態の不確実性に比べ、新しい測定が相対的に多くの情報を持つという意味です。
フィルターは未来観測を使いますか?
フィルタリングは使わず、平滑化は使います。歴史分析では両者を区別します。
常に最適ですか?
厳密な条件付き平均の最適性には指定した線形 Gaussian 状態空間モデルの正しさが必要です。
出典
この考え方を自分のオプションに当てはめる
契約と目標を選ぶと、価格・時間・ボラティリティの前提を一つの分析で確認できます
自分のオプションを分析する(英語)