Avaliação dos fatores oficiais de correção da produção até o primeiro controle leiteiro em vacas Holandesas

Evaluation of correction factors used to estimate milk production from calving to the first monthly control

Resumos

Para o cálculo oficial do desempenho de vacas Holandesas no Brasil, com base em controles leiteiros mensais, são utilizados os fatores de correção propostos por Everett & Carter (1968) para estimar a produção parcial do parto ao primeiro controle diário de uma dada lactação. A adequacidade destes fatores foi avaliada em 262 lactações oriundas de dois rebanhos que adotavam controle diário da produção de leite por vaca. Os dados foram divididos em classes com base na idade ao parto (CIP): CIP 1 (até 3 anos), CIP 2 (3,1 a 4 anos), CIP 3 (4,1 a 5 anos), CIP 4 (5,1 a 6 anos) e CIP 5 (acima de 6,1 anos) com 104, 56, 44, 18 e 40 lactações, respectivamente. A produção de leite do parto até um dia pré-definido da lactação foi estimada com os fatores oficiais e assumindo que o primeiro controle da produção naquela lactação ocorreu nos dias 6º, 20º, 34º, 48º, 62º ou 75º. Os dados estimados foram comparados aos dados reais calculados com a pesagem diária do leite. Em todas as CIP, os fatores oficiais superestimaram a produção acumulada até o primeiro controle e o erro relativo de superestimação decresceu com o avançar do dia do primeiro controle. A quantidade superestimada de leite no período aumentou com o avançar do dia do primeiro controle e foi 13,9 e 87,1 kg na CIP 1 e 20,6 e 196,3 kg na CIP 2 para os dias 6º e 75º, respectivamente. As CIPs 3, 4 e 5 apresentaram valores intermediários. Novos fatores para dias do primeiro controle variando de 6 a 75 e nas 5 CIPs foram gerados pela mesma metodologia descrita por Everett & Carter (1968). Os fatores gerados foram menores que os importados e rotineiramente utilizados para cálculo oficial de lactações no Brasil.

vacas Holandesas; Brasil; produção de leite

MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Avaliação dos fatores oficiais de correção da produção até o primeiro controle leiteiro em vacas Holandesas

¹

Evaluation of correction factors used to estimate milk production from calving to the first monthly control

Nilson Nunes Morais Júnior^I; Rilke Tadeu Fonseca de Freitas^II; Marcos Neves Pereira^II; Antonio Ilson Gomes de Oliveira^III; Idalmo Garcia Pereira^IV; Júlio Sílvio de Souza Bueno Filho^II; Tarcísio de Moraes Gonçalves^V

^IM.Sc. em Zootecnia, Professor da Escola Agrotécnica Federal de Colatina-ES, BR 259, km 17, C.P. 256, CEP 29709-910. E.mail: [e-mail]

^IIProfessor Adjunto da UFLA. E.mail: [e-mail]

^IIIProfessor Titular aposentado da UFLA. E.mail: [e-mail]

^IVProfessor Visitante da UFLA. E.mail: [e-mail]

^VProfessor Assistente da UFLA. E.mail: [e-mail]

RESUMO

Para o cálculo oficial do desempenho de vacas Holandesas no Brasil, com base em controles leiteiros mensais, são utilizados os fatores de correção propostos por Everett & Carter (1968) para estimar a produção parcial do parto ao primeiro controle diário de uma dada lactação. A adequacidade destes fatores foi avaliada em 262 lactações oriundas de dois rebanhos que adotavam controle diário da produção de leite por vaca. Os dados foram divididos em classes com base na idade ao parto (CIP): CIP 1 (até 3 anos), CIP 2 (3,1 a 4 anos), CIP 3 (4,1 a 5 anos), CIP 4 (5,1 a 6 anos) e CIP 5 (acima de 6,1 anos) com 104, 56, 44, 18 e 40 lactações, respectivamente. A produção de leite do parto até um dia pré-definido da lactação foi estimada com os fatores oficiais e assumindo que o primeiro controle da produção naquela lactação ocorreu nos dias 6^o, 20^o, 34^o, 48^o, 62^o ou 75^o. Os dados estimados foram comparados aos dados reais calculados com a pesagem diária do leite. Em todas as CIP, os fatores oficiais superestimaram a produção acumulada até o primeiro controle e o erro relativo de superestimação decresceu com o avançar do dia do primeiro controle. A quantidade superestimada de leite no período aumentou com o avançar do dia do primeiro controle e foi 13,9 e 87,1 kg na CIP 1 e 20,6 e 196,3 kg na CIP 2 para os dias 6^o e 75^o, respectivamente. As CIPs 3, 4 e 5 apresentaram valores intermediários. Novos fatores para dias do primeiro controle variando de 6 a 75 e nas 5 CIPs foram gerados pela mesma metodologia descrita por Everett & Carter (1968). Os fatores gerados foram menores que os importados e rotineiramente utilizados para cálculo oficial de lactações no Brasil.

