[論文レビュー] An algebra of Stein operators

Robert E. Gaunt, Guillaume Mijoule+1 more|arXiv (Cornell University)|Apr 22, 2016

Random Matrices and Applications参考文献 52被引用数 4

ひとこと要約

本稿では、独立した確率変数の積として得られる確率変数を対象として、多項式スティーン作用素を構築するための新規な代数的枠組みを導入する。微分作用素と多項式係数の組み合わせを用い、双対性の原則を活用することで、正規分布、ベータ分布、分散ガンマ分布、一般化ガンマ分布など広範なクラスの分布に対して、明示的なk次スティーン作用素を導出する。特に、k個の独立した対称分散ガンマ分布に従う確率変数の積の密度に関する新しい公式を提供する。この作用素は、ターゲット密度が閉形式で与えられない場合でも、計算が容易であることを保証する。

ABSTRACT

We build upon recent advances on the distributional aspect of Stein's method to propose a novel and flexible technique for computing Stein operators for random variables that can be written as products of independent random variables. We show that our results are valid for a wide class of distributions including normal, beta, variance-gamma, generalized gamma and many more. Our operators are $k$th degree differential operators with polynomial coefficients; they are straightforward to obtain even when the target density bears no explicit handle. As an application, we derive a new formula for the density of the product of $k$ independent symmetric variance-gamma distributed random variables.

研究の動機と目的

• 独立した確率変数の積分布に対して、系統的かつ柔軟な代数的枠組みを構築し、スティーン作用素を導出すること。
• 密度関数が明示的に与えられていない、あるいは特殊関数によって定義される分布に対しても、スティーン法の適用範囲を拡張すること。
• 広範なターゲット分布クラスに対して、多項式係数を有するk次微分スティーン作用素を統一的に構築するアプローチを提供すること。
• k個の独立した対称分散ガンマ分布に従う確率変数の積の密度関数に対する明示的公式を導出すること。
• メイジャーG関数やクーマーU関数によって定義される古典的および特殊関数に基づく分布への応用を通じて、提案手法の計算可能性と一般性を示すこと。

提案手法

• 著者らは、乗法作用素Mと微分作用素Dの単項式の有限線形結合として定義される多項式スティーン作用素のクラスを導入する。すなわち、A = ∑ a_ij M^i D^j である。
• 密度関数が線形常微分方程式を満たす場合に特に有効な、フォッカー・プランク方程式と部分積分に基づく双対性の議論を用いて、密度関数の微分方程式から作用素を導出する。
• 特殊関数（例：クーマーU関数、メイジャーG関数）によって定義される分布に対しては、これらの関数が満たす既知の微分方程式を用いて、対応するスティーン作用素を構築する。
• 独立成分の積分布への一般化は、同時密度が因数分解可能であり、作用素構造が独立性を保つことを利用することで実現される。
• 正規化およびスケーリング戦略を含む：正規分布および分散ガンマ分布では、作用素に右からMを乗じる。PRR分布では、M²を乗じる。
• 古典的分布（例：正規分布、ガンマ分布、ベータ分布、F分布、逆ガンマ分布、ステューデントt分布、PRR分布、分散ガンマ分布、一般化ガンマ分布）に対して既知および新規のスティーン作用素を導出し、既存の結果と整合することを確認することで、フレームワークの妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 ターゲット密度が明示的に与えられていない場合に、独立した確率変数の積分布に対して、多項式スティーン作用素を体系的かつ一貫して構築する方法は何か？
• RQ2 広範な分布クラスに対して多項式係数を有するk次微分スティーン作用素を構築する背後にある代数的構造は何か？
• RQ3 フォッカー・プランク方程式とスティーン作用素の双対性を活用することで、特殊関数によって定義される分布に対しても、計算可能な作用素を導出できるか？
• RQ4 k個の独立した対称分散ガンマ分布に従う確率変数の積に対するスティーン作用素の形は何か？
• RQ5 第一および第二の作用素に既知のアプローチを一般化し、高次の場合に適用するにはどうすればよいか？

主な発見

• 本稿では、k個の独立した対称分散ガンマ分布に従う確率変数の積に対して、k次多項式スティーン作用素を成功裏に構築し、その密度関数に対する新しい明示的公式を提供した。
• 本手法により、正規分布、ベータ分布、ガンマ分布、F分布、逆ガンマ分布、スチューデントt分布、PRR分布、分散ガンマ分布、一般化ガンマ分布など、広範な分布に対して一貫性があり体系的な方法でスティーン作用素を導出でき、すべて多項式係数を有する微分作用素として明示的に表現された。
• μ=0 の分散ガンマ分布に対しては、スティーン作用素が σ²T₁Tᵣ + 2θMTᵣ/₂ − M² で与えられ、Tᵢ は f ↦ x^i f’(x) を表す。
• 既知の結果の一般化が可能である：例えば、古典的な正規スティーン作用素 Af(x) = f’(x) − xf(x) は、一般枠組みの特殊ケースとして回復される。
• 特殊関数で定義される分布に対しても本手法は頑健である：PRR分布では、クーマーU関数の微分方程式から導かれる作用素 sT₁T₂ − M²T₂ₛ が得られる。
• ターゲット密度が明示的に与えられていなくても、線形微分方程式を満たす限り、本フレームワークによりスティーン作用素を導出可能である。これは、標準的なパラメトリック分布族を超える広範な適用可能性を示している。