Na raça Guzerá, BARBOSA et al. (1986), estudando 517 lactações de 142 vacas, encontraram efeito significativo do mês do parto sobre a produção de leite, porém observaram que as maiores produções foram de vacas que tiveram parto nos primeiros meses da estação seca do ano.

BIANCHINI SOBRINHO e DUARTE (1988), por intermédio da função linear hiperbólica, mostraram que a forma da curva de lactação das vacas Gir foi alterada pela estação de parto e pelo número de ordenhas. No entanto, LOPES et al. (1994 e1996) e MACMANUS et al. (1997) não confirmaram o efeito (P>0,05) da estação de parição sobre os parâmetros das funções exponencial parabólica, gama e quadrática, respectivamente.

QUEIROZ et al. (1991) encontraram efeito (P>0,05) da idade da vaca ao parto sobre a produção de leite de 672 vacas da raça Holandesa, distribuídas em sete rebanhos. Os autores relataram que o modelo polinomial de segundo grau é o que melhor explica este efeito. Na raça Gir, SOUZA et al. (1996) observaram que a produção de leite foi influenciada pela idade da vaca ao parto. JUNQUEIRA et al. (1997), com 2036 lactações de vacas da raça Holandesa, verificaram que todos os parâmetros da função gama incompleta foram influenciados pela idade da vaca ao parto.

O presente trabalho objetivou identificar, entre as diversas funções matemáticas citadas na literatura, as que melhor descrevem a produção de leite das vacas da raça Guzerá; estudar os efeitos dos fatores de ambiente sobre a forma da curva de lactação; e estimar os parâmetros genéticos para a produção de leite total e os parâmetros da curva de lactação que melhor se ajustou aos dados de produção de leite das vacas da raça Guzerá.

O arquivo inicial continha 1584 lactações e, após a exclusão dos registros com informações anormais que poderiam comprometer o presente estudo, restaram os registros de 1130 lactações de 583 vacas Guzerá, filhas de 165 touros, que tiveram suas produções observadas no período de 1983 a 1997. Estes dados são provenientes de nove rebanhos, distribuídos nas regiões Nordeste e Sudeste do Brasil, os quais fazem parte do Arquivo Zootécnico Nacional, mantido pela Embrapa Gado de Leite, Juiz de Fora - MG, e são gerados em parceria com a Associação Brasileira dos Criadores de Zebu (ABCZ).

Os parâmetros das 22 funções matemáticas utilizadas neste estudo (Tabela 1) foram estimados por meio de regressão não-linear da produção de leite, em função do número de dias do parto até a avaliação do controle leiteiro, utilizando-se o método de Gauss Newton. Deve-se considerar, também, que, embora os símbolos usados para representar os parâmetros das diferentes funções sejam os mesmos, eles não têm necessariamente a mesma interpretação. Após o ajuste das funções, procederam-se às análises para a escolha dos modelos que melhor se ajustaram às produções de leite das vacas. A escolha baseou-se no percentual de lactações em que os parâmetros das funções foram significativos e no número de curvas atípicas encontradas para cada modelo. Nas etapas posteriores, foram avaliadas as influências dos fatores não-genéticos sobre a produção de leite total e os parâmetros da função que melhor descreveu as curvas de lactações. Em seguida, foram estimados os coeficientes de herdabilidade e repetibilidade para as características produção de leite total, produção inicial de leite "a" e taxa de declínio da produção de leite "c".

O modelo estatístico usado para verificar quais fatores de ambiente poderiam influenciar a forma da curva de lactação foi:

em que

Para estimar os componentes de variância e covariância necessários para a estimativa dos parâmetros genéticos para as características em estudo, foi utilizado o método da Máxima Verossimilhança Restrita (REML), por intermédio do programa MTDFREML (Multiple Trait Derivate Free Restricted Maximum Likelihood - BOLDMAN et al., 1995), utilizando-se o seguinte modelo:

em que é o vetor das variáveis dependentes (produção de leite total, produção inicial e taxa de declínio da produção); X, a matriz de incidência de efeitos fixos; , o vetor de efeitos fixos (rebanho-ano-estação de parto, idade da vaca ao parto com termos linear e quadrático); Z, a matriz de incidência dos efeitos genéticos diretos; , o vetor dos efeitos aleatórios dos valores genéticos diretos do animal; W, a matriz de incidência dos efeitos permanentes de ambiente, o vetor dos efeitos permanentes de ambiente; e e, o vetor dos erros aleatórios associados a cada observação.

A média de produção de leite por lactação (período médio de 290 dias) foi de 2359 kg, com coeficiente de variação de 29,51%. O formato da curva de lactação apresentou-se linear, com pico de produção entre o primeiro e o segundo mês de lactação (Figura 2). Formato similar foi obtido por MADALENA et al. (1979) e GONÇALVES et al. (1996), para as lactações de vacas mestiças e puras da raça Gir, respectivamente.