Palavras-chave: vacas Holandesas, Brasil, produção de leite

ABSTRACT

Official Holstein cow performance data in Brazil, based on monthly milk weights, is calculated using the correction factors proposed by Everett & Carter (1968) to estimate the partial production from calving to the first daily control in a given lactation. The adequacy of these factors was evaluated using 262 lactations from two herds with daily control of milk production per cow. Data were divided into classes based on age at calving (CIP): CIP 1 (up to 3 years), CIP 2 (3.1 to 4 years), CIP 3 (4.1 to 5 years), CIP 4 (5.1 to 6 years) and CIP 5 (above 6.1 years) which included 104, 56, 44, 18 and 40 lactations, respectively. Milk yield from calving to a set day of lactation was estimated with the official factors and assuming that the first control of production of that given lactation occurred on days 6^th, 20^th, 34^th, 48^th, 62^nd or 75^th. Estimated data was compared with real data calculated from daily milk weights. For all CIPs, official factors overestimated the production from calving to first control and the relative overestimation error decreased with advance in the first control day. The overestimated amount of milk in that period increased with advance in the first control day and was 13.9 and 87.1 kg in CIP 1 and 20.6 and 196.3 kg in CIP 2 for days 6^th and 75^th, respectively. CIPs 3, 4, and 5 had intermediate values. New factors for first control days ranging from 6 to 75 and in all 5 CIPs were generated using the same methodology described by Everett & Carter (1968). The generated factors were smaller than the ones imported and routinely used to calculate official lactation records in Brazil.

Key words: Holstein cows, Brazil, milk production

Introdução

O controle da produção de leite em animais da raça Holandesa vem sendo supervisionado desde 1946 pela Associação Brasileira de Criadores - ABC e, mais recentemente, pela Associação Brasileira de Criadores de Bovinos da Raça Holandesa - ABCBRH, ambas oficialmente reconhecidas pelo Ministério da Agricultura (Machado, 1997). Esse controle tem como finalidades nacionais e institucionais: aquilatar o potencial produtivo, coletar dados estatísticos para teste de progênie e melhoramento genético e divulgar os resultados atualizados para a comunidade.

No Brasil, o controle leiteiro oficial é realizado mensalmente e as pesagens são efetuadas em intervalos regulares, entre controles consecutivos, sendo permitidos intervalos de 15 a 45 dias. O primeiro controle da lactação é realizado oficialmente entre o 6^o e 75^o dia após o parto.

O total de leite produzido em cada lactação é estimado por uma equação matemática contendo fatores propostos por Everett & Carter (1968). Como na maioria das vezes não é possível iniciar o controle logo após o parto, nessa equação são utilizados fatores de correção que permitem estimar a produção de leite da vaca do parto até o dia do 1o controle. Dessa forma, a produção de leite no primeiro controle é multiplicada pelo número de dias que o separa do parto e pelo fator de ajuste para o mesmo. Esses fatores de correção foram estimados por Everett & Carter (1968), baseando-se em controles diários de 403 lactações completas, realizadas entre junho de 1964 e outubro de 1967, no rebanho de gado Holandês da Universidade de Cornell -EUA.

Considerando-se as condições climáticas brasileiras e as características zootécnicas próprias de nossos rebanhos no que se refere aos sistemas de produção, alimentação e manejo, níveis de produção e composição genética, é possível que o emprego desses fatores de ajuste em controles leiteiros no Brasil possam estar causando erros no cálculo da produção total na lactação das vacas controladas.

