[論文レビュー] Rate of Convergence of the Expected Spectral Distribution Function to the Marchenko -- Pastur Law

Friedrich Götze, А. Н. Тихомиров|arXiv (Cornell University)|Dec 19, 2014

Random Matrices and Applications参考文献 17被引用数 6

ひとこと要約

この論文は、行列要素のモーメントおよび有界性条件の下で、標本共分散行列の期待スペクトル分布関数がMarchenko-Pastur分布に収束するレートを確立する。平均0、分散1、有限な4次モーメント、および要素ごとの有界性 $ Dn^{1/4} $ の仮定の下で、期待スペクトル分布関数とMarchenko-Pastur分布との間のコルモゴロフ距離が $ O(n^{-1}) $ であることが示され、非漸近的収束レートが得られる。

ABSTRACT

Let $\mathbf X=(X_{jk})$ denote a $n imes p$ random matrix with entries $X_{jk}$, which are independent for $1\le j\le n, 1\le k\le p$. Let $n,p$ tend to infinity such that $\frac np=y+O(n^{-1})\in(0,1]$. For those values of $n,p$ we investigate the rate of convergence of the expected spectral distribution function of the matrix $\mathbf W=\frac1{ p}\mathbf X\mathbf X^*$ to the Marchenko-Pastur law with parameter $y$. Assuming the conditions $\mathbf E X_{jk}=0$, $\mathbf E X_{jk}^2=1$ and $ \quad \quad \quad \quad \quad \quad \quad \sup_{n,p\ge1}\sup_{1\le j\le n,1\le k\le p}\mathbf E |X_{jk}|^4=: μ_4

研究の動機と目的

• 標本共分散行列の期待スペクトル分布関数がMarchenko-Pastur分布に非漸近的収束するレートを確立すること。
• Wigner行列に関する先行研究を、標本共分散行列へと拡張するために、期待スペクトル分布と極限となるMarchenko-Pastur分布との間のコルモゴロフ距離を分析すること。
• 最小限のモーメントおよび要素ごとの有界性仮定の下で、コルモゴロフ距離に鋭い $ O(n^{-1}) $ 界を与えること。
• パrameter範囲 $ y \in (0,1] $ 全体にわたって一様に成り立つ定量的収束レートを提供すること。このレートは $ D $, $ \mu_4 $, $ y $, および $ c_y $ に明示的な依存関係を示す。

提案手法

• スペクトル分布のステュルツェス変換とMarchenko-Pastur分布のステュルツェス変換を比較することで収束を分析する。
• 解体法（resolvent）技術と、$ \mathbf{X} $ および $ \mathbf{X}^* $ から構成されるヘルミート行列 $ \mathbf{V} $ を用いて、固有値の振る舞いを解体法の成分に関連付ける。
• 集中不等式およびモーメントの有界性を用いて、解体法展開における誤差項を制御し、特に $ m_n(z) $ と $ s_y(z) $ の差に注目する。
• 条件付きモーメント推定とデカップリング技術を用いて、解体法恒等式における誤差項の4次モーメントを有界化する。これには $ \mu_4 $ および要素ごとの有界性に関する仮定を活用する。
• 複数段階の誤差分解を実行し、$ \beta_{j1}, \beta_{j2}, \beta_{j3} $ を用いて行列構造の異なる部分からの寄与を分離・制御する。
• 解体法成分の虚部に対する事前有界性を適用し、グリーン関数比較法を制御的に用いて、スペクトルパラメータ $ z $ 全体にわたる一様推定を得る。

実験結果

リサーチクエスチョン

• RQ1 標本共分散行列の期待スペクトル分布関数がMarchenko-Pastur分布に非漸近的収束するレートは何か？
• RQ2 収束レートは行列要素の4次モーメントおよび要素ごとの有界性にどのように依存するか？
• RQ3 最小限のモーメントおよび有界性仮定の下で、鋭い $ O(n^{-1}) $ レートを確立できるか？
• RQ4 比 $ n/p \to y \in (0,1] $ のとき、$ |n/p - y| = O(n^{-1}) $ ならば、コルモゴロフ距離はどのように振る舞うか？

主な発見

• 条件 $ \mathbb{E}[X_{jk}] = 0 $, $ \mathbb{E}[X_{jk}^2] = 1 $, $ \sup \mathbb{E}[|X_{jk}|^4] = \mu_4 < \infty $, および $ \sup |X_{jk}| \leq Dn^{1/4} $ の下で、コルモゴロフ距離 $ \Delta_n = \sup_x |F_n(x) - G_y(x)| $ は $ C n^{-1} $ で有界である。
• 有界性における定数 $ C $ は $ D $, $ \mu_4 $, $ y $, および $ c_y $ のみに依存し、$ n $ や $ p $ とは無関係である。
• 結果は $ y \in (0,1] $ 全体にわたって一様に成り立ち、与えられたモーメントおよび有界性仮定の下で収束レートは鋭いものである。
• ステュルツェス変換および解体法誤差項の詳細な分析を通じて、$ O(n^{-1}) $ レートが得られる。この分析にはモーメントの有界性および条件付き期待値技術が用いられる。
• 論文は、4次モーメントが一様に有界であれば、要素が同一分布でなくても、期待スペクトル分布とMarchenko-Pastur分布との差が $ n^{-1} $ 以上で減少することを示している。