[論文レビュー] On the eigenvalues of the spatial sign covariance matrix in more than two dimensions

Alexander Dürre, David E. Tyler+1 more|arXiv (Cornell University)|Dec 9, 2015

Morphological variations and asymmetry参考文献 22被引用数 21

ひとこと要約

この論文は、楕円分布における空間符号共分散行列（SSCM）の固有値が、形状行列の固有値の1対1関数であり、元の固有値よりもより緊密に集中していることを確立している。また、SSCM固有値の1次元積分表現を導出し、効率的な数値計算を可能にするとともに、p > 2次元における楕円対称性のもとで、固有値の単射性、順序の保存、固有値の集中性を証明している。

ABSTRACT

We gather several results on the eigenvalues of the spatial sign covariance matrix of an elliptical distribution. It is shown that the eigenvalues are a one-to-one function of the eigenvalues of the shape matrix and that they are closer together than the latter. We further provide a one-dimensional integral representation of the eigenvalues, which facilitates their numerical computation.

研究の動機と目的

• 高次元の楕円分布における空間符号共分散行列（SSCM）の固有値と形状行列の関係を調査すること。
• SSCM固有値が形状行列固有値の1対1関数であることを確立し、同定可能性を保証すること。
• SSCM固有値が形状行列固有値よりもより緊密に集中している（近接している）ことを示すこと。
• SSCM固有値の1次元積分表現を導出し、効率的な数値計算を可能にすること。
• 高次元における空間符号法を用いたロバストな多変量解析の理論的基盤を提供すること。

提案手法

• 著者らはトレース正規化された形状行列 $ V_0 = O\Lambda O^\top $ の固有値分解を用い、ここで $ \Lambda = \text{diag}(\lambda_1, \dots, \lambda_p) $ であり、$ \lambda_1 \geq \cdots \geq \lambda_p \geq 0 $ かつ $ \sum \lambda_j = 1 $ を満たす。
• 母集団のSSCMを $ S = O\Delta O^\top $ と表現し、$ \Delta = \text{diag}(\delta_1, \dots, \delta_p) $ とし、$ \delta_i = \mathbb{E}\left[ \lambda_i Y_i^2 \left( \sum_{j=1}^p \lambda_j Y_j^2 \right)^{-1} \right] $ とする。ここで $ Y = RU $ は球面分布に従う。
• 解析的取り扱いの容易さを図るため、著者らは $ Y $ の球対称性を活用し、便利な分布（単位球内の一様分布）を選択して期待値積分を簡略化する。
• GradshteynとRyzhikの公式4.646を適用し、固有値 $ \delta_i $、$ \eta_{ii} $、$ \eta_{ij} $ の1次元積分表現を導出し、数値評価を可能にする。
• 単射性の証明は背理法を用い、球対称性の仮定のもとで、異なる $ \lambda $ ベクトルが異なる $ \delta $ ベクトルをもたらすことを示した。
• 固有値の集中性の証明は、重み付き調和平均に基づく不等式を用い、$ \rho_{(1)} < \rho_{(p)} $ の厳密な不等式を適用することで、$ \delta_i $ が $ \lambda_i $ よりも厳密に等しくなることを示した。

実験結果

リサーチクエスチョン

• RQ1 3次元以上において、形状行列固有値からSSCM固有値への写像は単射的か？
• RQ2 SSCM固有値の分散は形状行列固有値の分散と比べてどうか？
• RQ3 SSCM固有値を1次元積分で表現できるか？ これにより効率的な数値計算が可能か？
• RQ4 楕円対称性下で、SSCM固有値と形状行列固有値の関数的関係は何か？
• RQ5 SSCM固有値ベクトルは形状行列固有値の順序を保存するか？

主な発見

• 形状行列固有値 $ \lambda $ からSSCM固有値 $ \delta $ への写像は単射的であり、$ \phi(\lambda) = \phi(\tilde{\lambda}) $ ならば $ \lambda = \tilde{\lambda} $ が成り立つ。
• SSCM固有値 $ \delta_1, \dots, \delta_p $ は形状行列固有値 $ \lambda_1, \dots, \lambda_p $ よりもより緊密に集中しており、すなわち $ \delta_1 - \delta_p < \lambda_1 - \lambda_p $ が成り立つ。
• SSCM固有値 $ \delta_i $ は形状行列固有値の1対1関数であり、順序を保存する：$ \delta_1 \geq \cdots \geq \delta_p $。
• 単位球内の一様分布とGradshteynとRyzhikの公式4.646を用いて、$ \delta_i $ の1次元積分表現を導出し、数値計算を可能にした。
• SSCM推定量の漸近的共分散行列は $ \delta_i $ と $ \eta_{ij} $ で表現され、$ \eta_{ij} $ に対しても1次元積分表現が与えられている。
• 固有値 $ \delta_i $ は $ \delta_i = \frac{p\Gamma(p/2)}{2\pi^{p/2}} \int_{S_{1,p}} \frac{\lambda_i z_i}{\sum_{j=1}^p \lambda_j z_j} \prod_{j=1}^p z_j^{-1/2} dz $ と表現され、ここで $ S_{1,p} = \{ z \in \mathbb{R}^p_+ \mid \sum z_j \leq 1 \} $ である。