Estimation of lactation curves of Gyr cattle and some associated production parameters in the Colombian low tropic


  • Daniela Ferro Pedagogical and Technological University of Colombia
  • Jhon Gil Pedagogical and Technological University of Colombia
  • Ariel Jiménez Colombian Association of Zebu Cattle Breeders
  • Carlos Manrique National University of Colombia
  • Carlos A. Martínez Colombian Agricultural Research Corporation



Bos indicus, bovines, cows, lactation, lactation function shape, lactation parameters, linear models, milk production, milk yield, model comparison, nonlinear models, statistical modelling, tropical dairy operations


Background: The Gyr breed is widely used in Colombian low tropic dairy production systems. During the last 10 years, the Asociación Colombiana de Criadores de Ganado Cebú† - ASOCEBU, has been leading a dairy milk control program which led to the creation of a dataset that permits to carry out the first analysis of milk yield in Gyr cattle in the country using records from several herds. Objectives: To study milk production dynamics of Gyr cattle in the Colombian low tropic through the estimation of lactation curves and four derived production parameters: total milk yield between 5 and 305 days (TMY305), peak milk yield (PMY), days at peak (DP) and persistency (P). Methods: 13,798 daily milk yield records from 1,510 cows performing in 103 herds were used; the total number of lactations was 2,480. Four models were considered: Wood, Wiltmink, Papajcsik & Bordero, and a second-degree polynomial. Mean square error, mean absolute error, mean square error of prediction, Akaike and Bayesian information criteria were used to select the model better describing each lactation using the majority rule, that is, the model selected by most criteria was the chosen one. The shape of each fitted lactation curve was checked using basic results from calculus which permitted the classification of the estimated curves into two groups: typical and atypical; only typical functions were used to compute the four aforementioned production parameters. Results: The second-order polynomial was the model most frequently selected, while the Papajcsik & Bordero model had the lowest frequency. Average TMY305, PMY, DP and P were 3,489.86 kg, 17.28 kg, 57.17 days, and 0.83, respectively, with coefficients of variation: 0.27, 0.21, 0.41, and 0.16. Conclusions: This study permitted to identify individuals with outstanding phenotypic performance. To the best of our knowledge, this is the first study of this kind involving thousands of lactations from Gyr cows performing in several regions of Colombian low tropic.

= 673 veces | PDF
= 290 veces|


Download data is not yet available.


Metrics Loading ...

Author Biographies

Daniela Ferro, Pedagogical and Technological University of Colombia
School of Mathematics and Statistics, Pedagogical and Technological University of Colombia, Duitama Sectional Faculty, Colombia.

Jhon Gil, Pedagogical and Technological University of Colombia
School of Mathematics and Statistics, Pedagogical and Technological University of Colombia, Duitama Sectional Faculty, Colombia.

Ariel Jiménez, Colombian Association of Zebu Cattle Breeders
Colombian Association of Zebu Cattle Breeders (ASOCEBU), Bogotá, Colombia.

Carlos Manrique, National University of Colombia
Department of Animal Production, National University of Colombia, Bogotá, Colombia.

Carlos A. Martínez, Colombian Agricultural Research Corporation
Research and Development Directorate, Colombian Agricultural Research Corporation (AGROSAVIA), Headquarters, Mosquera, Colombia.


Akaike H. A New Look at the Statistical Model Identification. IEEE Transactions on Automatic Control 1974; 1 (6). DOI:

Analla M. Model Validation Through the Linear Regression Fit to Actual Versus Predicted Values. Agric Syst 1998; 57: 115–119. DOI:

ASOCEBU. La raza Gyr especializada en leche. [accessed date: May 2019] URL:

Brody S, Turner C. W, Ragsdale A. C. The Relation Between The Initial Rise And The Subsequent Decline Of Milk Secretion Following Parturition. J Gen Physiol 1924; 6: 541–545. DOI:

Butler Farms. Gyr Cattle. n.d. [accessed date: May 2018] URL:

Cañas A, Cerón-Muñoz M, Juan Corrales A. Modelación y parámetros genéticos de curvas de lactancia en bovinos Holstein en Colombia. Revista MVZ Córdoba 2012; 17(2): 2998–3003. DOI:

