Skip to content


  • Research
  • Open Access

Non-linear modelling to describe lactation curve in Gir crossbred cows

Journal of Animal Science and Technology201759:3

  • Received: 27 October 2016
  • Accepted: 3 February 2017
  • Published:



The modelling of lactation curve provides guidelines in formulating farm managerial practices in dairy cows. The aim of the present study was to determine the suitable non-linear model which most accurately fitted to lactation curves of five lactations in 134 Gir crossbred cows reared in Research-Cum-Development Project (RCDP) on Cattle farm, MPKV (Maharashtra). Four models viz. gamma-type function, quadratic model, mixed log function and Wilmink model were fitted to each lactation separately and then compared on the basis of goodness of fit measures viz. adjusted R2, root mean square error (RMSE), Akaike’s Informaion Criteria (AIC) and Bayesian Information Criteria (BIC).


In general, highest milk yield was observed in fourth lactation whereas it was lowest in first lactation. Among the models investigated, mixed log function and gamma-type function provided best fit of the lactation curve of first and remaining lactations, respectively. Quadratic model gave least fit to lactation curve in almost all lactations. Peak yield was observed as highest and lowest in fourth and first lactation, respectively. Further, first lactation showed highest persistency but relatively higher time to achieve peak yield than other lactations.


Lactation curve modelling using gamma-type function may be helpful to setting the management strategies at farm level, however, modelling must be optimized regularly before implementing them to enhance productivity in Gir crossbred cows.


  • Gir crossbred
  • Lactation curve
  • Non-linear modelling
  • Gamma-type function


The lactation curve in dairy cow provides graphical representation of milk production over the course of lactation. It shows that the milk production increases with high rate up to the point of peak yield and then after milk production started decreasing with lower rate [1, 2]. The shape of lactation curve could be useful in preparation of farm management strategy regarding health, breeding, nutrition and risk of potential diseases in dairy cows [36].

The lactation curve cannot be fitted using linear model since the curve don’t have linear trend with time. To explain the flow of milk production over the course of lactation in dairy cows, various mathematical models have been developed [79]. The lactation curve is influenced by several animal level and herd level factors such breed, age, parity, season of calving and managemental practices [10]. Thus, the parameters estimated from the different mathematical models may influences due to animal-level and environmental-level factors, and leads to variation in fitting to a typical lactation curve. A suitable model which can gave best fit to the lactation curve among the various models may only be efficient to predict the peak yield or phase of potential production, and to adopt the improved managemental practices.

Therefore, the aim of the present study was to determine the most appropriate non-linear model to describe lactation curves in Gir crossbred cows.



The present study was conducted at Research-Cum-Development Project on Cattle (RCDP), Mahatma Phule Krishi Vidyapeeth (MPKV), Rahuri, Maharashtra (India). This Project maintains Gir crossbred cows which have inheritance either of Jersey, Holstein or both. The data of half-bred and triple crossbred cows were considered as whole for analysis due to sample size in the study. The animals were kept in loose housing system. The separate buyers were provided for young calves, heifers, milch cows, dry cows, pregnant cows, breeding bulls etc. The feeding practices included lucerne, berseem, cowpea, maize, jowar and oat as green fodder and jowar kadabi as dry fodder. The concentrate mixture used for feeding different categories of animals was as per their nutritional requirement. Hand and machine milking was practiced at project and cows were milked twice a day. The production details regarding monthly milk yields of 134 of crossbred cows were collected for the period of April 2012 to March 2013. The monthly milk yields obtained from monthly milk registers for first 10 months of five lactations completed by Gir crossbred cows. The data were corrected for different date of calving by correction factor.

Statistical analysis

The descriptive statistics were obtained for monthly lactation milk yield in five lactations. One way ANOVA was used to determine if there was significant (p < 0.05) difference between milk yields of five lactation at each month of lactation. Pairwise comparisons between five lactation were done using by Tukey test.

