 Research article
 Open Access
 Published:
Bayesian analysis of herdlevel risk factors for bovine digital dermatitis in New Zealand dairy herds
BMC Veterinary Research volume 15, Article number: 125 (2019)
Abstract
Background
Bovine digital dermatitis (BDD) is considered the most important infectious cause of lameness in dairy cattle worldwide, but has only recently been observed in New Zealand. Although many studies have investigated the risk factors for BDD in confined dairy systems, information on risk factors in pasturebased system is limited. Therefore a crosssectional study including 59,849 animals from 127 dairy herds in four regions of New Zealand was conducted to identify the herdlevel factors associated with the probability of a herd being BDDlesion positive and with withinherd BDD prevalence.
Results
Purchasing heifers was associated with increased odds of a herd being BDDlesion positive (odds ratio [OR]: 2.33, 95% probability interval [PI]: 1.26–4.42) and a cow being BDD affected (OR: 3.76, 95%PI: 1.73–8.38), respectively. Higher odds of a herd being BDDlesion positive (OR: 2.06, 95%PI: 1.17–3.62) and a cow being BDD affected (OR: 2.87, 95%PI: 1.43–5.94) were also seen in herds where heifers cograzed with cattle from other properties. In addition, using outside staff to treat lameness was associated with higher odds of a cow being BDD affected (OR: 2.18, 95%PI: 0.96–4.98).
Conclusion
This study highlighted that movements of heifers are significantly associated with the spread of BDD within and between dairy herds in New Zealand. To minimise the risk of disease introductions in herds where moving heifers cannot be avoided, it is best to purchase heifers only from herds where BDDfreedom has been confirmed and, if heifers have to grazeoff a farm, they should be reared as a single biosecure management group, especially since animals may be BDDinfected without having clinically obvious lesions.
Background
Bovine digital dermatitis (BDD) has been found throughout the world in both confined and pasturebased dairy systems [1, 2]. In many countries, BDD appears to be endemic in dairy herds [3] and is commonly considered as the most important infectious cause of cattle lameness [4]. Clinically, BDD lesions progress or regress through different morphological stages, commonly described using M scores [5, 6]. A rapid BDD lesion detection method such as visual examination during milking is widely used in many studies [7]; although interpretation of such diagnostic outcome is subjective, which usually requires additional validation studies to assess the agreement across the examiners [8].
Multiple studies have evaluated the risk factors associated with BDD prevalence within herds in confined dairy systems. These studies have identified a widerange of potential risk factors including type of housing [9], using outside staff to trim hooves [10], footbath regimen [11] and access to pasture [12]. In contrast, very few studies [13,14,15] have been undertaken in cattle that are principally pasturebased with no or very limited use of housing, where many of the risk factors identified in confined animals are irrelevant. Specific research in such systems is essential as there can be large variation between pasturebased dairy herds in the prevalence of BDD [16].
In New Zealand, one previous study has evaluated herdlevel risk factors for BDD, but that was undertaken in only one region [14]. In that study we used a Bayesian hurdle model to explore the associations between risk factors and BDD prevalence at both the herd and animal levels. The initial separation of the herds into BDDlesionfree and BDDlesion positive was based on whether BDD lesions were observed; i.e. a herd with ≥1 lesion was defined as being BDDlesion positive, otherwise it was defined as being BDDlesion free [14]. However, simply basing herd status on the presence/absence of visible lesions probably leads to loss of information regarding probability of a herd having BDD and may introduce misclassification bias at the herd level, as there is a chance that a herd where BDD lesions are truly present could be wrongly classified as being BDDlesionfree due to a combination of limited diagnostic sensitivity and low cowlevel prevalence [17].
One method for overcoming this limitation is by using a Bayesian latent class model, which estimates the mean probability of a herd being BDDlesion positive conditional on the number of test positive animals, the total number of animals tested, and the test characteristics [17]. Thus, the mean probability contains more precise information than the simple dichotomised outcome and increases the power of the study to determine the impact of risk factors on the likelihood of a herd being BDDlesion positive.
The aim of this study was to use Bayesian methods to investigate the impact of farm management practices on pasturebased dairy herds across New Zealand on 1) the probability of a herd being BDDlesion positive obtained from a previous Bayesian latent class analysis [16] and 2) the withinherd BDD prevalence, namely the probability of a cow within a herd having BDD lesions.
Methods
Target and source population
The target population was the pasturebased dairy herds in New Zealand and the source population was the herds in the four regions across New Zealand: Waikato and Manawatu in the North Island and the West Coast and Canterbury in the South Island. These regions encompass most of the dairy systems (all grass fed and selfcontained; feed imported, either supplement or grazingoff and feed imported to extend lactation) used in New Zealand [16].
Data collection
The dataset was collected as described by Yang et al. [16]. Briefly, the data collection started in the Waikato and moved south following the seasonal pattern of calving to ensure that the great majority of the herds were milking at the herd examinations. In the first phase, half the sampled herds were visited in each region before moving on to the next. In the second phase, the order was reversed, starting in Canterbury and going back north. Within each herd, visual assessment was performed on cows’ rear feet in the milking parlour after hosing the feet gently [7].
The farm management practices undertaken in the previous 12 months were collected alongside the visual inspection for BDD using a questionnaire given to the owners or managers of the study herds. The questionnaire was modified by the authors from that used in Yang et al. [14]. The questionnaires were answered after the herd inspection while the first author was still on the farm, so that if the owners or managers were unsure of a question, the first author could explain the intent of the question. The categorical management predictors collected via the questionnaire are shown in Table 1 and a copy of the questionnaire is provided as an additional file (see Additional file 1).
Data processing
The data were imported into Stata 13.1 for cleaning and analysis (StataCorp, USA). Oneway tables were used to examine the frequency of responses for each level within the categorical variables. Levels with low frequencies were combined with adjacent levels where biologically plausible. As hoof trimmers were rarely used to trim cows or treat lame cows, this level was combined with using vets to treat lame cows, to create a new dichotomous variable of whether or not the farm had outside staff trimming cows or treating lame cows. Since few farmers reported chemically disinfecting hoof trimming equipment, “chemical disinfection” and “washed by water” were combined to create a new variable of whether or not trimming equipment was cleaned between animals. Since few farmers reported purchasing dairy heifers or cows from saleyards, “saleyards” was combined with “other farms”, to create new variables of whether heifers and cows were purchased from outside. Similarly, cow houses were rarely used, so “cow house” was combined with “standoff pad” and a new variable was created to describe whether the cows were permanently pasturebased (except for milking) or not. All categorical variables included in the final analysis had at least two levels and each level had at least ≥15% of the total responses for the question.
Evaluating seasonal variation
As this was a crosssectional study, the impact of season on BDDlesion status and lesion prevalence was not of primary interest. However, as all the data were collected by the first author, it was not possible to complete the data collection in a short time frame and therefore herds were sampled at different points throughout the 2016/2017 lactation season. To confirm that BDDlesion status and lesion prevalence did not vary significantly between different months, generalized estimating equations [18] and beta regression models were respectively used to examine whether the average cowlevel prevalence or probabilities of a herd being BDDlesion positive differed significantly between months in Waikato and Manawatu regions. This process was not applied to Canterbury since 18/19 herds were visited in the same month.
Univariable models
For the two outcome variables (withinherd prevalence and probability of a herd being BDDlesion positive), univariable logistic regression models and univariable beta regression models in the frequentist framework were respectively used to select predictors for fitting in the multivariable models. Any predictors with pvalue ≤0.2 were included in the further analyses.
Multivariable model 1
This analysis was designed to quantify the strength of associations between farm management practices and withinherd prevalence. A Bayesian binomial model was constructed. The model was built using a forward stepwise strategy. The predictors were retained in the model when the 90% probability intervals of their corresponding regression coefficients did not overlap 0. If inclusion of a predictor altered the coefficient of any of the existing predictors by > 15%; the newly included predictor was considered as a confounder and was forced into the model regardless its 90% probability interval [19]. Once the preliminary main effect model was constructed, twoway interactions between all predictors in the model were created. An interaction term was retained if its 95% probability interval excluded 0. The final model structure is presented below:
where y_{j} was the number of the cows with visible BDD lesions in the j^{th} herd of all the regions, which was modelled using a binomial distribution with the parameters: the proportion of cows with visible lesions (p_{j}) and number of cows being examined (n_{j}); β_{0} was the intercept, β_{1}, β_{2} and β_{3} were the regression coefficients for the predictors x_{j}, g_{j} and h_{j} which represented whether heifers were purchased from outside sources, whether heifers were cograzed with heifers from other farms and whether outside personnel treated lame cows. Finally, U_{region(j)}, W_{j} were the random effects at regional and herd level, respectively and modelled using two independent normal distributions with zero means and standard deviations σ_{U} and σ_{W}.
The choice of prior distributions contributes to the posterior distributions, thus utilising informative priors results in better inferences compared to “vague” priors [20]. It is difficult to place informative priors for the regression coefficients. However, such priors can be indirectly induced to define probabilities for different combinations of predictors. Partially informative priors were assigned to β_{0}, β_{1} and β_{2}. First, the proportion of cows with visible lesions (prevalence) in a “typical” closed herd was defined as \( {\overset{\sim }{p}}_0 \). This meant x = 0, g = 0 and h = 0. Therefore, according to Eq. 1, \( {\beta}_0=\mathrm{logit}\ \left({\overset{\sim }{p}}_0\right) \). Second, specify \( {\overset{\sim }{p}}_1 \) as the prevalence of a herd where some heifers were purchased from outside, in this case, x = 1, but g = 0 and h = 0. Thus, \( {\beta}_1=\mathrm{logit}\ \left({\overset{\sim }{p}}_1\right){\beta}_0 \). Finally, let \( {\overset{\sim }{p}}_2 \) denote the prevalence of a herd which contained purchased heifers, and, at the same time, sent its own heifers to cograze with heifers from other farms (x = 1 and g = 1, but h = 0). This gave \( {\beta}_2=\mathrm{logit}\ \left({\overset{\sim }{p}}_2\right)\mathrm{logit}\ \left({\overset{\sim }{p}}_1\right) \). Logitnormal distribution was used for these prevalence priors. Below we use \( {\overset{\sim }{p}}_0 \) as an example to illustrate the way to convert a prevalence estimate to its corresponding logitnormal distribution such as \( \mathrm{logit}\ \left({\overset{\sim }{p}}_0\right)\sim N\left({\mu}_{\beta 0},{\tau}_{\beta 0}\right) \), where τ is the precision term defined as the reciprocal of the variance. Our best estimate of the prevalence in a closed herd was m_{0} and we were 95% confident that it was less than l_{0}; then μ_{β0} = logit(m_{0}). The standard deviation σ_{β0} = [logit(l_{0}) − logit(m_{0})]/1.645, eventually \( {\tau}_{\beta 0}=1/{\sigma}_{\beta 0}^2 \).
The best estimates of \( {\overset{\sim }{p}}_0 \), \( {\overset{\sim }{p}}_1 \) and \( {\overset{\sim }{p}}_2 \) came from the previous analyses of BDD data in Taranaki and the authors’ expert opinion. One important observation was that in contrast to previous studies of housed cattle, the apparent cowlevel prevalence of BDD was very low (mean = 1.2%, Yang et al. [1]) with 26.8% of herds having fewer than1% of cows with observed lesions. In Canterbury region, where median herd size was 840 and > 21% of herds had ≥1000 cows; we were able to detect BDD at an apparent withinherd prevalence of 0.1% (i.e. one cow with lesions in a 1000cow herd). Thus to reflect our belief that a closed herd was likely to have no or extremely rare BDD lesions, we took 0.05% as our “best point estimate” for \( {\overset{\sim }{p}}_0 \). Furthermore we were also 95% confident that it was less than 0.35%, i.e., one cow with BDD lesion(s) in a 300cow herd. Based on the method described in the last paragraph, this led to μ_{β0} = − 7.6 and τ_{β0} = 0.71. Table 2 summarises our “best estimates” for \( {\overset{\sim }{p}}_0 \), \( {\overset{\sim }{p}}_1 \) and \( {\overset{\sim }{p}}_2 \). Uniform priors (0, 3) and (0, 2) were set for σ_{U} and σ_{W}, respectively. This reflected our belief that the variability of herdlevel prevalence across regions was bigger than the variability across herds. However, the parameter values assigned to the uniform priors were considered to be nonspecific as we did not know the standard deviations of the two random effects.
Under the partially informative priors, the fit of the model to the data was evaluated using posterior predictive checks which compared the observed outcome data to the data simulated/predicted by the posterior predictive distribution [21]. The Bayesian Pvalue quantifies the probability that the discrepancy between the predicted and observed values. A Bayesian Pvalue close to 0.50 indicates adequate model fit, although a value between 0.20 and 0.80 is also accepted [22].
Sensitivity analysis was used to assess the sensitivity of the posteriors to the priors. Table 3 summarises the distributions of the model priors and the priors used for the sensitivity analysis. The model was developed using OpenBUGS [23]. Posterior inferences were obtained using Markov chain Monte Carlo (MCMC) approximation. The posterior distribution of each parameter was reported using median and 95% probability interval (PI). After discarding the first 10,000 iterations as burnin period, the model was further run for 100,000 iterations. Convergence was assessed using BGRplots by running three chains with different sets of initial values [24].
Multivariable model 2
This analysis was designed to assess the associations between farm management predictors and the probability of herd being BDDlesion positive (PP). This analysis did not include herds on the West Coast as the region was determined to be free of the disease [16].
The data were modelled using a Bayesian beta model [25]. π_{k} was used to denote the PP_{k} for k^{th} herd . The variable “region” was initially modelled as a random effect \( {V}_{region(k)}\sim N\left(0,{\frac{1}{\sqrt{\tau}}}_V\right) \), where τ_{V} was the precision term. The model was constructed as follows:
where z_{k} was the predictor vector, γ denoted the regression coefficient vector and μ_{k} the mean and φ a measure of variability, with a larger value of φ indicating less variability [26].
Diffuse normal distributions (mean = 0, precision = 0.01) were set for all the regression coefficients, and a vague gamma distribution (1, 1) was set for φ and τ_{V}. The model was built using a forward stepwise strategy. Predictors were retained if the 90% probability interval for the regression coefficients excluded 0. Confounders were assessed using the method described as per Multivariable model 1. Twoway interactions between all predictors in model were investigated after building the main effect model. Inclusion criteria for an interaction term were the same as for Multivariable model 1. In this model, the linear predictor was the logodds. The odds were defined as the probability of a herd being BDDlesion positive divided by the probability of a herd being BDDlesionnegative at each level of a predictor. The model was therefore able to identify any farm management practice associated with higher odds of being BDDlesion positive for a randomly selected herd in any BDDaffected region.
Multivariable model 3
Although the mixed effects beta model modelled the overall variability of the probability in different regions; it was not able to describe the difference between particular regions, therefore we also built a model which treated “region” as a fixed effect. Assuming the model had in total t farm management practices, Eq. (4) was changed to:
with \( {V}_{region(k)}\sim N\left(0,{\frac{1}{\sqrt{\tau}}}_V\right)\ \mathrm{dropped} \). Here, z_{ck} and z_{wk} were the dummy variables for the regions Canterbury and Waikato (level “Manawatu” was treated as reference level). This fixed effects model can be used to predict the probability of a herd being BDDlesion positive with different covariates in any particular region.
The deviance information criteria (DIC) of both beta models were compared. In addition, a global measure of variation explained by each of the beta models was obtained by computing pseudoR^{2} defined as the squared correlation between the linear predictor and the logittransformed outcome variable [27]. Both beta models were developed using OpenBUGS [23]. After discarding the first 5000 iterations as the burnin period, the model was further run for 100,000 iterations. Convergence was assessed using BGRplots by running three chains with different sets of initial values [24]. The OpenBUGS code for Multivariable model 1, 2 and 3 is provided as an additional file (see Additional file 2).
Results
There was no evidence to support seasonal differences in any of the outcome variables. In Waikato region, the average cowlevel prevalences in September 2016 and in January 2017 were not significantly different (P = 0.94). The probabilities of BDDlesion positive also did not differ significantly between these two months (P = 0.65). In Manawatu, the average cowlevel prevalences in September (P = 0.46) and November (P = 0.22) were not significantly different to that in December. Similarly, significant differences in the probabilities of BDDlesion positive in September (P = 0.86), November (P = 0.28) and December were not evident. These findings ruled out the potential seasonal impact on BDD prevalences in this study. Table 4 displays the total herds and animals sampled as well as the proportions of herds/animals having BDD lesions in each region during the data collecting period.
The outputs from the Bayesian binomial model (Multivariable model 1) with our partially informative priors are shown in Table 5. Lack of model fit was not evident (Bayesian Pvalue = 0.5). The posteriors for β_{1}, β_{2} and β_{3} were robust in all sensitivity analysis scenarios. The posterior median for β_{0} increased slightly (− 8.05 vs. − 7.67) and its 95%PI was also wider (− 10.65, − 5.707) given the diffuse prior N (− 5, 0.001) rather than the informative prior. The posterior for σ_{W} was not sensitive to its prior, although the posterior for σ_{U} was sensitive to its prior. The posterior median of σ_{U} increased from 2.4 to 3.2 when the prior changed from Uniform (0, 3) to Uniform (0, 5). It further increased to 4 if the prior changed to Uniform (0, 9). Nevertheless, there was no impact on the posteriors for the regression coefficients. The results of the sensitivity analyses are provided as an additional file (see Additional file 3).
Based on Multivariable model 1, cattle in a herd which purchased heifers from outside were more likely to have BDD lesions than cattle in a herd that did not purchase heifers (OR: 3.76, 95%PI: 1.73–8.38). Being in a herd which cograzed heifers with animals from other properties also increased the odds of a cow having BDD lesions (OR: 2.87, 95%PI: 1.43–5.94). The use of outside staff to treat lameness was found to be associated with the increased withinherd prevalence (OR: 2.18, 95%PI: 0.96–4.98).
Except for the intercepts, the posteriors for the parameters reported by the Bayesian mixed effects beta model and fixed effects beta model were nearly identical. Table 6 summarises the models’ outputs. The DIC for each model was also very similar, − 533.3 for the fixed effects model and − 533.5 for the mixed effects model. Two farm management practices were identified as being significantly associated with the odds of a herd being BDDlesion positive. Based on the mixed effects beta model, the odds of a herd being BDDlesion positive was 2.33 times (95%PI: 1.26–4.42) higher in a herd with purchased heifers compared to one without, and 2.06 times (95%PI: 1.17–3.62) higher if heifers cograzed with cattle from other properties. The predicted probabilities from the fixed effects model of a herd being BDDlesion positive conditional on different farm management practices and region are displayed in Fig. 1. However, 74% of the variation in the probability of a herd being BDDlesion positive remained unexplained by either beta model (pseudoR^{2} = 0.26).
Discussion
This study found that both cograzing with heifers from other properties and purchasing heifers from other farms were associated with an increased probability of a herd being BDDlesion positive as well as increased withinherd prevalence. Our previous study [14] also found that that youngstock movement between farms considerably increased the probability of a farm having at least one visible lesion (OR: 4.15, 95%PI: 1.39–15.27). Compared to Yang et al. [14], our current study evaluated youngstock movement in a more detailed way by dividing such movements into heifer purchasing and heifer cograzing. Unlike Yang et al. [14] who reported that youngstock movement affected only a herd’s probability of having at least one cow with BDD lesions but not withinherd prevalence, this broaderscale study (along with the more detailed way of recording the predictors) found that youngstock movement increased both the probability of a herd being BDDlesion positive and the withinherd prevalence. The likely inference is that heifers act as a reservoir for BDD transmission between dairy herds and between cows within herds in New Zealand [28]. In contrast, both this analysis and Yang et al. [14] found no effect of purchasing adult cattle on BDD risk. This lack of effect is most likely due to the much smaller numbers of purchased adult cows compared to the numbers of heifers purchased for replacement [13, 14].
Yang et al. [14] reported that herds with a rotary platform were more likely to have at least one cow with BDD lesions than herds with a herringbone (OR: 3.19, 95%PI: 1.31–8.51), though as with heifer movement, no effect was seen on withinherd prevalence. This may have been due to the ease of finding at least one lesion in herds with rotary platforms rather than being an actual risk factor [14]. Our current analysis did not include parlour type in the final model as the analysis found it to be neither statistically significant nor a confounder.
Two New Zealand studies [14, 15] reported that on BDDpositive farms, the withinherd prevalence was higher on farms where the outside staff came for hoof trimming (prevalence ratio [PR]: 3.13, 95%PI: 1.25–7.29 and risk ratio [RR]: 2.06, 95% confidence interval [CI]: 1.05–4.06). Although our current analysis did not confirm this finding, the calculated OR 2.18 was still in the realm considered to be biologically important [10] and the 95%PI: 0.96–4.98 only just included 1. It is not entirely clear how these effects could be mediated under New Zealand conditions. The use of outside staff for lame cows is typically unrelated to BDD since BDD rarely causes lameness in New Zealand dairy cattle. However, failure to clean trimming equipment properly between cows and between herds could represent a mechanism for spread [14]. To confirm this hypothesis, our current study included whether trimming equipment was cleaned between cows as a potential risk factor. However, no effect of cleaning/washing equipment between cows was found; this suggests that if there is an effect of outside staff on the withinherd prevalence of BDD that it is not mediated via dirty equipment. Further research is required to better estimate the impact of using outside staff to treat lame cows on BDD prevalence and to investigate potential pathways by which such an effect could be mediated.
The only other study of risk factors for BDD in pasturebased cows is that by RodriguezLainz et al. [13]. However, of the 22 farms in that study only 2 kept their cattle at pasture all year round, with 13/22 keeping cattle in an open corral or loose yard for at least part of the year, whereas in this study, all 127 farms grazed their cattle throughout the year. As such many of the factors analysed by RodriguezLainz et al. [13] (e.g. housing type and season of calving) are not directly relevant to the New Zealand situation and thus not included in our analysis. Although it is difficult to directly compare the study findings, RodriguezLainz et al. [13] did find that there was an effect of purchasing replacement heifers on withinherd prevalence of BDD (OR: 3.16, 95%CI: 1.61–6.21), but not purchasing adult cows (OR: 1.31, 95% CI: 0.72–2.38). The data from RodriguezLainz et al. [13] provided no evidence as to whether, in pasturebased cattle, using outside staff to trim feet increases the withinherd prevalence of BDD as in that study all cattle were treated or trimmed by farm staff. However, using hoof trimmers who operated on multiple farms was found to be significantly associated with higher BDD withinherd incidence in housed cattle (OR: 2.8, 95% CI: 1.9–4.2) [10].
Many studies on dairy cattle from intensive housing systems in the northern hemisphere have also identified herdlevel risk factors for this disease. Decreasing the access to pasture was found to increase the risk of BDD [12, 29, 30]. The type of housing for animals was also associated with BDD prevalence, i.e. cows that housed in cubicles had higher BDD prevalence and more severe BDD cases [31] than cows in straw yards, which also agreed with Onyiro et al. [9]. In cubicles houses, the size of cubicles was linked to the risk of BDD [29]. This is because cows tend to spend longer time standing in shorter and narrower cubicles; therefore the contact between heels and slurry was increased [32]. However, these factors tend not to be an issue in New Zealand pasturebased systems and were therefore not included in the current study. It could be interesting in future studies to evaluate cleanliness of legs in cattle since higher prevalence of BDD had been found in cows with dirty legs [33]. This is one possible explanation why no BDD lesions were seen on the West Coast where the cows’ feet were generally much cleaner compared to other regions.
The results of this study show that even in New Zealand where BDD prevalence is very low, heifers are the most likely source of disease spread between and within herds. Particular care should be taken when purchasing heifers as replacement animals and ideally, replacement heifers should only be purchased from herds where BDDfreedom has been confirmed. The latter may be difficult in New Zealand since many heifers are purchased in lateautumn when cows are not being milked and it is therefore not possible to observe the milking herd for BDD. In such cases, visual inspection of the whole heifer group (not just the heifers for purchase) is a potential alternative to increase the probability of finding at least one animal with BDD lesions. If any of the heifers have visible lesion(s), then the entire group should not be purchased as animals can still be infected with BDD in the absence of visible lesions [34]. Where heifers are cograzing with animals from multiple herds, it becomes much more difficult to ensure that cograzing heifers will not come in contact with BDD infected cattle, although little is currently known about the transmission dynamics of BDD in grazed dry stock. Thus, the only reliable method to ensure that heifers grazed away from the farm do not become infected with BDD is to require that they are grazed as a single biosecure management group. This is important to prevent the spread of many infectious diseases as well as BDD.
Bayesian methods were adopted as the analytical approach in this paper. Bayesian analyses incorporate previous scientific understanding, e.g. such as the likely association between a farm management practice and BDD withinherd prevalence, into analysis (see Multivariable model 1), so that the inference (i.e. the posterior distribution) is based on both the data and our prior information. This is in contrast to other methods which typically ignore such previous understanding [35]. Furthermore even if previous information of a research question is not available, the Bayesian methods still has significant advantages such as being able to directly compare the relative probabilities of two or more hypotheses rather than simply using the probability of the data given the null hypothesis to determine whether an alternative hypothesis was plausible.
Multivariable model 2 and 3 used uniform priors, as this was the first use of beta models to study risk factors on the herdlevel BDD outcome estimated from a previous Bayesian latent class analysis. This use of the outcome from the latent class analysis reduced the likelihood of misclassification errors at the herd level, as the effect of diagnostic sensitivity and specificity on the herd level diagnosis was factored into the latent class model [16].
Although misclassification bias has been adjusted at the herd level in Multivariable model 2 and 3, our Multivariable model 1 did not account for animal level misclassifications. This could potentially have influenced the analysis of risk factors affecting withinherd prevalence. Misclassification at the individual level, as at the herd level, can be minimised by incorporating the known sensitivity and specificity of a diagnostic method [36]. However, when the impact of specificity and sensitivity on the diagnosis of BDD in the individual animal was included during the modelling process, it resulted in nonconvergence of the Markov chains. This may be related to the model being nonidentifiable. Using a more sensitive detection method inspecting lifted cows’ feet in the trimming chute, would have decreased any potential impact but would have been cost prohibitive [7].
The other limitation was that Multivariable models 2 and 3 explained only 26% of the variation in the probability of a herd being BDDpositive. This indicates that further investigation of more factors which could potentially affect the probability of herd being BDDlesion positive was required.
Conclusions
Our study investigated potential risk factors for BDD across New Zealand and identified that purchasing replacement heifers and cograzing heifers with animals from other herds were significantly associated with a higher probability of a herd being BDDlesion positive and higher withinherd prevalence of BDD. This is consistent with previous findings from pasturebased systems. However, the identified risk factors only explained a small proportion of the variation in probability of a herd being BDDlesion positive. Our study also found that using outside staff for trimming had a large effect on withinherd prevalence (doubling the odds of an individual cow having BDD). Given that we can’t rule out the possibility of contaminated hoof trimming equipment contributing to the betweenherd spread of BDD, it would be advisable for farms to maintain their own set of equipment. Further research should be undertaken to better estimate the impact of this factor on BDD and how it can be mediated through different biosecurity interventions.
Abbreviations
 BDD:

