Page 74 - Linear Models for the Prediction of Animal Breeding Values 3rd Edition
P. 74

−1
         Using the procedure outlined above, A  for the example data is:
               1       2      3      4      5      6       7     8      9     10

          1    1.83   0.50   0.00  −0.67    0.00  −1.00   0.00  0.00  −0.50  −0.17
          2    0.50   2.00   0.50   0.00   −1.00  −1.00   0.00  0.00  −0.50  −0.50
          3    0.00   0.50   2.00   0.00   −1.00   0.50   0.00  −1.00  −0.50  −0.50
          4   −0.67   0.00   0.00   1.83    0.50   0.00  −1.00  0.00   0.00  −0.67
          5    0.00  −1.00  −1.00   0.50    2.50   0.00  −1.00  0.00   0.00   0.00
          6   −1.00  −1.00   0.50   0.00    0.00   2.50   0.00  −1.00  0.00   0.00
          7    0.00   0.00   0.00  −1.00   −1.00   0.00   2.00  0.00   0.00   0.00
          8    0.00   0.00  −1.00   0.00    0.00  −1.00   0.00  2.00   0.00   0.00
          9   −0.50  −0.50  −0.50   0.00    0.00   0.00   0.00  0.00   0.75   0.75
         10   −0.17  −0.50  −0.50  −0.67    0.00   0.00   0.00  0.00   0.75   1.08

                                                           −1
         and A a is easily obtained by multiplying every element of A  by 2, the value of a. The
              −1
                −1
         matrix A a is added to equations for animal and group to obtain the MME, which are:
           ⎡  b 1 ⎤  ⎡ 3.000 0.000  0.000  0.000  0.000  1.000  0.000 0  0.000  1.000  1.000  0.000  0.000⎤  −1  ⎡ 13.0 ⎤
            ˆ
           ⎢  ˆ  ⎥ ⎥  ⎢                                                    ⎥  ⎢  ⎥
           ⎢  b 2 ⎥  ⎢ 0.000 2.000  0.000  0.000  0.000  0.000  1..000  1.000  0.000  0.000  0.000  0.000 ⎥ ⎥  ⎢  6.8  ⎥
           ⎢  a ˆ 1 ⎥  ⎢ 0.000 0.000  3.667  1.000  0.000 − 1.3333  0.000 − 2.000  0.000  0.000 − 1.000 − 0.333 ⎥  ⎢  0.00 ⎥
           ⎢  ⎥  ⎢                                                         ⎥  ⎢  ⎥
           ⎢  a ˆ 2⎥  ⎢ 0.000 0.000  1.000  4.000  1..000  0.000 − 2.000 − 2.000  0.000  0.000 − 1.000 − 1.000 ⎥  ⎢  0.00 ⎥
           ⎢  ⎥  ⎢ 0.000 0.000  0.000 0  1.000  4.000  0.000 − 2.000  1.000  0.000 − 2.000 − 1.000 − 1.000 ⎥  ⎢  0.00 ⎥
           ⎢  a ˆ 3 ⎥  ⎢                                                   ⎥  ⎢  ⎥
           ⎢  a ˆ 4 ⎥  ⎢ 1.000 0.0000 − 1.333  0.000  0.000  4.667  1.000  0.000 − 2.000  0.000  0.000 − 1.333⎥  ⎢  4.50 ⎥
           ⎢  ⎥ = ⎢                                                        ⎥  ⎢  ⎥
           ⎢  a ˆ 5 ⎥  ⎢ 0..000 1.000  0.000 − 2.000 − 2.000  1.000  6.000  0.000 − 2.000  0.000  0.000  0 0.000 ⎥  ⎢  2.90 ⎥
           ⎢  ⎥  ⎢ 0.000 1.000 − 2.000 − 2.000  1.000  0.000  0.000  6.000  0.000 − 2.000 0  0.000  0.000 ⎥  ⎢  3.90 ⎥
           ⎢  a ˆ 6 ⎥  ⎢                                                   ⎥  ⎢  ⎥
           ⎢  a ˆ 7 ⎥  ⎢ 1.000 0.000  0.000  0.000  0.000 − 2.000 − 2.000  0.000  5.000 0  0.000  0.000  0.000⎥  ⎢  3.50 ⎥
           ⎢  ⎥  ⎢                                                         ⎥  ⎢  ⎥
           ⎢  a ˆ 8 ⎥  ⎢ 1.000 0.000  0.000  0.000 − 2.000  0.000  0.000 − 2.000 0  0.000  5.000  0.000  0.000 ⎥  ⎢  5.00 ⎥
           ⎢  1 ˆ g ⎥  ⎢ 0.000 0.000 − 1.000 − 1.000 − 1.000  0.000  0.0000  0.000  0.000  0.000  1.500  1.500 ⎥  ⎢  0 0.00 ⎥
           ⎢  ⎥  ⎢                                                         ⎥  ⎢  ⎥
           ⎣ ⎢  g ˆ ⎥  ⎣ ⎢ 0.000 0.000 − 0.333 − 1.000 − 1.000 − 1..333  0.000  0.000  0.000  0.000  1.500  2.167 ⎦ ⎥  ⎣ ⎢  0.00 ⎦ ⎥
             2⎦
            There is dependency in the equations, that is, all effects cannot be estimated;
         therefore, the equation for the first group has been set to zero to obtain the following
         solutions:
                                 Sex effects
                                   Males          5.474
                                   Females        4.327
                                 Animals
                                   1             −0.780
                                   2             −0.936
                                   3             −0.977
                                   4             −1.287
                                   5             −1.113
                                   6             −0.741
                                   7             −1.354
                                   8             −0.782
                                 Groups
                                   9              0.000
                                  10             −1.795


          58                                                              Chapter 3
   69   70   71   72   73   74   75   76   77   78   79