A comparison of 4 different machine learning algorithms to predict lactoferrin content in bovine milk from mid-infrared spectra.

Lactoferrin (LF) is a glycoprotein naturally present in milk. Its content varies throughout lactation, but also with mastitis; therefore it is a potential additional indicator of udder health beyond somatic cell count. Condequently, there is an interest in quantifying this biomolecule routinely. First prediction equations proposed in the literature to predict the content in milk using milk mid-infrared spectrometry were built using partial least square regression (PLSR) due to the limited size of the data set. Thanks to a large data set, the current study aimed to test 4 different machine learning algorithms using a large data set comprising 6,619 records collected across different herds, breeds, and countries. The first algorithm was a PLSR, as used in past investigations. The second and third algorithms used partial least square (PLS) factors combined with a linear and polynomial support vector regression (PLS + SVR). The fourth algorithm also used PLS factors, but included in an artificial neural network with 1 hidden layer (PLS + ANN). The training and validation sets comprised 5,541 and 836 records, respectively. Even if the calibration prediction performances were the best for PLS + polynomial SVR, their validation prediction performances were the worst. The 3 other algorithms had similar validation performances. Indeed, the validation root mean squared error (RMSE) ranged between 162.17 and 166.75 mg/L of milk. However, the lower standard deviation of cross-validation RMSE and the better normality of the residual distribution observed for PLS + ANN suggest that this modeling was more suitable to predict the LF content in milk from milk mid-infrared spectra (R2v = 0.60 and validation RMSE = 162.17 mg/L of milk). This PLS +ANN model was then applied to almost 6 million spectral records. The predicted LF showed the expected relationships with milk yield, somatic cell score, somatic cell count, and stage of lactation. The model tended to underestimate high LF values (higher than 600 mg/L of milk). However, if the prediction threshold was set to 500 mg/L, 82% of samples from the validation having a content of LF higher than 600 mg/L were detected. Future research should aim to increase the number of those extremely high LF records in the calibration set.

Predicting milk mid-infrared spectra from first-parity Holstein cows using a test-day mixed model with the perspective of herd management.

The use of test-day models to model milk mid-infrared (MIR) spectra for genetic purposes has already been explored; however, little attention has been given to their use to predict milk MIR spectra for management purposes. The aim of this paper was to study the ability of a test-day mixed model to predict milk MIR spectra for management purposes. A data set containing 467,496 test-day observations from 53,781 Holstein dairy cows in first lactation was used for model building. Principal component analysis was implemented on the selected 311 MIR spectral wavenumbers to reduce the number of traits for modeling; 12 principal components (PC) were retained, explaining approximately 96% of the total spectral variation. Each of the retained PC was modeled using a single trait test-day mixed model. The model solutions were used to compute the predicted scores of each PC, followed by a back-transformation to obtain the 311 predicted MIR spectral wavenumbers. Four new data sets, containing altogether 122,032 records, were used to test the ability of the model to predict milk MIR spectra in 4 distinct scenarios with different levels of information about the cows. The average correlation between observed and predicted values of each spectral wavenumber was 0.85 for the modeling data set and ranged from 0.36 to 0.62 for the scenarios. Correlations between milk fat, protein, and lactose contents predicted from the observed spectra and from the modeled spectra ranged from 0.83 to 0.89 for the modeling set and from 0.32 to 0.73 for the scenarios. Our results demonstrated a moderate but promising ability to predict milk MIR spectra using a test-day mixed model. Current and future MIR traits prediction equations could be applied on the modeled spectra to predict all MIR traits in different situations instead of developing one test-day model separately for each trait. Modeling MIR spectra would benefit farmers for cow and herd management, for instance through prediction of future records or comparison between observed and expected wavenumbers or MIR traits for the detection of health and management problems. Potential resulting tools could be incorporated into milk recording systems.

Diagnosing the pregnancy status of dairy cows: How useful is milk mid-infrared spectroscopy?

