Решение Коллегии Евразийской экономической комиссии от 24.11.2020 N 149
КОЛЛЕГИЯ ЕВРАЗИЙСКОЙ ЭКОНОМИЧЕСКОЙ КОМИССИИ
РЕШЕНИЕ
от 24 ноября 2020 г. N 149
ОБ УТВЕРЖДЕНИИ
МЕТОДИК ОЦЕНКИ ПЛЕМЕННОЙ ЦЕННОСТИ СЕЛЬСКОХОЗЯЙСТВЕННЫХ
ЖИВОТНЫХ В ГОСУДАРСТВАХ - ЧЛЕНАХ ЕВРАЗИЙСКОГО
ЭКОНОМИЧЕСКОГО СОЮЗА
В целях реализации подпункта 12 пункта 7 статьи 95 Договора о Евразийском экономическом союзе от 29 мая 2014 года и в соответствии с частью второй статьи 3 Соглашения о мерах, направленных на унификацию проведения селекционно-племенной работы с сельскохозяйственными животными в рамках Евразийского экономического союза, от 25 октября 2019 года (далее - Соглашение) Коллегия Евразийской экономической комиссии решила:
1. Утвердить прилагаемые:
Методику оценки племенной ценности крупного рогатого скота молочного направления продуктивности;
Методику оценки племенной ценности крупного рогатого скота мясного направления продуктивности;
Методику оценки племенной ценности свиней.
2. В целях ведения реестров учета племенных животных государств - членов Евразийского экономического союза (далее - реестры, государства-члены) и перерасчета значений племенной ценности крупного рогатого скота молочного и мясного направления продуктивности и свиней (далее - животные), включенных в реестры, в соответствии с методиками, утвержденными настоящим Решением (далее - методики), а также в целях внедрения методик в процессы производства племенной продукции установить переходный период сроком на 2 года с даты вступления настоящего Решения в силу.
3. Установить, что:
а) значения племенной ценности животных, включенных в реестры до даты вступления настоящего Решения в силу, пересчитываются уполномоченными органами (операторами) государств-членов, ответственными за ведение реестров, до окончания переходного периода в соответствии с методиками;
б) учет продуктивных показателей, оценка селекционируемых признаков, расчет значений племенной ценности животных, вновь включаемых в реестры, во время переходного периода допускается осуществлять в соответствии с законодательством государств-членов. При этом значения племенной ценности таких животных, рассчитанные в соответствии с законодательством государств-членов, должны быть пересчитаны уполномоченными органами (операторами) государств-членов, ответственными за ведение реестров, в соответствии с методиками до окончания переходного периода;
в) при отсутствии возможности осуществить перерасчет значений племенной ценности животных, указанных в подпунктах "а" и "б" настоящего пункта, допускается использование таких животных в процессе воспроизводства с оценкой племенной ценности, полученной ранее в соответствии с законодательством государств-членов. При этом значения оценки племенной ценности, не пересчитанные в соответствии с методиками, не учитываются в расчете значений племенной ценности животных, рожденных после завершения переходного периода;
г) племенная ценность животных, включаемых в реестры по истечении переходного периода, оценивается в соответствии с методиками.
4. Настоящее Решение вступает в силу по истечении 30 календарных дней с даты его официального опубликования, но не ранее даты вступления в силу Соглашения.
Председатель Коллегии
Евразийской экономической комиссии
М.МЯСНИКОВИЧ
Утверждена
Решением Коллегии
Евразийской экономической комиссии
от 24 ноября 2020 г. N 149
МЕТОДИКА
ОЦЕНКИ ПЛЕМЕННОЙ ЦЕННОСТИ КРУПНОГО РОГАТОГО СКОТА
МОЛОЧНОГО НАПРАВЛЕНИЯ ПРОДУКТИВНОСТИ
I. Общие положения
1. Настоящая Методика разработана в целях реализации подпункта 12 пункта 7 статьи 95 Договора о Евразийском экономическом союзе от 29 мая 2014 года, в соответствии со статьей 3 Соглашения о мерах, направленных на унификацию проведения селекционно-племенной работы с сельскохозяйственными животными в рамках Евразийского экономического союза, от 25 октября 2019 года и устанавливает порядок оценки, определения продуктивности и расчета племенной ценности крупного рогатого скота молочного направления продуктивности.
2. Настоящая Методика предназначена для применения на территориях государств - членов Евразийского экономического союза (далее - государства-члены) в племенных организациях, хозяйствах, осуществляющих выращивание и (или) реализацию племенного крупного рогатого скота молочного направления продуктивности, а также в сервисных организациях, информационно-аналитических, селекционных, селекционно-генетических центрах, союзах, ассоциациях (палатах), научных организациях, осуществляющих деятельность в области племенного молочного скотоводства.
3. Оценке племенной ценности подлежат коровы и быки молочного направления продуктивности, зарегистрированные в качестве племенных животных в соответствии с законодательством государств-членов.
4. Сведения о результатах оценки племенной ценности животных в соответствии с настоящей Методикой вносятся в реестр учета племенных животных (племенную книгу) государства-члена и в племенные свидетельства (паспорта, сертификаты).
5. Для целей настоящей Методики используются понятия, которые означают следующее:
"база данных" - структурированный набор данных о племенных животных, вовлеченных в селекционный процесс;
"биометрическая модель животного" (animal model, AM) - математическая форма описания взаимосвязи наблюдаемых фенотипических характеристик животного и влияния на них внешних факторов наряду с происхождением;
"индекс племенной ценности" - результат прогноза племенной ценности животного по комплексу селекционируемых признаков согласно их значимости для селекции;
"комплексный селекционный индекс" - индекс, включающий в себя частные селекционные индексы с весовыми коэффициентами согласно целям селекции;
"молочная продуктивность" - количество и качество молока, получаемого от животного;
"наилучший линейный несмещенный прогноз" (best linear unbiased prediction, BLUP) - статистический метод прогнозирования племенной ценности животного по селекционируемому признаку на основе биометрической модели животного линейного типа;
"племенная ценность" (estimated breeding value, EBV) - прогнозируемая племенная ценность животного по конкретному селекционируемому признаку, рассчитанная на основе метода BLUP AM;
"племенное животное" - сельскохозяйственное животное, используемое для разведения, зарегистрированное в реестре учета племенных животных в порядке, установленном законодательством государства-члена в области племенного животноводства, и имеющее в случае его реализации племенное свидетельство (паспорт, сертификат);
"племенное свидетельство (паспорт, сертификат)" - документ установленного образца, подтверждающий происхождение, племенную ценность и иные качества племенного животного (племенного стада);
"реестр учета племенных животных" - база данных, которая содержит сведения о племенных животных и племенных стадах и ведется в государстве-члене;
"селекционируемые признаки" - количественные и качественные показатели животных, по которым проводится целенаправленная селекция;
"сельскохозяйственные животные" - животные, разводимые в целях получения животноводческой продукции;
"частный селекционный индекс" - значение племенной ценности, выраженное в долях стандартного отклонения конкретного селекционируемого признака.
II. Оценка племенной ценности коров и быков молочного
направления продуктивности
6. Племенная ценность коров и быков молочного направления продуктивности определяется по удою, молочному жиру и молочному белку и рассчитывается:
а) у коров - по окончании лактации;
б) у быков - при наличии данных о продуктивности дочерей за 305 дней лактации.
7. Сведения о племенной ценности быка публикуются при появлении информации об окончании лактации у дочерей (с указанием количества его дочерей, стад, в которых они находятся, степени достоверности). Племенная ценность коров и быков молочного направления продуктивности пересчитывается не реже 1 раза в год.
8. Оценка коров по молочной продуктивности проводится согласно приложению N 1.
9. Расчет племенной ценности (EBV) коров и быков молочного направления продуктивности проводится на основе метода BLUP AM согласно приложению N 2.
10. Комплексные и частные селекционные индексы рассчитываются с учетом:
а) племенной ценности (EBV) по удою, молочному жиру и молочному белку;
б) весовых коэффициентов по удою, молочному жиру и молочному белку в соответствии с методиками, применяемыми в селекционно-племенной работе в государствах-членах;
в) информации о среднепопуляционных значениях и среднеквадратичных отклонениях селекционных показателей, размещенной на сайтах уполномоченных органов государств-членов в области племенного животноводства.
11. Результаты расчета племенной ценности (EBV) коров и быков молочного направления продуктивности вносятся в реестр учета племенных животных (племенную книгу) государства-члена, племенные свидетельства (паспорта, сертификаты).
Приложение N 1
к Методике оценки
племенной ценности крупного
рогатого скота молочного
направления продуктивности
ОЦЕНКА
КОРОВ ПО МОЛОЧНОЙ ПРОДУКТИВНОСТИ
1. Оценка уровня продуктивности коров и качества молока за лактацию или другой период производится путем обобщения результатов регулярно проводимых контрольных доек.
2. Контрольная дойка осуществляется сотрудниками, ответственными за проведение данного селекционного мероприятия.
3. Контрольная дойка проводится в течение суток одновременно у всех животных, содержащихся в одном помещении, за исключением сухостойных коров и новотельных коров до вечера 4-го дня после отела.
4. Качественный анализ контрольной пробы молока должен проводиться только в лаборатории по определению качества молока, аккредитованной в порядке, установленном законодательством государства - члена Евразийского экономического союза.
5. Днем начала лактации считается следующий день после отела. Окончанием лактации считается начало сухостойного периода. При отсутствии сухостойного периода у коровы днем окончания лактации считается день перед следующим отелом.
6. Для определения количества надоенного молока используются технические средства: весы с погрешностью взвешивания не более 0,1 кг, мерные емкости, молокомеры и электронные автоматические приборы.
7. Технические средства, используемые для определения количества надоенного молока, подвергаются проверке на точность показаний в порядке, установленном законодательством государства - члена Евразийского экономического союза.
8. Количество надоенного молока за контрольные сутки определяется путем сложения всех удоев, последовательно полученных в течение суток контрольного доения, с точностью до 0,1 кг. Удой за контрольный период рассчитывается с точностью до 1 кг.
9. Расчет количества молока, молочного жира и молочного белка за лактацию производится в соответствии с методиками, рекомендованными Международным комитетом по учету животных (ICAR).
10. Средний процент молочного жира и молочного белка за лактацию определяется путем деления количества однопроцентного молока на удой за соответствующую лактацию.
11. Массовая доля молочного жира и молочного белка за контрольные сутки и контрольный период определяется соответственно с точностью до 0,01 процента.
12. Количество молочного жира и молочного белка рассчитывается соответственно с точностью до 0,1 кг.
13. При расчете племенной ценности коров по молочной продуктивности используются данные за 305 дней лактации.
Приложение N 2
к Методике оценки
племенной ценности крупного
рогатого скота молочного
направления продуктивности
РАСЧЕТ
ПЛЕМЕННОЙ ЦЕННОСТИ КОРОВ И БЫКОВ МОЛОЧНОГО НАПРАВЛЕНИЯ
ПРОДУКТИВНОСТИ НА ОСНОВЕ МЕТОДА BLUP AM
1. Племенная ценность (EBV) коров и быков по молочной продуктивности рассчитывается на основе метода BLUP AM.
2. Расчет комплексных селекционных индексов племенной ценности коров и быков молочного направления продуктивности на основе метода BLUP AM состоит из следующих этапов:
а) разработка оптимальных статистических моделей, значимо описывающих развитие селекционируемых признаков в оцениваемой популяции;
б) расчет селекционно-генетических параметров оцениваемой популяции по оптимальным статистическим моделям (наследуемость, изменчивость (вариансы));
в) расчет прогнозных значений племенной ценности (EBV) на основе метода BLUP AM, надежности (точности) прогноза (REL, r2) и стандартизация прогнозных значений племенной ценности;
г) разработка комплексных селекционных индексов племенной ценности коров и быков молочного направления продуктивности на основе теории селекционного индекса и их расчет.
3. Для разработки статистических моделей развития селекционируемых признаков в популяции используются модели смешанного типа:
yij = hi + aij + eij,
где:
yij - показатель признака j-го животного в i-х условиях среды;
hi - эффекты условий среды (фиксированные);
aij - аддитивный генетический эффект j-го животного в i-х условиях среды (племенная ценность, EBV) (рандомизированный);
eij - эффект не учтенных в модели факторов (рандомизированный).
Для выбора оптимальной статистической модели используются информационный критерий Акаике (AIC) и Байесовский информационный критерий (BIC).
При использовании информационного критерия Акаике (AIC) выбирается модель, минимизирующая значение статистики:
,
где:
- остаточная сумма квадратов, деленная на количество наблюдений;
n - число наблюдений;
r - число оцененных параметров модели.
Байесовский информационный критерий (BIC) рассчитывается по формуле:
.
Лучшая статистическая модель соответствует минимальному значению критерия.
4. Для расчета прогнозных значений племенной ценности коров и быков молочного направления продуктивности по разработанным оптимальным статистическим моделям применяется метод BLUP AM.
Скалярная форма уравнения BLUP имеет вид:
y = Xb + Za + e,
где:
y = n x 1 - вектор наблюдений (оценок) (n - число записей);
b = p x 1 - вектор фиксированных эффектов (p - число уровней фиксированных эффектов);
a = q x 1 - вектор случайных эффектов пробандов (q - число уровней случайных эффектов);
e = n x 1 - вектор случайных эффектов;
X - матрица порядка n x p, которая связывает оценку животных с фиксированными эффектами;
Z - матрица порядка n x q, которая связывает оценку животных со случайными эффектами.
Матрицы X и Z называются матрицами случаев. Предполагается, что математическое ожидание (E) переменных:
E(y) = Xb;
E(a) = E(e) = 0.
Главная цель уравнения смешанной линейной модели - предсказать линейную функцию и b (EBV) от y.
Для вычисления и b необходимо решить уравнения смешанной линейной модели (MME) для вычисления значений b (фиксированных эффектов) и предсказать решения для значений (случайных эффектов). Формула для биометрической модели животного (AM) в матричном виде имеет вид:
.
Коэффициент рассчитывается по формуле:
,
отсюда искомые коэффициенты равны:
.
Таким образом, - лучшая линейная оценка фиксированных эффектов модели; - лучший линейный несмещенный прогноз (BLUP) племенной ценности (EBV) животного.
5. Матрица аддитивных генетических связей (A), соответствующая матрице числителей коэффициентов родства, рассчитывается по следующему рекурсивному алгоритму:
а) животные в родословной кодируются от 1 до n (n - число животных) и упорядочиваются таким образом, что родители предшествуют потомкам.
Если оба родителя (s и d) животного i известны, используются формулы:
aji = aij = 0,5 (ajs + ajd);
j = 1 до (i - 1);
aii = 1 + 0,5 (asd).
Если только один из родителей (s) известен и предполагается, что он не связан родством с другим, используются формулы:
aji = aij = 0,5 (ajs);
j = 1 до (i - 1);
aii = 1.
Если оба родителя неизвестны, используются формулы:
aji = aij = 0;
j = 1 до (i - 1);
aii = 1;
б) произведение матрицы A и аддитивной генетической вариансы дает описание вариационно-ковариационной структуры аддитивных генетических ценностей оцениваемых животных.
6. Для прогнозирования племенной ценности используются обратная матрица родства A-1, метод расчета A-1 без применения матрицы A.
Первоначально элементы матрицы родства A-1 задаются нулями, и применяются следующие правила.
Диагональные элементы задаются как 2, или 4/3, или 1 для животных с двумя известными, одним известным и с неизвестными родителями соответственно.
Если известны оба родителя i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i), (i, s), (d, i) и (i, d);
ai/4 - к элементам (s, s), (s, d), (d, s) и (d, d).
Если известен один из родителей i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i) и (i, s);
ai/4 - к элементу (s, s).
Если неизвестны оба родителя, добавляется ai к элементу (i, i).
При применении метода BLUP AM смешанной модели (MME) вида
матрица коэффициентов имеет вид
.
При этом обобщенная обратная матрица коэффициентов имеет вид
.
7. Вариансы ошибки прогноза (доля аддитивной генетической вариансы, не учитываемая прогнозом) (prediction error variance, PEV) рассчитываются по формуле:
,
где:
PEV - доля аддитивной генетической вариансы, не учитываемая прогнозом;
r2 - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности.
Для расчета PEV необходимы диагональные элементы матрицы коэффициентов уравнений животных.
8. Точность прогноза (r) - корреляция между истинными и прогнозируемыми оценками племенной ценности. Однако при оценке точность обычно выражается как надежность - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности (r2). Для расчета r или r2 требуются диагональные элементы инвертированной смешанной модели (MME).
Корень квадратный из PEV дает стандартную ошибку прогноза (standard error prediction, SEP):
.
9. Для снижения ошибки прогноза необходимо использовать такие доступные методы, которые максимизируют r при имеющемся количестве информации.
Надежность оценки (reliability, REL) рассчитывается по формуле:
.
Утверждена
Решением Коллегии
Евразийской экономической комиссии
от 24 ноября 2020 г. N 149
МЕТОДИКА
ОЦЕНКИ ПЛЕМЕННОЙ ЦЕННОСТИ КРУПНОГО РОГАТОГО СКОТА
МЯСНОГО НАПРАВЛЕНИЯ ПРОДУКТИВНОСТИ
I. Общие положения
1. Настоящая Методика разработана в целях реализации подпункта 12 пункта 7 статьи 95 Договора о Евразийском экономическом союзе от 29 мая 2014 года, в соответствии со статьей 3 Соглашения о мерах, направленных на унификацию проведения селекционно-племенной работы с сельскохозяйственными животными в рамках Евразийского экономического союза, от 25 октября 2019 года и устанавливает порядок оценки, определения продуктивности и расчета племенной ценности крупного рогатого скота мясного направления продуктивности.
2. Настоящая Методика предназначена для применения на территориях государств - членов Евразийского экономического союза (далее - государства-члены) в племенных организациях, хозяйствах, осуществляющих выращивание и (или) реализацию племенного крупного рогатого скота мясного направления продуктивности, а также в сервисных организациях, информационно-аналитических, селекционных, селекционно-генетических центрах, союзах, ассоциациях (палатах), научных организациях, осуществляющих деятельность в области племенного мясного скотоводства.
3. Оценке племенной ценности животных подлежат особи всех половозрастных групп крупного рогатого скота мясного направления продуктивности (в том числе телки и ремонтные бычки, коровы, быки-производители), зарегистрированные в качестве племенных животных в соответствии с законодательством государств-членов.
4. Сведения о результатах оценки племенной ценности животных в соответствии с настоящей Методикой вносятся в реестр учета племенных животных (племенную книгу) государства-члена и в племенные свидетельства (паспорта, сертификаты).
5. Для целей настоящей Методики используются понятия, которые означают следующее:
"база данных" - структурированный набор данных о племенных животных, вовлеченных в селекционный процесс;
"биометрическая модель животного" (animal model, AM) - математическая форма описания взаимосвязи наблюдаемых фенотипических характеристик животного и влияния на них внешних факторов наряду с происхождением;
"наилучший линейный несмещенный прогноз" (best linear unbiased prediction, BLUP) - статистический метод прогнозирования племенной ценности животного по селекционируемому признаку на основе биометрической модели животного линейного типа;
"племенная ценность" (estimated breeding value, EBV) - прогнозируемая племенная ценность животного по конкретному селекционируемому признаку, рассчитанная методом BLUP AM;
"племенное животное" - сельскохозяйственное животное, используемое для разведения, зарегистрированное в реестре учета племенных животных в порядке, установленном законодательством государства-члена в области племенного животноводства, и имеющее в случае его реализации племенное свидетельство (паспорт, сертификат);
"племенное свидетельство (паспорт, сертификат)" - документ установленного образца, подтверждающий происхождение, племенную ценность и иные качества племенного животного (племенного стада);
"продуктивность" - совокупность хозяйственно полезных признаков племенного животного, включая качество получаемой от него продукции;
"реестр учета племенных животных" - база данных, которая содержит сведения о племенных животных и племенных стадах и ведется в государстве-члене;
"селекционируемые признаки" - количественные и качественные показатели животных, по которым проводится целенаправленная селекция;
"селекционно-племенная работа" - комплекс мероприятий, направленных на совершенствование племенных и продуктивных качеств сельскохозяйственных животных;
"сельскохозяйственные животные" - животные, разводимые в целях получения животноводческой продукции;
"частный селекционный индекс" - значение племенной ценности, выраженное в долях стандартного отклонения конкретного селекционируемого признака.
II. Оценка племенной ценности телок, ремонтных бычков, коров
и быков-производителей мясного направления продуктивности
6. Племенная ценность телок, ремонтных бычков (далее - молодняк), коров и быков-производителей мясного направления продуктивности определяется:
а) у молодняка - по фактической живой массе при рождении и скорректированной на 205 дней (в соответствии с пунктом 8 настоящей Методики) и 365 дней (в соответствии с пунктом 9 настоящей Методики);
б) у коров - по скорректированной живой массе на 365 дней, легкости отела, молочности;
в) у быков-производителей - по легкости отела дочерей, молочности дочерей, фактической живой массе потомков при рождении, скорректированной на 205 дней и 365 дней.
7. Информация о племенной ценности рассчитывается и публикуется не реже 1 раза в год.
8. Скорректированная живая масса при отъеме животного в возрасте 205 дней рассчитывается по формуле:
,
где:
CMо - скорректированная живая масса при отъеме (кг);
Mо - фактическая живая масса при отъеме (кг);
Mр - фактическая живая масса при рождении (кг);
Bм - возраст животного на момент отъема (дней).
9. Скорректированная на 365 дней (в диапазоне 300 - 430 дней) живая масса рассчитывается по формуле:
,
где:
CMг - скорректированная живая масса в годовалом возрасте (кг);
Mг - фактическая живая масса в годовалом возрасте (кг);
Mо - фактическая живая масса при отъеме (кг);
Bг - возраст животного при взвешивании в годовалом возрасте (дней);
Bм - возраст животного при взвешивании на момент отъема (дней);
160 - числовой показатель разницы между 365 днями (годовалый возраст) и 205 днями (скорректированная живая масса при отъеме);
CMо - скорректированная живая масса при отъеме (кг).
10. Расчет племенной ценности молодняка по фактической живой массе при рождении, скорректированной на 205 дней и 365 дней, а также коров и быков-производителей проводится на основе метода BLUP AM согласно приложению N 1.
11. Оценка коров по легкости отела проводится по шкале согласно приложению N 2 по среднему значению (по всем отелам).
12. Оценка коров по молочности проводится по весу потомка при отъеме в пересчете на 205 дней по среднему значению (по всем отелам).
13. Частные селекционные индексы рассчитываются на основании племенной ценности (EBV) по фактической живой массе при рождении, скорректированной на 205 дней и 365 дней, по легкости отела, молочности.
Комплексный селекционный индекс рассчитывается на основании частных селекционных индексов с учетом их весовых коэффициентов в соответствии с методиками, применяемыми в государствах-членах при проведении селекционно-племенной работы.
14. Результаты расчета племенной ценности (EBV) крупного рогатого скота мясного направления продуктивности вносятся в реестр учета племенных животных (племенную книгу) государства-члена, племенные свидетельства (паспорта, сертификаты).
Информация о среднепопуляционных значениях и среднеквадратичных отклонениях селекционных показателей размещается на официальных сайтах уполномоченных органов государств-членов в области племенного животноводства в информационно-телекоммуникационной сети "Интернет".
Приложение N 1
к Методике оценки
племенной ценности крупного
рогатого скота мясного
направления продуктивности
РАСЧЕТ
ПЛЕМЕННОЙ ЦЕННОСТИ КРУПНОГО РОГАТОГО СКОТА МЯСНОГО
НАПРАВЛЕНИЯ ПРОДУКТИВНОСТИ НА ОСНОВЕ МЕТОДА BLUP AM
1. Племенная ценность (EBV) крупного рогатого скота мясного направления продуктивности определяется по живой массе при рождении, на 205-й день, на 365-й день, по легкости отела, молочности и рассчитывается на основе метода BLUP AM.
2. Расчет комплексных селекционных индексов племенной ценности крупного рогатого скота мясного направления продуктивности на основе метода BLUP AM состоит из следующих этапов:
а) разработка оптимальных статистических моделей, значимо описывающих развитие селекционируемых признаков в оцениваемой популяции;
б) расчет селекционно-генетических параметров оцениваемой популяции по оптимальным статистическим моделям (наследуемость, изменчивость (вариансы));
в) расчет прогнозных значений племенной ценности (EBV) на основе метода BLUP AM, надежности (точности) прогноза (REL, r2) и стандартизация прогнозных значений племенной ценности;
г) разработка комплексных селекционных индексов племенной ценности крупного рогатого скота мясного направления продуктивности на основе теории селекционного индекса и их расчет.
3. Для разработки статистических моделей развития селекционируемых признаков в популяции используются модели смешанного типа:
yij = hi + aij + eij,
где:
yij - показатель признака j-го животного в i-х условиях среды;
hi - эффекты условий среды (фиксированные);
aij - аддитивный генетический эффект j-го животного в i-х условиях среды (племенная ценность, EBV) (рандомизированный);
eij - эффект не учтенных в модели факторов (рандомизированный).
Для выбора оптимальной статистической модели используются информационный критерий Акаике (AIC) и Байесовский информационный критерий (BIC).
При использовании информационного критерия Акаике (AIC) выбирается модель, минимизирующая значение статистики:
,
где:
- остаточная сумма квадратов, деленная на количество наблюдений;
n - число наблюдений;
r - число оцененных параметров модели.
Байесовский информационный критерий (BIC) рассчитывается по формуле:
.
Лучшая статистическая модель соответствует минимальному значению критерия.
4. Для расчета прогнозных значений племенной ценности крупного рогатого скота мясного направления продуктивности по разработанным оптимальным статистическим моделям применяется метод BLUP AM.
Скалярная форма уравнения BLUP имеет вид:
y = Xb + Za + e,
где:
y = n x 1 - вектор наблюдений (оценок) (n - число записей);
b = p x 1 - вектор фиксированных эффектов (p - число уровней фиксированных эффектов);
a = q x 1 - вектор случайных эффектов пробанда (q - число уровней случайных эффектов);
e = n x 1 - вектор случайных эффектов;
X - матрица порядка n x p, которая связывает оценку животных с фиксированными эффектами;
Z - матрица порядка n x q, которая связывает оценку животных со случайными эффектами.
Матрицы X и Z называются матрицами случаев. Предполагается, что математическое ожидание (E) переменных:
E(y) = Xb;
E(a) = E(e) = 0.
Главная цель уравнения смешанной линейной модели - предсказать линейную функцию и b (EBV) от y.
Для вычисления и b необходимо решить уравнения смешанной линейной модели (MME) для вычисления значений b (фиксированных эффектов) и предсказать решения для значений (случайных эффектов). Формула для биометрической модели животного (AM) в матричном виде имеет вид:
.
Коэффициент рассчитывается по формуле:
,
отсюда искомые коэффициенты равны:
.
Таким образом, - лучшая линейная оценка фиксированных эффектов модели, - лучший линейный несмещенный прогноз (BLUP) племенной ценности (EBV) животного.
5. Матрица аддитивных генетических связей (A), соответствующая матрице числителей коэффициентов родства, рассчитывается по следующему рекурсивному алгоритму:
а) животные в родословной кодируются от 1 до n (n - число животных) и упорядочиваются таким образом, что родители предшествуют потомкам.
Если оба родителя (s и d) животного i известны:
aji = aij = 0,5 (ajs + ajd);
j = 1 до (i - 1);
aii = 1 + 0,5 (asd).
Если только один из родителей (s) известен и предполагается, что он не связан родством с другим:
aji = aij = 0,5 (ajs);
j = 1 до (i - 1);
aii = 1.
Если оба родителя неизвестны:
aji = aij = 0;
j = 1 до (i - 1);
aii = 1;
б) произведение матрицы A и аддитивной генетической вариансы дает описание вариационно-ковариационной структуры аддитивных генетических ценностей оцениваемых животных;
в) для прогнозирования племенной ценности используются обратная матрица родства A-1, метод расчета A-1 (без применения матрицы A).
Первоначально элементы матрицы родства A-1 задаются нулями, и применяются следующие правила.
Диагональные элементы задаются как 2, или 4/3, или 1 для животных с двумя известными, одним известным и с неизвестными родителями соответственно.
Если известны оба родителя i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i), (i, s), (d, i) и (i, d);
ai/4 - к элементам (s, s), (s, d), (d, s) и (d, d).
Если известен один из родителей i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i) и (i, s);
ai/4 - к элементу (s, s).
Если неизвестны оба родителя, добавляется ai к элементу (i, i).
При применении метода BLUP AM смешанной модели (MME) вида
матрица коэффициентов имеет вид:
.
При этом обобщенная обратная матрица коэффициентов имеет вид:
.
Вариансы ошибки прогноза (доля аддитивной генетической вариансы, не учитываемая прогнозом) (prediction error variance, PEV) рассчитываются по формуле:
,
где:
PEV - доля аддитивной генетической вариансы, не учитываемая прогнозом;
r2 - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности.
Для расчета PEV необходимы диагональные элементы матрицы коэффициентов уравнений животных.
Точность прогноза (r) - корреляция между истинными и прогнозируемыми оценками племенной ценности. Однако при оценке точность обычно выражается как надежность - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности (r2). Для расчета r или r2 требуются диагональные элементы инвертированной смешанной модели (MME).
Корень квадратный из PEV дает стандартную ошибку прогноза (standard error prediction, SEP):
.
Для снижения ошибки прогноза необходимо использовать такие доступные методы, которые максимизируют r при имеющемся количестве информации.
Надежность оценки (reliability, REL) рассчитывается по формуле:
.
Приложение N 2
к Методике оценки
племенной ценности крупного
рогатого скота мясного
направления продуктивности
ШКАЛА
ОЦЕНКИ ЛЕГКОСТИ ОТЕЛА КОРОВ
Балл (код)
Характеристика легкости отела
Описание
1
Самостоятельный отел
корова (первотелка) отелилась без посторонней помощи
2
Легкое родовспоможение
без применения специализированного инструмента
3
Тяжелый отел
с применением специализированного инструмента
4
Неправильное предлежание плода
требуется помощь при отеле
5
Хирургическое вмешательство
требуется хирургическое вмешательство
Утверждена
Решением Коллегии
Евразийской экономической комиссии
от 24 ноября 2020 г. N 149
МЕТОДИКА
ОЦЕНКИ ПЛЕМЕННОЙ ЦЕННОСТИ СВИНЕЙ
I. Общие положения
1. Настоящая Методика разработана в целях реализации подпункта 12 пункта 7 статьи 95 Договора о Евразийском экономическом союзе от 29 мая 2014 года, в соответствии со статьей 3 Соглашения о мерах, направленных на унификацию проведения селекционно-племенной работы с сельскохозяйственными животными в рамках Евразийского экономического союза, от 25 октября 2019 года и устанавливает порядок оценки, определения продуктивности и расчета племенной ценности свиней.
2. Настоящая Методика предназначена для применения на территориях государств - членов Евразийского экономического союза (далее - государства-члены) в племенных организациях, хозяйствах, осуществляющих выращивание и (или) реализацию племенных свиней, а также в сервисных организациях, информационно-аналитических, селекционных, селекционно-генетических центрах, союзах, ассоциациях (палатах), научных организациях, осуществляющих деятельность в области племенного свиноводства.
3. Оценке племенной ценности подлежат особи всех половозрастных групп свиней (в том числе ремонтные свинки и хрячки, свиноматки, хряки), зарегистрированные в качестве племенных животных в соответствии с законодательством государств-членов.
4. Результаты оценки племенной ценности свиней в соответствии с настоящей Методикой вносятся в реестр учета племенных животных (племенную книгу) государства-члена, племенные свидетельства (паспорта, сертификаты).
5. Для целей настоящей Методики используются понятия, которые означают следующее:
"база данных" - структурированный набор данных о племенных животных, вовлеченных в селекционный процесс;
"биометрическая модель животного" (animal model, AM) - математическая форма описания взаимосвязи наблюдаемых фенотипических характеристик животного и влияния на них внешних факторов наряду с происхождением;
"выращивание" - период содержания животных от перевода в группу ремонтного молодняка до проведения оценки собственной продуктивности и развития;
"доращивание" - период содержания поросят с даты отъема от подсосных свиноматок до перевода их в группу ремонтного молодняка или на откорм;
"индексная оценка" - метод определения племенной ценности животного по комплексу селекционируемых признаков согласно их значимости для селекции;
"комплексный селекционный индекс" - индекс, включающий в себя частные селекционные индексы с весовыми коэффициентами согласно целям селекции;
"линейная оценка экстерьера" - метод экспертной оценки статей экстерьера животных с помощью количественной шкалы;
"матрица" - математический объект, записываемый в виде прямоугольной таблицы элементов чисел, представляющей собой совокупность строк и столбцов, на пересечении которых находятся ее элементы. Количество строк и столбцов задает размер матрицы;
"наилучший линейный несмещенный прогноз" (best linear unbiased prediction, BLUP) - статистический метод прогнозирования племенной ценности животного по селекционируемому признаку на основе биометрической модели животного линейного типа;
"обратная матрица" - матрица A-1, при умножении на которую исходная матрица A дает в результате единичную матрицу (матрицу, все диагональные элементы которой равны 1);
"племенная ценность" (estimated breeding value, EBV) - прогнозируемая племенная ценность животного по конкретному селекционируемому признаку, рассчитанная на основе метода BLUP AM;
"племенное животное" - сельскохозяйственное животное, используемое для разведения, зарегистрированное в реестре учета племенных животных (в племенной книге) в порядке, установленном законодательством государства-члена в области племенного животноводства, и имеющее в случае его реализации племенное свидетельство (паспорт, сертификат);
"племенное свидетельство (паспорт, сертификат)" - документ установленного образца, подтверждающий происхождение, племенную ценность и иные качества племенного животного (племенного стада);
"племенное стадо" - группа племенных животных определенного вида и породы, используемых в селекционных целях;
"популяция" - совокупность особей животных определенного вида, в пределах которой происходит размножение;
"порода" - группа животных общего происхождения, созданная человеком, обладающая генетически обусловленными биологическими и морфологическими хозяйственно полезными свойствами, специфичными для данной группы животных, позволяющими отличить ее от других пород этого вида и устойчиво передающимися по наследству;
"продуктивность" - совокупность хозяйственно полезных признаков племенного животного, включая качество получаемой от него продукции;
"реестр учета племенных животных" - база данных, которая содержит сведения о племенных животных и племенных стадах и ведется в государстве-члене;
"селекционируемые признаки" - количественные и качественные показатели животных, по которым проводится целенаправленная селекция;
"селекционно-племенная работа" - комплекс мероприятий, направленных на совершенствование племенных и продуктивных качеств сельскохозяйственных животных;
"сельскохозяйственные животные" - животные, разводимые в целях получения животноводческой продукции;
"статистическая модель" - математическая форма описания взаимосвязи фенотипических характеристик животного и факторов паратипического и генетического влияния на них;
"частный селекционный индекс" - значение племенной ценности, выраженное в долях стандартного отклонения конкретного селекционируемого признака;
"экстерьер животного" - внешний вид животного, его наружные формы в целом, а также внешние особенности и развитость статей, характеризующие тип телосложения.
II. Основные требования к определению
племенной ценности свиней
6. Племенная ценность свиней оценивается по комплексу признаков на основе метода BLUP AM, расчет селекционных индексов проводится согласно теории селекционного индекса на основе экономических весов признаков.
7. Для применения метода BLUP обязательным является наличие достоверной базы данных.
8. При расчете индексов племенной ценности свиней на основе метода BLUP для каждой особи осуществляется расчет селекционного индекса с учетом следующих факторов:
а) информация обо всех внесенных в базу данных родственниках животного (родителях, прародителях, боковых родственниках, потомках) с учетом степени родства, что расширяет сведения о его генетической ценности;
б) отклонения в показателях продуктивности животного, которые корректируются в зависимости от влияния условий среды;
в) продуктивность, переданная потомству и скорректированная по уровню спаривания;
г) генетическая и фенотипическая корреляция между признаками (учет генетической конкуренции, уровня спаривания). Показатели племенной ценности корректируются по отношению друг к другу. В биометрической модели животного учитывается корреляция между ними (например, между скоростью роста и толщиной шпика). При этом племенная ценность по каждому критерию умножается на степень ее достоверности.
III. Оценка племенной ценности ремонтных
свинок и хрячков
9. Первый отбор поросят для целей селекции проводится по результатам визуальной оценки в день отъема от свиноматки (при переводе на доращивание) с учетом индексов племенной ценности матери и отца и их продуктивности. Индекс родословной рассчитывается по формуле:
Ир = (Ио + Им) x 0,5,
где:
Ир - индекс родословной;
Ио - селекционный индекс отца (рассчитанный согласно разделу VI настоящей Методики);
Им - селекционный индекс матери (рассчитанный согласно разделу VI настоящей Методики).
Поросята должны быть здоровы, без каких-либо пороков, хорошо развиты. К отбору не допускаются животные с количеством сосков менее 12 (6/6).
10. При переводе на выращивание (в группу ремонтного молодняка) ремонтный молодняк подлежит индивидуальному взвешиванию, сведения о результатах взвешиваний вносятся в базу данных.
11. Молодняк в конце периода выращивания периодически взвешивается и осматривается, при этом особое внимание уделяется наличию пороков экстерьера. При достижении веса 90 - 110 кг весь ремонтный молодняк оценивается по показателям собственной продуктивности, указанным в пункте 13 настоящей Методики. Сведения о результатах взвешиваний и измерений вносятся в базу данных.
12. По результатам оценки фенотипических показателей собственной продуктивности, индексной оценки и линейной оценки экстерьера для ремонта собственного стада отбираются свинки и хрячки, показавшие лучший результат. Остальное поголовье (за исключением животных больных, с пороками, отстающих в росте и развитии) реализуется.
13. К показателям собственной продуктивности молодняка относятся следующие селекционируемые признаки:
а) возраст достижения живой массы 100 кг (дней);
б) среднесуточный прирост живой массы от рождения до достижения живой массы 100 кг (г);
в) среднесуточный прирост живой массы на выращивании (г);
г) затраты корма на 1 кг прироста живой массы от 30 до 100 кг (кг) - для хозяйств, оснащенных станциями контрольного выращивания;
д) длина туловища (см);
е) толщина шпика (мм);
ж) высота (глубина) длиннейшей мышцы спины (мм).
14. Пересчет фактических показателей собственной продуктивности на 100 кг осуществляется при достижении живой массы 90 - 110 кг. При живой массе менее 90 и более 110 кг пересчет не производится и измерения не используются в расчете племенной ценности.
15. Для определения возраста достижения живой массы 100 кг осуществляется взвешивание свиней на весах с пределом взвешивания до 500 кг и погрешностью взвешивания не более 0,5 кг.
16. Возраст достижения живой массы 100 кг рассчитывается по формуле:
,
где:
X - возраст достижения живой массы 100 кг (дней);
В - фактический возраст в день последнего взвешивания (дней);
М - фактическая живая масса животного в день последнего взвешивания (кг);
П - среднесуточный прирост живой массы на выращивании (кг).
Полученный результат округляется до целого числа.
17. Среднесуточный прирост живой массы от рождения до достижения живой массы 100 кг определяется у ремонтного молодняка (свинок и хрячков) путем периодического взвешивания животных до достижения живой массы 90 - 110 кг и рассчитывается по формуле:
,
где:
C - среднесуточный прирост живой массы от рождения до достижения живой массы 100 кг (г);
m2 - живая масса животного при последнем взвешивании (кг) (от 90 до 110 кг);
n - возраст животного при последнем взвешивании (дней);
1000 - коэффициент пересчета в граммы.
18. Среднесуточный прирост живой массы на выращивании с точностью до грамма рассчитывается по формуле:
П = (m2 - m1) / (n2 - n1),
где:
П - среднесуточный прирост живой массы на выращивании (г);
m2 - живая масса животного в день проведения оценки (от 90 до 110 кг) (кг);
m1 - живая масса животного в начале постановки на выращивание (кг);
n2 - возраст животного в день проведения оценки (от 90 до 110 кг) (дней);
n1 - возраст животного в начале постановки на выращивание (дней).
19. Затраты корма на 1 кг прироста живой массы от 30 до 100 кг определяются от первого до предпоследнего дня контроля путем ежедневного взвешивания сухого корма, предназначенного для кормления животных в период контроля, на весах с пределом взвешивания до 10 кг и погрешностью взвешивания не более 0,05 кг. Кормление животных осуществляется не реже 2 раз в сутки, при этом не допускаются остатки и потери корма.
20. Затраты корма на 1 кг прироста живой массы за период контроля рассчитываются по формуле:
X2 = K / X1,
где:
X2 - затраты корма на 1 кг прироста живой массы от 30 до 100 кг (кг);
K - масса съеденного сухого корма за период контроля (кг);
X1 - прирост живой массы за период контроля (кг).
Полученный результат фиксируется с точностью до 0,01 кг.
21. Длина туловища измеряется по средней линии спины от затылочного гребня до корня хвоста с помощью стальной мерной ленты с ценой деления 1 см.
22. Длина туловища рассчитывается с учетом поправочного коэффициента на 1 кг живой массы, уменьшающего или увеличивающего фактическую длину туловища в зависимости от отклонения фактической живой массы от стандартной величины 100 кг, по формуле:
Д100 = ДФ + 0,35 x (100 - М),
где:
Д100 - длина туловища (см);
ДФ - фактическая длина туловища (см);
М - фактическая живая масса животного в день последнего взвешивания (кг);
0,35 - поправочный коэффициент.
23. Толщина шпика определяется на живых свиньях с помощью прибора ультразвукового исследования мясных качеств свиней с погрешностью не более 1 мм (Piglog-105 или аналог) в 2 точках (Р1, Р2) по схеме согласно приложению N 1.
24. Толщина шпика в точках Р1 и Р2 рассчитывается с учетом поправочного коэффициента на 1 кг живой массы, уменьшающего или увеличивающего фактическую толщину шпика в зависимости от отклонения фактической живой массы от стандартной величины 100 кг, по формуле:
ТР1, Р2 = ТФ Р1, Р2 + 0,15 x (100 - М),
где:
ТР1, Р2 - толщина шпика в точках Р1 и Р2 (мм);
ТФ Р1, Р2 - фактическая толщина шпика в точках Р1 и Р2 (мм);
М - фактическая живая масса животного в день последнего взвешивания (кг);
0,15 - поправочный коэффициент.
25. Высота (глубина) длиннейшей мышцы спины определяется на живых свиньях с помощью прибора ультразвукового исследования мясных качеств свиней с погрешностью не более 1 мм (Piglog-105 или аналог) в точке Р2, предусмотренной приложением N 1 к настоящей Методике.
26. Высота (глубина) длиннейшей мышцы спины рассчитывается с учетом поправочного коэффициента на 1 кг живой массы, уменьшающего или увеличивающего фактическую высоту (глубину) длиннейшей мышцы спины в зависимости от отклонения фактической живой массы от стандартной величины 100 кг, по формуле:
В = ВФ + 0,25 x (100 - М),
где:
В - высота (глубина) длиннейшей мышцы спины (мм);
ВФ - фактическая высота (глубина) длиннейшей мышцы спины (мм);
М - фактическая живая масса животного в день последнего взвешивания (кг);
0,25 - поправочный коэффициент.
27. Экстерьер ремонтных свинок и хрячков оценивается визуально. Животные, имеющие кратерные соски, менее 12 (6/6) сосков, пороки (сильную иксобразность передних ног, резкий перехват за лопатками или в пояснице, провислую спину, мопсовидность, криворылость, неправильный прикус), оценке не подлежат и выбраковываются из стада.
28. С учетом сведений, содержащихся в базе данных, рассчитываются прогноз племенной ценности (EBV) ремонтных свинок и хрячков на основе метода BLUP AM и комплексные селекционные индексы согласно разделу VI настоящей Методики.
IV. Оценка племенной ценности свиноматок
29. Свиноматки оцениваются по:
а) собственной продуктивности (оценка, полученная на стадии ремонтной свинки согласно разделу III настоящей Методики);
б) воспроизводительным качествам.
30. Воспроизводительные качества проверяемых свиноматок оцениваются по первому опоросу, а основных - в среднем по всем опоросам (включая первый) по следующим показателям:
а) многоплодие (голов);
б) количество поросят при отъеме (голов);
в) масса гнезда при отъеме в 30 дней (кг).
31. Многоплодие определяется количеством родившихся живых поросят (в том числе слабых).
32. Масса гнезда при отъеме в 30 дней определяется путем корректировки фактической массы гнезда на 21-й - 45-й день после опороса (в зависимости от используемой технологии) с применением поправочных коэффициентов согласно приложению N 2 и фиксируется с точностью до 1 кг.
33. С учетом сведений, содержащихся в базе данных, рассчитываются прогноз племенной ценности (EBV) свиноматок на основе метода BLUP AM и комплексные селекционные индексы согласно разделу VI настоящей Методики.
V. Оценка племенной ценности хряков
34. Хряки оцениваются по:
а) собственной продуктивности (оценка, полученная на стадии ремонтного хрячка согласно разделу III настоящей Методики);
б) воспроизводительным качествам (оценка, полученная на основе данных о всех женских родственниках (матери, сестрах, дочерях), определенных согласно разделу IV настоящей Методики);
в) оплодотворяющей способности.
35. Оплодотворяющая способность хряка рассчитывается по формуле:
,
где:
ОСх - оплодотворяющая способность хряка;
Оп - количество опоросов свиноматок, осемененных хряком (спермопродукцией хряка);
А - количество абортов свиноматок, осемененных хряком (спермопродукцией хряка);
В - количество свиноматок, осемененных хряком (спермопродукцией хряка), выбракованных во второй период супоросности;
Ос - количество свиноматок, осемененных хряком (спермопродукцией хряка).
36. С учетом сведений, содержащихся в базе данных, рассчитываются прогноз племенной ценности (EBV) хряков на основе метода BLUP AM и комплексные селекционные индексы согласно разделу VI настоящей Методики.
VI. Расчет прогноза племенной ценности свиней на основе
метода BLUP AM и комплексных селекционных индексов
37. Расчет комплексных селекционных индексов на основе метода BLUP AM состоит из следующих этапов:
а) разработка оптимальных статистических моделей, значимо описывающих развитие селекционируемых признаков в оцениваемой популяции;
б) расчет селекционно-генетических параметров оцениваемой популяции по оптимальным статистическим моделям (наследуемость, изменчивость (вариансы));
в) расчет прогнозных значений племенной ценности (EBV) свиней на основе метода BLUP AM, определение надежности (точности) прогноза (REL, r2) и стандартизация прогнозных значений племенной ценности;
г) разработка комплексных селекционных индексов на основе теории селекционного индекса, их расчет и стандартизация.
38. Для разработки статистических моделей развития селекционируемых признаков в популяции используются модели смешанного типа:
yij = hi + aij + eij,
где:
yij - показатель признака j-го животного в i-х условиях среды;
hi - эффекты условий среды (фиксированные);
aij - аддитивный генетический эффект j-го животного в i-х условиях среды (племенная ценность, EBV) (рандомизированный);
eij - эффект не учтенных в модели факторов (рандомизированный).
39. Для выбора оптимальной статистической модели используются информационный критерий Акаике (AIC) и Байесовский информационный критерий (BIC).
При использовании информационного критерия Акаике (AIC) выбирается модель, минимизирующая значение статистики:
,
где:
- остаточная сумма квадратов, деленная на количество наблюдений;
n - число наблюдений;
r - число оцененных параметров модели.
Байесовский информационный критерий (BIC) рассчитывается по формуле:
.
Лучшая статистическая модель соответствует минимальному значению критерия.
40. Коэффициенты наследуемости селекционируемых признаков в оцениваемой популяции рассчитываются с помощью дисперсионного анализа по формуле:
,
где:
h2 - коэффициент наследуемости селекционируемого признака;
- дисперсия (варианса), обусловленная генотипом животного;
- дисперсия (варианса), обусловленная неучтенными (случайными) эффектами.
Дисперсия (вариансы) рассчитывается методом ограниченного максимального правдоподобия (restricted maximum likelihood, REML).
41. Для расчета прогнозных значений племенной ценности свиней по разработанным оптимальным статистическим моделям применяется метод BLUP AM. Скалярная форма уравнения BLUP имеет вид:
y = Xb + Za + e,
где:
y - вектор наблюдений (оценок), y = n x 1 (n - число записей);
b - вектор фиксированных эффектов, b = p x 1 (p - число уровней фиксированных эффектов);
a - вектор случайных эффектов пробанда, a = q x 1 (q - число уровней случайных эффектов);
e - вектор случайных эффектов, e = n x 1 (n - число записей);
X - матрица порядка n x p, которая связывает оценку животных с фиксированными эффектами;
Z - матрица порядка n x q, которая связывает оценку животных со случайными эффектами.
42. Матрицы X и Z называются матрицами случаев, предполагается, что математическое ожидание (E) переменных имеет вид:
E(y) = Xb; E(a) = E(e) = 0.
Главная цель уравнения смешанной линейной модели - предсказать линейную функцию b и a (EBV) от y.
43. Для вычисления a и b необходимо решить уравнения смешанной линейной модели (MME) для вычисления значений b (фиксированных эффектов) и предсказать решения для значений a (случайных эффектов). Формула для биометрической модели животного (AM) в матричном виде имеет вид:
.
44. В свиноводстве прогноз племенной ценности производится по признакам собственной продуктивности, поэтому коэффициент рассчитывается по формуле:
.
Отсюда искомые коэффициенты равны:
.
Таким образом, - лучшая линейная оценка фиксированных факторов модели, - наилучший несмещенный прогноз племенной ценности (EBV) животного.
45. Матрица аддитивных генетических связей (A), соответствующая матрице числителей коэффициентов родства, рассчитывается по следующему рекурсивному алгоритму:
а) животные в родословной кодируются от 1 до n (n - число животных) и упорядочиваются таким образом, что родители предшествуют потомкам.
Если оба родителя (s и d) i-го животного известны, используются формулы:
aji = aij = 0,5 (ajs + ajd),
j = 1 до (i - 1),
aii = 1 + 0,5 (asd).
Если только один из родителей (s) известен и не связан родством с другим, используются формулы:
aji = aij = 0,5 (ajs),
j = 1 до (i - 1),
aii = 1.
Если оба родителя неизвестны, используются формулы:
aji = aij = 0,
j = 1 до (i - 1),
aii = 1;
б) произведение матрицы аддитивных генетических связей (A) и аддитивной генетической вариансы дает описание вариационно-ковариационной структуры аддитивных генетических ценностей оцениваемых животных;
в) для прогнозирования племенной ценности используются обратная матрица родства A-1, метод расчета A-1 (без применения матрицы аддитивных генетических связей A).
Первоначально элементы матрицы родства A-1 задаются нулями, и применяются следующие правила.
Диагональные элементы задаются как 2, или 4/3, или 1 для животных с 2 известными, 1 известным и с неизвестными родителями соответственно.
Если известны оба родителя i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i), (i, s), (d, i) и (i, d);
ai/4 - к элементам (s, s), (s, d), (d, s) и (d, d).
Если известен один из родителей i-го животного, добавляются:
ai - к элементу (i, i);
- ai/2 - к элементам (s, i) и (i, s);
ai/4 - к элементу (s, s).
Если неизвестны оба родителя, добавляется ai к элементу (i, i).
При применении скалярной формы уравнения BLUP AM смешанной модели (MME) вида
матрица коэффициентов имеет вид:
.
При этом обобщенная обратная матрица коэффициентов имеет вид:
;
г) вариансы ошибки прогноза (доля аддитивной генетической вариансы, не учитываемая прогнозом) (prediction error variance, PEV) рассчитываются по формуле:
,
где:
PEV - доля аддитивной генетической вариансы, не учитываемая прогнозом;
r2 - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности.
Точность прогноза (r) - корреляция между истинными и прогнозируемыми оценками племенной ценности. Однако при оценке точность обычно выражается как надежность - квадрат коэффициента корреляции между истинными и прогнозируемыми оценками племенной ценности (r2). Для расчета r или r2 требуются диагональные элементы инвертированной смешанной модели (MME).
Корень квадратный из PEV дает стандартную ошибку прогноза (standard error prediction, SEP):
.
Для снижения ошибки прогноза необходимо использовать доступные методы, которые максимизируют r при имеющемся количестве информации.
Надежность оценки для животного (reliability, REL) рассчитывается по формуле:
.
46. Для того чтобы значения племенной ценности по селекционируемым признакам можно было легче интерпретировать, проводится их стандартизация.
47. Среднее значение племенной ценности всех животных в оцениваемой популяции принимается за 100, а значение стандартного отклонения приравнивается к 12. Полученное значение является стандартизированной племенной ценностью животного (ПЦ, частный селекционный индекс), рассчитываемой по формуле:
,
где:
EBV - прогнозируемое значение племенной ценности селекционируемого признака, полученное на основе метода BLUP AM;
- стандартное отклонение значения прогнозируемой племенной ценности селекционируемого признака, полученного на основе метода BLUP AM;
12 - коэффициент, определяющий одну двенадцатую часть стандартного отклонения как 1 балл;
100 - коэффициент, определяющий уровень среднего значения.
Альтернативным методом стандартизации значений племенной ценности свиней по количественным селекционируемым признакам является спектральный метод. Стандартизация таких значений с использованием спектрального метода осуществляется в порядке согласно приложению N 3. Данный метод позволяет с помощью компактной записи (спектральной оценки) в деталях продемонстрировать прогнозируемую племенную ценность животного по каждому оцениваемому признаку.
48. Из стандартизированных значений племенной ценности для каждого селекционируемого признака с учетом весовых коэффициентов формируется комплексный селекционный индекс животного на основе метода BLUP AM по формуле:
И = К1 x ПЦ1 + К2 x ПЦ2 +... Кn x ПЦn,
где:
И - комплексный селекционный индекс;
К1 - весовой коэффициент первого селекционируемого признака;
ПЦ1 - частный селекционный индекс первого селекционируемого признака;
Кn - весовой коэффициент последнего селекционируемого признака;
ПЦn - частный селекционный индекс последнего селекционируемого признака;
n - количество селекционируемых признаков.
49. Расчет весовых коэффициентов селекционируемых признаков осуществляется согласно приложению N 4.
50. Стандартизация значений комплексных селекционных индексов (ИС) осуществляется по формуле:
,
где:
ИС - стандартизованный комплексный селекционный индекс;
Ри - значение комплексного селекционного индекса оцениваемого животного;
- среднее значение комплексного селекционного индекса в оцениваемой популяции;
- стандартное отклонение комплексного селекционного индекса в оцениваемой популяции;
12 - коэффициент, определяющий одну двенадцатую часть стандартного отклонения как 1 балл;
100 - коэффициент, определяющий уровень среднего значения.
51. В случае использования животного из другой популяции значение его племенной ценности приравнивается к среднему значению популяции, равному 100, а значения специфической племенной ценности будут равны нулю до тех пор, пока в данной популяции от него не будет получено потомство и определена продуктивность этого потомства.
Приложение N 1
к Методике оценки
племенной ценности свиней
АНАТОМИЧЕСКАЯ СХЕМА
РАСПОЛОЖЕНИЯ ТОЧЕК ДЛЯ УЛЬТРАЗВУКОВОГО ИЗМЕРЕНИЯ ТОЛЩИНЫ
ШПИКА И ВЫСОТЫ (ГЛУБИНЫ) ДЛИННЕЙШЕЙ МЫШЦЫ СПИНЫ
Р1 - точка находится между третьим и четвертым позвонками (с дорсальной стороны туловища) поясничного отдела позвоночника в 7 см от средней линии спины. В данной точке измеряют только толщину шпика.
Р2 - точка находится на уровне третьего-четвертого ребра (с дорсальной стороны туловища) в 7 см от средней линии спины. В данной точке последовательно измеряют толщину шпика и высоту (глубину) длиннейшей мышцы спины.
Приложение N 2
к Методике оценки
племенной ценности свиней
ПОПРАВОЧНЫЕ КОЭФФИЦИЕНТЫ
ДЛЯ ПЕРЕРАСЧЕТА МАССЫ ГНЕЗДА ПРИ ОТЪЕМЕ
НА МАССУ ГНЕЗДА В 30 ДНЕЙ
Возраст при взвешивании, дней
Коэффициент
Возраст при взвешивании, дней
Коэффициент
Возраст при взвешивании, дней
Коэффициент
21
1,47
35
0,86
49
0,54
22
1,40
36
0,82
50
0,52
23
1,32
37
0,79
51
0,51
24
1,26
38
0,76
52
0,50
25
1,20
39
0,73
53
0,48
26
1,15
40
0,70
54
0,47
27
1,14
41
0,68
55
0,46
28
1,07
42
0,66
56
0,45
29
1,04
43
0,64
57
0,44
30
1,00
44
0,62
58
0,42
31
0,97
45
0,60
59
0,41
32
0,94
46
0,58
60
0,40
33
0,91
47
0,57
61
0,39
34
0,88
48
0,55
62
0,38
Приложение N 3
к Методике оценки
племенной ценности свиней
ПОРЯДОК
СТАНДАРТИЗАЦИИ ЗНАЧЕНИЙ ПЛЕМЕННОЙ ЦЕННОСТИ СВИНЕЙ
ПО КОЛИЧЕСТВЕННЫМ СЕЛЕКЦИОНИРУЕМЫМ ПРИЗНАКАМ
С ИСПОЛЬЗОВАНИЕМ СПЕКТРАЛЬНОГО МЕТОДА
При стандартизации значений племенной ценности (EBV) свиней по количественным селекционируемым признакам с использованием спектрального метода племенная ценность (EBV), выражаемая в фактических единицах измерения селекционируемого признака, переводится в безразмерные единицы и выражается в баллах от 0 до 9, что позволяет сравнивать племенную ценность свиней по различным селекционируемым признакам в единой системе отсчета.
Данный метод оценки предполагает линейную зависимость между племенной ценностью свиней и оцениваемым количественным селекционируемым признаком. Баллы рассчитываются с точностью, достаточной для их обратного пересчета в исходные значения племенной ценности (EBV) без потери точности, по формуле:
для прямого показателя селекционируемого признака свиней (при увеличении фактического значения которого племенная ценность животного увеличивается):
;
для обратного показателя селекционируемого признака свиней (при увеличении фактического значения которого племенная ценность животного уменьшается):
,
где:
Б - балл;
Пn - фактическое значение племенной ценности по прямому показателю селекционируемого признака свиней;
По - фактическое значение племенной ценности по обратному показателю селекционируемого признака свиней;
- среднее арифметическое значение;
- стандартное отклонение.
Максимальное значение племенной ценности по анализируемому селекционируемому признаку в данной системе нормирования соответствует 9 баллам, минимальное - 0 баллов, а среднеарифметическое значение и стандартное отклонение племенной ценности по селекционируемому признаку - 3 баллам. С учетом этого 9 баллам будет соответствовать значение племенной ценности, равное , а 0 баллов - значение, равное . Таким образом, за пределами располагается 2,1% всех измерений анализируемого селекционируемого признака, относящихся к лучшим животным, оцененным по данному признаку, в то время как за пределами располагается 15,8% измерений, относящихся к худшим животным по данному признаку, как указано на рисунке.
Короткая оценочная запись, содержащая присвоенные баллы по каждому количественному селекционируемому признаку, округленные до целой части, в которой баллы за воспроизводительные показатели отделяются от баллов за мясооткормочные показатели двоеточием, называется спектральной оценкой животного. При этом рассчитанным значениям баллов меньше 0 соответствует балл, равный 0, а значениям больше 9 - балл, равный 9.
На основании точных (неокругленных) значений баллов по каждому количественному селекционируемому признаку рассчитываются комплексный балл за воспроизводительные показатели и комплексный балл за мясооткормочные показатели по формуле:
И = (К1 x Б1 + К2 x Б2 +... Кn x Бn) / (К1 + К2 +... Кn),
где:
И - комплексный балл;
Б1 - балл за первый анализируемый селекционируемый признак;
Бn - балл за n-й анализируемый селекционируемый признак;
К1 - весовой коэффициент первого анализируемого селекционируемого признака;
Кn - весовой коэффициент n-го анализируемого селекционируемого признака.
Приложение N 4
к Методике оценки
племенной ценности свиней
РАСЧЕТ
ВЕСОВЫХ КОЭФФИЦИЕНТОВ СЕЛЕКЦИОНИРУЕМЫХ ПРИЗНАКОВ
1. Экономическая значимость (Wi) селекционируемых признаков, входящих в комплексный селекционный индекс, определяется с учетом суммы денежного выражения всех селекционируемых признаков при денежном выражении каждого из них (V1, V2, V3... Vn - денежное выражение первого, второго, третьего... n-го признаков) и с учетом стандартного отклонения по селекционируемым признакам.
2. Экономическая значимость (Wi) по всем селекционируемым признакам рассчитывается по формуле:
,
где:
Wi - экономическая значимость i-го признака;
Vi - денежное выражение i-го селекционируемого признака;
Si - стандартное отклонение i-го селекционируемого признака;
n - количество селекционируемых признаков в комплексном селекционном индексе;
i - порядковый номер селекционируемого признака в комплексном селекционном индексе.
3. Весовые коэффициенты селекционируемых признаков (Ki) комплексной оценки племенной ценности (всех половозрастных групп свиней) рассчитываются согласно теории селекционного индекса (в матричной форме):
Pb = Gw,
где:
P - матрица фенотипических корреляций (ковариаций) m x m селекционируемых признаков, входящих в комплексный селекционный индекс;
b - столбец вектора m весовых коэффициентов селекционируемых признаков (Ki) комплексного селекционного индекса;
G - матрица m x n генетических корреляций (ковариаций) по m селекционируемых признаков, входящих в комплексный индекс, и n селекционируемых признаков, входящих в агрегатный генотип (совокупность признаков, для которых рассчитаны весовые коэффициенты экономической значимости (идеальное с экономической точки зрения животное));
w - вектор целей разведения (весовых коэффициентов экономической значимости Wi) признаков в агрегатном генотипе.
Отсюда весовые коэффициенты селекционируемых признаков комплексных селекционных индексов равны:
b = P-1Gw.
Элементы P и G находим:
- фенотипическая варианса (диагональные элементы матрицы P);
- фенотипическая коварианса (внедиагональные элементы матрицы P);
- генетическая варианса (диагональные элементы матрицы G);
- генетическая коварианса (внедиагональные элементы матрицы G).
4. Весовые коэффициенты селекционируемых признаков (Ki) рассчитываются так, чтобы отбор индивидуумов по комплексному селекционному индексу (И) максимизировал ответ по агрегатному генотипу (Wi).