[論文レビュー] Generalized inverse Gaussian distributions and the time of first level crossing

V. K. Malinovskiǐ|arXiv (Cornell University)|Aug 29, 2017

Stochastic processes and statistical mechanics参考文献 20被引用数 5

ひとこと要約

本稿では、高水準 $u$ を超える複合再生過程 $V_s - cs$ の最初の通過時刻分布に対する、一般化逆ガウス（GIG）分布を活用した新規近似を提案する。この手法は、最初の通過時刻分布をGIG関連の積分および不完全ベッセル関数を用いて表現することで、特に臨界ドリフト $c = c^*$ の近傍で、リスク理論および再生過程における既存の漸近近似よりも顕著に精度を向上させる。

ABSTRACT

We propose a new approximation for the distribution of the time of the first crossing of a high level $u$ by random process $\homV{s}-cs$, where $\homV{s}$, $s>0$, is compound renewal process and $c>0$. It significantly outperforms the existing approximations, particularly in the region around the critical point $c=\cS$ which separates processes with positive and negative drifts. This approximation is tightly related to generalized inverse Gaussian distributions.

研究の動機と目的

• 複合再生過程における最初の通過時刻 $\Upsilon = \inf\{s > 0 : V_s - cs < -u\}$ の分布に対するより正確な近似を構築すること。
• 特に臨界ドリフト $c^*$ の近傍で従来の手法が破綻する問題を解消すること。
• 最初の通過時刻分布と一般化逆ガウス（GIG）分布との間のきめ細やかな数学的関係を確立すること。
• GIG関数および不完全ベッセル関数を用いて、破綻確率 $\mathsf{P}(\Upsilon \leq t)$ の計算可能かつ解析的に取り扱いやすい表現を提供すること。

提案手法

• 最初の通過時刻分布を、正規近似に基づく条件付き分布からの導出による三成分の混合 $\mathcal{M}_t(u,c,v)$, $\mathcal{F}_t(u,c,v)$, および $\mathcal{S}_t(u,c,v)$ でモデル化する。
• 最初の到着間隔時間 $T_1 = v$ における条件付きアプローチを採用し、生存確率を即時破綻と遅延破綻の二部分に分解する。
• 複合再生過程の漸近的正規性に依拠し、平均 $cM(1+x)$ および分散 $c^2D^2(1+x)/(u+cv)$ を用いる。ここで $M = \mathbb{E}[T]/\mathbb{E}[Y]$ であり、$D^2$ は分散尺度を表す。
• 主な数学的枠組みは、第二種変形ベッセル関数の $\frac{1}{2}$ 階数の不完全関数であり、ハンケルの級数展開を通じて逆ガウス分布と関連づけられる。
• 導出にはビネの積分恒等式を用い、変数変換と対称性を活用して最初の通過時刻積分を標準正規分布関数で表現可能な形に変換する。
• 最終的な近似式は $\mathcal{E}_t(u,c,v)$ として表され、主項、補正項 $C_{\mathcal{F}}$、歪度補正 $C_{\mathcal{S}}$ を含み、$c = c^*$ の近傍での精度向上を目的としてキャリブレーションされる。

実験結果

リサーチクエスチョン

• RQ1 複合再生過程の最初の通過時刻分布を、特に臨界ドリフト $c^*$ の近傍で、より正確に近似する方法は何か？
• RQ2 最初の通過時刻分布と一般化逆ガウス（GIG）分布との間の明確な数学的関係は何か？
• RQ3 不完全ベッセル関数と正規近似を用いることで、従来の漸近的手法よりもより強固で正確な近似が可能か？
• RQ4 補正項 $C_{\mathcal{F}}$ と $C_{\mathcal{S}}$ は、$c^*$ の周辺の臨界領域で性能をどのように向上させるか？

主な発見

• 提案された近似 $\mathcal{E}_t(u,c,v)$ は、特に臨界ドリフト $c^*$ の近傍で、従来手法を著しく上回り、$\mathsf{P}(\Upsilon \leq t)$ の推定において優れた性能を示す。
• 最初の通過時刻分布は、第二種変形ベッセル関数の $\frac{1}{2}$ 階数の不完全関数を通じて、一般化逆ガウス（GIG）分布と密接に結びついている。
• 逆ガウス分布の累積分布関数は $K_{1/2}(z)$ を用いて明示的に表現可能であり、これによりGIGフレームワークと接続され、解析的取り扱いが可能になる。
• ビネの積分恒等式を用いて最初の通過時刻積分を、標準正規分布関数を含む積分の和に変換し、効率的な数値計算を可能にする。
• 古典的近似に欠けている歪度およびスケール補正（$C_{\mathcal{F}}$, $C_{\mathcal{S}}$）を組み込むことで、精度が向上している。
• 理論的結果により、$u \to \infty$ であっても近似が高精度を維持することが確認され、特に破綻確率が急激に変化する $c = c^*$ の周辺領域で顕著な優位性を示す。