[論文レビュー] A CLT in Stein's distance for generalized Wishart matrices and higher order tensors

Dan Mikulincer|arXiv (Cornell University)|Feb 25, 2020

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

ひとこと要約

本稿は、一般化ワイシャールト行列および高次テンソルに対するスティンの距離における中心極限定理（CLT）を確立し、i.i.d. な $p$-次のテンソル累乗の和が $n^{2p-1} \ll d$ の下でガウス分布に収束することを示している。正則性条件のもとで、本研究の主な革新は、テンソル累乗の低次元構造を活用するため、スティン法における最適輸送の新規応用であり、標準的なスティン法に比べて収束バインディングをより厳密に得ている。

ABSTRACT

We study the central limit theorem for sums of independent tensor powers, $\frac{1}{\sqrt{d}}\sum\limits_{i=1}^d X_i^{\otimes p}$. We focus on the high-dimensional regime where $X_i \in \mathbb{R}^n$ and $n$ may scale with $d$. Our main result is a proposed threshold for convergence. Specifically, we show that, under some regularity assumption, if $n^{2p-1}\gg d$, then the normalized sum converges to a Gaussian. The results apply, among others, to symmetric uniform log-concave measures and to product measures. This generalizes several results found in the literature. Our main technique is a novel application of optimal transport to Stein's method which accounts for the low dimensional structure which is inherent in $X_i^{\otimes p}$.

研究の動機と目的

• 高次元CLTを、$\mathbb{R}^n$ 内のi.i.d. ランダムベクトルの $p$-次テンソル累乗からなる一般化ワイシャールト行列へ拡張すること。
• 積測度または対数凹型測度を仮定する従来の制約を緩和し、より広いクラスの等方的測度を許容すること。
• スティンの距離における収束の鋭い閾値を確立し、$p=2$ の既知の結果を $p>2$ へ一般化すること。
• テンソル累乗に内在する低次元構造を考慮する手法を開発し、収束速度バインディングを改善すること。

提案手法

• スティン法の枠組み内で最適輸送理論を新規に応用し、テンソル累乗の和とガウス分布との距離を分析すること。
• 交換可能なペアと勾配に基づくバインディングを用いたスティン法により、スティン散らばんを制御すること。
• テンソル累乗 $X_i^{\otimes p}$ の低次元構造を組み込む修正されたスティン作用素の導入により、問題の有効次元を低減すること。
• 輸送写像の条件付き期待値のヒルベルト＝シュミットノルムによる推定を用いて、スティン解の2次モーメントのバインディングを導出すること。
• タラグランの輸送-エントロピー不等式およびピンスカーリンの不等式を適用し、相対エントロピーとワーサーテイン距離を関連させること。
• 一様対数凹型および積測度への特殊化を行い、$n$、$d$、$p$ の関数としての収束速度の明示的バインディングを導出すること。

実験結果

リサーチクエスチョン

• RQ1 高次元において $n$ が $d$ とスケーリングするとき、正規化された $p$-次テンソル累乗のガウス分布への収束の鋭い閾値は何か？
• RQ2 ワイシャールト行列のCLTを、積測度および対数凹型測度を越えて、より一般の等方的分布へ拡張できるか？
• RQ3 最適輸送をスティン法に統合することで、低次元構造を持つランダムテンソルの収束バインディングをどのように改善できるか？
• RQ4 テンソル和に対角成分を含めると、CLTの閾値に影響を与えるか？もし与えるならば、その影響は何か？
• RQ5 定常ガウス過程のような、従属するエントリを有する基本的なランダム行列に対しても、この手法を適合可能か？

主な発見

• 本稿は、$n^{2p-1} \ll d$ ならば、測度 $\mu$ に対してやや強い正則性仮定のもとで、正規化された和 $\frac{1}{\sqrt{d}}\sum_{i=1}^d X_i^{\otimes p}$ がスティンの距離においてガウス分布に収束することを確立している。
• 閾値 $n^{2p-1} \ll d$ は、$p=2$ のときの既知の $n^3 \ll d$ の閾値を一般化し、高次テンソルへCLTを拡張している。
• 等方的 $L$-一様対数凹型かつ無条件測度に対しては、スティン散らばんのバインディングが $\frac{C_p}{L^4} n^{2p-1} \frac{\|\alpha\|_4^4}{\|\alpha\|_2^4}$ と表され、$C_p$ は $p$ のみに依存する。
• テンソル累乗 $X_i^{\otimes p}$ が次元 $\binom{n+p-1}{p}$ の空間に存在するが、$\mathbb{R}^n$ からのランダムネスを引き継ぐという低次元構造を活用することで、標準的なスティン法よりもより厳密なバインディングが達成されている。
• 基礎となる行列 $\mathbb{X}$ が従属する列を有する場合（例：定常ガウス過程）、直交的対角化による変数変換により、結果は依然として成り立つ。
• フレームワークは非一様和に対しても適用可能であり、一般の $p$-次テンソル和のためのスティン解の2次モーメントのバインディングを提供し、定量的収束速度を可能にしている。