Limiting spectral distribution of large-dimensional sample covariance matrices generated by VARMA

The existence of a limiting spectral distribution (LSD) for a large-dimensional sample covariance matrix generated by the vector autoregressive moving average (VARMA) model is established. In particular, we obtain explicit forms of the LSDs for random matrices generated by a first-order vector autoregressive (VAR(1)) model and a first-order vector moving average (VMA(1)) model, as well as random coefficients for VAR(1) and VMA(1). The parameters for these explicit forms are also estimated. Finally, simulations demonstrate that the results are effective.     

The convergence on spectrum of sample covariance matrices for information-plus-noise type data

In this paper,we consider the limiting spectral distribution of the information-plusnoise type sample covariance matrices Cn =1/N (Rn + σXn) (Rn + σXn)*,under the assumption that the entries of Xn are ...     

Spectra of upper triangular operator matrices

Let be given Banach spaces. For and , let be the operator defined on by . We give sufficient conditions on to get where runs over a large class of spectra. We also discuss the case of some spectra for which the latter equality fails.

     

Some limit theorems on the eigenvectors of large dimensional sample covariance matrices

Let {v_ij} i,j = 1, 2,…, be i.i.d. standardized random variables. For each n, let V_n = (v_ij) i = 1, 2,…, n; j = 1, 2,…, s = s(n), where $\text{(}\text{n}\text{s}\text{) → y > 0}$ as n → ∞, and let $\text{M}{}_{\mathrm{n}}\text{= (}\text{1}\text{s}\text{)V}{}_{\mathrm{n}}\text{V}{}_{\mathrm{n}}{}^{\mathrm{T}}$. Previous results [7, 8] have shown the eigenvectors of M_n to display behavior, for n large, similar to those of the corresponding Wishart matrix. A certain stochastic process X_n on [0, 1], constructed from the eigenvectors of M_n, is known to converge weakly, as n → ∞, on D[0, 1] to Brownian bridge when v₁₁ is N(0, 1), but it is not known whether this property holds for any other distribution. The present paper provides evidence that this property may hold in the non-Wishart case in the form of limit theorems on the convergence in distribution of random variables constructed from integrating analytic function w.r.t. X_n(F_n(x)), where F_n is the empirical distribution function of the eigenvalues of M_n. The theorems assume certain conditions on the moments of v₁₁ including E(v₁₁⁴) = 3, the latter being necessary for the theorems to hold.     

Randomly generating portfolio-selection covariance matrices with specified distributional characteristics

In portfolio selection, there is often the need for procedures to generate “realistic” covariance matrices for security returns, for example to test and benchmark optimization algorithms. For application in portfolio optimization, such a procedure should allow the entries in the matrices to have distributional characteristics which we would consider “realistic” for security returns. Deriving motivation from the fact that a covariance matrix can be viewed as stemming from a matrix of factor loadings, a procedure is developed for the random generation of covariance matrices (a) whose off-diagonal (covariance) entries possess a pre-specified expected value and standard deviation and (b) whose main diagonal (variance) entries possess a likely different pre-specified expected value and standard deviation. The paper concludes with a discussion about the futility one would likely encounter if one simply tried to invent a valid covariance matrix in the absence of a procedure such as in this paper.     

A note on the largest eigenvalue of a large dimensional sample covariance matrix

Let {v_ij; i, J = 1, 2, …} be a family of i.i.d. random variables with E(v₁₁⁴) = ∞. For positive integers p, n with p = p(n) and p/n → y > 0 as n → ∞, let M_n = (1/n) V_n V_n^T , where V_n = (v_ij)_{1 ≤ i ≤ p, 1 ≤ j ≤ n}, and let λ_max(n) denote the largest eigenvalue of M_n. It is shown that a.s. This result verifies the boundedness of E(v₁₁⁴) to be the weakest condition known to assure the almost sure convergence of λ_max(n) for a class of sample covariance matrices.     

Two step-down tests for equality of covariance matrices

The classical problem of testing the equality of the covariance matrices from k ? 2 p-dimensional normal populations is reexamined. The likelihood ratio (LR) statistic, also called Bartlett’s statistic, can be decomposed in two ways, corresponding to two distinct component-wise decompositions of the null hypothesis in terms of the covariance matrices or precision matrices, respectively. The factors of the LR statistic that appear in these two decompositions can be interpreted as conditional and unconditional LR statistics for the component-wise null hypotheses, and their mutual independence under the null hypothesis allows the determination of the overall significance level.     

基于样本协方差矩阵的多维随机数生成方法

对于概率模型未知的多维数据样本容量扩充问题,根据主成分分析原理以及多维正态分布的性质,讨论并给出了与已知多维样本数据有相同协方差结构的模拟数据生成算法,并在此基础上给出了变量的离散化处理方法。实现了在小样本数据基础上不改变变量间协方差结构的样本容量扩充,为小样本条件下的数学建模、检验和分析提供样本数据支撑。     

The application of spectral distribution of product of two random matrices in the factor analysis

In the factor analysis model with large cross-section and time-series dimensions,we pro- pose a new method to estimate the number of factors.Specially if the idiosyncratic terms satisfy a linear time series model,the estimators of the parameters can be obtained in the time series model. The theoretical properties of the estimators are also explored.A simulation study and an empirical analysis are conducted.     

An asymptotic representation for products of random matrices

Let {X_k} be a stationary ergodic sequence of nonnegative matrices. It is shown in this paper that, under mild additional conditions, the logarithm of the i, jth element of X_t···X₁ is well approximated by a sum of t random variables from a stationary ergodic sequence. This representation is very useful for the study of limit behaviour of products of random matrices. An iterated logarithm result and an estimation result of use in the theory of demographic population projections are derived as corollaries.     