Various lactation curve models, i.e., gamma-type function, quadratic model, mixed log function and Wilmink model was fitted to describe the lactation curve of Gir crossbred cows. These models fitted to average milk yields (kg/day) as follows:

Gamma-type function [7]: Y t  = at b e − ct

Quadratic model: Y t  = a + bt – ct 2

Mixed Log [8]: Y t  = a + bt 0.5 + c log t

Wilmink [9]: Y t  = a + be − kt  + ct

Where, Yt is the average milk yield on the tth month, a is the initial milk yield of lactation, b and c are the ascending slope parameter up to the peak yield and the descending slope parameter after peak yield, respectively. The constant k in Wilmink model was considered as 0.61 instead of 0.05 due to superiority of goodness of fit of model due to 0.61 in preliminary analysis.

The models were tested for goodness of fit (quality of prediction) using adjusted coefficient of determination (R adj 2 ), root mean square error (RMSE), Akaike’s Informaion Criteria (AIC) and Bayesian Information Criteria (BIC). The details of goodness fir criteria used in this study are as follows.
$$ \begin{array}{cc}\kern1em {R}_{adj}^2=1-\left[\frac{\left( n-1\right)}{\left( n- p\right)}\right]\left(1-{R}^2\right)\kern1em & \kern1em RMSE=\sqrt{\frac{ R SS}{n- p-1}}\kern1em \\ {}\kern1em AIC= n\ \ln (RSS)+2 p\kern1em & \kern1em BIC= n\ \ln \left(\frac{ R SS}{n}\right)+ p\ \ln ( n)\kern1em \end{array} $$

Where R2 = 1- (RSS/TSS), RSS and TSS represents residual and total sum of square, respectively; n and p are the number of observations (data points) and parameters in the model, respectively. ln represents natural logarithm.

The value of R adj 2 close to 1 indicates satisfactory fitting due to model. A smaller value of RMSE, AIC and BIC indicates a better fit of models. Therefore, the model which has highest adjusted R 2 value and lowest value of RMSE, AIC and BIC value was considered most appropriate for describing lactation curve of Gir crossbred cows. Statistical analyses was done was using SAS 9.3 Version [11]. SAS software uses Gauss-Newton method for iteration and provides least squares estimates. Further, peak yield, persistency and months in milk at peak yield were calculated using gamma-type function as described by Tekerli et al. [10] as Peak yield = a * (b/c)e‐ b, Persistency = ‐ (b + 1) * ln(c) and Months in milk at peak yield = b/c.

Results and discussion

The average milk yield across lactation in first five lactations in Gir crossbred cows is presented in Table 1. The highest milk yield was observed in second month of lactation in first (10.08 kg), second (13.19), third (14.92 kg), fourth (16.59 kg) and fifth (14.41 kg) lactation. One way analysis of variance revealed that there was significant (p < 0.05) difference between milk yields of five lactations up to four months only. Overall there was significantly higher milk yield in fourth lactation followed by third, fifth and second lactation than first lactation of Gir crossbred cows. Maximum yield in fourth lactation than other lactations was in accordance with reports of Jingar et al. [12] in Karan Fries (crossbred) cows, but contradicted with reports of Rekik et al. [13] who reported maximum milk yield in third lactation. Because of this significant difference in five lactations, the non-linear modelling of lactation curve was undertaken lactation-wise separately to describe model fitting precisely.
Table 1

Monthly average milk yield (kg/day) in first five lactations of Gir crossbred cows








7.41 ± 0.51a (49)

11.46 ± 0.71b (26)

11.94 ± 0.69b (24)

13.20 ± 1.07b (18)

10.50 ± 1.15ab (17)


10.08 ± 0.50a (48)

13.19 ± 0.82ab (26)

14.92 ± 0.76bc (24)

16.59 ± 1.03c (18)

14.41 ± 1.15bc (17)


9.81 ± 0.49a (48)

10.95 ± 0.83ab (26)

13.29 ± 0.83bc (24)

14.82 ± 0.85c (18)

13.60 ± 1.15bc (17)