Constam da Tabela 1 as médias dos coeficientes de determinação e o percentual médio de lactações em que os parâmetros dos modelos foram significativos (P<0,05). As estimativas dos diferentes parâmetros, assim como o percentual de curvas atípicas encontradas para os modelos, são apresentadas na Tabela 2. Os percentuais de curvas atípicas obtidos neste trabalho foram inferiores aos encontrados na literatura, à exceção dos resultados de SHANKS et al. (1981) e LOPES et al. (1996), que obtiveram, respetivamente, 0,73 e 1,92% de lactações com formas atípicas. Essas são caracterizadas por apresentarem estimativas negativas para os parâmetros da curva de lactação. Assim, estimativas negativas para o parâmetro "a" implicam que a produção inicial de leite estimada seja menor que zero.

Com base nos resultados apresentados nas Tabelas 1 e 2, conclui-se que o modelo 19 descreveu melhor a curva de lactação das vacas da raça Guzerá, seguido, na ordem, pelos modelos 22, 8 e 1. A representação gráfica da curva de lactação estimada para esses modelos é apresentada na Figura 1.

A análise de variância da produção de leite e dos parâmetros da curva de lactação encontra-se na Tabela 3. Observa-se que o efeito de rebanho, ano e idade da vaca ao parto influenciou (P<0,05) a produção de leite e as estimativas dos parâmetros da curva de lactação. Não houve variação (P>0,05) da produção inicial e produção de leite total nas diferentes estações de parto. Isto pode ter ocorrido em virtude do mesmo manejo e alimentação nas duas estações. Resultados semelhantes foram observados por LOPES et al. (1994), MACMANUS et al. (1997) e JUNQUEIRA et al. (1997), ajustando as lactações de vacas da raça Holandês, por meio das funções exponencial parabólica, gama incompleta e quadrática, respectivamente.

A curva de lactação média observada para as duas estações de parto estão representadas na Figura 2. Na Tabela 4, são apresentadas as médias da produção de leite observada e estimativas da produção inicial e taxa de declínio da produção para os diferentes rebanhos e estações de parto. A variação destas estimativas, segundo o ano de parto, é representada na Figura 3. As oscilações observadas podem ter como causa as variações climáticas, de manejo e da alimentação dos animais, durante esses anos. As variações desses parâmetros, no decorrer dos anos, foram observadas por LOPES et al. (1994), JUNQUEIRA et al. (1997) e GONÇALVES et al. (1997).

Nas Figuras 4, 5 e 6, são apresentadas, respectivamente, a variação da produção de leite, produção inicial e taxa de declínio da produção, em função da idade da vaca ao parto. As variações que ocorrem com o avançar da idade da vaca são, principalmente, causadas por fatores fisiológicos e proporcionam desempenhos máximos com a maturidade do animal. Observa-se, nas Figuras 5 e 6, que as vacas mais jovens apresentaram menores produção inicial e taxa de declínio da produção, quando comparadas com as mais velhas.

Na Tabela 5, são apresentadas as estimativas dos coeficientes de herdabilidade e repetibilidade para as características produção de leite total, produção inicial e taxa de declínio da produção obtidas por meio de análises com característica única. Nas análises com duas características, obtiveram-se duas estimativas de herdabilidade e repetibilidade para a característica produção de leite total, em alguns casos com valores diferentes, tendo sido calculada a média destas estimativas (Tabela 6). As estimativas de herdabilidade obtidas para os quatro modelos selecionados indicam que a produção inicial de leite tem maior potencial para seleção que a taxa de declínio da produção.

Em estudos que visam ajustar equações de regressão aos dados de produções de leite de vacas da raça Guzerá, indica-se utilizar uma das seguintes funções (y= a n e ^-cn), (y= a - cn + ln(n)), (y= a -cn) ou (y= a e^{-c n}).

Em conseqüência da baixa herdabilidade da taxa de declínio da produção de leite, o formato da curva de lactação das vacas da raça Guzerá não pode ser alterado pela seleção com base nesse componente da curva de lactação.

BARBOSA, S.B.P., PEREIRA, C.S., FONSECA, C.G. et al. 1986. Estudo de características produtivas de um rebanho Guzerá na região agreste de Pernambuco. I - Produção de leite. Arq. Esc. Med. Vet. Zootec., 38(3):381-389.

BATRA, T.R., LIN, C.Y., McALLISTER, A.J. et al. 1987. Multitrait estimation of genetic parameters of lactation curves in Holstein heifers. J. Dairy Sci., 70(10):2104-2111.

BIANCHINI SOBRINHO, E. Estudo da curva de lactação de vacas da raça Gir. Ribeirão Preto,SP: FMVRP/USP, 1984, 88p. Dissertação (Doutorado em Genética) Faculdade de Medicina Veterinária de Ribeirão Preto, Universidade de São Paulo, 1984.

BIANCHINI SOBRINHO, E., DUARTE, F.A.M. 1988. Genetic and environmental aspects of the linear hyperbolic lactation curve. Rev. Bras. Genet., 11(3):671-678.