Pregnancy diagnosis is an essential part of successful breeding programs on dairy farms. Milk composition alters with pregnancy, and this is well documented. Fourier-transform mid-infrared (MIR) spectroscopy is a rapid and cost-effective method for providing milk spectra that reflect the detailed composition of milk samples. Therefore, the aim of this study was to assess the ability of MIR spectroscopy to predict the pregnancy status of dairy cows. The MIR spectra and insemination records were available from 8,064 Holstein cows of 19 commercial dairy farms in Australia. Three strategies were studied to classify cows as open or pregnant using partial least squares discriminant analysis models with random cow-independent 10-fold cross-validation and external validation on a cow-independent test set. The first strategy considered 6,754 MIR spectra after insemination used as independent variables in the model. The results showed little ability to detect the pregnancy status as the area under the receiver operating characteristic curve was 0.63 and 0.65 for cross-validation and testing, respectively. The second strategy, involving 1,664 records, aimed to reduce noise in the MIR spectra used as predictors by subtracting a spectrum before insemination (i.e., open spectrum) from the spectrum after insemination. The accuracy was comparable with the first approach, showing no superiority of the method. Given the limited results for these models when using combined data from all stages after insemination, the third strategy explored separate models at 7 stages after insemination comprising 348 to 1,566 records each (i.e., progressively greater gestation) with single MIR spectra after insemination as predictors. The models developed using data recorded after 150 d of pregnancy showed promising prediction accuracy with the average value of area under the receiver operating characteristic curve of 0.78 and 0.76 obtained through cross-validation and testing, respectively. If this can be confirmed on a larger data set and extended to somewhat earlier stages after insemination, the model could be used as a complementary tool to detect fetal abortion.

Relationships between milk mid-IR predicted gastro-enteric methane production and the technical and financial performance of commercial dairy herds.

Considering economic and environmental issues is important in ensuring the sustainability of dairy farms. The objective of this study was to investigate univariate relationships between lactating dairy cow gastro-enteric methane (CH4) production predicted from milk mid-IR (MIR) spectra and technico-economic variables by the use of large scale and on-farm data. A total of 525 697 individual CH4 predictions from milk MIR spectra (MIR-CH4 (g/day)) of milk samples collected on 206 farms during the Walloon milk recording scheme were used to create a MIR-CH4 prediction for each herd and year (HYMIR-CH4). These predictions were merged with dairy herd accounting data. This allowed a simultaneous study of HYMIR-CH4 and 42 technical and economic variables for 1024 herd and year records from 2007 to 2014. Pearson correlation coefficients (r) were used to assess significant relationships (P<0.05). Low HYMIR-CH4 was significantly associated with, amongst others, lower fat and protein corrected milk (FPCM) yield (r=0.18), lower milk fat and protein content (r=0.38 and 0.33, respectively), lower quantity of milk produced from forages (r=0.12) and suboptimal reproduction and health performance (e.g. longer calving interval (r=-0.21) and higher culling rate (r=-0.15)). Concerning economic results, low HYMIR-CH4 was significantly associated with lower gross margin per cow (r=0.19) and per litre FPCM (r=0.09). To conclude, this study suggested that low lactating dairy cow gastro-enteric CH4 production tended to be associated with more extensive or suboptimal management practices, which could lead to lower profitability. The observed low correlations suggest complex interactions between variables due to the use of on-farm data with large variability in technical and management practices.

The feeding system impacts relationships between calving interval and economic results of dairy farms.

The calving interval (CI) can potentially impact the economic results of dairy farms. This study highlighted the most profitable CI and innovated by describing this optimum as a function of the feeding system of the farm. On-farm data were used to represent real farm conditions. A total of 1832 accounts of farms recorded from 2007 to 2014 provided economic, technical and feeding information per herd and per year. A multiple correspondence analysis created four feeding groups: extensive, low intensive, intensive and very intensive herds. The gross margin and some of its components were corrected to account for the effect of factors external to the farm, such as the market, biological status, etc. Then the corrected gross margin (cGMc) and its components were modelled by CI parameters in each feeding system by use of GLM. The relationship between cGMc and the proportion of cows with CI<380 days in each feeding group showed that keeping most of the cows in the herd with CI near to 1 year was not profitable for most farms (for the very intensive farms there was no effect of the proportion). Moreover, a low proportion of cows (0% to 20%) with a near-to-1-year CI was not profitable for the extensive and low intensive farms. Extending the proportion of cows with CI beyond 459 days until 635 days (i.e. data limitation) caused no significant economic loss for the extensive and low intensive farms, but was not profitable for the intensive and very intensive farms. Variations of the milk and feeding components explained mainly these significant differences of gross margin. A link between the feeding system and persistency, perceptible in the milk production and CI shown by the herd, could explain the different relationships observed between the extent of CI and the economic results in the feeding groups. This herd-level study tended to show different economic optima of CI as a function of the feeding system. A cow-level study would specify these tendencies to give CI objectives to dairy breeders as a function of their farm characteristics.