9.02 ± 0.43a (45)

9.82 ± 0.85ab (26)

11.87 ± 0.84ab (24)

12.80 ± 0.92b (17)

12.28 ± 0.97b (17)


8.34 ± 0.45 (44)

8.74 ± 0.89 (25)

10.62 ± 0.68 (23)

10.95 ± 0.79 (17)

10.16 ± 0.94 (17)


7.46 ± 0.38 (40)

7.38 ± 0.80 (24)

9.02 ± 0.68 (23)

9.74 ± 0.84 (17)

8.61 ± 1.03 (17)


7.19 ± 0.39 (37)

6.94 ± 0.92 (19)

8.11 ± 0.72 (21)

7.86 ± 0.88 (17)

7.74 ± 1.03 (16)


6.63 ± 0.34 (34)

6.41 ± 0.94 (15)

7.47 ± 0.68 (19)

7.21 ± 0.99 (12)

6.53 ± 1.08 (14)


5.64 ± 0.37 (30)

6.29 ± 0.92 (13)

6.44 ± 0.63 (18)

6.96 ± 0.61 (10)

6.94 ± 1.20 (10)


5.82 ± 0.35 (23)

6.01 ± 0.76 (10)

5.55 ± 0.65 (14)

6.49 ± 0.31 (9)

6.17 ± 1.02 (9)

Figure in parenthesis indicates number of observations. Different superscript (a, b, c) differ significantly (p < 0.05) in same row

The parameters estimates (along with standard error) due to various lactation curve models fitted average milk yield for different lactations are presented in Table 2. The graphical presentation of lactation curve for 5 lactations due to observed and predicted values is shown in Figs. 1, 2, 3, 4, and 5. In general, the estimates of initial production (parameter A) varied between five lactations and it was lowest in first lactation. This finding was similar to reports of Madalena et al. [1] in Holtein-Frisian and Holtein-Friesian × Gir cows. The parameter b and c were also found to be wide-ranging for different lactations. These findings were in accordance with reports of Boujenane [14] in Moroccan Holstein‐Friesian dairy cows. The estimates parameters due to gamma-type functions fitted to different lactations were similar to findings reported by Jingar et al. [12] in Karan Fries (crossbred) cows and Rekik et al. [13] in Holstein–Friesian cows.
Table 2

Lactation-wise non-linear modelling to average milk yield (kg/day) in Gir crossbred cows



Parameters estimates (Standard error)

Goodness of fit




R adj 2






9.38 (0.45)

0.47 (0.10)

0.17 (0.03)







17.53 (1.22)

−9.69 (1.44)

7.94 (1.39)







12.93 (0.88)

−7.67 (2.03)

−0.78 (0.12)







8.54 (1.03)

0.32 (0.43)

0.07 (0.04)







13.83 (0.72)

0.14 (0.12)

0.13 (0.03)







18.49 (2.05)

−6.25 (2.41)

2.82 (2.33)







12.83 (1.34)

0.60 (3.10)

−0.76 (0.18)







14.03 (0.99)

−1.27 (0.41)

−0.04 (0.04)







15.45 (0.59)

0.43 (0.09)

0.21 (0.02)







26.06 (1.82)

−13.42 (2.14)

9.28 (2.07)







17.76 (1.31)

−6.64 (3.02)

−1.28 (0.18)







14.25 (1.3)

−0.48 (0.54)

0.04 (0.05)







17.42 (0.94)

0.42 (0.12)

0.22 (0.04)







27.97 (2.85)

−13.74 (3.36)

8.97 (3.24)







18.8 (1.94)

−5.42 (4.49)

−1.36 (0.26)







16.60 (1.72)

−1.04 (0.72)

0.01 (0.06)







14.33 (0.95)

0.55 (0.15)

0.23 (0.04)







26.06 (2.81)

−14.62 (3.31)

10.92 (3.20)







17.86 (1.91)

−9.00 (4.41)

−1.29 (0.26)







