Multiple Trait Evaluation

Mutiple trait genetic evaluation differs from a single trait evaluation in that several phenotypic traits are evaluated at the same time. Reasons for doing a multiple trait genetic evaluation include greater prediction accuracy and reduced prediction bias.

As with single trait genetic evaluation most traits use BLUP to obtain EPD by solving the mixed model equations. For single trait with breeding values as the only random effect the mixed model equations take the form

${\begin{pmatrix}\mathbf {X} '\mathbf {X} &\mathbf {X} '\mathbf {Z} \\\mathbf {Z} '\mathbf {X} &\mathbf {Z} '\mathbf {Z} +\mathbf {A} ^{-1}\lambda \end{pmatrix}}{\begin{pmatrix}{\hat {\boldsymbol {\beta }}}\\{\hat {\mathbf {u} }}\end{pmatrix}}={\begin{pmatrix}\mathbf {X} '\mathbf {y} \\\mathbf {Z} '\mathbf {y} \end{pmatrix}}$

where $\mathbf {X}$ and $\mathbf {Z}$ are the incidence matrices for the fixed and random effects, $\mathbf {A}$ is the numerator relationship matrix, ${\hat {\boldsymbol {\beta }}}$ are estimated fixed effects, ${\hat {\mathbf {u} }}$ the predicted breeding values, $\mathbf {y}$ the vector of observed phenotypes, and $\lambda =\sigma _{e}^{2}/\sigma _{g}^{2}$ is the ratio of the environmental variance $\sigma _{e}^{2}$ and the additive genetic variance $\sigma _{g}^{2}$.

In a multiple trait genetic evaluation we have covariances in addition to addition to variances associated with the random effects. In addition, different animals could have observed phenotypes for different subsets of traits. In the case where we have two traits and with breeding values as the only random effect he mixed model equations take the form

${\begin{pmatrix}\mathbf {X} '\mathbf {R} ^{-1}\mathbf {X} &\mathbf {X} '\mathbf {R} ^{-1}\mathbf {Z} \\\mathbf {Z} '\mathbf {R} ^{-1}\mathbf {X} &\mathbf {Z} '\mathbf {R} ^{-1}\mathbf {Z} +\mathbf {G} ^{-1}\end{pmatrix}}{\begin{pmatrix}{\hat {\boldsymbol {\beta }}}\\{\hat {\mathbf {u} }}\end{pmatrix}}={\begin{pmatrix}\mathbf {X} '\mathbf {R} ^{-1}\mathbf {y} \\\mathbf {Z} '\mathbf {R} ^{-1}\mathbf {y} \end{pmatrix}}$

where $G$ is the genetic covariance matrix and is a function of the numerator relationship matrix $\mathbf {A}$ and genetic variances of the two traits and the genetic covariance between the two traits, and $R$ is the environmental covariance matrix and is a function of the environmental variances of the two traits and the environmenatl covariance between the two traits.