BOLDMAN, K.G., KRIESE, L.A., VAN VLECK, L.D. 1995. A Manual for use of MTDFREML; a set programs to obtain estimates of variances and covariances [DRAFT]. Lincoln: Department of Agriculture, Agricultural Research Service. 115p.

BRODY, S., RAGSDALE, A.C., TURNER, C.W. 1923. The rate of decline of milk secretion with the advance of the period of lactation. J. Gen. Physiol., 5:441-444.

BRODY, S., TURNER, C.W., RAGSDALE, A.C. 1924. The relation between the initial rise and the subsequent decline of milk secretion following parturition. J. Gen. Physiol., 6:541-545.

COBBY, J.M., LE DU, Y.L.P. 1978. On fitting curves to lactation data. Anim. Prod., 26(2):127-133.

DAVE, B.K. 1971. First lactation curve of Indian water buffalo. JNKVV Res. J., 5:93.

GALDINI, C.H., FARO, L.E, FREITAS, M.A.R. Seleção para produção de leite auxiliada pela curva de lactação de vacas mestiças. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 35, 1998, Botucatu, SP. Anais... Botucatu: SBZ, 1998, p.291.

GONÇALVES, T.M., MARTINEZ, M.L., MILAGRES, J.C. 1996. Curvas de lactação na raça Gir. I. Escolha do modelo de melhor ajuste. R. Soc. Bras. Zootec., 25(4):616-626.

GONÇALVES, T.M., MARTINEZ, M.L., MILAGRES, J.C. 1997. Curvas de lactação na raça Gir. 2. Influência dos fatores de meio ambiente, estimativas de repetibilidade e herdabilidade para os parâmetros da curva de lactação quadrática logarítmica. R. Soc. Bras. Zootec., 26(1):88-97.

GROSSMAN, M., KOOPS, W.J. 1988. Multiphasic analysis of lactation curves in dairy cattle. J. Dairy Sci., 71(6):1598-1608.

JUNQUEIRA, L.V., NEIVA, R.S., VEIGA, R.D. et al. 1997. 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. R. Soc. Bras. Zootec., 26(6):1109-1118.

LOPES, M.A., NEIVA, R.S., VALENTE, J. et al. Aplicação da função exponencial parabólica no estudo da curva de lactação de vacas da raça Holandesa variedade preto e branco. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 31, 1994, Maringá, PR. Anais... Maringá: SBZ, 1994. p.204.

LOPES, M.A., NEIVA, R.S., VALENTE, J. et al. 1996. Aplicação da função tipo gama incompleta no estudo da curva de lactação de vacas da raça Holandesa, variedade preta-branca, mantidas em sistema intensivo de produção. R. Soc. Bras. Zootec., 25(6):1086-1101.

MADALENA, F.E., MARTINEZ, M.L., FREITAS, A.F. 1979. Lactation curves of Hostein-Friesian and Holstein-Friesian x Gir cows. Anim. Prod., 29:101-107.

MCMANUS, C., GUTH, T.L.F., SAUERESSIG, M.G. Curvas de lactação em gado holandês em confinamento total no DF. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 34, 1997, Juiz de Fora, MG. Anais... Juiz de Fora: SBZ, 1997. p.74.

MOLINA, J.R., BOSCHINI, C. 1979. Adjustment of the dairy curve of a Holstein herd with a linear modal model. Agron. Costarric., 3(2):167.

NELDER, J.A. 1966. Inverse polynomials a useful group of multi-factor response functions. Biometrics, 22(1):128-141.

PAPAJCSIK, I.A., BODERO, J. 1988. Modeling lactation curves of friesian cow in a subtropical climate. Anim. Prod., 47(2):201-207.

QUEIROZ, S.A., ALBUQUERQUE, L.G, FREITAS, M.A.R. et al. 1991. Fatores genético e de meio que influenciam os componentes da curva de lactação de bovinos da raça Holandesa. Arq. Bras. Med. Vet. Zootec., 43(4):357-370.

SCHNEEBERGER, M. 1981. Inheritance of lactation curve in Swiss brown cattle. J. Dairy Sci., 64(3):475-483.

SIKKA, L.C. 1950. A study of lactation as effected by buredity and environment. J. Dairy Res., 17(3):231-252.

SINGH, R.P., GOPAL, R. 1982. Lactation curves analysis of buffaloes maintained under village conditions. Ind. J. Anim. Sci., 52(12):1157-1160.

SOUZA, E.M., MILAGRES, J.C., REGAZZI, A.J. et al. 1996. Efeitos de fatores genéticos e de meio ambiente sobre a produção de leite em rabanhos de Gir leiteiro. R. Soc. Bras. Zootec., 25(4):889-901.

WOOD, P.D.P. 1967. Algebric model of lactation curve in cattle. Nature, 216(5111):164-165.

² Zootecnista, MS. Bolsista da CAPES.

³ Professor da Universidade Federal de Viçosa.

⁴ Pesquisador do CNPGL/EMBRAPA, Juiz de Fora, MG.

⁵ Professor da Universidade Federal de Minas Gerais.