[論文レビュー] On the consistency of Sobol indices with respect to stochastic ordering of model parameters

Areski Cousin, Alexandre Janon+2 more|arXiv (Cornell University)|Jul 21, 2014

Probabilistic and Robust Engineering Design参考文献 29被引用数 5

ひとこと要約

本稿は、ソボル感度インデックスと確率的順序理論の理論的整合性を確立し、モデル出力に凸性および単調性条件が成り立つ場合、ソボルインデックスが超過財産順序および分散順序を満たすことを示している。主な貢献は、確率的順序の意味での入力の不確実性が増加すると、それに応じてソボルインデックスも増加することを証明したことであり、これは、不確実性評価およびリスク管理におけるソボルインデックスの使用を正当化するものである。

ABSTRACT

In the past decade, Sobol's variance decomposition have been used as a tool - among others - in risk management. We show some links between global sensitivity analysis and stochastic ordering theories. This gives an argument in favor of using Sobol's indices in uncertainty quantification, as one indicator among others.

研究の動機と目的

• ソボルインデックスがモデル入力の確率的順序に関して一貫性を示すかどうかを調査すること。
• 入力の不確実性（確率的順序の意味で）が増加する場合に、それがソボルインデックスの増加を引き起こす理論的条件を確立すること。
• 確率的順序理論と既存の理論を結びつけることで、ソボルインデックスが不確実性評価およびリスク管理における信頼性の高いツールとして使用可能であることを検証すること。
• 入力分布の変更が、特にヘストンおよびヴァシチェクモデルのような金融モデルにおけるソボルインデックスの順位に与える影響を調査すること。
• ソボルインデックスを、入力分布のクラス全体にわたる普遍的な上限として使用する理論的基盤を提供すること。

提案手法

• 主な確率的順序（標準的確率的順序、凸順序、分散順序、超過財産順序、スターパターン順序、ローレンツ順序）を定義し、適用する。
• 分散分解を用いて、ソボルインデックスを $ S_X = \frac{\mathrm{Var}(\mathbb{E}(Y|X))}{\mathrm{Var}(Y)} $ として表現し、入力が出力分散に与える影響を定量化する。
• 出力関数 $ f $ が凸または単調である場合に、ソボルインデックスが超過財産順序および分散順序と整合することを証明する。
• 理論的結果を金融モデル（ヴァシチェクおよびヘストンモデル）に適用し、現実的な入力不確実性下での整合性を検証する。
• モンテカルロサンプリングおよび高速フーリエ変換を用いて、オプション価格および全ソボルインデックスを計算する数値実験を実施する。
• 異なる入力分布（例：一様分布対確率的順序に従う分布）におけるソボルインデックスを比較し、分布仮定への感度を評価する。

実験結果

リサーチクエスチョン

• RQ1 出力関数 $ f $ にどのような条件下で、ソボルインデックスが超過財産順序や分散順序といった確率的順序と整合するのか。
• RQ2 入力パラメータの分布が変化した場合（例：$ r^* \leq_{\text{disp}} r $）に、ソボルインデックスの順位および大きさにどのような影響が生じるか。
• RQ3 入力不確実性が確率的順序の意味で増加する場合、ソボルインデックスがどの程度まで整合性を保つのか。
• RQ4 入力分布が曖昧または正確に特定できない状況下でも、ソボルインデックスをリスク管理に信頼性を持って使用できるか。
• RQ5 確率的順序との理論的整合性を活用して、入力分布のクラス全体にわたるソボルインデックスの普遍的上限を導出できるか。

主な発見

• 出力関数 $ f $ が凸または単調である場合、定理2.1および2.3で示されるように、ソボルインデックスは超過財産順序および分散順序と整合する。
• ヘストンモデルにおいて、金利分布を $ r^* \sim \mathcal{U}[0,1] $ から $ r \sim \mathcal{U}[0,2] $ に変更した（$ r^* \leq_{\text{disp}} r $ を満たす）と、全ソボルインデックスは0.32から0.73に増加し、理論的整合性が確認された。
• 金利 $ r $ が一様分布 [0,1] に従うとき、配当利回り $ q $ は2番目に影響力の強いパラメータであるが、$ r $ の分布が拡大すると、$ r $ の影響力が $ q $ を上回るようになる。
• パrameter $ \kappa, \theta, \sigma, \rho $ は、両方の分布仮定下でも、0.01未満の非常に小さいソボルインデックスを示しており、オプション価格の不確実性にほとんど影響しないことが示された。
• 本研究は、ソボルインデックスの順位付けが入力分布の仮定に極めて敏感であることを示しており、堅牢な分布モデリングの重要性を強調している。
• 確率的順序との理論的整合性は、ソボルインデックスを、入力分布のクラス全体にわたる普遍的な上限として導出する可能性を示しており、不確実性評価における信頼性を高めるとともに、応用の幅を広げるものである。