Two kinds of variance/covariance estimates in linear mixed models

For longitudinal data, the within-subject covariance matrix plays an important role in statistical inference and it is of great interest to investigate this. In the paper, two kinds of estimators are investigated for the random effect covariance matrix D ₁ and the error variance σ ² in linear mixed models. One is to estimate D ₁ first and then to estimate σ ²; the other kind is to estimate σ ² first and then for D ₁. Both kinds of estimators are consistent. The covariance matrices of these covariance estimators and the variances of these two error variance estimators are calculated. In particular, the mean square errors of these estimators are also derived for one dimensional random effects. Besides, a simulation study is conducted to investigate the performances of these estimators.