[論文レビュー] Generalized Four Moment Theorem and an Application to CLT for Spiked Eigenvalues of Large-dimensional Covariance Matrices

Dandan Jiang, Zhidong Bai|arXiv (Cornell University)|Aug 16, 2018

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

ひとこと要約

本稿は、高次元スパiked共分散行列のモーメントマッチング条件を緩和する一般化された4次モーメント定理（G4MT）を導入する。G4MTは4次モーメントの要件を、尾確率条件に置き換えることで、スパiked固有値の極限分布の普遍性を確立する。これにより、有界な4次モーメントや対角ブロックの独立性を仮定しなくても、これらの固有値に対する中心極限定理（CLT）が成り立つようになり、BaiとYao（2012）の結果を、より広範かつ現実的である設定へと拡張する。

ABSTRACT

We consider a more generalized spiked covariance matrix $Σ$, which is a general non-definite matrix with the spiked eigenvalues scattered into a few bulks and the largest ones allowed to tend to infinity. By relaxing the matching of the 4th moment to a tail probability decay, a {\it Generalized Four Moment Theorem} (G4MT) is proposed to show the universality of the asymptotic law for the local spectral statistics of generalized spiked covariance matrices, which implies the limiting distribution of the spiked eigenvalues of the generalized spiked covariance matrix is independent of the actual distributions of the samples satisfying our relaxed assumptions. Moreover, by applying it to the Central Limit Theorem (CLT) for the spiked eigenvalues of the generalized spiked covariance matrix, we also extend the result of Bai and Yao (2012) to a general form of the population covariance matrix, where the 4th moment is not necessarily required to exist and the spiked eigenvalues are allowed to be dependent on the non-spiked ones, thus meeting the actual cases better.

研究の動機と目的

• 厳密な4次モーメントマッチングではなく、緩和されたモーメント条件の下で、一般化されたスパiked共分散行列の局所的スペクトル統計の漸近法則の普遍性を確立すること。
• BaiとYao（2012）がスパiked固有値のCLTに課していた、厳密な4次モーメントおよび対角ブロック独立性の制約を排除すること。
• 正確な4次モーメントマッチングを尾確率の減衰条件に置き換えることで、4次モーメント定理を一般化すること。
• 新しい仮定の下で、スパiked固有値の極限分布が分布に依存しないことを示すこと。
• 高次元統計的推論における共分散構造の分析に、より実用的かつ広く適用可能なフレームワークを提供すること。

提案手法

• 正確な4次モーメントマッチングの代わりに、尾確率条件（仮定B）に基づく一般化された4次モーメント定理（G4MT）を提唱する。
• 高次元のランダム行列における高階微分係数の推定を避けるために、Lindeberg型戦略を用いて非ガウス型要素をガウス型に置き換える。
• G4MTを応用して、一般化されたスパiked共分散行列の正規化されたスパiked固有値に対するCLTを導出する。
• C₀条項（一様な指数的尾部減衰）を必要としない低次元の行列固有値解析に、普遍性問題を還元する。
• 緩和されたモーメントおよび尾確率条件を満たす分布のクラスにおいて、スパiked固有値の極限分布が普遍的であることを確立する。
• ガウス分布および二項分布の仮定の下でシミュレーションを実施し、標本サイズの増加に伴うスパイクの位置と大きさの推定が一貫していることを確認する。

実験結果

リサーチクエスチョン

• RQ1 正確な4次モーメントマッチングよりも緩いモーメント条件の下でも、スパiked固有値の極限分布の普遍性を確立できるか？
• RQ2 母共分散行列に対角ブロックの独立性が存在しない場合、スパiked固有値のCLTに影響を与えるか？
• RQ3 C₀条項（指数的尾部減衰）を緩和しても、スペクトル統計の普遍性は維持されるか？
• RQ4 4次モーメント条件を、4次モーメントの代わりに尾確率の減衰条件に置き換えると、スパiked固有値の漸近的正規性は失われるのか？
• RQ5 提案されたG4MTは、重尾または非i.i.d.成分を含む実世界データに対して、高次元統計におけるCLTの適用性をどの程度向上させるか？

主な発見

• G4MTは、尾確率条件（仮定B）の下で、スパiked固有値の極限分布の普遍性を確立し、正確な4次モーメントマッチングの必要性を排除する。
• 有界な4次モーメントや対角ブロック独立性を仮定しなくても、スパiked固有値のCLTが成り立つ。これは、BaiとYao（2012）の結果を一般化された母共分散行列構造へと拡張する。
• ガウス分布および二項分布の仮定の下でのシミュレーションでは、スパイクの位置と大きさが一貫して推定された。p=200、n=1000の場合、推定されたスパイクの位置は(1,2,3,198,199,200)、大きさはそれぞれ≈(3.949, 3.217, 2.887, 0.231, 0.188, 0.102)であった。
• p=400、n=1000の場合、推定されたスパイクの位置は(1,2,3,198,199,200)のままであり、大きさは≈(4.123, 3.211, 3.001, 0.216, 0.195, 0.096)であった。これは次元の増加に対しても安定であることを示している。
• 慢性腎臓病データ（n=114、p=24）の実データ解析では、9つのスパイクが検出され、位置は(1,2,18,19,20,21,22,23,24)、大きさは≈(10.818, 2.143, 0.219, ..., 0.009)であった。これにより、データ内に3つの明確な要因群が存在することが明らかになった。
• C₀条項および高階微分係数の推定を回避することで、重尾または非サブガウス型データに対しても、よりロバストで適用可能な手法が得られた。