The objective of this research was to identify the mathematical function that best fits the lactation curve of Guzerá and Sindi cattle. 840 lactations records (467-Guzerá, 373-Sindi) were used from animals raised at Carnauba ranch, in the town of Taperoá, Paraíba. Milk production control were taken every in 35 days, between years 1987 and 2004. Lactation records used contained a minimum of five and a maximum of ten controls. Six mathematical functions were fitted to the mean curve and individual lactations: polynomial inverse (PI), hyperbolic lineal (LH), gamma incomplete (GI), quadratic logarithmic (QL), lineal (L) and quadratic (Q) functions were fitted using an iteration process. Criteria used to verify fitting quality were: determination coefficient (Ra²), percent of deviation between the total observed estimated productions, and percent of typical curves. Residual distribution graph were used only to evaluate fitting of mean curve. Ra² for mean curve fitting were always larger than 0.93. Good fits were also obtained (Ra²>0.80), in 61.4% and 66.7% of the lactations fitted by PI and GI functions, respectively, for Guzerá cattle. For Sindi cattle, these values were 57.2% and 50.0%. PI and GI, respectively, were closest to the production oscillations throughout lactation, even though estimated deviations between total productions were observed and estimated that were similar to the other functions. All studied functions fitted well means curve data. For individual lactations, however, PI and GI give better estimates of total milk yields.
O objetivo deste trabalho foi comparar algumas funções matemáticas citadas na literatura, identificar a que melhor descreve a produção de leite e estimar alguns componentes da curva de lactação de vacas Guzerá e Sindi criadas na região semi-árida do estado da Paraíba. Foram utilizadas 840 lactações (467-Guzerá; 373-Sindi) de animais pertencentes à Fazenda Carnaúba, município de Taperoá, Paraíba. O controle leiteiro foi feito em intervalos de 35 dias, no período de 1987 a 2004. As lactações utilizadas continham no mínimo cinco e no máximo doze controles. O ajuste foi feito para a curva média e para lactações individuais, utilizando-se seis funções matemáticas. A função polinomial inversa (PI), linear hiperbólica (LHE), gama incompleta (GI), quadrática logarítmica (QL), linear (L) e quadrática (Q) foram ajustadas usan-do-se processos iterativos por meio de regressão não-linear. Os critérios utilizados para verificar a qualidade do ajuste para cada função foram: coeficiente de determinação ajustado (Ra²), percentagem de desvios entre as produções totais observadas e estimadas e percentagem de curvas típicas. Os gráficos de distribuição de resíduos foram utilizados apenas para avaliar o ajuste para a curva média. Para a curva média, os valores de Ra² foram superiores a 0,93 para todas as funções. Bons ajustes baseados nos Ra²>0,80 foram obtidos, respectivamente, para 61,4% e 66,7% das lactações ajustadas pelas funções PI e GI para vacas Guzerá. Para vacas Sindi, esses valores foram de 57,2% e 55,0%, indicando boa qualidade no ajuste. As funções PI e GI acompanharam mais de perto as oscilações da produção ao longo da lactação. Para curva média, todas as funções poderiam ser utilizadas, pois os desvios foram muito próximos. Para as lactações individuais, melhores estimativas dos componentes da curva de lactação poderiam ser obtidos pelas funções PI e GI.