Artículo original

Estimación de la curva de lactancia mediante modelos matemáticos lineales y no lineales en búfalas colombianas^*

Estimation of the lactation curve by the use of mathematic lineal and non-lineal models in Colombian female buffalos

Estimação da curva de latência mediante modelos matemáticos lineares e não lineares em bufa-as colombianas

Juan Carlos Quintero Vélez¹, Jorge Serna Gallo¹, Mario Cerón Muñoz², Naudin Hurtado Lugo³, Divier Antonio Agudelo Gómez⁴

^* Proyecto financiado por Vicerrectoría de Investigación de la Universidad de Antioquia.
¹ Investigador de GRICA (Grupo de Investigación en Ciencias Agrarias) de la Facultad de Ciencias Agrarias de la Universidad de Antioquia.
² Grupo de Investigación en Genética y Mejoramiento Animal de la Universidad de Antioquia, Medellín, Colombia. Profesor de la Facultad de Ciencias Agrarias de la Universidad de Antioquia, Correspondencia: e-mail: mceronm@agronica.udea.edu.co.
³ Grupo de Investigación en Genética y Mejoramiento Animal de la Universidad de Antioquia. Profesor, Facultad de Ciencias Agrarias de la Universidad de Antioquia. Estudiante de Maestría, FCAV/UNESP, São Paulo, Brasil.
⁴ Grupo de Investigación en Genética y Mejoramiento Animal de la Universidad de Antioquia y Grupo de Investigación en Producción, Desarrollo y Transformación Agropecuaria de la Corporación Universitaria Lasallista.

Recibido: 16/10/2007; fecha de aprobación: 05/02/2008

Resumen

Introducción. La curva de lactancia de los mamíferos puede ser explicada por funciones matemáticas que permiten relacionar la producción de leche de una hembra a través del tiempo, teniendo en cuenta la fisiología de la hembra a lo largo de la lactancia. Objetivo. Identificar las funciones matemáticas que mejor expliquen la curva de producción de leche en Búfalas del corregimiento de la sierra, municipio de Puerto Nare, departamento de Antioquia-Colombia. Materiales y métodos. Fueron utilizadas 3580 lactancias las cuales fueron agrupadas de acuerdo al número de partos y época de parto de las búfalas. Se utilizaron 9 modelos matemáticos de curvas de lactancia entre lineales y no lineales. Los criterios de selección para el mejor modelo matemático, fue el que presentó el mayor porcentaje de curvas significativas y la media del coeficiente de determinación. Resultados. Los resultados mostraron que el modelo que mejor ajustó la curva de lactancia en búfalas fue el propuesto por Papajcsik y Bordero, en 1988, el cual ajusto un 100 % las curvas en todos los números de partos con todas las épocas de parto. Conclusión. La metodología para la obtención del tiempo al pico de producción y la producción máxima puede ser una herramienta útil para trabajos futuros.

Palabras clave: Análisis longitudinales. Medidas repetidas. Modelos no lineales. Modelos polinomiales. Búfalas.

Abstract

Introduction.The lactation curve of mammals can be explained by mathematic functions which allow relating milk production of a female individual through time, taking into account her physiology throughout the lactation period. Objective. To identify the mathematic functions that can explain best the curve of milk production in a population of female buffalos from La Sierra, a village in the Puerto Nare town, Antioquia, Colombia. Materials and methods. 3850 lactations were used for the study, they were divided in groups according to their number and time of births the buffalos had. Nine mathematic models of lactation curves were used, both lineal and non-lineal. The selection criteria for choosing the best model were the highest percentage of significant curves and the media of the determination coefficient. Results. The model that best adjusted the lactation curve in female buffalos was that proposed by Papajcsik and Bordero, in 1988, which adjusted the curves in all the numbers of births with all the birth époques a 100%. Conclusion. The methodology to obtain the time at the peak of the production and the production at its maximum can be a useful tool for future works.

Key words: Longitudinal analysis. Repeated measures. Non- lineal models. Polinomial models. Female buffalos.

Resumo

