An iterative procedure is presented for the estimation of variance component s in the mixed model with a numerator relationship matrix from non-balanced data .