Covariance estimation

Esp Gaussian

November 17, 2014 β€” April 26, 2023

Figure 1

Estimating the thing that is given to you by oracles in statistics homework assignments: the covariance matrix. Or, if the data is indexed by some parameter we might consider the covariance kernel. We are especially interested in this in Gaussian processes, where the covariance structure characterises the process up to its mean.

I am not introducing a complete theory of covariance estimation here, merely some notes.

Two big data problems problems can arise here: large \(p\) (ambient dimension) and large \(n\) (sample size). Large \(p\) is a problem because the covariance matrix is a \(p \times p\) matrix and frequently we need to invert it to calculate some target estimand.

Often life can be made not too bad for large \(n\) with Gaussian structure because, essentially, the problem has nice nearly low rank structure.

1 Bayesian

Inverse Wishart priors. πŸ— Other?

2 Precision estimation

The workhorse of learning graphical models under linearity and Gaussianity. See precision estimation for a more complete treatment.

3 Continuous

See kernel learning.

4 Parametric

4.1 Cholesky methods

(Huang et al. 2006; Wu and Pourahmadi 2003).

4.2 on a lattice

Estimating a stationary covariance function on a regular lattice? That is a whole field of its own. Useful keywords include circulant embedding. Although strictly more general than Gaussian processes on a lattice, it is often used in that context and some extra results are on that page for now.

5 Unordered

Thanks to Rothman (2010) I now think about covariance estimates as being different in ordered versus exchangeable data.

6 Sandwich estimators

Figure 2

For robust covariances of vector data. AKA Heteroskedasticity-consistent covariance estimators. Incorporating Eicker-Huber-White sandwich estimator, Andrews kernel HAC estimator, Newey-West and others. For an intro see Achim Zeileis, Open-Source Econometric Computing in R.

7 Incoming

8 Bounding by harmonic and arithmetic means

There are some known bounds for the univariate case. Wikipedia says, in Relations with the harmonic and arithmetic means that it has been shown (Mercer 2000) that for a sample \(\left\{y_i\right\}\) of positive real numbers, \[ \sigma_y^2 \leq 2 y_{\max }(A-H) \] where \(y_{\max }\) is the maximum of the sample, \(A\) is the arithmetic mean, \(H\) is the harmonic mean of the sample and \(\sigma_y^2\) is the (biased) variance of the sample. This bound has been improved, and it is known that variance is bounded by \[ \begin{gathered} \sigma_y^2 \leq \frac{y_{\max }(A-H)\left(y_{\max }-A\right)}{y_{\max }-H}, \\ \sigma_y^2 \geq \frac{y_{\min }(A-H)\left(A-y_{\min }\right)}{H-y_{\min }}, \end{gathered} \] where \(y_{\min }\) is the minimum of the sample (Sharma 2008).

Mond and Pec̆arić (1996) says

Let us define the arithmetic mean of \(A\) with weight \(w\) as \[ A_n(A ; w)=\sum_{i=1}^n w_i A_i \] and the harmonic mean of \(A\) with weight \(w\) as \[ H_n(A ; w)=\left(\sum_{i=1}^n w_i A_i^{-1}\right)^{-1} \] It is well known \([2,5]\) that \[ H_n(A ; w) \leqslant A_n(A ; w) \] Moreover, if \(A_{i j}(i, j=1, \ldots, n)\) are positive definite matrices from \(H_m\), then the following inequality is also valid [1]: \[ \frac{1}{n} \sum_{j=1}^n\left(\frac{1}{n} \sum_{i=1}^n A_{i j}^{-1}\right)^{-1} \leqslant\left[\frac{1}{n} \sum_{i=1}^n\left(\frac{1}{n} \sum_{j=1}^n A_{i j}\right)^{-1}\right]^{-1} \]

For multivariate covariance we are interested in the PSD matrix version of this.

9 References

