Multivariate dynamic linear models for estimating the effect of experimental interventions in an evolutionary operations setup in dairy herds.

Evolutionary operations is a method to exploit the association of often small changes in process variables, planned during systematic experimentation and occurring during the normal production flow, to production characteristics to find a way to alter the production process to be more efficient. The objective of this study was to construct a tool to assess the intervention effect on milk production in an evolutionary operations setup. The method used for this purpose was a dynamic linear model (DLM) with Kalman filtering. The DLM consisted of parameters describing milk yield in a herd, individual cows from a herd, and an intervention effect on a given day. The model was constructed to handle any number of cows, experimental interventions, different data sources, or presence of control groups. In this study, data from 2 commercial Danish herds were used. In herd 1, data on 98,046 and 12,133 milkings registered from an automatic milking system (AMS) were used for model building and testing, respectively. In herd 2, data on 3,689 milkings on test days were used for estimating the initial model parameters. For model testing, data from both bulk tank milk yield (85 observations) and test-day milkings (1,471) were used. In herd 1, the manager wanted to explore the possibility of reducing the amount of concentrate provided to the cows in an AMS. In herd 2, the manager wanted to know if the milk yield could be increased by elevating the energy level provided to the cows in a total mixed ration. The experiment conducted in herd 1 was designed with a treatment and a control group, whereas in herd 2 we used a pretest/posttest design. The constructed tool provided estimates (mean and confidence intervals) for each of 3 interventions carried out in both herds. In herd 1, we concluded that the reduction in concentrate amount provided in the AMS had no negative influence on milk yield. For herd 2, the increased level of energy had a significant positive effect on milk yield but only for the first intervention. In this herd, the effect of intervention was also evaluated for cows in the first lactation and without bulk tank records. The presented model proved to be a flexible and dynamic tool, and it was successfully applied for systematic experimentation in dairy herds. The model can serve as a decision support tool for on-farm process optimization exploiting planned changes in process variables and the response of production characteristics.

[1]  Anders Ringgaard Kristensen,et al.  Optimal slaughter pig marketing with emphasis on information from on-line live weight assessment , 2012 .

[2]  Nils Toft,et al.  Estimation of farm specific parameters in a longitudinal model for litter size with variance components and random dropout , 2002 .

[3]  George E. P. Box,et al.  Evolutionary Operation (EVOP) , 2006 .

[4]  A R Kristensen,et al.  Monitoring growth in finishers by weighing selected groups of pigs - A dynamic approach. , 2016, Journal of animal science.

[5]  M. West,et al.  Bayesian forecasting and dynamic models , 1989 .

[6]  R. Veerkamp,et al.  Associations of breed and feeding management with milk production curves at herd level using a random regression test-day model. , 2010, Journal of dairy science.

[7]  W Junge,et al.  Monitoring daily milk yields with a recursive test day repeatability model (Kalman filter). , 1999, Journal of dairy science.

[8]  S. Thamsborg,et al.  Effect of mastitis treatment and somatic cell counts on milk yield in Danish organic dairy cows. , 2003, Journal of dairy science.

[9]  N. Gengler,et al.  Prediction of daily milk, fat, and protein production by a random regression test-day model. , 2004, Journal of dairy science.

[10]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[11]  V. L. Anderson Evolutionary Operation : A Method for Increasing Industrial Productivity , 1970 .

[12]  A. Kristensen,et al.  Use of inline measures of l-lactate dehydrogenase for classification of posttreatment mammary Staphylococcus aureus infection status in dairy cows. , 2016, Journal of dairy science.

[14]  A. Bach,et al.  Associations between nondietary factors and dairy herd performance. , 2008, Journal of dairy science.

[15]  L R Nielsen,et al.  Optimal replacement policies for dairy cows based on daily yield measurements. , 2010, Journal of dairy science.

[16]  L. Pedersen,et al.  Using Evolutionary Operation technique to evaluate different management initiatives at herd level , 2016 .

[17]  S. Calsamiglia,et al.  Effect of amount of concentrate offered in automatic milking systems on milking frequency, feeding behavior, and milk production of dairy cattle consuming high amounts of corn silage. , 2007, Journal of dairy science.

[18]  A De Vries,et al.  Application of statistical process control charts to monitor changes in animal production systems. , 2010, Journal of animal science.

[19]  Douglas C. Montgomery,et al.  Response Surface Methodology: Process and Product Optimization Using Designed Experiments , 1995 .

[20]  J. Pedersen,et al.  A Nordic proposal for an energy corrected milk (ECM) formula , 1991 .

[21]  I. Halachmi,et al.  Feeding soyhulls to high-yielding dairy cows increased milk production, but not milking frequency, in an automatic milking system. , 2009, Journal of dairy science.

[22]  Iver Thysen,et al.  Monitoring Bulk Tank Somatic Cell Counts by a Multi-Process Kalman Filter , 1993 .

[23]  Henk Hogeveen,et al.  Bayesian integration of sensor information and a multivariate dynamic linear model for prediction of dairy cow mastitis. , 2016, Journal of dairy science.

[24]  J M Siegford,et al.  Invited review: The impact of automatic milking systems on dairy cow management, behavior, health, and welfare. , 2012, Journal of dairy science.

[25]  H. Volden,et al.  Norfor - the nordic feed evaluation system. , 2011 .

[26]  G. H. Kroeze,et al.  Description of a detection model for oestrus and diseases in dairy cattle based on time series analysis combined with a Kalman filter , 1999 .