[論文レビュー] Stein characterizations for linear combinations of gamma random variables

Benjamin Arras, Ehsan Azmoodeh+2 more|arXiv (Cornell University)|Aug 31, 2017

Random Matrices and Applications参考文献 16被引用数 3

ひとこと要約

本稿では、特性関数が単純な常微分方程式（ODE）を満たすことを活用して、ガンマ確率変数の線形結合に対するステイン作用素を明示的に導出する新規で明示的な手法を提案する。この手法は、第二ウィエンル・コーン上でのマリヤン・カレンスラに接続され、独立または非独立なガンマ変数の和に対する閉形式のステイン作用素を導く。マケイ型I分布への応用も含む。

ABSTRACT

In this paper we propose a new, simple and explicit mechanism allowing to derive Stein operators for random variables whose characteristic function satisfies a simple ODE. We apply this to study random variables which can be represented as linear combinations of (non necessarily independent) gamma distributed random variables. The connection with Malliavin calculus for random variables in the second Wiener chaos is detailed. An application to McKay Type I random variables is also outlined.

研究の動機と目的

• 特性関数が単純なODEを満たす確率変数に対して、一般的で明示的なステイン作用素の導出メカニズムを構築すること。
• 従来の手法が十分にカバーしていない非独立なガンマ確率変数の線形結合に対し、ステイン法を拡張すること。
• 第二ウィエンル・コーンに属する確率変数に対して、ステイン作用素とマリヤン・カレンスラの明確な数学的関係を確立すること。
• マケイ型I分布および関連するガンマ混合分布に対して明示的なステイン作用素を提供すること。
• ODEに基づく特徴付けを用いて、ステイン法の適用範囲をガウス分布およびポアソン近似を超えて、複雑な非楕円型分布へと一般化すること。

提案手法

• 目的分布の特性関数が満たすODEを分析することで、ステイン作用素を導出するためのフーリエ解析的フレームワークを提案する。
• 特性関数の導関数の構造を用いて、ステイン作用素の微分作用素形式を同定する。
• ヴィエートの定理を用いて、微分作用素の係数をガンマパラメータの基本対称多項式で表現する。
• マリヤン・カレンスラを活用し、ステイン核τF(F)を発散作用素およびオルンシュタイン＝ウーレンバック過程と関連づけ、分散の上限評価を可能にする。
• 作用素のフーリエ記号における多項式係数の一致を用いて、ステイン作用素の明示的表現を導出する。
• フーリエ領域におけるコーシー＝リプシッツの定理を用いて、特性関数の解の一意性を証明し、分布の同値性を確立する。

実験結果

リサーチクエスチョン

• RQ1 特性関数が単純なODEを満たす分布に対して、体系的で明示的なステイン作用素の導出法を開発可能か？
• RQ2 第二ウィエンル・コーンに属する確率変数に対して、ステイン作用素とマリヤン・カレンスラの関係をどのように形式化できるか？
• RQ3 非独立または独立なガンマ確率変数の線形結合に対して、ステイン作用素の明示的形は何か？
• RQ4 このフレームワークを用いて、マケイ型I分布をステイン法で特徴づけられるか？
• RQ5 ガンマ混合分布のパラメータに対して、それが一次のステイン作用素を許容するための必要十分条件は何か？

主な発見

• 特性関数が多項式係数をもつ線形ODEを満たす任意の確率変数に対して、明示的なステイン作用素を導出するための新規で一般的なメカニズムが確立された。
• パラメータ(λj, αj, cj)を持つガンマ確率変数の線形結合に対して、ステイン作用素は、λj cjの逆数の基本対称多項式を係数に含むd次線形微分作用素として明示的に与えられる。
• マケイ型I分布（カイ二乗分布の混合）に対しても、対応する特性関数を同定し、関連ODEを解くことで、閉形式のステイン作用素が得られる。
• マリヤン・カレンスラとの関係が形式化された：ステイン核τF(F)は⟨DF, −DL⁻¹F⟩Hとして表現され、分散の上限評価にマリヤン・カレンスラの道具が利用可能になる。
• フーリエ解析とコーシー＝リプシッツの定理を用いた証明技法により、解の一意性が保証され、導出された作用素が目標分布を一意に特徴付けることが保証される。
• このフレームワークは独立または非独立なガンマ成分の両方に対応でき、ステイン法の適用範囲を第二ウィエンル・コーン内での非i.i.d.ガンマ混合にまで拡張する。