По результатам линейной оценки потомства селекционер может выбрать быка для улучшения основных признаков, представляющих для него интерес, а затем обратить пристальное внимание на описательные признаки, чтобы иметь уверенность, что данный бык будет удовлетворять специфическим нуждам конкретной высокопродуктивной коровы.
Цель оценки быков по потомству заключается в том, чтобы по имеющимся данным получить как можно более точный прогноз их генотипа (племенной ценности). Для этого необходимо устранить влияние на продуктивность дочерей многочисленных систематических факторов внешней среды (уровня кормления, условий содержания дочерей быков в разных хозяйствах, года и сезона отела, продолжительности лактации и т. д.). Опыт зарубежных стран свидетельствует, что получение более достоверных оценок племенной ценности быков становится возможным при использовании метода наилучшего линейного несмещенного прогноза (метод BLUP ).
Для апробации метода BLUP в Республике Беларусь использовали материалы зоотехнического и племенного учета по быкам-производителям и маточному поголовью Несвижского и Минского районов Минской области. Общее число первотелок составило 891 гол. Они являлись дочерьми 29 быков-производителей. Средняя молочная продуктивность первотелок, занесенных в базу данных, и ее изменчивость представлена в табл. 51. Корректурфакторы представлены в табл. 52
Таблица 51. Средняя продуктивность первотелок по хозяйствам
| Хозяйство | Количество | Удой, кг | Жир, % | Жир, кг | ||||
|---|---|---|---|---|---|---|---|---|
| быков | дочерей | X | Сv,% | X | Сv,% | X | Сv,% | |
| МРУП «Агро- комбинат «Ждановичи» | 3 | 116 | 5692 | 16,2 | 3,81 | 7,2 | 216,9 | 14,6 |
| СПК «Городея» | 2 | 62 | 5521 | 17,4 | 3,64 | 6,4 | 201,0 | 13,7 |
| ОАО «Гастел- ловское» | 2 | 59 | 3984 | 18,1 | 3,71 | 6,7 | 147,8 | 14,9 |
| СПК «Агроком- бинат «Снов» | 1 | 21 | 4960 | 16,5 | 3,48 | 7,1 | 172,6 | 15,1 |
| ОАО «Грицкевичи» | 3 | 98 | 4518 | 14,9 | 3,52 | 6,2 | 159,0 | 14,2 |
| ОАО «Новая жизнь» | 2 | 52 | 4180 | 19,6 | 3,49 | 5,8 | 145,9 | 13,6 |
| ЗАО «1 Мая» | 2 | 57 | 4618 | 20,4 | 3,61 | 6,1 | 166,7 | 15,7 |
| Хозяйство | Количество | Удой, кг | Жир, % | Жир, кг | ||||
|---|---|---|---|---|---|---|---|---|
| быков | дочерей | X | Сv,% | X | Сv,% | X | Сv,% | |
| ОАО «Юшевичи» | 1 | 23 | 3680 | 17,6 | 3,52 | 7,4 | 129,5 | 14,9 |
| ОАО «Сейловичи» | 3 | 104 | 3460 | 18,4 | 3,46 | 6,3 | 119,7 | 15,3 |
| ОАО «17 сентября» | 2 | 76 | 3890 | 17,6 | 3,61 | 5,9 | 140,4 | 14,6 |
| СПК «Карцевичи» | 1 | 18 | 3910 | 18,7 | 3,56 | 6,7 | 139,2 | 16,1 |
| КПСУП «Свекловичная» | 1 | 24 | 3770 | 19,8 | 3,44 | 7,2 | 129,7 | 16,3 |
| СПК «Несвижский» | 3 | 88 | 3190 | 18,7 | 3,63 | 6,9 | 115,8 | 16,7 |
| КПД «Беларусь» | 3 | 93 | 3740 | 16,8 | 3,61 | 6,5 | 135,0 | 15,7 |
| В среднем | 29 | 891 | 4220 | 18,1 | 3,58 | 6,8 | 176,6 | 14,9 |
Таблица 52. Аддитивные корректур-факторы для корректировки удоя и количества молочного жира на продолжительность лактации, кг
| Группа | Количе- ство дойных дней | Отелы 1996–1999 гг. Отелы 2000–2003 гг. | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Уровень продуктивности | Уровень продуктивности | ||||||||||||
| высокий | средний | низкий | высокий | средний | низкий | ||||||||
| удой | жир | удой | жир | удой | жир | удой | жир | удой | жир | удой | жир | ||
| I | 230–240 | 1109 | 43 | 835 | 33 | 793 | 27 | 902 | 35 | 724 | 26 | 650 | 21 |
| II | 241–250 | 990 | 38 | 740 | 28 | 682 | 23 | 658 | 27 | 687 | 25 | 600 | 20 |
| III | 251–260 | 739 | 30 | 603 | 22 | 521 | 19 | 899 | 31 | 610 | 24 | 525 | 19 |
| IV | 261–270 | 720 | 27 | 530 | 20 | 245 | 16 | 609 | 24 | 285 | 18 | 470 | 17 |
| V | 271–280 | 510 | 19 | 438 | 17 | 345 | 14 | 625 | 23 | 390 | 16 | 350 | 15 |
| VI | 281–290 | 292 | 14 | 365 | 15 | 251 | 9 | 390 | 18 | 310 | 13 | 213 | 9 |
| VII | 291–300 | 260 | 10 | 320 | 11 | 118 | 6 | 271 | 13 | 46 | 4 | 155 | 5 |
| VIII | 301–305 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Чтобы учесть влияние продолжительности лактации на продуктивность, была проведена предварительная корректировка данных на число дойных дней. Корректур-факторы рассчитывали методом наименьших
квадратов с учетом периода отела первотелок (1996–1999 и 2000–2003 гг.) и уровня продуктивности стада (высокий, средний, низкий). Посредством аддитивных корректур-факторов продуктивность каждой первотелки была приведена к продуктивности коров, имеющих продолжительность лактации 301…305 дн.:

