The variance components were analyzed using an iterative version of the Laplace approximation for integration of the conditional likelihood function and SE for the estimates of variance and covariance components were obtained from the approximated observed information matrix implemented in the DMU software package, version 6, release 4.7 (Madsen and Jensen, 2006; Madsen et al., 2010).