Application of Householder`s transformations and the QL algorithm to REML estimation of variance components
Restricted maximum likelihood (REML) is widely regarded as the preferred procedure for estimating variance components in animal breeding problems. The size of the coefficient matrix, however, often leads to computational difficulties and many simplified algorithms, including diagonalization, have been proposed.…