Objetivou-se, com este estudo, avaliar os efeitos do emprego direto dos fatores de correção oficiais na estimação da produção de leite do parto até o dia do primeiro controle, comparada a produções reais obtidas de dados diários de produção de vacas da raça Holandesa sob controle leiteiro no Brasil.

Material e Métodos

Os dados utilizados neste estudo são originários de controles diários de produção provenientes de dois rebanhos da raça Holandesa, que participam do controle leiteiro oficial junto à Associação de Criadores de Gado Holandês de Minas Gerais.

O rebanho denominado "Cova da Onça", localizado no município de Pedro do Rio - RJ, participou neste estudo com 163 lactações oriundas do período de 04/94 a 12/99. Do segundo rebanho, denominado "EPAMIG", localizado no município de Três Pontas - MG, foram utilizadas 99 lactações correspondentes ao período de 04/97 a 12/2000.

O estudo foi conduzido no Departamento de Zootecnia da Universidade Federal de Lavras. A fase inicial constou da verificação dos dados diários de produção, eliminando-se lactações com falhas ou alterações drásticas de produção nos primeiros 75 dias por motivos conhecidos (mastite graves ou dados de produção incompletos). Eliminaram-se também lactações com duração inferior a 90 dias. Foram, então, eliminadas 44 lactações, sendo 29 do rebanho "Cova da Onça" e 15 do rebanho "Epamig". Após essa etapa de consistência, foram gerados cinco arquivos para análise, conforme a classe de idade ao parto (CIP), mostrada na Tabela 1.

Com base nos dados diários reais dos primeiros 75 dias de lactação de cada arquivo, calculou-se a produção acumulada real dia-a-dia para esse período.

A produção acumulada estimada com o uso dos fatores foi obtida conforme método oficial (BRASIL, 1986) definido com base no trabalho de Everett & Carter (1968). Utilizou-se parte da fórmula do cálculo da quantidade total de leite produzida em uma lactação, em que a quantidade de leite produzida ao primeiro controle é obtida pelas fórmulas:

em que: PLT1ºc é a produção total de leite acumulada até o primeiro controle (em kg); C₁ é a quantidade de leite (kg) obtida no primeiro controle da lactação; E₁ é o intervalo em dias, entre a data do parto e a data do primeiro controle da lactação, podendo variar de 6 a 75 dias; F₁ é o fator de correção para o primeiro controle da lactação; µ, b₁, b₂, b₃ e a são parâmetros do modelo, conforme a idade ao parto (Tabela 2).

Como cada CIP utiliza diferentes fatores de correção para o primeiro controle, foram utilizados, em cada arquivo de dados, os respectivos fatores (F1) para o cálculo da produção acumulada, conforme a CIP. Calculou-se, assim, a produção acumulada até o 6o dia após o parto, até o 7^o dia e, assim, sucessivamente, até uma produção acumulada ao 75^o dia de lactação para todas as vacas em cada arquivo.

Para avaliar os fatores de Everett & Carter (1968), criou-se uma estrutura de análise pelo PROC MIXED (Littell et al.,1996), analisando-se, como medidas repetidas, as produções acumuladas de cada lactação ao longo dos dias. Foram utilizados nessa análise os dias de controle 6, 20, 34, 48, 62 e 75 e comparados os valores de produção de leite reais acumulados e os valores acumulados, utilizando-se os fatores de correção para todas as vacas em cada uma das classes de idade ao parto.

A estrutura de covariância utilizada foi aquela com maior valor para o critério de informação de Akaike, considerando as estruturas CS (compound symmetry), UN (unstructured) e AR(1) (first-order autoregressive). O modelo utilizado para a produção de leite acumulada foi:

em que: Y_ijk é a produção de leite acumulada da vaca i pelo método j, até o dia do controle k; µ é a média geral; V_i é o efeito da vaca i( i = 1, 2, 3,..., n)1; M_j é o efeito do método de cálculo da produção acumulada j ( j = 1 e 2); VM_ij é o efeito da interação da vaca i com o método j, sendo usado como erro para testar o efeito do método; D_k é o efeito do dia do primeiro controle k (k = 6, 20, 34, 48, 62 e 75); MD_jk é o efeito da interação do método j com o dia do primeiro controle k; e_ijk é o erro associado a cada observação.

