Parameter expansion for estimation of reduced rank covariance matrices.
Parameter expanded and standard expectation maximisation algorithms are described for reduced rank estimation of covariance matrices by restricted maximum likelihood, fitting the leading principal components only. Convergence behaviour of these algorithms is examined for several examples and contrasted to that of the average information algorithm, and implications for practical analyses are discussed. It is shown that expectation maximisation type algorithms are readily adapted to reduced rank estimation and converge reliably. However, as is well known for the full rank case, the convergence is linear and thus slow. Hence, these algorithms are most useful in combination with the quadratically convergent average information algorithm, in particular in the initial stages of an iterative solution scheme.
Meyer
Karin
K
Animal Genetics and Breeding Unit, University of New England, Armidale NSW 2351, Australia. kmeyer@didgeridoo.une.edu.au
D000465:Algorithms; D000704:Analysis of Variance; D003627:Data Interpretation, Statistical; D016013:Likelihood Functions; D008957:Models, Genetic; D008962:Models, Theoretical; D025341:Principal Component Analysis