[論文レビュー] An Exposition of Götze's Estimation of the Rate of Convergence in the Multivariate Central Limit Theorem

Rabi Bhattacharya, Susan Holmes|arXiv (Cornell University)|Mar 22, 2010

Random Matrices and Applications参考文献 11被引用数 19

ひとこと要約

本稿は、Götze が Stein の方法を用いて、凸 Borel 集合上での多変量中心極限定理に対する Berry-Esseen 型の上限を導出する応用を詳細に解説している。収束速度は $ O(k^{5/2} ⑶_3) $ であり、$ ⑶_3 $ は3次の絶対モーメントノルムを表す。本稿は Götze の元々の議論を完全な導出とともに再検討し、以前に報告されたものよりも高次元的な依存関係を特定した。その後、Ball の不等式を用いて $ O(k^{3/2}) $ に改善されたが、Bentkus の $ O(k^{1/4}) $ は別法により依然として最適である。

ABSTRACT

We provide an explanation of the main ideas underlying Götze's main result in using Stein's method. We also provide detailed derivations of various intermediate estimates. Curiously, we are led to a different dimensional dependence of the constant than that given Götze's paper. We would like to dedicate this to Charles Stein on the occasion of his 90th birthday.

研究の動機と目的

• Stein の方法を用いた Götze の多変量中心極限定理の誤差項の境界に関する明確でアクセス可能な説明を提供すること。
• 確率論者や研究者にとって技術的枠組みをより明確にするために、中間的な推定をすべて詳細に導出すること。
• 誤差項の次元依存性を再表現し、Götze の元々の $ O(k) $ 評価との不一致を特定すること。
• Ball の不等式を用いて境界を改善し、$ O(k^{5/2}) $ から $ O(k^{3/2}) $ に次元依存性を低減すること。また、既知の最適結果と比較すること。

提案手法

• Ornstein-Uhlenbeck過程の生成子を Stein 演算子として用い、分布距離を微分方程式の解に結びつける Stein の方法を用いる。
• 平滑化作用素 $ T_t $ を凸集合の特性関数に作用させ、O-U過程の遷移密度を介して滑らかな関数に変換する。
• Stein 方程式に基づく摂動論法を用い、期待値の差を $ E[Lg_0 - L_αg_α] $ の形で表現する。ここで $ L $ は生成子で、$ g $ は Stein 方程式の解である。
• 3階微分係数を含む積分剰余項を伴う3階テイラー展開による誤差の表現を導出する。関数 $ \psi_t $ の3階微分が関与する。
• 同分布独立な確率的ベクトルの和の正規分布近似の分布的性質を用いて、剰余項の条件付き期待値を推定する。
• ガウス積分を $ Q_{(n-1),j} - \Phi $ と $ \Phi $ に分解し、前者を $ \delta_{n-1} $ で上界評価し、後者を変数変換と行列ノルムを用いて評価する。

実験結果

リサーチクエスチョン

• RQ1 Stein の方法を用いた Götze の多変量 CLT 評価における誤差項の正しい次元依存性は何か？
• RQ2 本稿で導出された境界が Götze の元々の $ O(k) $ 評価と異なるのはなぜで、真のオーダーは何か？
• RQ3 既知の不等式（例：Ball の不等式）を用いて誤差項の境界を改善できるか？
• RQ4 高次元における収束速度に、3次の絶対モーメント $ \beta_3 $ がどのように影響するか？
• RQ5 この手法は従属する確率的ベクトルへと拡張可能か。その制限は何か？

主な発見

• 本稿では収束速度 $ O(k^{5/2} \beta_3) $ を導出し、Götze の元々の $ O(k) $ の主張とは矛盾する。これは次元に高次の依存性があることを示唆している。
• Ball の不等式を用いることで、境界を $ O(k^{3/2} \beta_3) $ に改善でき、次元の成長を顕著に低減できる。
• $ \delta_n \leq c k \gamma_3 $ という境界が得られ、ここで $ \gamma_3 = \sum_{j=1}^n E\left(\sum_{i=1}^k |X_j^{(i)}|\right)^3 $ である。$ \gamma_3 \ll k^{3/2} \beta_3 $ のとき、$ O(k^{5/2} \beta_3) $ よりもタイトな推定が可能である。
• 解析により、$ S_n - X_j $ の共分散の逆行列 $ N_j $ を用いた正規化が、剰余項の上界評価において中心的な役割を果たすことが明らかになった。
• ガウス積分の細かな分解と、行列 $ A_j = e^{-s} \sqrt{1 - e^{-2s}} N_j $ を用いた変数変換が、条件付き期待値における摂動効果を捉える上で不可欠である。
• 著者らは、最適でないが、Bentkus の手法よりも透明性が高く、従属設定への拡張にも適していると結論づけている。