[論文レビュー] Stein's method of exchangeable pairs for absolutely continuous, univariate distributions with applications to the Polya urn model

Christian Döbler|arXiv (Cornell University)|Jul 2, 2012

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

ひとこと要約

本稿は、交換可能な対を用いたスティーンの方法の枠組みを、密度アプローチに焦点を当て、一変量の絶対連続分布に拡張する。特にベータ分布に適用し、ポリアのウールモデルにおける赤玉の相対数の収束速度 $ O(1/n) $ を確立する。密度の対数微分による新しい特徴付けと、スティーン方程式の有界解を用いることで達成される。

ABSTRACT

We propose a way of finding a Stein type characterization of a given absolutely continuous distribution $μ$ on $\R$ which is motivated by a regression property satisfied by an exchangeable pair $(W,W')$ where $\calL(W)$ is supposed or known to be close to $μ$. We also develop the exchangeable pairs approach within this setting. This general procedure is then specialized to the class of Beta distributions and as an application, a convergence rate for the relative number of drawn red balls among the first $n$ drawings from a Polya urn is computed.

研究の動機と目的

• 絶対連続な一変量分布の一般化に対して、密度アプローチを用いた交換可能な対のスティーンの方法を拡張すること。
• 密度の対数微分に基づく、このような分布のためのスティーン特徴付けを導出すること。
• 分布の密度およびその導関数を用いて、スティーン方程式の解のバインドを構築する枠組みを開発すること。
• 本手法をポリアのウールモデルに適用し、赤玉の割合の収束速度を分析すること。
• 交換可能な対の構成と解の正則性を用いて、ウォッシャー形式およびコルモゴロフ距離に対する定量的バインドを確立すること。

提案手法

• 密度 $ p $ を持つ分布 $ \mu $ に対して、対数微分 $ \psi(x) = \frac{d}{dx} \log p(x) $ を用いたスティーン恒等式を提案し、適切なクラスのすべての $ f $ に対して $ \mathbb{E}[f'(X) + \psi(X)f(X)] = 0 $ が成り立つようにする。
• スティーン方程式 $ f'(x) + \psi(x)f(x) = h(x) - \mu(h) $ を導入し、解として $ f_h(x) = \frac{1}{p(x)} \int_{-\infty}^x (h(y) - \mu(h)) p(y) dy $ を与える。これは台上でほとんど everywhere で有効である。
• 交換可能な対 $ (W, W') $ を用い、線形回帰性 $ \mathbb{E}[W' - W \mid W] = -\lambda W + R $ を満たし、$ R $ が $ \lambda $ に対して小さいものとすることで、分布距離をバインドする。
• 区間 $ (-1,1) $ 上のベータ分布にこの枠組みを適用し、関数 $ H(x), G(x), I(x), \eta(x), \gamma(x) $ を用いて、スティーン解の導関数に対する明示的バインドを導出する。
• 端点 $ \pm 1 $ における比関数 $ S(x) = \frac{H(x)G(x)}{I(x)\eta(x)} $ の振る舞いを解析し、ロピタルの定理を用いて、$ f_h' $ の有界性を確立する。
• スティーン方程式の解 $ f_h $ が有界であることは、$ \mathbb{E}[u(Z)] = 0 $ であることに同値であり、誤差バインドの適用可能性を保証する解の正則性を確保する。

実験結果

リサーチクエスチョン

• RQ1 交換可能な対のスティーンの方法を、密度アプローチを用いて一般の絶対連続一変量分布にどのように拡張できるか？
• RQ2 ベータ分布におけるスティーン方程式の構造とその解は何か？ また、その導関数をどのようにバインドできるか？
• RQ3 本フレームワークを用いて、ポリアのウールモデルにおける赤玉の割合の収束速度をどの程度確立できるか？
• RQ4 スティーン方程式の解が有界であるための条件は何か？ また、これは分布のモーメントとどのように関係するか？
• RQ5 ウォッシャー形式およびコルモゴロフ距離における誤差バインドは、交換可能な対の構成と残差 $ R $ の大きさにどのように依存するか？

主な発見

• 本稿では、密度の対数微分を用いた新しいスティーン特徴付けを導出し、正規分布を超えるスティーンの方法の体系的かつ一貫したアプローチを可能にする。
• 区間 $ (-1,1) $ 上のベータ分布に対して、スティーン解の導関数は $ \frac{2\|h'\|_\infty}{\alpha + \beta + 2} \cdot S(x) $ でバインド可能であり、ここで $ S(x) $ は区間上で一様有界な積分と密度の比である。
• 端点 $ \pm 1 $ における極限解析により、$ \sup_{x \in (-1,1)} S(x) < \infty $ が示され、$ -1 $ では明示的な有限極限 $ \frac{\alpha + \beta + 2}{2\beta + 4} $、$ 1 $ では $ \frac{\alpha + \beta + 2}{2\alpha + 4} $ をとる。これにより解の正則性が保証される。
• スティーン方程式の解 $ f_h $ が有界であることは、$ \mathbb{E}[u(Z)] = 0 $ であることに同値であり、誤差バインドの適用を保証する。
• ポリアのウールモデルにおける赤玉の相対数の収束速度は $ O(1/n) $ であることが示され、解の導関数の有界性と残差 $ R $ からの誤差バインドが導出される。
• 残差 $ R $ の大きさを制御し、解の導関数の有界性を活用することで、ウォッシャー形式およびコルモゴロフ距離において明示的なバインドが得られる。