Use of a computer-aided design to develop a stress simulation model for lactating dairy sheep

Introducing a sinusoidal equation to describe lactation curves for cumulative milk yield and composition in Holstein cows

The aim of the work reported here was to investigate the appropriateness of a sinusoidal function by applying it to model the cumulative lactation curves for milk yield and composition in primiparous Holstein cows, and to compare it with three conventional growth models (linear, Richards and Morgan). Data used in this study were 911 144 test-day records for milk, fat and protein yields, which were recorded on 834 dairy herds from 2000 to 2011 by the Animal Breeding Centre and Promotion of Animal Products of Iran. Each function was fitted to the test-day production records using appropriate procedures in SAS (PROC REG for the linear model and PROC NLIN for the Richards, Morgan and sinusoidal equations) and the parameters were estimated. The models were tested for goodness of fit using adjusted coefficient of determination $\lpar {R_{{\rm adj}}^2 } \rpar $, root mean square error (RMSE), Akaike's information criterion (AIC) and the Bayesian information criterion (BIC). $R_{{\rm adj}}^2 $ values were generally high (>0.999), implying suitable fits to the data, and showed little differences among the models for cumulative yields. The sinusoidal equation provided the lowest values of RMSE, AIC and BIC, and therefore the best fit to the lactation curve for cumulative milk, fat and protein yields. The linear model gave the poorest fit to the cumulative lactation curve for all production traits. The current results show that classical growth functions can be fitted accurately to cumulative lactation curves for production traits, but the new sinusoidal equation introduced herein, by providing best goodness of fit, can be considered a useful alternative to conventional models in dairy research.

Modelling lactation curve for milk fat to protein ratio in Iranian buffaloes (Bubalus bubalis) using non-linear mixed models

The aim of this study was to compare seven non-linear mathematical models (Brody, Wood, Dhanoa, Sikka, Nelder, Rook and Dijkstra) to examine their efficiency in describing the lactation curves for milk fat to protein ratio (FPR) in Iranian buffaloes. Data were 43 818 test-day records for FPR from the first three lactations of Iranian buffaloes which were collected on 523 dairy herds in the period from 1996 to 2012 by the Animal Breeding Center of Iran. Each model was fitted to monthly FPR records of buffaloes using the non-linear mixed model procedure (PROC NLMIXED) in SAS and the parameters were estimated. The models were tested for goodness of fit using Akaike's information criterion (AIC), Bayesian information criterion (BIC) and log maximum likelihood (-2 Log L). The Nelder and Sikka mixed models provided the best fit of lactation curve for FPR in the first and second lactations of Iranian buffaloes, respectively. However, Wood, Dhanoa and Sikka mixed models provided the best fit of lactation curve for FPR in the third parity buffaloes. Evaluation of first, second and third lactation features showed that all models, except for Dijkstra model in the third lactation, under-predicted test time at which daily FPR was minimum. On the other hand, minimum FPR was over-predicted by all equations. Evaluation of the different models used in this study indicated that non-linear mixed models were sufficient for fitting test-day FPR records of Iranian buffaloes.