y 1 = y 1 + k 1 ,
где y 1 – скорректированная продуктивность первотелки;
y 1– фактическая продуктивность первотелки;
k 1 – корректировочный коэффициент для 1-й группы продолжительности лактации.
Проведена группировка и кодировка данных по месяцу отела на три сезона:
I сезон – январь, февраль, март, апрель;
II сезон – май, июнь, июль, август;
III сезон – сентябрь, октябрь, ноябрь, декабрь.
Этот паратипический фактор соединен с годом отела и номером хозяйства, в котором лактировала первотелка. В результате такого объединения был сформирован один комплексный паратипический фактор: стадо–год–сезон ( HVS ), который учитывал влияние как каждого фактора в отдельности, так и всевозможные взаимодействия между ними.
В зависимости от кровности по голштинской породе все быки отнесены к пяти генетическим группам: чистопородные черно-пестрые, чистопородные голштинские, быки с кровностью по голштинской породе 50; 75; 87,5 %.
Система линейных уравнений смешанной модели (ММЕ) имеет вид:

X X ′ X Z ′ b ⌢ˆ X ′ Y
Z X ′ Z Z ′ + λ s ⌢ˆ = Z ′ Y ,
где λ – дисперсия по отцам (равна 0,25 σ 2 s ), λ = σ e 2 /σ s 2 = (4 − h 2 )/ h 2 ;
h 2 – коэффициент наследуемости признака;

/ – единичная матрица.
Решение системы ММЕ обеспечивает несмещенную оценку эффекта генетической группы, к которой относится оцениваемый бык, и наилучший несмещенный прогноз ½ аддитивной генетической ценности быка, не зависящий от влияния содержания и кормления дочерей на ферме, года и сезона отела, а также всех взаимодействий между этими факторами.
Племенная ценность ( BV ) была выражена как удвоенная оценка суммы эффекта быка-отца и эффекта генетической группы:
