Ajustamento para heterogeneidade de variância da produção de leite de vacas da Raça Holandesa no Est (página 2)

Os dados usados fazem parte dos arquivos da Associação de Criadores de Gado Holandês de Minas Gerais. Duas estações de parto foram consideradas: outubro a março e abril a setembro. Na edição dos dados, as classes rebanho-ano-estação com menos de dez lactações foram eliminadas. As produções de leite até 305 dias de 25.688 vacas com partos entre 1988 e 1999 foram pré-ajustadas para idade e estação de parto por meio de fatores multiplicativos, conforme Durães et al. (1999).

As variâncias das produções podem variar dentro de rebanho e entre anos para um determinado rebanho. A variação entre os anos ocorre porque as produções tendem a aumentar ao longo dos anos, e espera-se que as variâncias aumentem com a produção. Estimaram-se variâncias para cada subclasse rebanho-ano, por meio de:

em que: = variância estimada dentro de rebanho-ano ij; y_ijk = produção do animal k no rebanho-ano ij; n_ij = número de produções no rebanho-ano ij;

Em uma análise preliminar ajustou-se o modelo linear fixo seguinte para identificar fontes de variação dos desvios-padrão fenotípicos dentro de cada subclasse rebanho-ano:

em que: D_ijklmno = desvio-padrão fenotípico da produção de leite até 305 dias ajustada para idade e estação no ano i, rebanho j, mês k, grupo genético l, classe de número de observação m, classe de produção média de leite n; Ano_i = efeito fixo de ano de parto i ( i = 1988, ..., 1999); Reb_j = efeito fixo de rebanho j ( j = 1,..., 385); Mês_k = efeito fixo de mês de parto k ( k = 1,...,12); CR_l = efeito fixo de grupamento genético l ( l = 1,...,6); CLASNOB_m = efeito fixo de classe de número de observações m ( m = 1,..., 6); CLASLEI_n = efeito fixo de classe de produção média de leite n ( n = 1,..., 8); e_ijklmno = efeito aleatório de erro.

Uma vez que o número de lactações em cada subclasse rebanho-ano pode ser pequeno, os erros de amostragem das estimativas das variâncias podem ser grandes. Foram então estimadas variâncias ponderadas dentro de rebanho-ano , com menor erro de amostragem, combinando-se informação de anos adjacentes, isto é, a variância dentro de ano e as de rebanho-anos adjacentes conforme procedimento sugerido por Wiggans & Van Raden (1991) e também usado por Van Der Werf et al. (1994). Os pesos para a ponderação das diferentes fontes de informação foram os mesmos de Wiggans & Van Raden (1991). As produções foram, em seguida, padronizadas por meio de:

As produções padronizadas foram, então, expressas na forma original usando-se como um valor constante a média geral das variâncias dentro de rebanho-ano. Para se obterem os valores genéticos, utilizou-se um modelo linear misto (modelo animal) contemplando os efeitos fixos de idade ao parto, rebanho, ano e estação de parto, grupos genéticos (PO e PC) e os efeitos aleatórios genético aditivo de animal e de meio ambiente permanente. O sistema MTDFREML (Boldman et al., 1995) foi usado.

Critérios para verificação de efeitos do ajustamento para heterogeneidade da variância foram: 1) correlações de ordem para touros e vacas; 2) mudanças de vacas elite entre classes de desvio-padrão; 3) tendência genética na população de vacas.

Na Tabela 1 encontram-se número de subclasses rebanho-ano, número de lactações, médias das produções e dos desvios-padrão dentro das subclasses por ano e mês de parto e grupo genético. Houve aumento significativo da média dos desvios-padrão de 1988 até 1997, fato ocorrido provavelmente devido ao aumento na produção no decorrer do período. As mudanças dos desvios-padrão com o mês de parto não foram significativas. À medida que os animais passaram a ser mais controlados quanto aos seus pedigrees, constatou-se maior produção de leite e aumento da média dos desvios-padrão, concordando com Weigel et al. (1993), que verificaram variâncias residuais para características produtivas maiores em rebanhos registrados. Observou-se aumento de 100% nas médias dos desvios-padrão entre os níveis extremos de produção (Tabela 2), indicando ser o nível de produção o principal fator de variabilidade dos rebanhos. Houve, também, aumento significativo na média dos desvios-padrão com o número de animais do rebanho, concordando com Weigel et al. (1993), os quais encontraram variâncias residuais maiores para produção de leite em rebanhos maiores.

Pela Tabela 3 e Figura 1, observam-se mudanças nas Capacidades Previstas de Transmissão (PTAs) dos touros e das vacas com o ajustamento para heterogeneidade da variância. As correlações de ordem entre PTAs obtidas antes e após o ajustamento para heterogeneidade, para todos os animais avaliados foram altas, respectivamente, 0,950 e 0,976 para touros e vacas. Quando os touros apresentavam PTA com maior confiabilidade e possuíam filhas em muitos rebanhos não foram afetados pela correção para heterogeneidade da variância. Entretanto, quando possuíam filhas em poucos rebanhos, as diferenças de ordem foram mais pronunciadas. Como constatado por Hill (1984), Wiggans & Van Raden (1991), o ajustamento apresentou efeitos sobre os valores genéticos das vacas elite. As correlações foram muito menores até negativa para as mesmas, indicando a ocorrência de um importante reordenamento à medida em que se reduziu a amostra. Na Tabela 3 tem-se também o número de novos animais entre os elite. Quando se ajustou para heterogeneidade, aproximadamente 40% das vacas elite foram substituídas. Como relatado também por Wiggans & Van Raden (1991), a origem das vacas elite foi também afetada (Tabela 4). Antes originavam-se de rebanhos com maior variabilidade e com o ajustamento houve uma troca de vacas destes para rebanhos de variabilidade média. Rebanhos com maior variabilidade contribuíram com 35 das 100 melhores vacas quando não houve correção para heterogeneidade e 16 quando houve correção.

Proporcionalmente, rebanhos com maior variabilidade mantiveram a contribuição dos rebanhos com baixa variabilidade, uma vez que alta variabilidade estava associada a um alto nível de produção e a um alto nível genético médio.

Como já foi sugerido (Winkelman & Schaeffer, 1988; Boldman & Freeman, 1990), o efeito do ajustamento para heterogeneidade da variância sobre os valores genéticos dos touros foi pequeno. Entretanto, foram significativas as mudanças nos valores genéticos de vacas, sugerindo que possivelmente as PTAs de vacas antes do ajustamento estivessem sujeitas a vícios, devido à heterogeneidade da variância. Estas mudanças eram esperadas, uma vez que as produções de uma vaca ocorrem, na sua maioria, em um mesmo rebanho, pelo que a sua PTA está mais sujeita a vícios devido a diferenças de variabilidade dentro do rebanho.

O efeito do ajustamento sobre a tendência genética para a população de vacas é ilustrado na Figura 2. No período de 1983 a 1997, o valor estimado da tendência aumentou de aproximadamente 12 kg/ano com o ajustamento para heterogeneidade.

Com o ajustamento para heterogeneidade de variância usado, constataram-se melhorias na avaliação genética. Mudanças de rebanhos que contribuíram com vacas elite para aqueles com variância média deverão possibilitar avaliações mais precisas e justas para identificação das melhores vacas. O procedimento deverá ser adotado nas futuras avaliações genéticas para produção de leite, na raça Holandesa, no Estado de Minas Gerais, até que novas pesquisas sugiram outro mais adequado.

BOLDMAN, K.G.; FREEMAN, A.E. Adjustment for heterogeneity of variances by herd production level in dairy cow and sire evaluation. Journal of Dairy Science, v.73, n.2, p.503-512, 1990.

BOLDMAN, K.G., L.A.; KRIESE, L., D.; VAN VLECK, C.P. et al. A manual for use of MTDFREML. A set of programs to obtain estimates of variances and covariances [DRAFT]. Beltsville: U.S. Department of Agriculture, Agricultural Research Service, 1995.

DONG, M.C.; MAO, I.L. Heterogeneity of (co)variance and heritability in different levels of intra-herd milk production variance and of herd average. Journal of Dairy Science, v.73, n.3, p.843-851, 1990.

DURÃES, M.C.; FREITAS, A.F.; TEIXEIRA, N.M. et al. Fatores de ajustamento para a produção de leite e gordura na raça Holandesa para rebanhos do Estado de Minas Gerais. Revista Brasileira de Zootecnia, v.28, n.6, p.1231-1238, 1999.

GIANOLA, D. On selection criteria and estimation of parameters when the variance is heterogeneous. Theoretical and Applied Genetics, v.72, p.671-677, 1986.

HENDERSON, C.R. Applications of linear models in animal breeding. Guelph: University Guelph, 1984. 439p.

HILL, W.G. On selection among groups with heterogeneous variance. Animal Production, v.39, n.3, p.473-477, 1984.

MEUWISSEN, T.H.E.; Van DER WERF, J.H.J. Impact of heterogeneous within herd variances on dairy cattle breeding schemes. Livestock Production Science, v.33, n.1, p.31-41, 1993.

MIRANDE, S.L.; Van VLECK, L.D. Trends in genetic and phenotypic variances for milk production. Journal of Dairy Science, v.68, n.9, p.2278-2286, 1985.

RAMOS, A.A.; VALENCIA, E.F.T.; WECHSLER, F.S. et al. Heterogeneidade da variância das características de produção de bovinos da raça Holandesa no trópico. I. estratificação por nível de produção de rebanho. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 33., 1996, Fortaleza. Anais... Fortaleza: Sociedade Brasileira de Zootecnia, 1996. p.71-73.

ROBERT-GRANIÉ, C.; BONAÏTI, B.; BOICHARD, D. et al. Accounting for variance heterogeneity in French dairy cattle genetic evaluation. Livestock Production Science, v.60, n.2, p.343-357, 1999.

TORRES, T.R.; BERGMAN, J.A.G.; COSTA, C.N. et al. Ajustamento para heterogeneidade de variância para produção de leite entre rebanhos da raça Holandesa no Brasil. Revista Brasileira de Zootecnia, v.28, n.2, p.295-303, 1999.

VALÊNCIA, E.F.T.; RAMOS, A.A.; WECHSLER, F.S. et al. Heterogeneidade da variância das características de produção de bovinos da raça Holandesa no trópico. II. Estratificação por nível de produção de vacas. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 33., 1996, Fortaleza. Anais... Fortaleza: Sociedade Brasileira de Zootecnia, 1996. p.68-70.

WEIGEL, K.A.; GIANOLA, D.; YANDELL, B.S. et al. Identification of factors causing heterogeneous within-herd variance components using a structural model for variances. Journal of Dairy Science, v.76, n.5, p.1466-1478, 1993.

WIGGANS, G.H.; Van RADEN, P.M. Method and effect of adjustment for heterogeneous variance. Journal of Dairy Science, v.74, n.2, p.4350-4357, 1991.

WINKELMAN, A.; SCHAEFFER, L.R. Effect of heterogeneity of variance on dairy sire evaluation. Journal of Dairy Science, v.71, n.11, p.3033-3039, 1988.