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.

+More

reduced rank estimation of covariance matrices expectation maximisation type algorithms full rank case restricted maximum likelihood fitting leading principal components quadratically convergent average information iterative solution scheme

Karin, Meyer,.Parameter expansion for estimation of reduced rank covariance matrices.. 40 (1),3-24.

Karin, Meyer,"Parameter expansion for estimation of reduced rank covariance matrices." 40

Karin, Meyer,"Parameter expansion for estimation of reduced rank covariance matrices."