Bovine digital dermatitis
 CI:

Confidence interval
 DIC:

Deviance information criteria
 MCMC:

Markov chain Monte Carlo
 OR:

Odds ratio
 PI:

Probability interval
 PP:

Probability of herd being BDDlesion positive
 PR:

Prevalence ratio
 RR:

Risk ratio
References
Yang DA, Heuer C, Laven R, Vink WD, Chesterton RN. Farm and cowlevel prevalence of bovine digital dermatitis on dairy farms in Taranaki, New Zealand. N Z Vet J. 2017;65(5):252–6.
Solano L, Barkema HW, Mason S, Pajor EA, LeBlanc SJ, Orsel K. Prevalence and distribution of foot lesions in dairy cattle in Alberta, Canada. J Dairy Sci. 2016;99(8):6828–41.
Orsel K, Plummer P, Shearer J, De Buck J, Carter S, Guatteo R, Barkema H. Missing pieces of the puzzle to effectively control digital dermatitis. Transbound Emerg Dis. 2018;65:186–98.
Laven R, Lawrence K. An evaluation of the seasonality of veterinary treatments for lameness in UK dairy cattle. J Dairy Sci. 2006;89(10):3858–65.
Döpfer D, Koopmans A, Meijer F, Szakall I, Schukken Y, Klee W, Bosma R, Cornelisse J, Van Asten A, Ter Huurne A. Histological and bacteriological evaluation of digital dermatitis in cattle, with special reference to spirochaetes and campylobacter faecalis. Vet Rec. 1997;140(24):620–3.
Berry SL, Read DH, Famula TR, Mongini A, Döpfer D. Longterm observations on the dynamics of bovine digital dermatitis lesions on a California dairy after topical treatment with lincomycin HCl. Vet J. 2012;193(3):654–8.
Yang DA, Laven RA. Detecting bovine digital dermatitis in the milking parlour: to wash or not to wash, a Bayesian superpopulation approach. Vet J. 2019;247:38–43.
Yang DA, Laven RA. Interobserver agreement between two observers for bovine digital dermatitis identification in New Zealand using digital photographs. N Z Vet J. 2019;67(3):143–7.
Onyiro O, Andrews L, Brotherstone S. Genetic parameters for digital dermatitis and correlations with locomotion, production, fertility traits, and longevity in HolsteinFriesian dairy cows. J Dairy Sci. 2008;91(10):4037–46.
Wells S, Garber L, Wagner B. Papillomatous digital dermatitis and associated risk factors in US dairy herds. Prev Vet Med. 1999;38(1):11–24.
Speijers M, Baird L, Finney G, McBride J, Kilpatrick D, Logue D, O’Connell N. Effectiveness of different footbath solutions in the treatment of digital dermatitis in dairy cows. J Dairy Sci. 2010;93(12):5782–91.
Holzhauer M, Brummelman B, Frankena K, Lam T. A longitudinal study into the effect of grazing on claw disorders in female calves and young dairy cows. Vet J. 2012;193(3):633–8.
RodriguezLainz A, MelendezRetamal P, Hird DW, Read DH, Walker RL. Farmand hostlevel risk factors for papillomatous digital dermatitis in Chilean dairy cattle. Prev Vet Med. 1999;42(2):87–97.
Yang DA, Laven RA, Heuer C, Vink WD, Chesterton RN. Farm level risk factors for bovine digital dermatitis in Taranaki, New Zealand: an analysis using a Bayesian hurdle model. Vet J. 2018;234:91–5.
Yang DA, Laven RA, Chesterton RN. Effects of climate and farm management practices on bovine digital dermatitis in springcalving pasturebased dairy farms in Taranaki, New Zealand. Vet J. 2019;247:75–80.
Yang DA, Johnson WO, Müller KR, Gates MC, Laven RA. Estimating the herd and cow level prevalence of bovine digital dermatitis on New Zealand dairy farms: a Bayesian superpopulation approach. Prev Vet Med. 2019;165:76–84.
Yang DA, Heuer C, Laven R, Vink WD, Chesterton RN. Estimating the true prevalence of bovine digital dermatitis in Taranaki, New Zealand using a Bayesian latent class model. Prev Vet Med. 2017;147:158–62.
Zeger SL, Liang KY, Albert PS. Models for longitudinal data: a generalized estimating equation approach. Biometrics. 1988;44:1049–60.
Dohoo I, Martin W, Stryhn H: Veterinary Epidemiologic Research, 2nd edn. Charlottetown, Canada: VER Inc; 2010.
Dunson DB. Commentary: practical advantages of Bayesian analysis of epidemiologic data. Am J Epidemiol. 2001;153(12):1222–6.
Gelman A, Meng XL, Stern H. Posterior predictive assessment of model fitness via realized discrepancies. Stat Sin. 1996:733–60.
Neelon BH, O’Malley AJ, Normand SLT: A Bayesian model for repeated measures zeroinflated count data with application to outpatient psychiatric service use. Stat Model 2010, 10(4):421–439.
Spiegelhalter D, Thomas A, Best N, Lunn D. OpenBUGS user manual, version 3.0. 2. Cambridge: MRC Biostatistics Unit; 2007.
Brooks SP, Gelman A. General methods for monitoring convergence of iterative simulations. J Comput Graph Stat. 1998;7(4):434–55.
Branscum AJ, Johnson WO, Thurmond MC. Bayesian beta regression: applications to household expenditure data and genetic distance between footandmouth disease viruses. Aust N Z J Stat. 2007;49(3):287–301.
Branscum A, Gardner I, Johnson W. Bayesian modeling of animaland herdlevel prevalences. Prev Vet Med. 2004;66(1):101–12.
Ferrari S, CribariNeto F. Beta regression for modelling rates and proportions. J Appl Stat. 2004;31(7):799–815.
Laven R, Logue D. The effect of precalving environment on the development of digital dermatitis in first lactation heifers. Vet J. 2007;174(2):310–5.
Somers J, Frankena K, NoordhuizenStassen E, Metz J. Risk factors for digital dermatitis in dairy cows kept in cubicle houses in the Netherlands. Prev Vet Med. 2005;71(1):11–21.
Read DH, Walker RL. Papillomatous digital dermatitis (footwarts) in California dairy cattle: clinical and gross pathologic findings. J Vet Diagn Investig. 1998;10(1):67–76.
Laven R. The environment and digital dermatitis. Cattle Practice. 1999;7:349–54.
Laven R: Determination of the factors affecting the cause, prevalence and severity of digital dermatitis as a major cause of lameness in dairy cows. Milk Dev Counc Study 2000, 95(May):1–5.
Relun A, Lehebel A, Bruggink M, Bareille N, Guatteo R. Estimation of the relative impact of treatment and herd management practices on prevention of digital dermatitis in French dairy herds. Prev Vet Med. 2013;110(3):558–62.
Vink W, Jones G, Johnson W, Brown J, Demirkan I, Carter S, French N. Diagnostic assessment without cutoffs: application of serology for the modelling of bovine digital dermatitis infection. Prev Vet Med. 2009;92(3):235–48.
Johnson WO. Comment: Bayesian statistics in the twenty first century. Am Stat. 2013;67(1):9–11.
McGlothlin A, Stamey JD, Seaman JW Jr. Binary regression with misclassified response and covariate subject to measurement error: a Bayesian approach. Biom J. 2008;50(1):123–34.
Acknowledgements
The authors thank all the farmers, veterinarians and vet technicians for their involvement in this study. We acknowledge the help for farm recruitment from the following companies, organisation and individual: VetEnt – Te Awamutu, Otorohanga and Ashburton clinics, Totally Vets – Palmerston North clinic, Southern Rangitikei Veterinary Service (SRVS) – Bulls clinic, West Coast Vets – Hokitika and Greymouth clinics, Farm Services, School of Veterinary Science, Massey University and Mr. Shove.
Funding
The first author is funded by the Massey University Doctoral Scholarship. The funding source had no role in the design of this study and will not have any role during its execution, analyses, interpretation of the data, or decision to submit results.
Availability of data and materials
The subset of the data was available in Additional file 2. The full data collected and analysed during the current study are not publicly available as they contain confidential information of the participated farmers. However, the datasets are available from the corresponding author on reasonable request.
Author information
Affiliations
Contributions
DAY, KRM and RAL participated in the study design and coordination. DAY collected the data and analysed the data. DAY, MCG and RAL contributed to the manuscript. All authors were involved in the manuscript preparation and approved the final manuscript.
Corresponding author
Ethics declarations
Authors’ information
Dr. Yang obtained his bachelor’s degree in veterinary medicine and a master’s degree in vet epidemiology; he is working towards his PhD in vet science.
Dr. Gates is a senior lecturer in veterinary epidemiology.
Dr. Müller (or Mueller as the Massey University website) is a senior veterinarian in dairy cattle health and production.
Prof Laven is a professor in production animal health and welfare.
Ethics approval and consent to participate
The New Zealand Animal Welfare Act (1999) states that if an animal is subject to a manipulation, it needs ethics approval.
Section three of this act defines a manipulation as follows:
In this Act, unless the context otherwise requires, the term manipulation, in relation to an animal, means, subject to subsections (1A) to (3), interfering with the normal physiological, behavioural, or anatomical integrity of the animal by deliberately—.
(a) subjecting it to a procedure which is unusual or abnormal when compared with that to which animals of that type would be subjected under normal management or practice and which involves—.
(i) exposing the animal to any parasite, microorganism, drug, chemical, biological product, radiation, electrical stimulation, or environmental condition; or.
(ii) enforced activity, restraint, nutrition, or surgical intervention; or.
(b) depriving the animal of usual care;—.
As washing of feet to observe BDD is normal management and the animals were observed during milking (so there is no restraint beyond normal) what was done does not meet the definition of a manipulation.
The farm owners were identified by local veterinary practices and verbal agreement was obtained prior to visiting farms.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1:
The questionnaire used to collect farm management practices. (DOCX 176 kb)
Additional file 2:
OpenBUGS code for the three Bayesian multivariable models used in this study. (DOCX 19 kb)
Additional file 3:
Results of the sensitivity analysis for Multivariable model 1. (DOCX 18 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Yang, D.A., Gates, M.C., Müller, K.R. et al. Bayesian analysis of herdlevel risk factors for bovine digital dermatitis in New Zealand dairy herds. BMC Vet Res 15, 125 (2019). https://doi.org/10.1186/s1291701918713
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1291701918713
Keywords
 Digital dermatitis
 Dairy cattle
 Lameness
 Risk factors
 Pastoral system
 Bayesian
 Multilevel modelling