Introdução. A curva de latência dos mamíferos pode ser explicada por funções matemáticas que permitem relacionar a produção de leite de uma fêmea através do tempo, tendo em conta a fisiologia da fêmea ao longo da latência. Objetivo. Identificar as funções matemáticas que melhor expliquem a curva de produção de leite em Bufa-as da cidade da serra, município de Porto Nare, departamento de Antioquia- Colômbia. Materiais e métodos. Foram utilizadas 3580 latências as quais foram agrupadas de acordo ao número de partos e época de parto das bufa-as. Utilizaram-se 9 modelos matemáticos de curvas de lactância entre lineares e não lineares. Os critérios de seleção para o melhor modelo matemático foi o que apresentou a maior percentagem de curvas significativas e a média do coeficiente de determinação. Resultados. Os resultados mostraram que o modelo que melhor ajustou a curva de latência em bufaas foi o proposto por Papajcsik e Bordero, em 1988, o qual ajusto um 100 % as curvas em todos os números de partos com todas as épocas de parto. Conclusão. A metodologia para a obtenção do tempo ao bico de produção e a produção máxima pode ser uma ferramenta útil para trabalhos futuros.

Palavras chaves: Análises longitudinais. Medidas repetidas. Modelos não lineares. Modelos polinomiais. Bufa-as.

Introducción

La curva de lactancia de los mamíferos puede ser explicada por funciones matemáticas que permiten relacionar la producción de leche de una hembra a través del tiempo, teniendo en cuenta la fisiología de la hembra a lo largo de la lactancia. En el caso de producción de leche de especies zootécnicas, el conocimiento de la curva de lactancia, permite predecir el desempeño futuro de los animales, la persistencia de la producción y tiempo en alcanzar el pico de producción, y ajustes por días en leche de hembras que no han terminado la lactancia.

Las curvas de lactancia permiten predecir el desempeño de las madres, información que puede utilizarse en la preselección de animales jóvenes destinados a las pruebas de progenie, también permiten entender mejor los sistemas de producción proporcionando información útil para el desarrollo de estrategias de selección, ayudando a determinar habilidades óptimas de producción y manejo con el fin de mejorar la eficiencia en la producción de leche^1,2.

Los modelos más utilizados para la estimación de curvas de lactancia son los modelos polinomiales y los modelos no lineales, entre los lineales se encuentran los modelos cuadráticos, cuadrático logarítmicos, polinomiales inversos, polinomios segmentados y el modelo lineal hiperbólico y entre los modelos no lineales están la función gamma incompleta, parabólica exponencial, modelo de Wilmink, regresión múltiple entre otros^1-5.

Fraga et al² (2003), realizaron un estudio de curvas de lactancia con 665 registros de producción de búfalas mestizas Murrah sometidas a un sistema de manejo con pastoreo de gramíneas y sin suplementación, Los modelos utilizados en el trabajo fueron: función gamma incompleta (Wood 1967)⁶, regresión múltiple (Ali y Shaeffer 1987)⁷, cuadrática logarítmica y lineal hiperbólica, concluyendo que el modelo cuadrático logarítmico fue el que mejor ajusto a la curva de lactancia en las búfalas.

Muñoz - Berrocal et al⁸ (2001), estudiaron el uso de polinomios ordinarios y segmentados en el ajuste de curvas de lactancia de búfalas Murrah y sus mestizos en Brasil. En este trabajo fueron utilizados 2868 lactancias correspondientes a 20959 controles lecheros. Los modelos utilizados fueron lineal-cuadrático, cúbico y polinomios segmentados. Los autores encontraron una mejor calidad de adherencia de la curva ajustada por el polinomio segmentado lineal-cuadrático.

Muñoz - Berrocal et al⁴ (2005), usaron modelos lineales y no lineales para el estudio de la curva de lactancia en búfalas Murrah en el estado de Sao Pablo Brasil, emplearon 3267 lactancias, referentes a 27600 controles de producción de leche. Los autores concluyeron que entre las funciones estudiadas, las funciones lineal cuadrática logarítmica y lineal hiperbólica fueron las que presentaron mejor ajuste de la curva media de lactancia.

Metry et al⁹ (1994), en estudios de curvas de lactancia en búfalas de Egipto concluyeron que los factores que afectan la forma de las curvas de lactancia de las búfalas son probablemente similares a aquellos del ganado vacuno. Un alto porcentaje de búfalas de primera lactancia en este estudio tienen lactancias cortas (37,3 % tiene < 150 días de lactancia).

