APPROXIMATE RESTRICTED MAXIMUM-LIKELIHOOD AND APPROXIMATE PREDICTION ERROR VARIANCE OF THE MENDELIAN SAMPLING EFFECT

In an Expectation-Maximization type Restricted Maximum Likelihood (REML) procedure, the estimation of a genetic (co-)variance component involves the trace of the product of the inverse of the coefficient matrix by the inverse of the relationship matrix. Computation of this trace is usually the limiting factor of this procedure. In this paper, a method is presented to approximate this trace in the case of an animal model, by using an equivalent model based on the Mendelian sampling effect and by simplifying its coefficient matrix and its inversion. This approximation appeared very accurate for low heritabilities but was downwards biased when the heritability was high. Implemented in a REML procedure, this approximation reduced dramatically the amount of computation, but provided downwards biased estimates of genetic variances. Several examples axe presented to illustrate the method.