In the case of variance component estimation, the original data set is replaced by a set of contrast calculated from the data, and the likelihood function is calculated from the probability distribution of these contrasts, according to the model for the complete data set. In particular, REML is used as a method for fitting linear mixed model. In contrast to the earlier likelihood function estimation, REML can produce unbiased estimates of variance and covariance parameters