Explicit invariant measures for products of random matrices

We construct explicit invariant measures for a family of infinite products of random, independent, identically-distributed elements of SL. The matrices in the product are such that one entry is gamma-distributed along a ray in the complex plane. When the ray is the positive real axis, the products are those associated with a continued fraction studied by Letac & Seshadri [Z. Wahr. Verw. Geb. 62 (1983) 485-489], who showed that the distribution of the continued fraction is a generalised inverse Gaussian. We extend this result by finding the distribution for an arbitrary ray in the complex right-half plane, and thus compute the corresponding Lyapunov exponent explicitly. When the ray lies on the imaginary axis, the matrices in the infinite product coincide with the transfer matrices associated with a one-dimensional discrete Schrödinger operator with a random, gamma-distributed potential. Hence, the explicit knowledge of the Lyapunov exponent may be used to estimate the (exponential) rate of localisation of the eigenstates.

     

Random matrices with complex Gaussian entries

In this paper we give new and purely analytical proofs of a number of classical results on the asymptotic behavior of large random matrices of complex Wigner type (the GUE-case) or of complex Wishart type: Wigner's semi-circle law, the Harer-Zagier recursion formula, the Marchenko-Pastur law, the Geman-Silverstein results on the largest and smallest eigenvalues and other related results. Our approach is based on the derivation of explicit formulae for the moment generating functions for random matrices of the two considered types.     

The rank of sparse random matrices

We determine the asymptotic normalized rank of a random matrix $\mathit{A}$ over an arbitrary field with prescribed numbers of nonzero entries in each row and column. As an application we obtain a formula for the rate of low-density parity check codes. This formula vindicates a conjecture of Lelarge (2013). The proofs are based on coupling arguments and a novel random perturbation, applicable to any matrix, that diminishes the number of short linear relations.     

Deformations of the infinite-dimensional Lie algebra L 3

For the nilpotent infinite-dimensional Lie algebra L ₃, we compute the second cohomology group H ²(L ₃, L ₃) with coefficients in the adjoint module. Nontrivial cocycles are found in closed form, and Massey powers are computed for them.     

Some limit theorems for the eigenvalues of a sample covariance matrix

Limit theorems are given for the eigenvalues of a sample covariance matrix when the dimension of the matrix as well as the sample size tend to infinity. The limit of the cumulative distribution function of the eigenvalues is determined by use of a method of moments. The proof is mainly combinatorial. By a variant of the method of moments it is shown that the sum of the eigenvalues, raised to k-th power, k = 1, 2,…, m is asymptotically normal. A limit theorem for the log sum of the eigenvalues is completed with estimates of expected value and variance and with bounds of Berry-Esseen type.     

A limit theorem for the eigenvalues of product of two random matrices

The existence of limit spectral distribution of the product of two independent random matrices is proved when the number of variables tends to infinity. One of the above matrices is the Wishart matrix and the other is a symmetric nonnegative definite matrix.     

Circular law for random matrices with exchangeable entries

An exchangeable random matrix is a random matrix with distribution invariant under any permutation of the entries. For such random matrices, we show, as the dimension tends to infinity, that the empirical spectral distribution tends to the uniform law on the unit disc. This is an instance of the universality phenomenon known as the circular law, for a model of random matrices with dependent entries, rows, and columns. It is also a non‐Hermitian counterpart of a result of Chatterjee on the semi‐circular law for random Hermitian matrices with exchangeable entries. The proof relies in particular on a reduction to a simpler model given by a random shuffle of a rigid deterministic matrix, on hermitization, and also on combinatorial concentration of measure and combinatorial Central Limit Theorem. A crucial step is a polynomial bound on the smallest singular value of exchangeable random matrices, which may be of independent interest. © 2015 Wiley Periodicals, Inc. Random Struct. Alg., 48, 454–479, 2016     

On the principal components of sample covariance matrices

We introduce a class of \(M \times M\) sample covariance matrices \({\mathcal {Q}}\) which subsumes and generalizes several previous models. The associated population covariance matrix \(\Sigma = \mathbb {E}{\mathcal {Q}}\) is assumed to differ from the identity by a matrix of bounded rank. All quantities except the rank of \(\Sigma - I_M\) may depend on \(M\) in an arbitrary fashion. We investigate the principal components, i.e. the top eigenvalues and eigenvectors, of \({\mathcal {Q}}\). We derive precise large deviation estimates on the generalized components \(\langle {\mathbf{{w}}} , {\varvec{\xi }_i}\rangle \) of the outlier and non-outlier eigenvectors \(\varvec{\xi }_i\). Our results also hold near the so-called BBP transition, where outliers are created or annihilated, and for degenerate or near-degenerate outliers. We believe the obtained rates of convergence to be optimal. In addition, we derive the asymptotic distribution of the generalized components of the non-outlier eigenvectors. A novel observation arising from our results is that, unlike the eigenvalues, the eigenvectors of the principal components contain information about the subcritical spikes of \(\Sigma \). The proofs use several results on the eigenvalues and eigenvectors of the uncorrelated matrix \({\mathcal {Q}}\), satisfying \(\mathbb {E}{\mathcal {Q}} = I_M\), as input: the isotropic local Marchenko–Pastur law established in Bloemendal et al. (Electron J Probab 19:1–53, 2014), level repulsion, and quantum unique ergodicity of the eigenvectors. The latter is a special case of a new universality result for the joint eigenvalue–eigenvector distribution.