Abrahamsen. 1997. β€œA Review of Gaussian Random Fields and Correlation Functions.”
Anderson. 2007. β€œExploring the Need for Localization in Ensemble Data Assimilation Using a Hierarchical Ensemble Filter.” Physica D: Nonlinear Phenomena, Data Assimilation,.
Azizyan, Krishnamurthy, and Singh. 2015. β€œExtreme Compressive Sampling for Covariance Estimation.” arXiv:1506.00898 [Cs, Math, Stat].
Baik, Arous, and PΓ©chΓ©. 2005. β€œPhase Transition of the Largest Eigenvalue for Nonnull Complex Sample Covariance Matrices.” The Annals of Probability.
Banerjee, Ghaoui, and d’Aspremont. 2008. β€œModel Selection Through Sparse Maximum Likelihood Estimation for Multivariate Gaussian or Binary Data.” Journal of Machine Learning Research.
Barnard, McCulloch, and Meng. 2000. β€œModeling Covariance Matrices in Terms of Standard Deviations and Correlations, with Application to Shrinkage.” Statistica Sinica.
Ben Arous, and PΓ©chΓ©. 2005. β€œUniversality of Local Eigenvalue Statistics for Some Sample Covariance Matrices.” Communications on Pure and Applied Mathematics.
Bickel, and Levina. 2008. β€œRegularized Estimation of Large Covariance Matrices.” The Annals of Statistics.
Bosq. 2002. β€œEstimation of Mean and Covariance Operator of Autoregressive Processes in Banach Spaces.” Statistical Inference for Stochastic Processes.
Cai, Zhang, and Zhou. 2010. β€œOptimal Rates of Convergence for Covariance Matrix Estimation.” The Annals of Statistics.
Chan, Golub, and Leveque. 1983. β€œAlgorithms for Computing the Sample Variance: Analysis and Recommendations.” The American Statistician.
Chen, Xiaohui, Xu, and Wu. 2013. β€œCovariance and Precision Matrix Estimation for High-Dimensional Time Series.” The Annals of Statistics.
Chen, Hao, Zheng, Al Kontar, et al. 2020. β€œStochastic Gradient Descent in Correlated Settings: A Study on Gaussian Processes.” In Proceedings of the 34th International Conference on Neural Information Processing Systems. NIPS’20.
Cook. 2018. β€œPrincipal Components, Sufficient Dimension Reduction, and Envelopes.” Annual Review of Statistics and Its Application.
Cunningham, Shenoy, and Sahani. 2008. β€œFast Gaussian Process Methods for Point Process Intensity Estimation.” In Proceedings of the 25th International Conference on Machine Learning. ICML ’08.
Damian, Sampson, and Guttorp. 2001. β€œBayesian Estimation of Semi-Parametric Non-Stationary Spatial Covariance Structures.” Environmetrics.
Daniels, and Pourahmadi. 2009. β€œModeling Covariance Matrices via Partial Autocorrelations.” Journal of Multivariate Analysis.
Dasgupta, and Hsu. 2007. β€œOn-Line Estimation with the Multivariate Gaussian Distribution.” In Learning Theory.
Efron. 2010. β€œCorrelated z-Values and the Accuracy of Large-Scale Statistical Estimates.” Journal of the American Statistical Association.
Fan, Liao, and Liu. 2016. β€œAn Overview of the Estimation of Large Covariance and Precision Matrices.” The Econometrics Journal.
Friedman, Hastie, and Tibshirani. 2008. β€œSparse Inverse Covariance Estimation with the Graphical Lasso.” Biostatistics.
Fuentes. 2006. β€œTesting for Separability of Spatial–Temporal Covariance Functions.” Journal of Statistical Planning and Inference.
Furrer, R., and Bengtsson. 2007. β€œEstimation of high-dimensional prior and posterior covariance matrices in Kalman filter variants.” Journal of Multivariate Analysis.
Furrer, Reinhard, Genton, and Nychka. 2006. β€œCovariance Tapering for Interpolation of Large Spatial Datasets.” Journal of Computational and Graphical Statistics.
Gneiting, Kleiber, and Schlather. 2010. β€œMatΓ©rn Cross-Covariance Functions for Multivariate Random Fields.” Journal of the American Statistical Association.
Goodman. 1960. β€œOn the Exact Variance of Products.” Journal of the American Statistical Association.
Hackbusch. 2015. Hierarchical Matrices: Algorithms and Analysis. Springer Series in Computational Mathematics 49.
Hansen. 2007. β€œGeneralized Least Squares Inference in Panel and Multilevel Models with Serial Correlation and Fixed Effects.” Journal of Econometrics.
Heinrich, and Podolskij. 2014. β€œOn Spectral Distribution of High Dimensional Covariation Matrices.” arXiv:1410.6764 [Math].
Huang, Liu, Pourahmadi, et al. 2006. β€œCovariance Matrix Selection and Estimation via Penalised Normal Likelihood.” Biometrika.
James, and Stein. 1961. β€œEstimation with Quadratic Loss.” In Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability.
JankovΓ‘, and van de Geer. 2015. β€œHonest Confidence Regions and Optimality in High-Dimensional Precision Matrix Estimation.” arXiv:1507.02061 [Math, Stat].
Kauermann, and Carroll. 2001. β€œA Note on the Efficiency of Sandwich Covariance Matrix Estimation.” Journal of the American Statistical Association.
Khoromskij, Litvinenko, and Matthies. 2009. β€œApplication of Hierarchical Matrices for Computing the Karhunen–LoΓ¨ve Expansion.” Computing.
Khoshgnauz. 2012. β€œLearning Markov Network Structure Using Brownian Distance Covariance.” arXiv:1206.6361 [Cs, Stat].
Kuismin, and SillanpÀÀ. 2017. β€œEstimation of Covariance and Precision Matrix, Network Structure, and a View Toward Systems Biology.” WIREs Computational Statistics.
Lam, and Fan. 2009. β€œSparsistency and Rates of Convergence in Large Covariance Matrix Estimation.” Annals of Statistics.
Ledoit, and Wolf. 2004. β€œA Well-Conditioned Estimator for Large-Dimensional Covariance Matrices.” Journal of Multivariate Analysis.
Ling. 1974. β€œComparison of Several Algorithms for Computing Sample Means and Variances.” Journal of the American Statistical Association.
Loh. 1991. β€œEstimating Covariance Matrices II.” Journal of Multivariate Analysis.
Mardia, and Marshall. 1984. β€œMaximum Likelihood Estimation of Models for Residual Covariance in Spatial Regression.” Biometrika.
Meier. 2018. β€œA matrix Gamma process and applications to Bayesian analysis of multivariate time series.”
Meier, Kirch, and Meyer. 2020. β€œBayesian Nonparametric Analysis of Multivariate Time Series: A Matrix Gamma Process Approach.” Journal of Multivariate Analysis.
Meinshausen, and BΓΌhlmann. 2006. β€œHigh-Dimensional Graphs and Variable Selection with the Lasso.” The Annals of Statistics.
Mercer. 2000. β€œBounds for A–G, A–H, G–H, and a Family of Inequalities of Ky Fan’s Type, Using a General Method.” Journal of Mathematical Analysis and Applications.
Minasny, and McBratney. 2005. β€œThe MatΓ©rn Function as a General Model for Soil Variograms.” Geoderma, Pedometrics 2003,.
Mond, and PecΜ†ariΔ‡. 1996. β€œA Mixed Arithmetic-Mean-Harmonic-Mean Matrix Inequality.” Linear Algebra and Its Applications, Linear Algebra and Statistics: In Celebration of C. R. Rao’s 75th Birthday (September 10, 1995),.
PΓ©bay. 2008. β€œFormulas for Robust, One-Pass Parallel Computation of Covariances and Arbitrary-Order Statistical Moments.” Sandia Report SAND2008-6212, Sandia National Laboratories.
Pleiss, Gardner, Weinberger, et al. 2018. β€œConstant-Time Predictive Distributions for Gaussian Processes.” In.
Pourahmadi. 2011. β€œCovariance Estimation: The GLM and Regularization Perspectives.” Statistical Science.
Prause, and Steland. 2018. β€œEstimation of the Asymptotic Variance of Univariate and Multivariate Random Fields and Statistical Inference.” Electronic Journal of Statistics.
Ramdas, and Wehbe. 2014. β€œStein Shrinkage for Cross-Covariance Operators and Kernel Independence Testing.” arXiv:1406.1922 [Stat].
Ravikumar, Wainwright, Raskutti, et al. 2011. β€œHigh-Dimensional Covariance Estimation by Minimizing β„“1-Penalized Log-Determinant Divergence.” Electronic Journal of Statistics.
Rigollet, and HΓΌtter. 2019. High Dimensional Statistics.
Rosenblatt. 1984. β€œAsymptotic Normality, Strong Mixing and Spectral Density Estimates.” The Annals of Probability.
Rothman. 2010. β€œSparse Estimation of High-Dimensional Covariance Matrices.”
Sampson, and Guttorp. 1992. β€œNonparametric Estimation of Nonstationary Spatial Covariance Structure.” Journal of the American Statistical Association.
SchΓ€fer, and Strimmer. 2005. β€œA shrinkage approach to large-scale covariance matrix estimation and implications for functional genomics.” Statistical Applications in Genetics and Molecular Biology.
Schmidt, and O’Hagan. 2003. β€œBayesian Inference for Non-Stationary Spatial Covariance Structure via Spatial Deformations.” Journal of the Royal Statistical Society: Series B (Statistical Methodology).
Shao, and Wu. 2007. β€œAsymptotic Spectral Theory for Nonlinear Time Series.” The Annals of Statistics.
Sharma. 2008. β€œSome More Inequalities for Arithmetic Mean, Harmonic Mean and Variance.” Journal of Mathematical Inequalities.
Shimotsu, and Phillips. 2004. β€œLocal Whittle Estimation in Nonstationary and Unit Root Cases.” The Annals of Statistics.
Stein. 2005. β€œSpace-Time Covariance Functions.” Journal of the American Statistical Association.
Sun, and Stein. 2016. β€œStatistically and Computationally Efficient Estimating Equations for Large Spatial Datasets.” Journal of Computational and Graphical Statistics.
Takemura. 1984. β€œAn Orthogonally Invariant Minimax Estimator of the Covariance Matrix of a Multivariate Normal Population.” Tsukuba Journal of Mathematics.
Warton. 2008. β€œPenalized Normal Likelihood and Ridge Regularization of Correlation and Covariance Matrices.” Journal of the American Statistical Association.
Whittle, Peter. 1952. β€œSome Results in Time Series Analysis.” Scandinavian Actuarial Journal.
Whittle, P. 1952. β€œTests of Fit in Time Series.” Biometrika.
β€”β€”β€”. 1953a. β€œThe Analysis of Multiple Stationary Time Series.” Journal of the Royal Statistical Society: Series B (Methodological).
β€”β€”β€”. 1953b. β€œEstimation and Information in Stationary Time Series.” Arkiv FΓΆr Matematik.
Wolter. 2007. Introduction to Variance Estimation. Statistics for Social and Behavioral Sciences.
Wu, and Pourahmadi. 2003. β€œNonparametric Estimation of Large Covariance Matrices of Longitudinal Data.” Biometrika.
Yuan, and Lin. 2007. β€œModel Selection and Estimation in the Gaussian Graphical Model.” Biometrika.
Zeileis. 2004. β€œEconometric Computing with HC and HAC Covariance Matrix Estimators.” Journal of Statistical Software.
β€”β€”β€”. 2006a. β€œImplementing a Class of Structural Change Tests: An Econometric Computing Approach.” Computational Statistics & Data Analysis.
β€”β€”β€”. 2006b. β€œObject-Oriented Computation of Sandwich Estimators.” Journal of Statistical Software.
Zhang, and Zou. 2014. β€œSparse Precision Matrix Estimation via Lasso Penalized D-Trace Loss.” Biometrika.