Cuellar M.A, Pereira R.J. Evaluación lechera en razas mestizas y gir lechero puro en condiciones sub-tropicales (período 1984 – 2001). Santa Cruz (Bolivia): Departamento de Santa Cruz Facultad de Ciencias Veterinarias, UAGRM; 2002. URL:

De-Souza R, Alcalde CR, Lopes-de-Oliveira CA, de-Labio-Molina BS, Fonseca-de-Macedo FDA, Gomes LC, Hygino B, Silva-Possamai AP. Lactation curves and economic results of Saanen goats fed increasing dietary energy levels obtained by the addition of calcium salts of fatty acids. Rev Bras Zootec 2014; 43(2): 73–79. DOI:

Faraway J. Linear Models with R. Boca Raton (FL): CRC Press, 2015.

Ferreira AG, Henrique S, Vieira RA, Maeda EM, Valotto AA. Fitting mathematical models to lactation curves from Holstein cows in the southwestern region of the state of Parana, Brazil. Acad Bras Cienc 2015; 1(87): 503–517. DOI:

Johansson I, Hansson A. Causes of Variation in Milk and Butterfat Yields of Dairy Cows. Kungliga Lantbruksakademiens Handlingar 1940; 79(62): 127. DOI:

López S, France J, Odongo NE, McBride RA, Kebreab E, AlZahal O, Dijkstra J. On the analysis of Canadian Holstein dairy cow lactation curves using standard growth functions. JDS 2015; 98(4): 2701–2712. DOI:

Nicolò P.P, Corrado D, Salvatore P.G, Roberto S, Giuseppe P. The Mathematical Description of Lactation Curves in Dairy Cattle. Italian Journal of Animal Science 2011; 10(51): 213–223. DOI:

Mellado M, Coronel F, Estrada A, Ríos F. Lactation Performance of Holstein and Holstein X Gyr Cattle Under Intensive Condition in a Subtropical Environment. Tropical and Subtropical Agroecosystems 2011; 14(3): 927–931. URL:

Ossa G, Torregroza L, Alvarado L. Determinación de la curva de lactancia en vacas mestizas de un hato de doble propósito en la Región Caribe de Colombia. Revista Corpoica 1997; 2(1): 54–57. DOI:

Wood P. D. P. Algebraic Model of the Lactation Curve in Cattle. Nature 1967; 216(5111), 164–165. DOI:

Pereira R.J. et al (). Milk yield persistency in Brazilian Gyr cattle based on a random regression model. Genetics and Molecular Research 2012; 11(2): 1599–1609. URL:

Papajcsik I, Bodero J. Modelling Lactation Curves Of Friesian Cows In A Subtropical Climate. Anim 1988; 47(2): 201–207. URL:

Quiroz, J, Granados L, and M. Barrón. Producción De Leche De Ganado Gyr Y F1 Y (Holstein X Gyr). AICA 2015; 6: 294–299. URL:

Rawlings J, Pantula S, Dickey. Applied regression analysis: a research tool. New York (NY): Springer-Verlag; 1998.

Rivera J. Producción bovina para leche. Pasto (Colombia): Universidad de Nariño; 1997.

Tedeschi L. Assessment of the Adequacy of Mathematical Models. Workshop on Mathematical Model Analysis and Evaluation. Sassari (Italy), 2004. URL:

Western Canadian Dairy Herd Improvement Services. DHI Lactation Curves. Western Diary Digest 2001, n.d. [accessed date: May 2019] URL:

Wiltmink J.B.M. Comparison of Different Methods of Predicting 305-Day Milk Yield Using Means Calculated from Within-Herd Lactation Curves. Livest Prod Sci 1987; 17: 1-17. DOI:

Yanar M, and Güler O. Factors Influencing the Shape of Lactation Curve and Persistency of Holstein Friesian Cows in High Altitude of Eastern Turkey. Journal of Applied Animal Research 2009; 35(1): 39-44. DOI:




How to Cite

Ferro, D., Gil, J., Jiménez, A., Manrique, C., & Martínez, C. A. (2021). Estimation of lactation curves of Gyr cattle and some associated production parameters in the Colombian low tropic. Revista Colombiana De Ciencias Pecuarias, 35(1), 3–13.



Original research articles

Similar Articles

> >> 

You may also start an advanced similarity search for this article.