Estimability analysis of variance and covariance components

Although variance and covariance components have been extensively investigated and a number of elegant formulae to compute them have been derived, nothing is known, without any ambiguity, about their estimability in the case of a fully unknown variance–covariance matrix. We prove that variance and covariance components in this case are not estimable, thus clarifying the ambiguity of the literature on the topic and correcting some erroneous statements in the literature. We also give a new theorem on the estimability of a linear function of variance and covariance components. Then we propose a new method to estimate the variance–covariance matrix with special structure, which can presumably be represented by, at most, r(r + 1)/2 independent parameters to guarantee its estimability in such a subspace, by directly implementing the positive definiteness of the matrix as constraint to the restricted maximum likelihood method, where r is the number of redundant measurements. Therefore, our estimates of the variance and covariance components always reconstruct a positive definite matrix and are always physically meaningful.