Las diferencias estaciónales sugieren que un manejo mejorado por medio de una adecuada alimentación para la producción de leche, prevención de enfermedades, control de parásitos, protección contra climas extremos y un ciclo reproductivo regular podrían ser útiles para mejorar la producción de leche⁹.

El objetivo de este trabajo fue comparar diferentes modelos matemáticos y determinar cuáles son los que mejor se ajustan la producción de leche, tiempo y producción al pico en un sistema de producción orgánica en búfalas colombianas, de diferentes partos que paren en diferentes épocas del año.

Materiales y Métodos

El estudio fue realizado en una lechería orgánica de trópico bajo localizada en el corregimiento de la sierra municipio de Puerto Nare (Departamento de Antioquia-Colombia), en selva húmeda tropical a 125 MSNM, la bio-temperatura promedio para la zona es de 27°C, con una humedad relativa del 75%, la precipitación es 2.016 mm anuales, con una distribución bianual del periodo lluvioso: marzo-mayo, septiembre-noviembre y del periodo seco: Diciembre-Febrero y Junio-Agosto.

La caracterización racial es de tipo bufalipso o búfalo colombiano, sometida en un sistema extensivo con pasturas naturales (Panicum maximum, Dichantum aristatum, Brachiaria Spp, Echynocloa polystachya, Paspalum virgatum), suplementados con melaza y sal mineralizada a voluntad, los animales debían de caminar del potrero hasta la sala de ordeño entre 2,0 y 2,5 km aproximadamente. Se realizaron dos ordeños manuales al día (3:00 AM y 3:00 PM) con apoyo del becerro (cría) y al finalizar el ordeño, las búfalas junto con las crías van hacia una sala post-ordeño, donde consumen leche por espacio de 30 minutos, el pesaje de la leche se realizó con una balanza análoga de 10 kg.

Fueron empleados 102.302 controles diarios de leche, de 3.580 lactancias, de 1.829 búfalas, los datos fueron recolectados durante el periodo de 1994 a 2004.

Se utilizaron modelos matemáticos polinomiales con los procedimientos PROC REG y para los modelos no lineales PROC NLIN, del paquete estadístico SAS®¹⁰ (2002).

Los registros de producción fueron analizados separadamente por orden de parto (1...,5); época de parto (Seca 1= Dic-Feb; Seca 2=Jun-Ago; Lluvia 1=Mar-May; Lluvia 2=Sep- Nov). Para la elección del mejor modelo se tuvo en cuenta el porcentaje de curvas significativas (PM) y el valor del coeficiente de determinación (R²). Los modelos matemáticos analizados se muestran en la tabla 1.

Las medias de producción de leche en las diferentes épocas de parto se muestran en la Tabla 2. Fue observado que los mejores modelos para hembras de 1 a 5 partos y en las diferentes épocas de parto, fueron Papajcsik y Bordero¹⁴ (1988) con valores de PM de 100% de las curvas en todos los partos diferentes épocas y Brody¹² (1923), presentaron valores de PM entre 98,5% y 100% de las curvas en todos los partos y diferentes épocas. Estos dos modelos presentaron valores de R² desde 94,2 ± 0,7 hasta 96,1 ± 2,5. De forma general en los modelos lineales, presentaron estimativas de PM bajas, variando de 26,9 hasta 76,3 (Tablas 2, 3, 4, 5, 6 y 7).

La función gamma incompleta (Wood 1967)⁶ presentó en todos los partos y en todas las diferentes épocas (tablas 3, 4, 5, 6 y 7), porcentajes de curvas que no alcanzaron el nivel de convergencia, para las época seca 1, seca 2, lluvia 1 y Lluvia 2 variaron de 0,61% - 4,35%, de 0,50% - 2,97%, 5,81% y de 0,52% -1,87%, respectivamente.

