-
где i – интенсивность селекции (в нашем случае 1,5899-1,7095);
-
а – экономическое значение признака;
-
σ – фенотипическое стандартное отклонение (0,13-0,32 и 0,10-0,30);
-
h2 – коэффициент наследуемости.
-
Δ g я.п. = 1,5899√(0,85 * 0,18 * 0,13)2 + (0,92 * 0,23 * 0,32)2 = 0,112;
-
Δ g у.г. = 1,7095√(0,85 * 0,25 * 0,10)2 + (0,92 * 0,21 * 0,30)2 = 0,152.
Полученные результаты (11,2% – улучшения за поколение селекции по ярославской породе и 15,2% по улучшенным генотипам) дают возможность реализации метода в виде компьютерной программы.
Следует отметить, что предлагаемый вариант соответствует следующим требованиям:
-
1) корреляция между суммарным генотипом и селекционным признаком, а также ожидаемое генетическое улучшение с точки зрения экономики максимальны;
-
2) среднее квадратов отклонений между оптимальным линейным индексом и суммарным генотипом (I-H2 ) имеет минимальное значение;
-
3) оптимальный линейный индекс максимизирует вероятность того, что из двух особей с неодинаковыми данными по продуктивности будет отобрано то животное, которое имеет более высокую племенную ценность.
С использованием предложенного нами алгоритма в отделе селекции и генетики сельскохозяйственных животных разработана программа для ЭВМ «СТРЕСС» (рег. № 15070.7822000013.06.8.002.2/011 в Федеральном агентстве по науке и инновациям). Программа является базовой надстройкой для ИВС селекционного профиля, может использоваться с «АРМС-W», «СЕЛЕКС» и «ИНСЕЛ». В основу метода оценки генотипа по стрессоустойчивости заложен принцип гнездовых классификаций, поэтому на уровне стад в корректировке на межгодовые различия нет необходимости.
Индивидуальные оценки дочерей (первичное звено) на стрессоустойчивость выполняются по формулам:
ИСУ1 = 0,95(Р12 - Р11) + 0,90(Р13 - Р12)……..
ОИСУ = (ИСУ + а) / а
Здесь и далее,
ИСУ… – абсолютное значение индекса;
ОИСУ… – относительное значение индекса;
а – взвешивающий компонент, равный максимальному значению индекса;
Р11, Р12… – изменение признака;
Сводные оценки по быкам-производителям получены суммированием внутриклассовых отклонений с учетом наследуемости и количества потомков. Рабочие формулы имеют вид:
ИСУБ = В(∑ИСУД / ПD),
где ∑ИСУД – суммарная продуктивность дочерей;
ПD – количество дочерей;
В = W / W + Q,
где В – коэффициент регрессии генотипа родителей на генотип потомков;
Q = 4 - h2 / h2 .
W – количество эффективных дочерей;
W = (ПР * ПD) / (ПР+ ПD),
где ПР – количество оцененных коров по стаду.
ОИСУБ = (ИСУБ + к) / к;
к – взвешивающий компонент, учитывающий максимальную величину индекса и снижение точности оценок.
Оценка перспективных ветвей по стрессоустойчивости учитывает взвешенные оценки всех относящихся к ней производителей.
Расчеты проводится по формулам:
ИСУВ = ∑(WI * ИСУБI / ПDI)/ ∑ WI,
где ∑(WI * ИСУБI / ПDI)=W1 * ИСУБ1 / ПD1 + W2 * ИСУБ2 / ПD2 +…+
- WП * ИСУБП / ПDП; ∑ WI = WБ1 + WБ2 +…+ WБП;
где WБ1 = (ПВ * ПБ1) / (ПВ + ПБ1) …… WБП = (ПВ * ПБП) / (ПВ + ПБП);
где ПВ – количество животных в ветви.
ОИСУВ = (ИСУВ + к) / к.
Оценка линий выполняется через взвешенные оценки ветвей по формулам:
ИСУЛ = ∑(WI * ИСУБI / ПDI) / ∑ WI,
где ∑(WI * ИСУБI / ПDI) = W1 * ИСУБ1 / ПD1 + W2 * ИСУБ2 / ПD2 +…+
- WП * ИСУБП / ПDП;
∑ WI = WБ1 + WБ2 +…+ WБП;
WБ1 = (ПЛ * ПБ1) / (ПЛ + ПБ1) …… WБП = (ПЛ * ПБП) / (ПЛ + ПБП),
где ПЛ – количество животных в линии.
ОИСУЛ = (ИСУЛ + к) / к.
При необходимости можно получить характеристики популяции. В силу снижения точности оценок результаты можно использовать только для теоретических выкладок и определения средних значений.
Рабочие формулы:
ИСУБП = ∑(WI * ИСУБI / ПDI) / ∑ WI,
где ∑(WI * ИСУБI / ПDI) = W11 * ИСУБ11 + W21 * ИСУБ21 +…+ WП * ИСУБТ1; ∑ WI = WБ11 + WБ21 +…+ WБТ1,
где WБ11 = (ПП * ПБ11) / (ПП + ПБ11)…… WБПП = (ПП * ПБП) / (ПП + ПБП), где ПП – количество животных в популяции.
ОИСУБП = (ИСУБП + к) / к.