Para o cálculo dos novos fatores, para cada classe de idade ao parto, utilizou-se a mesma metodologia de Everett & Carter (1968), sendo as produções acumuladas calculadas conforme a equação abaixo:

em que: Y_i é a média da produção de leite até o iéssimo dia da lactação; J é o número de dias do período a ser estimado (j = 6,...,75); Y_j é a média da produção de leite no jéssimo dia da lactação.

Cada fator é a proporção entre o somatório das produções médias até o dia j em relação à produção estimada pela multiplicação da produção num dado dia j por J.

Considerando-se esses fatores, estimaram-se os parâmetros por meio do método modificado de Gauss-Newton, disponível no procedimento NLIN SAS (1995). Foram utilizadas 70 variáveis por classe de idade de parto (6^o ao 75^o dia de produção de leite). Os parâmetros iniciais do modelo foram obtidos de Everett & Carter (1968) ou próximos a esses, utilizando-se a equação:

em que: F_n é o novo fator gerado; E₁ é o dia do primeiro controle (6, ..., 75); m, b₁, b₂, b₃ e a são parâmetros do novo modelo.

Resultados e Discussão

As vacas da classe de idade ao parto até 3 anos (CIP1) apresentaram produção média inicial de 7,2 (0,3) kg/dia e as vacas das demais classes apresentaram produções iniciais médias entre 9,3 e 11,4 kg de leite (Tabela 3).

Quanto à produção máxima ou produção no dia do pico, as vacas da CIP 1 apresentaram valores acima dos observados por Grossi (1999), em seu estudo com vacas de vários grupos genéticos Holandês-Zebu, e também superiores aos valores encontrados por Junqueira et al. (1997) em vacas Holandesas, com idade ao parto inferiores a 37 meses. Entretanto, estão bastante próximos aos 27,8 kg/dia estimados pela função gama incompleta obtidos por Lopes (1994) em animais da raça Holandesa de primeira ordem de lactação.

A classe de idade ao parto entre 4 e 5 anos (CIP3) apresentou a maior produção de leite no pico e a maior produção total nos primeiros 75 dias de lactação dentre as estudadas. Segundo Freitas et al. (2000), ocorre um efeito quadrático da idade da vaca ao parto sobre a produção total de leite na lactação, atingindo o máximo aos 78 meses (6,5 anos). Pires (1984) observou efeito semelhante, constatando produção máxima entre 72 e 90 meses de idade ao parto.

O dia do pico ocorreu entre 41,3 e 44,6 dias pós parto (Tabela 3), apresentando valores próximos ao limite mínimo observado por Junqueira et al. (1997), os quais, com base em 2.036 lactações completas no período entre 1989 e 1991, provenientes de registros da Associação de Criadores de Gado Holandês do Estado de Minas Gerais, em 122 rebanhos agrupados em cinco núcleos, concluíram que os picos de produção ocorreram entre 1,5 a 2,1 meses de lactação.

Os valores observados nesse período inicial (Figura 1) são bastante semelhantes aos encontrados na literatura, muito embora, não tenham sido comparados com dados diários de produção, o que demonstra que as lactações utilizadas têm comportamento lactacional próximo aos rebanhos de gado Holandês estudados no Brasil.

Para todas as cinco classes estudadas de idade ao parto, a produção acumulada real até o primeiro controle foi menor (P< 0,05) que a produção estimada pelos fatores de Everett & Carter (1968). Não foi encontrada interação (P>0,05) entre método de cálculo da produção acumulada x dia do primeiro controle, exceto para a CIP 2.

Por esses resultados infere-se que, nos rebanhos estudados, os fatores de ajuste para o primeiro controle utilizados atualmente superestimaram, em todas as classes de idade ao parto, a produção de leite acumulada. O erro relativo apresenta maiores valores no início da lactação (primeiro controle no 6o dia). Mas, como a produção acumulada é referente a apenas 6 dias, o erro é quantitativamente, pequeno (Tabela 4).

