This LaTeX document is available as postscript or asAdobe PDF.

L. R. Schaeffer, April 1999

In most animal breeding applications, only additive genetic effects are considered in the evaluation of animals. In those cases, an infinitesimal animal model is assumed, where animals have been assumed to be randomly mating. There have been models developed where a finite number of loci have been assumed. In the last few years there has been an increased interest in models that consider non-additive genetic effects. However, applications have been limited because dominance genetic relationships would be required (which can be calculated for inbred populations with the genomic relationship matrix), but also the inverse of the dominance genetic relationship matrix is needed. There has not been an easy way to calculate this inverse, as was the case with the additive genetic relationship matrix. Hoeschele and VanRaden (1991) presented a method of inverting a dominance relationship matrix by including an extra sire by dam subclass effect into the model for non-inbred populations, which has been used in several research studies. Lately, a method has been found such that the inverses of the covariance matrices of the non-additive genetic effects are not necessary.

If non-additive genetic effects are to be included in an animal model,
then the assumption of random mating is still required, otherwise
non-zero covariances can arise between additive and dominance genetic
effects, which complicates the model enormously. Many studies also
assume that there is no inbreeding because of the difficulty in
accounting for inbreeding in the covariance matrices of the non-additive
genetic effects. However, this assumption is not absolutely necessary.
Thus, the models
to be explored in these notes are approximations
(as is any model). Consider a simple
animal model with additive, dominance, and additive by dominance
genetic effects, and repeated observations per animal, i.e.,

where is the overall mean,

**Simulation of Data**

The desired data structure is given in the following table for four animals.

Animal | Number of Records |

1 | 3 |

2 | 2 |

3 | 1 |

4 | 4 |

Assume that

The additive genetic relationship matrix for the four animals is

The dominance genetic relationship matrix is

The additive by dominance genetic relationship matrix is the Hadamard product of and , which is the element by element product of matrices.

The Cholesky decomposition of each of these is necessary to simulate the separate genetic effects. The simulated genetic effects for the four animals are (with , , and being vectors of random normal deviates)

In the additive genetic animal model we were able to simulate base population animals and for progeny to average the additive genetic values of the parents and add a random Mendelian sampling effect to obtain the additive genetic value of the progeny. With non-additive genetic effects, such a simple process does not exist. The appropriate genetic relationship matrices are necessary and these need to be decomposed. The alternative is to determine the number of loci affecting the trait, and to generate genotypes for each animal after defining the loci with dominance genetic effects and those that have additive by dominance interactions. This might be the preferred method depending on the objectives of the study.

Let the permanent environmental effects for the four animals be

The observations on the four animals, after adding a new residual effect for each record, and letting , are given in the table below.

Animal | 1 | 2 | 3 | 4 | ||||

1 | 12.91 | 15.09 | -12.22 | 8.16 | 36.21 | 45.69 | 49.41 | |

2 | 13.28 | 5.32 | -1.32 | -8.05 | 9.14 | -14.10 | ||

3 | -10.15 | -17.74 | -4.30 | -1.67 | -20.74 | |||

4 | 38.60 | 3.89 | 5.76 | 15.12 | 24.13 | 83.09 | 64.67 | 50.13 |

**HMME**

Using the simulated data, the MME that need to be constructed
are as follows.

where

and

The solutions are

and .

The total genetic merit of an animal can be estimated by adding
together the solutions for the additive, dominance, and
additive by dominance genetic values,

On the practical side, it is not clear how the solutions for the individual dominance and additive by dominance solutions should be used in breeding programs. Dominance effects are generated by particular sire-dam matings, and thus, dominance genetic values could be used to determine which matings were better. However, it is not clear how additive by dominance genetic solutions can be utilized. Perhaps the main point is that if non-additive genetic effects are significant, then they should be removed through the model to obtain more accurate estimates of the additive genetic effects, assuming that these have a much larger effect than non-additive genetic effects.

**Computing Simplification**

Take the MME as shown earlier, i.e.

Now subtract the equation for dominance genetic effects from the equation for additive genetic effects, and similarly for the additive by dominance and permanent environmental effects, giving

Re-arranging terms, then

The only inverse that is needed is for , and the equations to solve are only as large as the usual animal model MME. The steps in the procedure would be iterative.

- 1.
- Adjust the observation vector for solutions to
,
,
and
(initially these would be zero) as

- 2.
- Solve the following equations:

- 3.
- Obtain solutions for
,
,
and
using

- 4.
- Go to step 1 and begin again until convergence is reached.

**Estimation of Variances**

Given the new computing algorithm, and using Gibbs sampling as
a tool the variances can be estimated. Notice from the above
formulas that

Again, the inverses of and are not needed. The necessary quadratic forms are then

and . Generate 4 random Chi-Square variates,

The residual variance would be estimated from

where

This LaTeX document is available as postscript or asAdobe PDF.