[論文レビュー] A note on universality of the distribution of the largest eigenvalues in certain sample covariance matrices

Alexander Soshnikov|ArXiv.org|Apr 10, 2001

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

ひとこと要約

この論文は、ガウス分布でない場合の標本共分散行列の最大固有値に関する普遍性結果を拡張する。Wigner行列理論からの組合せ的手法を用いて、一般の i.i.d. 要素（有限モーメントおよびサブガウス型尾部を有する）のもとで、上位固有値の連続分布がトレーシー・ウィドム分布に収束することを証明し、最大固有値のスケーリングに関する確実な上界 $(n^{1/2} + p^{1/2})^2 + O(p^{1/2} \log p)$ を確立する。

ABSTRACT

Recently Johansson and Johnstone proved that the distribution of the (properly rescaled) largest principal component of the complex (real) Wishart matrix $ X^* \* X (X^t \*X) $ converges to the Tracy-Widom law as $ n, p $ (the dimensions of $ X $) tend to $ \infty $ in some ratio $ n/p o γ>0. $ We extend these results in two directions. First of all, we prove that the joint distribution of the first, second, third, etc. eigenvalues of a Wishart matrix converges (after a proper rescaling) to the Tracy-Widom distribution. Second of all, we explain how the combinatorial machinery developed for Wigner matrices allows to extend the results by Johansson and Johnstone to the case of $ X $ with non-Gaussian entries, provided $ n-p =O(p^{1/3}) . $ We also prove that $ λ_{max} \leq (n^{1/2}+p^{1/2})^2 +O(p^{1/2}\*\log(p)) $ (a.e.) for general $ γ>0.$

研究の動機と目的

• ガウス分布でない場合の実および複素標本共分散行列の最大固有値に対するトレーシェ・ウィドム分布の普遍性を確立すること。
• ジョーハンソンおよびジョンストーンの最大固有値の極限分布に関する結果を、モーメントおよび尾部減衰条件の下で非ガウス型要素へ拡張すること。
• 適切なスケーリングの下で、最初の数個の最大固有値の連続収束をトレーシェ・ウィドム分布に示すこと。
• 一般の $\gamma = n/p > 0$ のもとで、$X^*X$ や $X^tX$ の最大固有値 $λ_{\max}$ に対する確実な上界を導出すること。

提案手法

• Wigner確率行列理論からの組合せ的道具を標本共分散行列の固有値モーメントの解析に適応する。
• 生成関数およびコーシーの積分公式を用いて、モーメント展開における関連パスの数を推定し、特に非閉路頂点に注目する。
• 複素解析を用いて、二次関数の平方根を含む積分を漸近的に評価し、パス数の $m^{-3/2}$ 減衰率を導出する。
• モーメントおよび尾部減衰条件を課す：$\mathbb{E}|x_{ij}|^{2m} \leq (\text{const} \cdot m)^m$ および対称分布を仮定して高次モーメントを制御する。
• チェビシェフの不等式およびボレル＝カンテリの補題を用いて、モーメント推定から確実な収束結果を導出する。
• 制約 $n - p = O(p^{1/3})$ の下でのパス数の漸近的挙動を解析し、この領域で普遍性が成立することを保証する。

実験結果

リサーチクエスチョン

• RQ1 標本共分散行列の上位 $k$ 個の固有値の連続分布は、非ガウス型 i.i.d. 要素のもとでトレーシェ・ウィドム分布に収束するか？
• RQ2 組合せ的モーメント手法を用いて、最大固有値に対するトレーシェ・ウィドム普遍性をガウス分布でない場合に拡張できるか？
• RQ3 一般の $\gamma > 0$ のもとで、$X^*X$ の最大固有値 $λ_{\max}$ に対する確実な上界は何か？
• RQ4 一般の要素分布のもとで、最大固有値の収束速度は標本および母集団の次元にどのように依存するか？

主な発見

• 非ガウス型 i.i.d. 要素を有する標本共分散行列の最初の $m$ 個の最大固有値の連続分布は、適切なスケーリングの下でトレーシェ・ウィドム分布に収束する。
• 条件 $n - p = O(p^{1/3})$ の下で収束が成立し、ガウス分布でない場合の普遍性が拡張される。
• 一般の $\gamma > 0$ に対して、最大固有値は $\lambda_{\max} \leq (n^{1/2} + p^{1/2})^2 + O(p^{1/2} \log p)$ almost surely を満たす。
• モーメント展開における関連パスの漸近的数は $\sim \frac{y^{1/4}(\sqrt{y}+1)}{2\sqrt{\pi}} \frac{(\sqrt{y}+1)^m}{m^{3/2}}$ と振る舞い、トレーシェ・ウィドムスケーリングに対応する。
• 各 $k \geq 2$ に対して $\sum_{k=2}^m k n_k > 0$ を満たすパスに対応するモーメントの部分和は $o\left(\frac{p \mu_{n,p}^m}{m^{3/2}}\right)$ である。これは、主要寄与は特定のパスクラスに由来することを確認する。
• 証明により、提示された条件下で、スケーリング済み最大固有値の $m$ 次モーメントがトレーシェ・ウィドム分布の $m$ 次モーメントに収束することが示された。