Journal
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION
Volume 97, Issue 460, Pages 1141-1153Publisher
AMER STATISTICAL ASSOC
DOI: 10.1198/016214502388618942
Keywords
Bayesian model; Cholesky decomposition; high-dimensional covariance matrix; Markov chain Monte Carlo; model averaging
Categories
Ask authors/readers for more resources
This article proposes a data-driven method to identify parsimony in the covariance matrix of longitudinal data and to exploit any such parsimony to produce a statistically efficient estimator of the covariance matrix. The approach parameterizes the covariance matrix through the Cholesky decomposition of its inverse. For longitudinal data this is a one-step-ahead predictive representation, and-the Cholesky factor is likely to have off-diagonal elements that are zero or close to zero. A hierarchical Bayesian model is used to identify any such zeros in the Cholesky factor, similar to approaches that have been successful in Bayesian variable selection. The model is estimated using a Markov chain Monte Carlo sampling scheme that is computationally efficient and can be applied to covariance matrices of high dimension. It is demonstrated through simulations that the proposed method compares favorably in terms of statistical efficiency with a highly regarded competing approach. The estimator is applied to three real examples in which the dimension of the covariance matrix is large relative to the sample size. The first two examples are from biometry and electricity demand modeling and are longitudinal. The third example is from finance-and highlights the potential of our method for estimating cross-sectional covariance matrices.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available