13.37 (1.88)

−0.34 (0.79)

0.05 (0.07)





LO Lactation order, GT Gamma-type function, QD Quadratic model, ML Mixed log function, WL Wilmink model, R adj 2 Adjusted coefficient of determination, RMSE Root mean square error, AIC Akaike’s information criteria, BIC Bayesian Information Criteria

Fig. 1
Fig. 1

Predicted milk yield due to non-linear modeling in 1st lactation of Gir Crossbred

Fig. 2
Fig. 2

Predicted milk yield due to non-linear modeling in 2nd lactation of Gir Crossbred

Fig. 3
Fig. 3

Predicted milk yield due to non-linear modeling in 3rd lactation of Gir Crossbred

Fig. 4
Fig. 4

Predicted milk yield due to non-linear modeling in 4th lactation of Gir Crossbred

Fig. 5
Fig. 5

Predicted milk yield due to non-linear modeling in 5th lactation of Gir Crossbred

In primiparous cows (first lactations), high value of adjusted R2 (0.679 to 0.893) indicated that non-linear modelling explained sufficient variability in shape of lactation curve, which was also reported by Boujenane [14]. However, lower estimates of adjusted coefficient of determination reported by Olori et al. [15]. The adjusted R2 was highest for mixed log function (0.893), followed by gamma-type function (0.891), wilmink (0.850) and least for Quadratic model (0.679). Further, The RMSE values ranged from 0.548 to 0.949 and mixed log function provided lowest value of RMSE than other models. Similarly, AIC and BIC values were also least for mixed model (11.889 and −10.229) followed by gamma-type function (12.027 and −10.091), wilmink (15.310 and −6.808) and quadratic model (22.862 and 0.744). Therefore, mixed log function was considered best model for fitting to lactation curve of primiparous Gir crossbred cows. The best fit due to mixed log function was also reported by Dongare et al. [16] while comparing gamma-type function, mixed log function and quadratic model in Sahiwal cows. The closeness of fitting between gamma-type function and mixed log function was observed, as also reported by Dongare et al. [16].

In second lactation, range of adjusted R2 (0.866 to 0.905) indicated better non-linear modelling in explaining the variability in lactation curve than first lactation. Gamma-type function had explained higher variation (adjusted R 2  = 0.905) and fitted better than other models. Further, RMSE, AIC and BIC values due to gamma-type function were lowest as 0.841, 20.462, −1.656 than that of other models. The superiority of gamma-type function for fitting of lactation curve among various models was also reported Cankaya et al. [17] in Jersey cattle. Whereas, least fitting was observed due to wilmink model with highest values of RMSE (0.994), AIC (23.797) and BIC (1.679). The trend of goodness of fit criteria for lactation curve of second lactation was in order (best first): gamma-type function < quadratic model < mixed log function < wilmink model.

Similarly, in case of third and higher lactations, highest adjusted R2 was observed due to gamma-type function. RMSE, AIC and BIC values due to gamma-type function were lowest for third, fourth and fifth lactations as compared to other models. The best fit due to gamma-type function model was reported by Boujenane [14] in Moroccan Holstein-Friesian dairy cows and Jingar et al. [12] in Karan Fries (crossbred) cows. The superiority in variability explained by gamma-type function in multiparous (second or more lactations) cows was in agreement with findings of previous studies [18, 19] but contrast with reports of Koçak and Ekiz [20] in Holstein cows and Dohare et al. [21] in Frieswal cows (62% Friesian and 38% Sahiwal inheritance). However, lowest adjusted R 2 value and highest values of RMSE, AIC and BIC were observed for quadratic model fitting, which in accordance with reports of Cilek and Keskin [22] who fitted gamma-type function, mixed log, quadratic model, cubic and exponential and polynomial regression model to lactation curve of Simmental cows. The trend of goodness of fit for third, fourth and fifth lactation was in order with best due to gamma-type function followed by mixed log function, Wilmink model and least with quadratic model.

