中文

基于扩展AIC型方法的主成分个数检测

统计方法学 2020-09-01 v7

摘要

估计主成分个数是信号处理(或尖峰协方差模型)等许多科学领域中的基本问题之一。在本文中,我们首先证明,对于固定的pp,任何形式为k(pk/2+1/2)Cnk'(p-k'/2+1/2)C_n的惩罚项在条件CnC_n\to\inftyCn/n0C_n/n\to0下可能导致渐近一致估计量。我们还将结果扩展到n,pn,p\to\inftyp/nc>0p/n\to c>0的情形。在此情形下,对于k=o(n13)k=o(n^{\frac{1}{3}}),我们首先在条件λk>1+c\lambda_k>1+\sqrt{c}下研究样本协方差矩阵SnS_n的前导特征值的极限定律。在低信噪比下,由于AIC倾向于低估信号数kk,AIC在此情形下应被重新定义。作为固定pp时AIC的自然扩展,我们提出扩展AIC(EAIC),即带调谐参数γ=φ(c)=1/2+1/clog(1+c)/c\gamma=\varphi(c)=1/2+\sqrt{1/c}-\log(1+\sqrt{c})/c的AIC型方法,并证明EAIC型方法(即带调谐参数γ>φ(c)\gamma>\varphi(c)的AIC型方法)能一致地选择信号数kk。在以下两种情形中,(1) pp固定,nn\to\infty,(2) n,pn,p\to\inftyp/n0p/n\to 0,若AIC定义为n,pn,p\to\inftyp/nc>0p/n\to c>0情形下EAIC的退化,即γ=limc0+0φ(c)=1\gamma=\lim_{c\rightarrow 0+0}\varphi(c)=1,则我们实质上已证明,为使上述两种情形中AIC型方法的一致性,需要γ>1\gamma>1。此外,我们表明EAIC型方法本质上是免调谐的,并且优于Kritchman和Nadler (2008)提出的著名KN估计量以及Bai、Choi和Fujikoshi (2018)提出的BCF估计量。数值研究表明所提方法表现良好。

关键词

引用

@article{arxiv.1708.07595,
  title  = {Detection of the number of principal components by extended AIC-type method},
  author = {Jianwei Hu and Jingfei Zhang and Ji Zhu},
  journal= {arXiv preprint arXiv:1708.07595},
  year   = {2020}
}

备注

56 pages; AIC is re-defined; for $k=o(n^{\frac{1}{3}})$, we investigate the limiting laws for the leading eigenvalues of the sample covariance matrix $S_n$ under the condition that $\lambda_k>1+\sqrt{c}$