Com o avanço do dia do primeiro controle (6^o, 20^o, 34^o, 48^o, 62^o e 75^o dia de produção de leite), ocorre uma queda no erro relativo de superestimação. Entretanto, o impacto do dia do primeiro controle mais tardio na quantidade de leite superestimada é maior (Figura 2).

A classe de idade ao parto 1 (até 3 anos) apresenta uma superestimativa entre 13,4 e 87,1 kg de leite, quando o primeiro controle foi ao 6^o ou o 75^o dia, respectivamente. Na classe de idade ao parto 2, esses valores variam de 20,6 e 196,3 kg de leite nos mesmos dias de primeiro controle.

As classes 3, 4 e 5 apresentaram valores intermediários aos das classes 1 e 2 e, em geral, seguem o mesmo comportamento dessas.

Como não foram feitas comparações estatísticas entre classes de idade ao parto pelo fato de utilizarem fatores de correção diferentes, não é possível concluir se as quantidades apresentadas indicam uma superestimação diferenciada entre classe de idade ao parto.

Os parâmetros das equações para estimar os novos fatores de correção para o primeiro controle (Fn) encontram-se na Tabela 5. Por meio desses parâmetros, podem-se calcular os fatores de ajuste para o primeiro controle (apresentados na Tabela 6), segundo a equação:

em que: F_n é o novo fator de correção para o primeiro controle da lactação; D₁ é o intervalo em dias, entre a data do parto e a data do primeiro controle da lactação, podendo variar de 6 a 75 dias; µ, b₁, b₂, b₃ e a são os parâmetros do modelo, conforme idade ao parto.

As equações geradas pelos novos parâmetros apresentam ajustes elevados e o R² dos modelos foram superiores a 0,97 para as cinco classes estudadas (Tabela 5). Por meio desses valores infere-se que o modelo é capaz de gerar bons fatores de ajuste para o primeiro período de lactação. Na verdade, o modelo, ou seja, a equação proposta por Everett & Carter (1968) é confiável, desde que os dados utilizados para estimar os parâmetros sejam compatíveis aos rebanhos onde se deseja estimar a produção. No caso do Brasil, essa avaliação não foi feita e os parâmetros do modelo foram simplesmente implementados sem uma avaliação do nível de produção e do comportamento no início da lactação do rebanho nacional.

Observa-se pela Figura 3 que os novos fatores gerados são inferiores aos fatores oficiais em todas as classes de idade. Diante desse fato, pode-se afirmar que os fatores atuais estão superestimando a produção de leite até o primeiro controle. Uma vez que os fatores reais estão abaixo dos oficiais, quanto maior for a produção da vaca no dia do primeiro controle, maior a superestimativa e quanto maior o número de dias, maior a superestimação, principalmente se o controle ocorrer próximo ao pico de lactação.

As diferenças numéricas entre os fatores nas diferentes CIPs seguem uma tendência bastante uniforme nas diferentes classes de idade. As maiores diferenças entre os fatores atuais e os novos ocorrem até o 20^o dia da lactação, oscilando pouco a partir desse ponto nas cinco classes de idade ao parto.

Embora os fatores gerados neste trabalho não apresentem valores definitivos, uma vez que a amostragem foi restrita a apenas dois rebanhos, os fatores utilizados atualmente também não deveriam ser utilizados para a raça Holandesa ou generalizados para outras raças leiteiras no Brasil, porque os mesmos são originados de apenas um rebanho dos Estados Unidos, com um número relativamente pequeno de lactações e com dados de mais de 30 anos atrás.

Os fatores novos não foram comparados a outros, uma vez que nenhum estudo relativo a fatores de ajuste do primeiro controle com dados de produção diária, na raça Holandesa ou outra, foi desenvolvido no Brasil nos últimos anos.

Considerando-se que os novos fatores tenham sido estimados com base nas produções médias diárias e médias de produção acumulada para cada classe de idade ao parto e que existem variações no nível de produção dentro das classes, os novos fatores permitem boas estimativas para os animais que se aproximam mais da média de produção da sua classe de idade ao parto.

Conclusões

