Quantifying inter-group variability in lactation curve shape and magnitude with the MilkBot(®) lactation model

Genetic selection programs have driven development of most lactation models, to estimate the magnitude of animals’ productive capacity from sampled milk production data. There has been less attention to management and research applications, where it may also be important to quantify the shape of lactation curves, and predict future daily milk production for incomplete lactations since residuals between predicted and actual daily production can be used to quantify the response to an intervention. A model may decrease the confounding effects of lactation stage, parity, breed, and possibly other factors depending on how the model is constructed and used, thus increasing the power of statistical analyses. Models with a mechanistic derivation may allow direct inference about biology from fitted production data. The MilkBot(®) lactation model is derived from abstract suppositions about growth of udder capacity. This permits inference about shape of the lactation curve directly from parameter values, but not direct conclusions about physiology. Individual parameters relate to the overall scale of the lactation, the ramp, or rate of growth around parturition, decay describing the senescence of productive capacity (inversely related to persistence), and the relatively insignificant time offset between calving and the physiological start of milk secretion. A proprietary algorithm was used to fit monthly test data from two parity groups in 21 randomly selected herds, and results displayed in box-and-whisker charts and Z-test tables. Fitted curves are constrained by the MilkBot(®) equation to a single peak that blends into an exponential decline in late lactation. This is seen as an abstraction of productive capacity, with actual daily production higher or lower due to random error plus short-term environmental effects. The four MilkBot(®) parameters, and metrics calculated directly from them including fitting error, peak milk and cumulative production, can be used to describe and compare individual lactations or groups of lactations. There is considerable intra-herd and inter-herd variability in scale, ramp, decay, RMSE, peak milk, and cumulative production, suggesting that management and environment have significant influence on both shape and magnitude of normal lactation curves.