Resilient cows are those that are minimally affected by environmental disturbances, such as pathogens or extreme weather, and that quickly recover if they are affected by these disturbances (Colditz and Hine, Reference Colditz and Hine2016; Berghof et al., Reference Berghof, Poppe and Mulder2019). Farms with good herd management could support less resilient cows, whereas limitations in herd management could potentially be compensated by the resilience of cows. Herd resilience is not simply the average of the individual cow resilience values in a herd, since it depends on the adaptive capacity of the animals together with management decisions that affect the performance of these animals and their environment (Blanc et al., Reference Blanc, Ollion, Puillet, Delaby, Ingrand, Tichit and Friggens2013). In the context of dairy herds, this means that resilient herds show less milk yield deviations at herd level and thus the herd as a whole is assumed to be less affected by disturbances. There are differences in herd resilience and these differences could partly be explained by herd management: Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021) have shown that resilient herds tended to have a lower somatic cell score (SCS), a lower proportion of cows with elevated somatic cell count (SCC), a higher survival to second lactation, a shorter calving interval, a lower proportion of cows with either a ketosis indication or a rumen acidosis indication and, finally, a lower age at first calving when compared to non-resilient herds.
The general idea of proxies for individual cow resilience (Elgersma et al., Reference Elgersma, de Jong, van der Linde and Mulder2018; Poppe et al., Reference Poppe, Veerkamp, van Pelt and Mulder2020) and dairy herd resilience (Poppe et al., Reference Poppe, Mulder, Kamphuis and Veerkamp2021) is based on fluctuations of the actual measured daily milk yield around the expected daily milk yield curves. These fluctuations can be used to compute variances or auto-correlations that express different aspects of resilience (Poppe et al., Reference Poppe, Veerkamp, van Pelt and Mulder2020). Daily milk yield recordings are, therefore, necessary to be able to compute proxies for resilience, and thus, resilience can only be estimated on farms that use automated milking systems (AMS), sometimes known as milking robots. In contrast to AMS, many milking parlours and other conventional milking systems (CMS) do not record daily milk yield. In the Netherlands the percentage of farms with an AMS has increased from approximately 10% in 2010 to 31% (4,574 farms) in 2022 (Stichting KOM, Reference Stichting2022). Worldwide, the estimated number of farms with an AMS increased from 8,000 in 2009 (de Koning, Reference de Koning2011) to approximately 25 000 farms in 2015 (Barkema et al., Reference Barkema, von Keyserlingk, Kastelic, Lam, Luby, Roy, LeBlanc, Keefe and Kelton2015). Even though the number of farms with an AMS is increasing, a large portion of farms still use CMS. These farms have less daily milk yield information and individual or herd resilience estimation becomes problematic to estimate. Furthermore, farms that install an AMS need to wait one full lactation before it is possible to estimate individual or herd resilience of their cows, since this estimation is based on one full lactation (Poppe et al., Reference Poppe, Veerkamp, van Pelt and Mulder2020). Even though individual cow resilience cannot be estimated for these CMS farms, it could be helpful to know their herd resilience. Farms with low herd resilience could try to improve the herd resilience via management decisions that affect the performance of their animals and their environment.
The aim of this study is to assess whether herd resilience can be predicted using herd performance data commonly available on CMS farms. To do so, we used data from AMS farms and used the methodology of Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021) to estimate herd resilience, which served as a reference value. We, subsequently, applied a Random Forest model to predict herd resilience using only variables that are also available on CMS farms.
Material and methods
To address the aim of this study the following information was needed: (1) a herd resilience indicator for each herd and (2) herd performance variables that could predict this indicator. Thus far, herd resilience was only possible to estimate using daily milk yield observation from the AMS and, therefore, this study used data from AMS farms. However, to predict herd resilience only data commonly available on CMS farms, i.e. herd performance variables, were used. First, the variable which is expected to represent general herd resilience is described. Second, a description is given of the variables that were used to predict this general herd resilience indicator and last, the analyses performed in this study are detailed.
General herd resilience indicator
The general herd resilience indicators used in this study were the herd-year (HY) effects, originating from Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021)'s model to estimate individual cow resilience. In that particular study, individual cow resilience was estimated by fitting an individual lactation curve using a quantile polynomial regression method. This lactation curve was used as the expected daily milk yield of each individual cow. Deviations between this expected curve and the actual observed milk yield were obtained and used for calculating the natural log-transformed variance (LnVar) from these deviations. Low values of this LnVar indicate that a cow would have a smaller deviation from this expected milk yield curve and thus is considered to be resilient. The mixed animal model, as used by Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021) to estimate individual cow resilience for primiparous cows, was as follows:
$${\boldsymbol y}_{{\boldsymbol ijk}} = {\boldsymbol H}{\boldsymbol Y}_{\boldsymbol i} + {\boldsymbol Y}{\boldsymbol S}_{\boldsymbol j} + {\boldsymbol \;}{\boldsymbol a}_{\boldsymbol k}{\boldsymbol \;} + {\boldsymbol \;}{\boldsymbol e}_{{\boldsymbol ijk}}$$where yijk was the ln-transformed variance from an individual expected lactation curve, i.e. the individual cow resilience of cow k in HY class i and year-season class j; HYi was the fixed effect of HY of calving i; YSj was the fixed effect of year-season of calving j; ak was the random genetic effect of animal k; and eijk was the random error term. A detailed description of this model is provided by Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021). The HY effects represent the effect on individual cow resilience in a given herd and a given year. These HY effects could be interpreted as the average resilience in a certain herd for all primiparous cows that calved in a certain year, corrected for their breeding values and general year-season effects. Low HY estimates would indicate less variation of deviations at herd level and thus the herd as a whole is assumed to be less affected by disturbances, indicating good resilience. Therefore, the HY effects originating from this mixed animal model will be used to represent herd resilience hereafter, where low values indicate good herd resilience.
Transformation of predictive variables
For this study, an existing dataset from 2,644 Dutch AMS farms between 2011 and 2017 was available (data previously described in Poppe et al., Reference Poppe, Mulder, Kamphuis and Veerkamp2021). The data consisted of herd performance variables, information obtained by AMS and the herd resilience indicator at HY level, i.e. one record per herd per year. To predict herd resilience, seventeen herd performance variables commonly available on farms that use milk recording were used: mean daily milk yield per herd (average of all cows), mean fat content, mean protein content, mean lactose content, mean urea content, mean somatic cell score (SCS), mean calving interval, mean age of first calving, proportion cows with elevated somatic cell count (PropSCC), proportion cows with rumen acidosis indication (PropACID), proportion cows with ketosis indication (PropKET), proportion cows that survived till second lactation (PropSURV), mean parity, mean age, number of cows in a herd, proportion of cows not 100% Holstein Friesian and proportion cows herd-book registered (Table 1). Mean SCS was derived from the SCC using (CRV, 2018):
$${\rm SCS} = 1, \;000 + 100 \times \;\left[{{\log }_2\left({\displaystyle{{{\rm SCC}} \over {1, \;000}}} \right)} \right]$$Herd performance variables, type of cows on which each herd performance variable was based, 5-year average (mean, min and max) and 5-year variance (mean, min and max) of 585 herds
a Milk = mean daily kilogram milk per farm (×10); fat = mean fat content; protein = mean protein content; lactose = mean lactose content; ureum = mean ureum content; SCS = mean somatic cell score; CI = mean calving interval from first to second lactation; AFC = mean age at first calving; PropSCC = proportion of cows with at least 1 elevated somatic cell count; PropACID = proportion of cows with at least 1 rumen acidosis indication; PropKET = proportion of cows with at least 1 ketosis indication; PropSURV = proportion of cows that survived to second lactation; Parity = mean parity; age = mean age; herd size = number of cows calved; PropNonHF = proportion of cows that are not 100% Holstein Friesian; PropREG = proportion of cows that are herd-book registered; daily milk = mean daily milk yield records from automated milking system; missing records = number of days with a missing milk yield record; records = number of days with a milk yield record.
b Primi = primiparous cows with a resilience indicator in the herd-year class; All = all cows in the herd-year class.
c Variables only available on farms with an automated milking system.
d Herd resilience indicator provided by Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021).
PropSCC was the proportion of cows with ≥1 cases of elevated SCC during lactation, where SCC was considered elevated if SCC >100 000 cells/ml. PropACID was the proportion of cows with ≥1 indications of rumen acidosis during lactation. A rumen acidosis indication was based on the fat and protein content; when fat content was lower than the protein content and below 4% the cow received a rumen acidosis indication. PropKET was the proportion of cows with a ketosis indication during lactation, where a ketosis indication was based on the fat–protein ratio and Fourier transformed infrared measurements of milk acetone and milk β-hydroxybutyric acid (Poppe et al., Reference Poppe, Veerkamp, van Pelt and Mulder2020). Incomplete observations, i.e. HY combinations with missing values on any of the aforementioned variables, were removed, resulting in a subset of 585 herds between 2012 and 2016 (2,925 HY combinations). To include possible fluctuations in herd size and health indicators such as the mean SCS, acidosis indication and ketosis indication, the variance of these HY combinations was computed (Table 1). Thus, in total 17 different 5-year averaged and 17 different 5-year variance predictive variables commonly available on dairy farms were used to predict the 5-year averaged herd resilience indicator. Predicting herd resilience for AMS farms using herd performance data has not been done before. Moreover, AMS farms have six additional herd performance variables specific for these types of farms (and thus unavailable for CMS farms). For reference purposes, we also performed the analysis including these additional six herd performance variables. These six additional variables are: both the 5-year average and the 5-year variance for (1) the daily herd milk yield, (2) the number of days with missing milk yield records, and (3) the number of days with records (Table 1). These additional six variables will be added to the 34 herd performance variables to study whether they improve the prediction of herd resilience for AMS farms specifically.
Analyses
Van der Heide et al. (Reference van der Heide, Veerkamp, van Pelt, Kamphuis, Athanasiadis and Ducro2019) and Ouweltjes et al. (Reference Ouweltjes, Spoelstra, Ducro, de Haas and Kamphuis2021) have shown previously that a Random Forest algorithm can produce reliable predictions with, respectively, the prediction of dairy cow survival till second lactation and the prediction of a lifetime resilience score in dairy cows. Lifetime resilience predicted by Ouweltjes et al. (Reference Ouweltjes, Spoelstra, Ducro, de Haas and Kamphuis2021) was defined as the cumulative results of a cow's ability to recalve and thus the ability to extend her productive lifespan (Adriaens et al., Reference Adriaens, Friggins, Ouweltjes, Scott, Aernouts and Statham2020). The current study also used the Random Forest model based on the algorithm of Breiman (Reference Breiman2001) using R package RandomForest (Liaw and Wiener, Reference Liaw and Wiener2002). Tuning of the model parameters was done using an extensive grid; number of generated trees were 500 (default), 1000, 1500, 2000 and 5000, minimum number of herds per branch were 5 (default), 10 and 15, and the random number of candidate variables at each split 1 till 34 (in total 34 predictive variables were available), where the default setting was the total number of variables divided by three. This extensive grid search resulted in 510 Random Forest models with different parameter settings. Each of these 510 Random Forest models was trained on 80% of the data and validated on the remaining 20% of the data using a 5-fold cross validation stratified to herd, i.e. each herd was used four times in the training data and once in the validation set. This grid was used on the dataset containing 34 predictive variables commonly available on CMS farms and the dataset containing 40 predictive variables (CMS variables extended with six AMS variables). To determine the best performing Random Forest model, the five validation datasets were combined. The best model was the model with the highest Pearson correlation between herd resilience and predicted herd resilience.
The importance of the predictive variables was assessed, and because a 5-fold cross validation was used, the order of importance variables could be different for each k-fold. For example, the top 10 predictive variables of k1 could be different from the top 10 of k3, because a different portion of the dataset was used as training and validation sets. To summarize this into one top 10 of most predictive variables, a ranking was made, where a value of 10 was assigned to the 1st important variable and a value of 1 to the 10th variable. For example, if the predictive variable ‘5-year average herd size’ was 2nd most important in k1, 2nd most important in k2, 3rd most important in k3, 2nd most important in k4 and 1st most important in k5, in total 5-year average herd size would be quantified with a value of 45 (9 + 9 + 8 + 9 + 10) in the total ranking of importance of predictive variables. Analyses were done using RStudio version 3.40.6 (R Studio Team, 2016). Data handling and visualization was done with, respectively R package dplyr version 0.8.5 (Wickham et al., Reference Wickham, François, Henry and Müller2020) and ggplot2 version 3.3.1 (Wickham, Reference Wickham2016).
Results
Prediction accuracy
Using the dataset containing 34 CMS features, the best performing Random Forest model was the one that generated 500 random trees, with a minimum of five herds per branch, and that randomly selected four predictive variables per split. The Pearson correlation between herd resilience and predicted herd resilience was 0.55 ± 0.06 and the Spearman rank correlation was 0.56 ± 0.06. The average for both the observed and the predicted herd resilience was 1.30. However, the range differed since the range of herd resilience was 0.70–1.86, whereas the range of the predicted herd resilience was 1.10–1.50. Nevertheless, herds with below average herd resilience tended to have below average predictions (197 herds) and herds with above average herd resilience tended to have above average predictions (212 herds; Figure 1). Thus with a 69.9% probability it was possible to predict whether a herd performs above average or below average.
Scatterplot of herd resilience vs. predicted herd resilience using 34 predictive variables commonly available on farms with conventional miking systems. The horizontal and vertical line represent the mean herd resilience (1.30). The number in each quartile represents the number of herds in that corresponding quartile.
Using the dataset containing 34 CMS features and 6 additional AMS features, the best performing Random Forest model had the same parameter settings as the one using the dataset containing only the 34 CMS features (500 generated trees, minimum of five herds per branch and randomly four selected variables per split). Adding these six AMS variables resulted in a slightly higher prediction performance, with a Pearson correlation of 0.58 ± 0.05 and a Spearman rank correlation of 0.58 ± 0.05. Consequentially, a few more herds with below average herd resilience also had below average predicted herd resilience (203 herds) and above average herds had above average predicted herd resilience (214 herds; Figure 2), increasing the probability to predict below or above average herd resilience to 71.3%. Furthermore, the Pearson and Spearman rank correlations of herd resilience predicted from either the dataset with 34 CMS or the dataset with 34 CMS and 6 AMS features were both 0.96. In other words, farms with high prediction in the first analysis also received high predictions in the second analysis.
Scatterplot of herd resilience vs. predicted herd resilience using 40 predictive variables commonly available on farms with conventional miking systems and available on farms that use an automated milking system. The horizontal and vertical line represent the mean herd resilience (1.30). The number in each quartile represents the number of herds in that corresponding quartile.
Importance of predictive variables
The importance of the predictive variables using the best performing model with 34 variables is quantified in Table 2. In each of the 5 folds the most important predictive variable was the 5-year average proportion of cows with at least 1 rumen acidosis indication (importance = 50). This was followed by the 5-year average proportion of cows with at least 1 elevated somatic cell count (SCC) and the 5-year variance of herd size. The Pearson correlations showed that farms with poor herd resilience tended to have a higher proportion of cows with a rumen acidosis indication (0.46), a higher proportion of cows with elevated SCC (0.29) and a more variable herd size over the years (0.17).
Top 10 predictive variables using 34 herd performance variables commonly available on farms with conventional milking systems
aPearson correlation with the herd resilience indicator.
bPropACID = proportion of cows with at least 1 rumen acidosis indication; PropSCC = proportion of cows with at least 1 elevated somatic cell count; Herd size = number of cows calved; SCS = mean somatic cell score; Milk = Mean kilogram milk (×10) from national database; Fat = mean fat content; PropREG = proportion of cows that are herd-book registered; CI = mean calving interval from first to second lactation.
The quantified top 10 of the best performing model using the 34 variables and additionally six AMS variables are presented in Table 3. Besides the proportion of cows with a rumen acidosis indication and a variable herd size over the years, the average number of missing records (rank 2) and variance of missing records over the years (rank 10) from the AMS tended to be important to predict herd resilience. Similar as the results from the dataset using 34 CMS variables, the Pearson correlations showed that farms with poor herd resilience tended to have a higher proportion of cows with a rumen acidosis indication (0.46), a more variable herd size over the years (0.17), and more missing daily milk yield records (0.15).
Top 10 predictive variables using 34 herd performance variables commonly available on farms with conventional milking systems and additionally six automated milking system variables (bold)
a Pearson correlation with the herd resilience indicator.
b PropACID = proportion of cows with at least 1 rumen acidosis indication; Missing records = Number of days with a missing milk yield record; Herd size = number of cows calved; PropSCC = proportion of cows with at least 1 elevated somatic cell count; SCS = mean somatic cell score; Daily milk = mean daily milk yield records from automated milking system; Milk = Mean kilogram milk (×10) from national database; Fat = mean fat content.
Discussion
This study investigated the possibility of predicting herd resilience using only herd performance variables commonly available on CMS farms. Low herd resilience values indicate less variation of milk yield deviations at herd level and thus the herd as a whole is assumed to be less affected by disturbances, indicating good resilience. This study showed that it is possible to predict whether a herd performs above average or below average with a 69.9% probability. We also showed that adding six additional variables that are only available on AMS farms increased this probability to 71.3%. Furthermore, the Pearson correlation between the analysis containing 34 CMS variables and adding six AMS variables was 0.96, indicating farms with a high prediction in the first analysis (34 variables) also received high predictions in the second analysis (40 variables). Thus, no difference was observed between the predictions using CMS or those using CMS + six additional AMS variables. This study also showed that the 5-year average proportion of cows with ≥1 rumen acidosis indication, the 5-year average proportion of cows with ≥1 elevated SCC and the 5-year variance of herd size are important variables to predict herd resilience. Three of the additional six AMS also appeared in the top 10 most important predictive variables (average number of missing records, daily milk yield observations and variance of missing records over the years). These results indicate that besides health parameters (acidosis indication and elevated SCC), the size of a farm and the information from the AMS is relevant to predict herd resilience, although this did not significantly increase the accuracy of prediction.
As previously suggested by Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021), herd resilience is the ability to control the number and severity of disturbances in a herd. Since we used HY effects from an animal model that is corrected for genetic effects as a herd resilience indicator, this indicator represents the ability to control the number and severity of disturbances in the herd through management decisions related to feed, health, housing and how quickly a farmer can respond to disturbances. For example, the provision of a well-balanced diet containing adequate amounts of vitamins and minerals could reduce the vulnerability to mastitis pathogens (Heinrichs et al., Reference Heinrichs, Costello and Jones2009) or exposure to extreme heat or extreme cold could be reduced by proper roof insulation (Fournel et al., Reference Fournel, Ouellet and Charbonneau2017), and implementing good hygiene protocols can reduce the exposure to pathogens and reduce risk of secondary infections (Deng et al., Reference Deng, Koop, Lam, van der Lans, Vernooij and Hogeveen2019).
As previously described, to estimate herd resilience, the daily milk yield observations from the AMS are needed. Thus, herd resilience of CMS farms is unknown and predicting herd resilience and simultaneously validating this prediction is not possible. Therefore, we needed to make the assumption that AMS farms also represent CMS farms. However, it is important to realize that CMS farms and AMS farms are actually different systems. For example, cows on an AMS farm can choose when they want to be milked. Hopster et al. (Reference Hopster, Bruckmaier, van der Werf, Korte, Macuhova, Korte-Bouws and van Reenen2002) has shown that social competition forces low-ranking cows to visit the AMS at different times than preferred and this potential irregularity in milking intervals could affect milk production (Ouweltjes, Reference Ouweltjes1998) and SCC (Mollenhorst et al., Reference Mollenhorst, Hidavat, van den Broek, Neijenhuis and Hogeveen2011). Both aspects may affect the results from this study as the resilience indicator itself is based on the variation of deviations from an expected lactation curve and the second most predictive variable was the proportion of cows with ≥1 elevated SCC event. Further research is needed to determine if these variables are indeed different between AMS and CMS farms and with what magnitude this affects the results of this study.
The dataset used in this study included 585 herds with data between 2012 and 2016 (after filtering). Herd performance variables were at herd level, meaning for each herd and each year an average of all cows was available. These herd level averages were averaged once more into five-year averages and to account for possible yearly fluctuations, the variance values of these five years were included. Results from Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021) showed a consistency in herd resilience estimates over multiple years (high correlation between years), indicating the fluctuation of herd resilience over years is limited. We observed a low 5-year variance (0.03) of herd resilience between 2012 and 2016 supporting these reported findings. Furthermore, predicting herd resilience of 2016 using HY variables of 2012 till 2015 resulted in a similar prediction accuracy and range (results not shown). This indicates that the abolishment of the Dutch milk quota in 2015, meaning farmers could increase herd size and total milk production per year, did not affect the results in this study.
In the analysis using 34 CMS variables we observed four predictors with a relative importance of >35. The proportion of cows with a rumen acidosis indication was the most important variable (importance 50) followed by the proportion of cows with ≥1 elevated SCC (importance 42), mean SCC in the herd (importance 36) and the variance in herd size over the five years (importance 36). Pearson correlations of these predictors with herd resilience were: 0.46, 0.29, 0.29 and 0.17, respectively. Similar correlations were observed in Poppe et al. (Reference Poppe, Veerkamp, van Pelt and Mulder2020) where a correlation between herd resilience and proportion acidosis (0.31), elevated SCC (0.20) and mean SCS (0.19) was observed. Both rumen acidosis (Krause and Oetzel, Reference Krause and Oetzel2006; Enemark, Reference Enemark2008; Abdela, Reference Abdela2016) and mastitis (Rajala-Schultz et al., Reference Rajala-Schultz, Gröhn, McCulloch and Guard1999; Gröhn et al., Reference Gröhn, Wilson, González, Hertl, Schulte, Bennett and Schukken2004; Halasa et al., Reference Halasa, Nielen, de Roos, van Hoorne, de Jong, Lam, van Werven and Hogeveen2009) can lead to an impaired milk production and thus more deviations from an expected lactation curve, with a consequential effect on the herd resilience indicator. Furthermore, Poppe et al. (Reference Poppe, Mulder, Kamphuis and Veerkamp2021) have shown that a general resilience indicator is affected by many factors and if a strong correlation between e.g. SCC and herd resilience was observed, it would rather be a mastitis indicator than a general resilience indicator. Also, the results showed that a consistent herd size (low variance of herd size over 5 years) is related to good herd resilience. A high variation indicates a large increase or a large reduction in size and we assume that major management changes are related to this, such as a new barn or different management strategies. Thus, a consistent herd size and consistent management practice indicate good resilience. Our results support the theory that the general herd resilience indicator is affected by a combination of multiple health traits and consequently, that herd management decisions related to food, health and herd size have an effect on herd resilience.
The maximum 5-year average of the two most important predictive variables was high with 0.55 and 0.96 for the proportion of cows with a rumen acidosis indication and the proportion of cows with an elevated SCC, respectively. These maximum values indicate that there is one farm where, on average, half of the cows had a rumen acidosis indication every year during 2012 and 2026, and another farm where, on average, almost all cows had a mastitis infection every year during 2012 and 2016. Comparing ten farms with high maximum 5-year averages for PropACID (range 0.42–0.55) with the ten lowest farms (range 0.02–0.04) revealed that farms with a high average of cows with a rumen acidosis indication are significantly less resilient (P < 0.001; results not shown). Furthermore, comparing ten farms with high maximum 5-year averages for PropSCC (range 0.87–0.96) with the ten lowest farms (range 0.30–0.42) revealed that also farms with a high average of cows with an elevated SCC are significantly less resilient (P < 0.001; results not shown). This indicates that herd resilience could be improved by reducing mastitis infections and rumen acidosis on farms and underlines the importance of proper herd management on farms, such as proper hygienic management practices and feeding adequate diets, to reduce the risk of mastitis infections and rumen acidosis. We can justifiably assume that poor herd resilience relates to poor health of the herd. Poor health is also related to poor fertility (Fourichon et al., Reference Fourichon, Seegers and Malher2000; Wolfenson et al., Reference Wolfenson, Leitner and Lavon2015) and cows show behavioral changes such as increased lying behavior, less feeding or physical activity if they are sick (Dittrich et al., Reference Dittrich, Gertz and Krieter2019). Therefore, it is likely that adding additional data such as fertility information or activity data would improve the prediction of herd resilience.
Table 3 shows that the AMS variable number of missing records as a mean over five years as well as the variance over five years is an important predictive variable when we include six AMS specific herd performance variables in the dataset. In practice, sick cows are separated from the herd and milked separately, meaning the number of missing records in the AMS would increase. Although low, the positive correlation of 0.27 between number of missing records and proportion of cows with an elevated SCC supports these results, meaning cows that might have mastitis are separated from the herd and not milked in the AMS and thus the number of missing records increased.
To conclude, until now, it was only possible to estimate herd resilience for farms with an AMS and only after one full lactation of a cow. However, results from the current study show that also for CMS farms it is possible to predict herd resilience. Moreover, with a 69.9% probability it is possible to predict whether a herd performs above or below average using only herd performance variables. Proportion of rumen acidosis, proportion of SCC in herds, mean SCS and fluctuation in herd size over five years are the better predictors of herd resilience. Results suggest that herd management decisions affect herd resilience, therefore, a lower predicted herd resilience could be an indication that, in general, tactical or strategic management changes could be taken to improve resilience of a herd.
Acknowledgements
We acknowledge European Union's Horizon 2020 research and innovation program (GenTORE) under grant agreement No. 727213 for their financial support. We acknowledge Prof. Dr Roel Veerkamp (Wageningen University & Research) for his contribution in getting funding for the GenTORE project. In addition, we acknowledge Cooperation CRV and CRV BV for providing the data. Finally, we acknowledge Dr Marieke Poppe for preparing the data used in this study and helping the authors to understand the data. The authors have not stated any conflicts of interest.
Financial support
None.
Open access