Os fatores de correção utilizados atualmente na fase inicial da lactação pela ABCBRH superestimam a produção de leite acumulada até o primeiro controle para todas as classes de idade ao parto nos rebanhos estudados.

A produção de leite em até 305 dias (P305) é afetada pelo dia do primeiro controle.

Os parâmetros da equação do fator de correção para o primeiro controle, estimados com base em dados de controle diário, apresentaram um bom ajuste do modelo.

Pelo que se observou, fica evidente a necessidade de mais estudos para a estimação desses fatores para as raças leiteiras no Brasil.

Agradecimento

Aos proprietários e responsáveis técnicos pelas fazendas "Cova da Onça" (Agropecuária Cova da Onça Ltda) e "EPAMIG" (Fazenda Experimental da Empresa de Pesquisa Agropecuária de Minas Gerais), pela concessão dos dados utilizados neste estudo.

Literatura Citada

Recebido em: 29/06/01

Aceito em: 09/01/02

• BRASIL. Ministério da Agricultura. Secretaria Nacional de Produção Agropecuária. Portaria nº 45, de 10 de outubro de 1986. Normas técnicas para execução do serviço de controle leiteiro em bovídeos. Diário Oficial [da] República Federativa do Brasil. Brasília, v.124, n.198, p.15532-15535, out. 1986. Seção 1.
• EVERETT, R.W.; CARTER, H.W. Accurracy of test interval method of calculating dairy herd improvement association records. Journal of Dairy Science, v.51, n.12, p.1936-1941, 1968.
• FREITAS, A.F.; DURÃES, M.C.; TEIXEIRA, N.M. Parâmetros genéticos da produção de leite de animais da raça Holandesa mantidos em sistema intensivo do tipo "free stall". Revista Brasileira de Zootecnia, v.9, n.6, p.2008-2012, 2000 (Suplemento 1).
• GROSSI, S.F. Eficiência reprodutiva e produtiva em rebanhos leiteiros monitorados por sistema de informação. Jaboticabal: Universidade Estadual Paulista, 1999. 87p. Dissertação (Mestrado em Zootecnia) - Universidade Estadual Paulista, 1999.
• JUNQUEIRA, L.V.; NEIVA, R.S.; VEIGA, R.D. et al. Estudo das curvas de lactação de vacas Holandesas de alguns rebanhos do Estado de Minas Gerais, por intermédio de uma função gama incompleta. Revista Brasileira de Zootecnia, v.26, n.6, p.1109-1118, 1997.
• LITTELL, R.C.G.A.; MILLIKEN, W.W.; WOLFINGER, R.D. SAS^® System for mixed models. Cary, North Caroline: SAS Institute, 1996. 633p.
• LOPES, M.A. Estudo da curva de lactação de um rebanho de gado Holandês variedade preto e branco mantido em um sistema intensivo de produção. Lavras: Escola Superior de Agricultura de Lavras, 1994. 42p. Dissertação (Mestrado em Zootecnia) - Escola Superior de Agricultura de Lavras, 1994.
• MACHADO, S.G. Parâmetros genéticos e de ambiente da produção de leite no dia do controle da primeira lactação de vacas da raça Holandesa. Jaboticabal: Universidade Estadual Paulista, 1997. 76p. Dissertação (Mestrado em Zootecnia) - Universidade Estadual Paulista, 1997.
• PIRES, A.V. Fatores genéticos e ambiente que afetam a produção de leite e duração da lactação de um rebanho com diferentes grupos genéticos Holandês:Guzerá. Piracicaba: Escola Superior de Agricultura "Luiz de Queiroz", 1984. 112p. Dissertação (Mestrado em Agronomia) - Escola Superior de Agricultura "Luiz de Queiroz", 1984.
• SAS INSTITUTE. SAS/STAT. User's guide. Version 6.12. 4.ed. Cary: 1995. v.2. 1686p.

1

Parte da Dissertação de Mestrado em Zootecnia.

Datas de Publicação

• Publicação nesta coleção
  08 Set 2004
• Data do Fascículo
  Abr 2002

Histórico

• Recebido
  29 Jun 2001
• Aceito
  09 Jan 2002