Bayesian spatio-temporal modelling of national milk-recording data of seasonal-calving New Zealand dairy herds

A spatio-temporal analysis was undertaken with the aim of identifying the dynamics of herd mean individual cow SCCs (MICSCC) in seasonally calving New Zealand dairy herds. Two datasets were extracted from the Livestock Improvement Corporation's extensive national dairy recording database: (1) milk-recording data aggregated at the herd-level and (2) sales questionnaire data containing information on the size, location and infrastructure of each farm. A Bayesian spatio-temporal modelling approach was applied to the analysis. The data were aggregated by 10 km² grid cells and linear regression models were developed with spatially structured and unstructured random effects, a linear temporal trend random effect and spatial-temporal interactions for log-transformed median MISCC (ln(median MISCC)).

Significant associations were found between ln(median MISCC) and milk yield, milk fat, milk protein, farm area and number of cups in the dairy. This led us to suggest that SCCs should be adjusted for volume and constituents prior to determining a threshold MISCC for identification of subclinical mastitis (SCM) problem herds. Part, or all, of the temporal trend in MISCC in the spatio-temporal model was accounted for by inclusion of yield and milk constituents as independent variables. This supports the hypothesis of a dilution effect with potential consequences for misdiagnosis of SCM, particularly in late lactation. Unmeasured covariates were similarly likely to be spatially structured and unstructured.