Peak yield, Persistency, and Months in milk at peak yield due to gamma-type function are presented in Table 3. The highest peak yield was observed in fourth lactation (15.02 kg) followed by third (13.68 kg), fifth (13.35 kg), second (12.15 kg) and least in first lactation (9.45 kg). Months in milk at peak were lowest for second lactation (1.08) and highest for first lactation (2.76). However, persistency was found to be high in first lactation (2.60 months) than other lactations, which was in accordance with reports of Rekik et al. [13].
Table 3

Peak yield, persistency and months in milk at peak in five lactations in Gir crossbred cows


Peak yield (kg)

Persistency (months)

Months in milk at peak yield





















The understanding of the lactation curve of Gir crossbred cows may be an efficient tool for adopting the feeding and management practices. The gamma-type function may be used for as leading model for achieving desire productivity in Gir crossbred cows. However, the accuracy in prediction due to nonlinear models varies with variability in lactation yield in herd structure over time. Therefore, it was suggested that the optimization of lactation curve models at regular interval is necessary before their implementation.


As there was significant difference in milk yield in different lactations, non-linear modelling showed varied fitting of lactation curve in first lactation and other lactations. Among the four models studied, mixed log function provided best fit of the lactation curve of primiparous cows, due to lower values of RMSE, AIC and BIC. However, in multiparous cows, gamma-type function described most appropriately the lactation curve as compared to other model. Quadratic model gave least fit to lactation curve in almost all lactations. Peak yield was highest in fourth lactation and least in first lactation. The persistency was higher in first lactation of Gir crossbred. It was suggested that lactation curve models may be helpful to setting the management strategies at farm level, however, modelling must be optimized regularly before implementing them to enhance productivity in Gir crossbred cows.



Akaike’s informaion criteria


Bayesian information criteria


p value


Root mean square error



The authors thank to Late Dr. B. Singh, Ex-Head, Division of LES & IT, Indian Veterinary Research Institute, Izatnagar (India) for valuable guidance during the study programme. The authors are highly thankful to Director of Research, and Head of Department, Division of Animal Husbandry and Dairy Science, M.P.K.V., Rahuri, (Maharashtra) for permission to use milk production data from RCDP on Cattle farm. The authors also express their gratitude to the learned referee and the Editor-in-Chief for their valuable comments on the original version of the paper.


This study was part of PhD research work which was facilitated by Indian Veterinary Research Institute, Izatnagar, Uttar Pradesh (India).

Availability of data and materials

Not applicable.

Authors’ contributions

YCB has collected data, performed analysis and wrote manuscript. MRV guided during the study and also corrected manuscript. Both authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval

Not required.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Network Project on Sheep (Deccani) Improvement, Mahatma Phule Krishi Vidyapeeth, Rahuri, Maharashtra, 413722, India
Division of Livestock Economics, Statistics and IT, Indian Veterinary Research Institute, Bareilly, Uttar Pradesh, 243122, India


  1. Madalena FE, Martinez ML, Freitas AF. Lactation curves of Holtein-Frisian and Holtein-Friesian × Gir cows. Anim Prod. 1979;29:101–7.View ArticleGoogle Scholar
  2. Val-Arreola D, Kebreab E, Dijkstra J, France J. Study of the lactation curve in dairy cattle on farms in central Mexico. J Dairy Sci. 2004;87:3789–99.View ArticlePubMedGoogle Scholar
  3. Pérochon L, Coulon JB, Lescourret F. Modelling lactation curves of dairy cows with emphasis on individual variability. Anim Sci. 1996;63:189–200.View ArticleGoogle Scholar
  4. Macciotta N, Vicario D, Cappio-Borlino A. Detection of different shapes of lactation curve for milk yield in dairy cattle by empirical mathematical models. J Dairy Sci. 2005;88:1178–91.View ArticlePubMedGoogle Scholar
  5. Macciotta N, Dimauro C, Rassu S, Steri R, Pulina G. The mathematical description of lactation curves in dairy cattle. Ital J Anim Sci. 2011;10(4):e51.View ArticleGoogle Scholar
  6. Sahoo SK, Singh A, Gupta AK, Chakravarty AK, Ambhore GS, Dash SK. Comparative evaluation of different lactation curve functions for prediction of bi-monthly test day yields in murrah buffaloes. Anim Sci Reporter. 2015;9(3):89–94.Google Scholar
  7. Wood PDP. Algebraic model of lactation curve in cattle. Nature. 1967;216:164–5.View ArticleGoogle Scholar
  8. Guo Z, Swalve HH. Modelling of the lactation curve as a sub-model in the evaluation of test day records. Proc. Interbull Mtg., Prague, INTERBULL Bull., No. 11; 1995.Google Scholar
  9. Wilmink JBM. Adjustment of test-day milk, fat and protein yield for age, season and stage of lactation. Livest Prod Sci. 1987;16:335–48.View ArticleGoogle Scholar
  10. Tekerli M, Akinci Z, Dogan J, Akcan A. Factors affecting the shape of lactation curves of Holstein cows from the Balikesir Province of Turkey. J Dairy Sci. 2000;83:1381–6.View ArticlePubMedGoogle Scholar
  11. SAS 9.3 software. SAS Institute Inc: SAS/STAT® 131 User’s Guide, Cary, NC. 2013.Google Scholar
  12. Jingar S, Mehla RK, Singh M, Roy AK. Lactation curve pattern and prediction of milk production performance in crossbred cows. J Vet Med. 2014. doi:10.1155/2014/814768.PubMedPubMed CentralGoogle Scholar
  13. Rekik B, Ben Gara A, Ben Hamouda M, Hammami H. Fitting lactation curves of dairy cattle in different types of herds in Tunisia. Livest Prod Sci. 2003;83(2–3):309–15.View ArticleGoogle Scholar
  14. Boujenane I. Comparison of different lactation curve models to describe lactation curve in Moroccan Holstein-Friesian dairy cows. Iran J Appl Anim Sci. 2013;3(4):817–22.Google Scholar
  15. Olori VE, Brotherstone S, Hill WG, Mc Guirk BJ. Fit of standard models of the lactation curve to weekly records of milk production of cows in a single herd. Livest Prod Sci. 1999;58:55–63.View ArticleGoogle Scholar
  16. Dongre V, Gandhi RS, Singh A. Comparison of different lactation curve models in Sahiwal cows. Turk J Vet Anim Sci. 2012;36:723–6.Google Scholar
  17. Cankaya S, Unalan A, Soydan E. Selection of a mathematical model to describe the lactation curves of Jersey cattle. Arch Tierz. 2011;54(1):27–35.Google Scholar
  18. Gantner V, Jovanovac S, Raguz N, Solic D, Kuterovac K. Nonlinear vs. Linear regression models in lactation curve prediction. Bulg J Agric Sci. 2010;16(6):794–800.Google Scholar
  19. Ferreira AGT, Henrique DS, Vieira RAM, Maeda EM, Valotto AA. Fitting mathematical models to lactation curves from holstein cows in the southwestern region of the state of Parana, Brazil. An Acad Bras Cienc. 2015;87(1):503–17.View ArticlePubMedGoogle Scholar
  20. Koçak O, Ekiz B. Comparison of different lactation curve models in Holstein cows raised on a farm in the south-eastern Anatolia region. Arch Tierz Dummerstorf. 2008;4:329–37.Google Scholar
  21. Dohare AK, Singh B, Verma MR, Perme B, Sharma VB, Gupta N, Kshandakar S. Comparison of standard lactation curve models using fortnightly milk records in Frieswal cattle. Vet World. 2014;7(10):831–4.View ArticleGoogle Scholar
  22. Cilek S, Keskin I. Comparison of six different mathematical models to the lactation curve of Simmental cows reared in Kazova state farm. J Anim Vet Advan. 2008;7:1316–9.Google Scholar


© The Author(s). 2017