El tiempo de producción máxima (pico de producción), se muestra en la Tabla 8, donde se observó el desempeño de las hembras según el número de parto y las diferentes épocas de parto. El promedio de días en alcanzar el pico fue 112,24, 102,98, 98,7, 100,8 y 101,12 para las hembras de primero, segundo, tercero, cuarto y quinto o más partos, respectivamente. Las búfalas de primer parto alcanzan el pico más tarde que las hembras multíparas, disminuyendo en el segundo y tercer parto, estabilizándose en el cuarto y más partos entre 100,8 y 101,12 días.

1. RAMÍREZ, VR. et al. Comparación de ecuaciones para estimar curvas de lactancia con diferentes estrategias de muestreo en bovinos angus, suizo y sus cruzas. En: Revista Veterinaria Mexicana. Vol 35 N° 3. (2004); p. 187-201.        [ Links ]

2. FRAGA, LM. et al. Estudio preliminar de las curvas de lactancia en las búfalas mestizas de Murrah. En: Revista Cubana Ciencias Agrícolas. Vol 37 (2003); p. 151-155.        [ Links ]

3. SHERCHAND, L. et al. Selection of mathematical model to generate lactation curves using daily milk yields of Holstein cows. En: Journal of Dairy Science. Vol 78 N° 11 (1995); p. 2507-2513.        [ Links ]

4. MUÑOZ-BERROCAL, M. et al. Uso de modelos lineales e nao lineares para o estudo da curva de lactacao em búfalos Murrah e sus mesticos en sistema de criacao extensivo, no estado de Sao Paulo. En: Archivo latinoamericano producción Animal. Vol 13 N°1 (2005); p. 7-12.        [ Links ]

5. FERNÁNDEZ, L. et al. Estimación de curvas de lactancia estándar de la raza siboney para su utilización en extensiones de lactancias. En: Revista Cubana Ciencias Agrícolas. 35 (2001); p. 99-104.        [ Links ]

6. WOOD, P. D. P. Algeraic model of the lactation curve in cattle. En: Nature (London). Vol. 216 (1967); p. 164 - 165.        [ Links ]

7. ALI, TE and SCHAEFFER, LR. Accounting for covariances among test day milk yields in dairy cows. En: Can. J. Anim. Sci. vol 67, (1987); p.637- 644.        [ Links ]

8. MUÑOZ-BERROCAL, M. et al. Uso de polinomios ordinarios y segmentados en el ajuste de curvas de lactancia de búfalas Murrah y sus mestizas en Brasil; The buffaloe an alternative for animal agricultural in the third Millenium Proceeding of the IV World buffalo congress Practical experiences. Vol. 2 (2001); 354 p.        [ Links ]

9. METRY, G.H. et al. Lactation curves for first lactation Egyptian buffalo. En: Journal of Dairy Science. 77 (1994); p. 1306-1314.        [ Links ]

10. SAS. System for regression SAS institute Inc, Cary, version 6 Edition Carv Nc North Caroline: SAS, 1995.        [ Links ]

11. WILMINK, J. B. M. Comporison of different methods of predicting 305 - day milk yield using means calculated from withinherd lactation curves. In: Livestock Production Science. Vol. 17 (1987); p.1 - 17.        [ Links ]

12. BRODY, S.; RAGSDALE, A.C. and TURNER, C.W. The rate of decline of milk secretion with the advance of period of lactation. In: Journal Genetic Physiology. Vol 5, (1923); p. 441-444.        [ Links ]

13. _. The relation between the initial rise and the subsequent decline of milk secretion following parturition. In: Journal Genetic Physiology. Vol 6, (1924); p.541-545.        [ Links ]

14. PAPAJCSIK, I.A. and BODERO, J. Modelling lactation curves of friesian cows in a subtropical climate. In: Animimal Production. Vol 47 N° 2 (1988); p. 201-207.        [ Links ]

15. COBUCI, J.A. et al. Curva de lactação na raça guzerá. En: Revista Brasilera de Zootecnia. Vol 29 N° 5 (2000); p. 1332-1339.        [ Links ]

16. MONTIEL, N.S; ANGULO, S.F. and PEROZO, F. Variaciones en los parámetros productivos en búfalas mestizas de primer y segundo parto en un bosque muy seco tropical de Venezuela. En: I Congreso internacional sobre mejoramiento animal CIMA La Habana Cuba.        [ Links ]