Maximum likelihood estimate of variance components

Using the orthogonal complement likehood function, an iterative procedure for the maximum likelihood estimates of the variance and covariance components is derived. It is shown that these estimates are identical with the reproducing estimates of the locally best invariant quadratic unbiased estimation of variance and covariance components. Successive approximations of the maximum likelihood estimates are given in addition.