The evaluation of variance component estimation software: generating benchmark problems by exact and approximate methods

The prediction of breeding values depends on the reliable estimation of variance components. This complex task leads to nonlinear minimization problems that have to be solved by numerical algorithms. In order to evaluate the reliability of these algorithms benchmark problems have to be constructed where the exact solution is a priori known. We develop techniques to construct such benchmark problems for mixed models including fixed and random effects, ANOVA, ML and REML predictors, balanced and unbalanced data for 1-way classification. Besides the construction of artificial data that produce the desired variance components we describe a projection method to construct benchmark data from simulated data. We discuss the cases where exact expressions for the projection can be given and where a numerical approximation procedure has to be used. Copyright Springer-Verlag Berlin Heidelberg 2013