[論文レビュー] The rank 1 real Wishart spiked model I. Finite N analysis

M. Y. Mo|arXiv (Cornell University)|Nov 24, 2010

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

ひとこと要約

本稿は、ランク1の実Wishartスプライクドモデルに対する有限N解析を確立し、直交多項式技法を用いて最大固有値分布のFredholm行列式表現を導出する。フェーズ転移点における漸近的最大小固有値分布の最初の正確な計算がなされ、確率的作用素法とは異なり、超楕円積分とO(N)上でのゾーン多項式展開を用いる。

ABSTRACT

This is the first part of a paper that studies the phase transition in the asymptotic limit of the rank 1 real Wishart spiked model. In this paper, we consider $N$-dimensional real Wishart matrices $S$ in the class $W_{\mathbb{R}}\left(Σ,M ight)$ in which all but one eigenvalues of $Σ$ is $1$. Let the non-trivial eigenvalue of $Σ$ be $1+τ$, then as $N$, $M ightarrow\infty$, with $N/M=γ^2$ finite and non-zero, the eigenvalue distribution of $S$ will converge into the Machenko-Pastur distribution inside a bulk region. As $τ$ increases from zero, one starts seeing stray eigenvalues of $S$ outside of the support of the Machenko-Pastur density. As the first of these stray eigenvalues leaves the bulk region, a phase transition will occur in the largest eigenvalue distribution of the Wishart matrix. In this paper will compute the asymptotics of the largest eigenvalue distribution when the phase transition occur. In the this first half of the paper, we will establish the results that are valid for all $N$ and $M$ and will use them to carry out the asymptotic analysis in the second half of the paper, which will follow shortly. In particular, we have derived a formula for the integral $\int_{O(N)}e^{- r(XgYg^T)}g^T\D g$ when $X$, $Y$ are symmetric and $Y$ is a rank 1 matrix. This allows us to write down a Fredholm determinant formula for the largest eigenvalue distribution and analyze it using orthogonal polynomial techniques. This approach is very different from a recent paper by Bloemendal and Virag, in which the largest eigenvalue distribution was obtained using stochastic operator method.

研究の動機と目的

• ランク1実Wishartスプライクドモデルにおける最大固有値分布の有限N挙動を分析すること。
• 直交多項式技法を用いて固有値の連関確率密度関数の正確な式を導出すること。
• スプライク固有値がバルク支持から脱出するフェーズ転移点における最大固有値分布の漸近的性質を計算すること。
• O(N)上でのゾーン多項式展開とハール積分評価を用いて、確率的作用素法の厳密な代替手法を提供すること。

提案手法

• Yがランク1対称行列であるとき、群積分 ∫_{O(N)} exp(−tr(XgYg^T)) g^T dg の閉形式表現を導出する。
• Zonal多項式 Z_p(X) を用いて、群積分を Z_{(k)}(X) と正規化定数 Z_{(k)}(I_N) を含む級数に表現する。
• Zonal多項式の母関数恒等式を適用し、積分を変数 t に関する複素積分に変換する。
• 無限遠点における留数を用いて複素積分を評価し、t^{-1} 項の係数を抽出することで、超楕円積分表現を獲得する。
• 導出された積分表現を用いて、最大固有値分布のFredholm行列式表現を確立する。
• 論文の後半で漸近的解析を実施し、N, M → ∞ かつ N/M = γ² を固定するフェーズ転移領域に焦点を当てる。

実験結果

リサーチクエスチョン

• RQ1 ランク1実Wishartスプライクドモデルにおける有限N状態での最大固有値分布はどのように振る舞うか？
• RQ2 ランク1の共分散摂動を伴う実Wishart行列の固有値の連関確率密度関数の正確な形は何か？
• RQ3 行列Yがランク1であるとき、O(N)上でのハール積分はどのように評価可能となり、固有値統計の計算が可能になるか？
• RQ4 スプライク固有値がバルクから脱出するフェーズ転移点における最大固有値分布の漸近的挙動は何か？
• RQ5 直交多項式アプローチは、同じ漸近的極限を導出する確率的作用素法と比較してどのように異なるか？

主な発見

• Yがランク1であるとき、積分 ∫_{O(N)} exp(−tr(XgYg^T)) g^T dg の正確な式が、超楕円積分として導出された。
• 最大固有値分布は、直交多項式とゾーン多項式展開を含むFredholm行列式として表現可能であることが示された。
• フェーズ転移点における漸近的解析は、有限Nの公式のおかげで可能となり、極限分布の計算が可能になった。
• 本手法は確率的作用素法を避け、代わりにゾーン多項式の恒等式と複素積分の留数計算に依存している。
• 本結果は、実Wishartスプライクドモデルのフェーズ転移漸近的性質の最初の正確な有限N導出を提供し、先行研究が残した空白を解消した。
• 導出結果は、最大固有値がMarchenko-Pastur分布のサポートから脱出する際にフェーズ転移が発生することを確認しており、普遍的スケーリング極限が得られることを示している。