We study the spectral properties of the Jacobian matrix in artificial neural networks and establish both its limiting spectral distribution and the second-order fluctuations of its spectrum. Our analysis uncovers fundamental phenomena, including universality properties with respect to the weight distribution and a phase transition in spectral behavior driven by the choice of activation function. Leveraging tools from random matrix theory, we first analyze a class of…
Zeqin Lin · Yiming Liu · Guangming Pan · Chi Yao · Jia Zhou·刊期 2026-07-06 · 卷 113(3)·Biometrika已编入正刊
识别谱尖峰的来源:因子还是混合?AI 解读
Abstract We consider the problem of identifying the pattern of latent variables in high-dimensional linear latent variable models, which can also be interpreted as determining the source of spiked singular values in the data matrix. Specifically, we test whether the latent variables are continuous or categorical, a distinction that is crucial for data interpretation, but challenging in the high-dimensional regime. To address this inference problem, we analyse the asymptotic…
Abstract We consider the problem of estimating the number of significant components in high-dimensional principal component analysis. We propose a new penalized approach using the explained variance ratio and the rigidity of the nonspiked sample eigenvalues of sample covariance matrices of $ p $ variables. Compared with methods in the existing literature, the consistency of the proposed estimator holds, not only for independent data, but also for some times series data when…
Consider a data matrix $Y = [\mathbf{y}_1, \cdots, \mathbf{y}_N]$ of size $M \times N$, where the columns are independent observations from a random vector $\mathbf{y}$ with zero mean and population covariance $Σ$. Let $\mathbf{u}_i$ and $\mathbf{v}_j$ denote the left and right singular vectors of $Y$, respectively. This study investigates the eigenvector/singular vector overlaps $\langle {\mathbf{u}_i, D_1 \mathbf{u}_j} \rangle$, $\langle {\mathbf{v}_i, D_2 \mathbf{v}_j} \rangle$…
Cross-sectional structures and temporal tendency are important features of high-dimensional time series. Based on eigen-analysis on sample covariance matrices, we propose a novel approach to identifying four popular structures of high-dimensional time series, which are grouped in terms of factor structures and stationarity. The proposed three-step method includes a ratio statistic of empirical eigenvalues;a projected Augmented Dickey-Fuller Test;a new unit-root test based on…
Huiqin Li · Guangming Pan · Yanqing Yin · Wang Zhou·刊期 2024-06-01 · 卷 52(3)·The Annals of Statistics已编入正刊
随机缺失观测下 Gram 矩阵的谱分析与统计推断AI 解读
Motivated by the statistical inference using the Gram matrix in the context of missing at random observations, this paper investigates the spectral properties of the random matrices Sn=1nZZ∗, where Z=D∘(Σ1/2X) represents a Hadamard random matrix with entries determined by independent Bernoulli variables D. Operating within the high-dimensional framework, we establish the convergence of the empirical spectral distribution of Sn to a well-defined limiting distribution. In…
We propose a new unsupervised learning method for clustering a large number of time series based on a latent factor structure.Each cluster is characterized by its own cluster-specific factors in addition to some common factors which impact on all the time series concerned.Our setting also offers the flexibility that some time series may not belong to any clusters.The consistency with explicit convergence rates is established for the estimation of the common factors, the…
Zhixiang Zhang · Shurong Zheng · Guangming Pan · Ping-Shou Zhong·刊期 2022-08-01 · 卷 50(4)·The Annals of Statistics已编入正刊
大样本协方差矩阵尖峰特征值与线性谱统计量的渐近独立性AI 解读
We consider general high-dimensional spiked sample covariance models and show that their leading sample spiked eigenvalues and their linear spectral statistics are asymptotically independent when the sample size and dimension are proportional to each other. As a byproduct, we also establish the central limit theorem of the leading sample spiked eigenvalues by removing the block diagonal assumption on the population covariance matrix, which is commonly needed in the…
This article considers the problem of estimating a change point in the covariance matrix in a sequence of high-dimensional vectors, where the dimension is substantially larger than the sample size. A two-stage approach is proposed to efficiently estimate the location of the change point. The first step consists of a reduction of the dimension to identify elements of the covariance matrices corresponding to significant changes. In a second step, we use the components after…
We study the asymptotic distributions of the spiked eigenvalues and the largest nonspiked eigenvalue of the sample covariance matrix under a general covariance model with divergent spiked eigenvalues, while the other eigenvalues are bounded but otherwise arbitrary. The limiting normal distribution for the spiked sample eigenvalues is established. It has distinct features that the asymptotic mean relies on not only the population spikes but also the nonspikes and that the…
Qingliang Fan · Xiao Han · Guangming Pan · Bibo Jiang·刊期 2020-06 · 卷 36(3)·Econometric Theory已编入正刊
大型似不相关回归系统:惩罚准极大似然估计AI 解读
In this article, using a shrinkage estimator, we propose a penalized quasi-maximum likelihood estimator (PQMLE) to estimate a large system of equations in seemingly unrelated regression models, where the number of equations is large relative to the sample size. We develop the asymptotic properties of the PQMLE for both the error covariance matrix and model coefficients. In particular, we derive the asymptotic distribution of the coefficient estimator and the convergence rate…
Zhigang Bao · Jiang Hu · Guangming Pan · Wang Zhou·刊期 2019-02-01 · 卷 47(1)·The Annals of Statistics已编入正刊
高维高斯向量的典型相关系数:有限秩情形AI 解读
Consider a Gaussian vector $\mathbf{z}=(\mathbf{x}',\mathbf{y}')'$, consisting of two sub-vectors $\mathbf{x}$ and $\mathbf{y}$ with dimensions $p$ and $q$, respectively. With $n$ independent observations of $\mathbf{z}$, we study the correlation between $\mathbf{x}$ and $\mathbf{y}$, from the perspective of the canonical correlation analysis. We investigate the high-dimensional case: both $p$ and $q$ are proportional to the sample size $n$. Denote by $\Sigma_{uv}$ the…
Let $\{Z_{ij}\}$ be independent and identically distributed (i.i.d.) random variables with $EZ_{ij}=0$, $E\vert Z_{ij}\vert^{2}=1$ and $E\vert Z_{ij}\vert^{4}<\infty$. Define linear processes $Y_{tj}=\sum_{k=0}^{\infty}b_{k}Z_{t-k,j}$ with $\sum_{i=0}^{\infty}\vert b_{i}\vert <\infty$. Consider a $p$-dimensional time series model of the form $\mathbf{x}_{t}=\boldsymbol{\Pi} \mathbf{x}_{t-1}+\Sigma^{1/2}\mathbf{y}_{t},\ 1\leq t\leq T$ with $\mathbf{y}_{t}=(Y_{t1},\ldots,Y_{tp})'$…
Summary Statistical inferences for sample correlation matrices are important in high dimensional data analysis. Motivated by this, the paper establishes a new central limit theorem for a linear spectral statistic of high dimensional sample correlation matrices for the case where the dimension p and the sample size n are comparable. This result is of independent interest in large dimensional random-matrix theory. We also further investigate the sample correlation matrices of…
This article considers testing equality of two population covariance matrices when the data dimension p diverges with the sample size n (p/n → c > 0). We propose a weighted test statistic that is data-driven and powerful in both faint alternatives (many small disturbances) and sparse alternatives (several large disturbances). Its asymptotic null distribution is derived by large random matrix theory without assuming the existence of a limiting cumulative distribution function…
The article is concerned with empirical Bayes shrinkage estimators for the heteroscedastic hierarchical normal model using Stein's unbiased estimate of risk (SURE). Recently, Xie, Kou, and Brown proposed a class of estimators for this type of problems and established their asymptotic optimality properties under the assumption of known but unequal variances. In this article, we consider this problem with unequal and unknown variances, which may be more appropriate in real…
Let ${\mathbf{A}}_{p}=\frac{{\mathbf{Y}}{\mathbf{Y}}^{*}}{m}$ and ${\mathbf{B}}_{p}=\frac{{\mathbf{X}}{\mathbf{X}}^{*}}{n}$ be two independent random matrices where ${\mathbf{X}}=(X_{ij})_{p\times n}$ and ${\mathbf{Y}}=(Y_{ij})_{p\times m}$ respectively consist of real (or complex) independent random variables with $\mathbb{E}X_{ij}=\mathbb{E}Y_{ij}=0$, $\mathbb{E}|X_{ij}|^{2}=\mathbb{E}|Y_{ij}|^{2}=1$. Denote by $\lambda_{1}$ the largest root of the determinantal equation $\det(\lambda{\mathbf{A}}_{p}-{\mathbf{B}}_{p})=0$…
Zhigang Bao · Liang-Ching Lin · Guangming Pan · Wang Zhou·刊期 2015-12-01 · 卷 43(6)·The Annals of Statistics已编入正刊
高维 Spearman 秩相关矩阵的谱统计及其应用AI 解读
Let $\mathbf{Q}=(Q_{1},\ldots,Q_{n})$ be a random vector drawn from the uniform distribution on the set of all $n!$ permutations of $\{1,2,\ldots,n\}$. Let $\mathbf{Z}=(Z_{1},\ldots,Z_{n})$, where $Z_{j}$ is the mean zero variance one random variable obtained by centralizing and normalizing $Q_{j}$, $j=1,\ldots,n$. Assume that $\mathbf{X}_{i},i=1,\ldots,p$ are i.i.d. copies of $\frac{1}{\sqrt{p}}\mathbf{Z}$ and $X=X_{p,n}$ is the $p\times n$ random matrix with $\mathbf{X}_{i}$…
This paper proposes a new statistic to test independence between two high dimensional random vectors $\mathbf{X}:p_{1}\times1$ and $\mathbf{Y}:p_{2}\times1$. The proposed statistic is based on the sum of regularized sample canonical correlation coefficients of $\mathbf{X}$ and $\mathbf{Y}$. The asymptotic distribution of the statistic under the null hypothesis is established as a corollary of general central limit theorems (CLT) for the linear statistics of classical and…
This paper is aimed at deriving the universality of the largest eigenvalue of a class of high-dimensional real or complex sample covariance matrices of the form $\mathcal{W}_{N}=\Sigma^{1/2}XX^{*}\Sigma^{1/2}$. Here, $X=(x_{ij})_{M,N}$ is an $M\times N$ random matrix with independent entries $x_{ij}$, $1\leq i\leq M$, $1\leq j\leq N$ such that $\mathbb{E}x_{ij}=0$, $\mathbb{E}|x_{ij}|^{2}=1/N$. On dimensionality, we assume that $M=M(N)$ and $N/M\rightarrow d\in(0,\infty)$ as…
Capturing dependence among a large number of high-dimensional random vectors is a very important and challenging problem. By arranging n random vectors of length p in the form of a matrix, we develop a linear spectral statistic of the constructed matrix to test whether the n random vectors are independent or not. Specifically, the proposed statistic can also be applied to n random vectors, each of whose elements can be written as either a linear stationary process or a…
Bing-Yi Jing · Guangming Pan · Qi-Man Shao · Wang Zhou·刊期 2010-12-01 · 卷 38(6)·The Annals of Statistics已编入正刊
样本协方差矩阵谱密度函数的非参数估计:初步研究AI 解读
The density function of the limiting spectral distribution of general sample covariance matrices is usually unknown. We propose to use kernel estimators which are proved to be consistent. A simulation study is also conducted to show the performance of the estimators.
Abstract We consider the problem of estimating the number of significant components in high-dimensional principal component analysis. We propose a new penalized approach using the explained variance ratio and the rigidity of the nonspiked sample eigenvalues of sample covariance matrices of p variables. Compared with methods in the existing literature, the consistency of the proposed estimator holds, not only for independent data, but also for some times series data when…
Consider a data matrix Y=[y1,⋯,yN] of size M×N, where the columns are independent observations from a random vector y with zero mean and population covariance Σ. Let ui and vj denote the left and right singular vectors of Y, respectively. This study investigates the eigenvector/singular vector overlaps ⟨ui,D1uj⟩, ⟨vi,D2vj⟩…
Consider a Gaussian vector z=(x′,y′)′, consisting of two sub-vectors x and y with dimensions p and q, respectively. With n independent observations of z, we study the correlation between x and y, from the perspective of the canonical correlation analysis. We investigate the high-dimensional case: both p and q are proportional to the sample size n. Denote by Σuv the…
Let {Zij} be independent and identically distributed (i.i.d.) random variables with EZij=0, E∣Zij∣2=1 and E∣Zij∣4<∞. Define linear processes Ytj=∑k=0∞bkZt−k,j with ∑i=0∞∣bi∣<∞. Consider a p-dimensional time series model of the form xt=Πxt−1+Σ1/2yt,1≤t≤T with yt=(Yt1,…,Ytp)′…
Let Ap=mYY∗ and Bp=nXX∗ be two independent random matrices where X=(Xij)p×n and Y=(Yij)p×m respectively consist of real (or complex) independent random variables with EXij=EYij=0, E∣Xij∣2=E∣Yij∣2=1. Denote by λ1 the largest root of the determinantal equation det(λAp−Bp)=0…
Let Q=(Q1,…,Qn) be a random vector drawn from the uniform distribution on the set of all n! permutations of {1,2,…,n}. Let Z=(Z1,…,Zn), where Zj is the mean zero variance one random variable obtained by centralizing and normalizing Qj, j=1,…,n. Assume that Xi,i=1,…,p are i.i.d. copies of p1Z and X=Xp,n is the p×n random matrix with Xi…
This paper proposes a new statistic to test independence between two high dimensional random vectors X:p1×1 and Y:p2×1. The proposed statistic is based on the sum of regularized sample canonical correlation coefficients of X and Y. The asymptotic distribution of the statistic under the null hypothesis is established as a corollary of general central limit theorems (CLT) for the linear statistics of classical and…
This paper is aimed at deriving the universality of the largest eigenvalue of a class of high-dimensional real or complex sample covariance matrices of the form WN=Σ1/2XX∗Σ1/2. Here, X=(xij)M,N is an M×N random matrix with independent entries xij, 1≤i≤M, 1≤j≤N such that Exij=0, E∣xij∣2=1/N. On dimensionality, we assume that M=M(N) and N/M→d∈(0,∞) as…