IngentaConnect Covariance matrix selection and estimation via penalised normal l...

In this: publication
By this: publisher
In this Subject: Biology ,&nbsp; Public Health
By this author: Huang, Jianhua Z. ;&nbsp; Liu, Naiping ;&nbsp; Pourahmadi, Mohsen ;&nbsp; Liu, Linxu

Authors: Huang, Jianhua Z.; Liu, Naiping; Pourahmadi, Mohsen; Liu, Linxu

Source: Biometrika, Volume 93, Number 1, March 2006 , pp. 85-98(14)

Purchase options

The full text electronic article is available for purchase. You will be able to download the full text electronic article after payment.

$33.35 plus tax Refund Policy

Abstract:

We propose a nonparametric method for identifying parsimony and for producing a statistically efficient estimator of a large covariance matrix. We reparameterise a covariance matrix through the modified Cholesky decomposition of its inverse or the one-step-ahead predictive representation of the vector of responses and reduce the nonintuitive task of modelling covariance matrices to the familiar task of model selection and estimation for a sequence of regression models. The Cholesky factor containing these regression coefficients is likely to have many off-diagonal elements that are zero or close to zero. Penalised normal likelihoods in this situation with L<inf>1</inf> and L<inf>2</inf> penalities are shown to be closely related to Tibshirani's (1996) LASSO approach and to ridge regression. Adding either penalty to the likelihood helps to produce more stable estimators by introducing shrinkage to the elements in the Cholesky factor, while, because of its singularity, the L<inf>1</inf> penalty will set some elements to zero and produce interpretable models. An algorithm is developed for computing the estimator and selecting the tuning parameter. The proposed maximum penalised likelihood estimator is illustrated using simulation and a real dataset involving estimation of a 102 × 102 covariance matrix.

Covariance matrix selection and estimation via penalised normal likelihood