|
|
||||||||





* Bou-Matic, LLC, Madison, WI 53708
Department of Mathematics and Statistics, University of Maryland at Baltimore County (UMBC), Baltimore 21250
Department of Animal and Food Sciences, University of Delaware, Newark 19717
Department of Mechanical Engineering, UMBC, Baltimore 21250
|| Agriculture Experiment Station (Central Maryland Research and Education Center), University of Maryland, Ellicott City 21042
# Department of Animal and Avian Science, University of Maryland, College Park 20742
1 Corresponding author: parimal{at}boumatic.com
| ABSTRACT |
|---|
|
|
|---|
Key Words: dairy cow lameness automated lameness detection gait score lesion score
| INTRODUCTION |
|---|
|
|
|---|
Methods of objective, automated scoring may enable earlier, more accurate detection of lameness. Early, accurate detection could lead to more efficacious practices that reduce lameness in commercial dairy units. Cattle caretakers, claw trimmers, and veterinarians normally detect lame cows by visual observation. One popular method for detecting lameness is to score cows from 1 (sound) to 5 (severely lame) based on subjective observations of postural and stride characteristics (Sprecher et al., 1997). To begin introducing automated, continuous, and objective methodology into lameness detection, the authors recently developed an automated system for detecting lameness based on vertical ground reaction force (GRF) measurements of individual limbs (Rajkondawar et al., 2002b; Tasch et al., 2004). The system is called the reaction force detection (RFD) system and consists of 2 parallel-force plates with leveling platforms before and after the plates. Side railings guide the cows so that the left limbs contact the left plate and right limbs the right plate. Vertical forces measured over time for each plate can be used to calculate a number of limb movement variables (LMV).
The step of cattle is divided into a weight-bearing (stance) and non-weight-bearing (swing) phase. The weight-bearing phase is initiated at the time the foot lands on the ground and is terminated after push-off when the foot is lifted off the ground. During this phase, vertical and horizontal GRF are transmitted by the foot to the ground. The vertical GRF increases until midstance and then decreases until the foot is lifted off the ground (Scott, 1988, 1989; Van der Tol et al., 2003). Lameness changes the maximal force applied by lame limbs, relative times of peak loading, and deacceleratory and acceleratory longitudinal forces in the lame as well as the nonlame contralateral limbs (Scott, 1989). Conceivably, these changes could provide the basis for the development of objective methods for lameness detection. In this regard, a pilot study demonstrated that 5 LMV generated from the RFD system can be used to detect lameness (Rajkondawar et al., 2002a).
The objective of the current study was to determine how best to use the 5 LMV to detect lameness. An additional objective was to develop statistical models that related measures of foot and hind limb health to LMV derived from the RFD system. The effect of including more days in model development also was investigated. The models were then tested using a previously unevaluated set of cows from the same herd.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Experimental Design and Procedures
Cows were fed a balanced TMR consisting primarily of corn and alfalfa silage; grass hay; soybeans; cottonseed; ground, shelled corn; and vitamins and minerals. They were housed in a free-stall facility with stalls containing rubber mats and sawdust, and hooves were trimmed by the herdsmen following an undefined procedure at the time of dry off. Floors of the alleyways were deeply grooved cement with crosshatched 6.4-cm squares. As cows were assigned to the study, the hooves were trimmed by the veterinarian according to the Dutch technique guidelines (van Amstel et al., 2000) and any lesions were noted and scored. Gaits were scored for all 4 limbs. Regardless of the number, score, and type of lesion, when all 4 limbs received a score of 1, cows were not considered lame.
Data were collected from 10 to 15 cows at a time. Cows were allowed to walk through an RFD system located in the return alley from the parlor during 2 to 3 h before milking. Each cow in the study was examined once weekly by the same veterinarian. The veterinarian assigned GS and LS, and treated cows as necessary. Cows with normal GS in all 4 limbs (score = 1) for 2 consecutive weeks were released from the study. The herdsman randomly selected replacement cows without regard for gait status from the low-production group (<16 kg of milk/d) or the dry cow lot. The modeling data were constructed during a 6-mo period using 21 Holstein cows of mixed age. A test data set was acquired using 10 cows during a subsequent 6-wk period (Table 1
).
LMV
Vertical GRF measured over time for individual hind and forelimbs were used to calculate 5 LMV (Rajkondawar et al., 2002b). The maximum of the sum of instantaneous GRF for a cow was used as an estimate of her BW. To enable between-cow comparisons, the LMV were scaled to dynamic BW.
GS
Once weekly, cows were first observed in a stance and then while walking in a straight line. To facilitate identification of lame limbs, cows also were circled to the right and to the left. All examinations were performed in an alleyway on concrete flooring. Gait observations were performed by observing the cow perpendicular, parallel, posterior, and anterior to the line of travel. During gait evaluation, cows and limbs were observed for freedom of motion, left- and right-sided stride length, length of the anterior and posterior swing phases, symmetry and arc of the foot flight, foot placement relative to body position and limb axis, foot rotation during weight bearing, symmetry of weight distribution at the walk and stance, and position of the top line at a walk and in a stance. A score of 1 to 5 was assigned for each limb, as described in Table 3
. The scoring procedure is an adaptation of other methods (Wells et al., 1993; Sprecher et al., 1997). To score 2 limbs on the same cow, the lamer of the 2 limbs was assigned the GS (Table 3
) normally assigned to the cow. The unscored second limb was assigned the same score or a lower score from Table 3
, depending on the relative level of gait changes displayed in the second limb.
|
Lesions were diagnosed as described elsewhere (Bergstem, 1997; Collick et al., 1997; Leach et al., 1998). All lesions were scored as summarized in Table 4
. A lesion score of 1 signified no visible lesion.
|
Modeling Approach
Analyses were based on data collected for individual hind limbs, and all weekly sets of measurements for individual limbs were considered to be statistically independent experimental trials. Several models of lameness were developed using GS, LS, and combined scores (CS). The CS was defined as the maximum of contemporaneous GS and LS, or max {GS, LS}, for a single limb. Gait and lesion scores for individual limbs were related to the 5 LMV for corresponding limbs using logistic regression (Hosmer and Lemeshow, 2000). Because of a very limited number of GS and LS greater than 3, raw scores of 4 and 5 were not included in the development of the lameness models. Therefore, the GS and LS used in the analysis had ordinal values of 1, 2, or 3. Logistic regression models were developed to generate predictive equations for the probability that a limb had been assigned a specific score as an S-shaped function of the 5 LMV. [Note: The mathematical form of an S-shaped function is f(x) = eß1x/1 + eß2x ]. The general form of these models is as follows:
![]() |
![]() |
![]() |
where P(A) is the probability that A is true, and f represents a weighted linear combination of the 5 LMV. The coefficients for these models were estimated using the maximum likelihood methods applied to the underlying multinomial distribution. For details, the reader is referred to Hosmer and Lemeshaw (2000) and Dyer et al. (2004). These cumulative logistic regression models are used widely in the statistical literature to model ordinal scores (Hosmer and Lemeshaw, 2000). Five LMV were calculated for each limb every day. Models were constructed based on 1 d (1D), 2 d (2D), and 3 d (3D) of LMV observation. Therefore, the 1D, 2D, and 3D models used 5, 10, and 15 LMV, respectively. Based on the above probabilities, a lameness index (LI) was calculated as a weighted sum:
![]() |
where "score" was GS, LS, or CS.
Model Validation
Model validity was evaluated using sums of squared errors and ROC (Hosmer and Lemeshow, 2000; Appendix A). The sum of squared error of the gait score (EGS) was calculated as
![]() |
where GS was the gait score and LIGS was the lameness index estimated by the model being tested. The summation extended over all N data points, and ELS and ECS were similarly defined.
The area under the ROC curve (AUC) represents the level of discrimination achieved by a model. The AUC values can be used to compare the performance of various models much as correlation coefficients are used to compare results of linear regression models. Interpretations of AUC values were as described in Hosmer and Lemeshow (2000): 1) 0.7
AUC < 0.8 (acceptable discrimination); 2) 0.8
AUC < 0.9 (very good discrimination); and 3) 0.9
AUC (excellent discrimination). A detailed description of ROC in terms of true and false positives, and true and false negatives is provided in Appendix A.
| RESULTS |
|---|
|
|
|---|
1. A scale parameter estimate of 1 corresponds to perfect independence. For GS, LS, and CS, the smallest values (1.27, 1.04, and 1.25, respectively) occurred for the 3D model. Thus, correlations among observations from the same limb were best accounted for by the inclusion of the LMV from the previous days, demonstrating an advantage of the 3D model. These overall small values provided justification for the assumption that all weekly sets of measurements for individual limbs could be considered to be statistically independent experimental trials.
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Other confounding factors may have had a potentially negative impact on the robustness of the lameness models reported herein. The nature of the lameness problem in the experimental population may have been problematic. Lameness was associated with single or multiple lesions within claws, single or multiple lesions between claws, and lameness in single or multiple limbs. It is unclear how single or multiple lesions between and within claws, or single and multiple limb lameness influenced the process of modeling. In addition, the models were generated with data restricted to one type of herd with one type of management, housing, ration, and flooring. Despite these potential limitations, the data were sufficiently robust that the general relations in terms of the relative importance and value of the 5 LMV, and the relative ability to associate the 5 LMV to GS and LS were of value for refining lameness detection models. Evidence of the robustness of the data was the finding that the effect of lack of independence attributable to repeated measures diminished when LMV from additional days were included in the models.
Cost is an important factor to consider when developing systems for agricultural use. To reduce the potential cost of a commercial system to detect lameness, the RFD system used in this study measured only vertical reaction forces. An earlier study demonstrated that 5 LMV derived from vertical GRF could be used to identify lame cows (Rajkondawar et al., 2002a). However, a variety of changes in the GRF in both the vertical and longitudinal components are associated with lameness (Scott, 1989; Clayton et al., 2000). Thus, measurement of all components of GRF could potentially provide a more robust approach to lameness modeling. The tradeoff would be a greatly increased cost of the measurement system.
The mean of the 5 LMV (Tables 5
and 6
) tended to decrease with increasing LS and GS. However, changes within any single LMV were not sufficient to distinguish sound from lame cows within either GS or LS. These findings are in accord with Scott (1989), who showed that bovine lameness was associated with a variety of inconsistent changes in GRF. Inconsistent GRF changes between and within lame cows would increase the variability of the 5 LMV and therefore decrease the discriminatory value of any single LMV in lameness detection. For equines, normal (nonlame) animals showed 5 to 8% variability in vertical GRF and 2 to 6% variability in stance time (Hodson et al., 2001; Khumsap et al., 2001). Lameness increased the variability 3- to 4-fold (Clayton et al., 2000). In our study, variability tended to increase with lameness, but did not approach the magnitude reported by Clayton et al. (2000).
Variability of the 5 LMV measured in our study generally exceeded that observed in equine studies for a few plausible reasons. It is presently unclear how much of the variability can be attributed to the RFD, the cows, or the methodology. Normal limbs in this study included all individual hind limbs scored sound as normal gait or normal gait with the absence of lesions. So-called normal limbs included limbs on completely sound cows or the single sound limb on a cow with unilateral lameness. This approach could be problematic particularly in the latter case. Compensatory changes in vertical forces applied to the contralateral sound limb of a unilaterally lame cow could be expected to generate vertical GRF that differ from GRF generated by normal limbs of a bilaterally sound cow with symmetrical weight bearing (Scott, 1989). Differences could be expected to increase LMV variability within the pool of normal limbs used in our modeling. Similar arguments could be raised about LMV for any limb contralateral to a lame limb. Second, LMV data used in the equine studies were selected from a larger pool of force signatures generated by multiple repeated passes. This selection process likely resulted in an underestimation of the variability in LMV in many equine studies. In this bovine study, all 5 LMV were derived from single observations. This method of data assimilation is considerably different from that used in equine studies and may present a more accurate assessment of the variability in the 5 LMV of sound and lame limbs.
Results showed that the lesion models identified cows with lesions better than the gait models identified them with abnormal gaits. The reason for this difference is not clear. One possibility may be that identification of lesions was less subjective than identifying changes in gait characteristics. Earlier studies showed that observer agreement tended to be good when gaits of moderately lame cows or limbs were scored. However, observer agreement was much poorer when cows and limbs with subtle lameness were scored (Wells et al., 1993; Mill and Ward, 1994; Keegan et al., 1998). In the current study, subtle lameness dominated the experimental group in that 34, 46, and 18% of the limbs in the investigation scored 1, 2 and 3, respectively. Although we approached the difficulties of scoring subtle lameness with a meticulously detailed lameness examination, the risk of subjectivity and low reproducibility in gait scoring may have remained problematic because of the large prevalence of mild lameness.
Accuracy of the CS model was no greater than the LS model because LS were generally
GS for most limbs. The magnitude of the maximum lesion severity score was
GS in 92.4% of the 656 limbs with LS and GS. Thus, modeling lameness by selecting the greater of the 2 scores for the CS model produced results that were nearly identical to that utilized in the LS model.
Another confounding problem is the temporal nature of the data collection. For each data record, the 5 LMV were measured daily and the corresponding GS and LS were obtained on the day after the last measurement of the 5 LMV. Conceivably, gait adjustment to accommodate pain, changes in lesion sensitivity, or movement of the nociceptive threshold (Whay et al., 1998) during the 24- to 96-h period of LI computation could have generated disparities between LI scores and GS or LS, or both. The AUC values for the 1D, 2D, and 3D lesion and gait models may support this contention. The AUC values for the GS models worsened with additional days of data. This could reflect the very fluid nature of pain-driven changes in locomotion with time. In contrast, AUC values for the LS models improved as additional days were incorporated into the model (Table 8
). This may reflect the inability of lesions (compared with gait) to change appearance during the same time period.
Opportunities also exist to improve the lesion scoring system. An improved lesion scoring system might allow construction of more accurate detection models. The scoring scheme used in this investigation (Table 4
) generated a degree of graded severity within a lesion, but did not present an accurate assessment of lesion severity across different types of lesions. For example, sole ulcers with a score of 2 and sole hemorrhages with a score of 2 were equivalent in terms of modeling, but could have very different effects on pain perception. Hemorrhage in horn tissue is regarded as historical evidence of corium damage (Greenough and Vermunt, 1991; Ossent et al., 1997), and therefore little correlation exists between lesions and gait manifestations in limbs bearing hemorrhagic stains of horn tissue. Temporal disassociation between the time of hemorrhage in the corium and the appearance of hemorrhage on weight-bearing surfaces likely created inaccuracies in the lesion model and underscores the need for weighting lesion severity scores for pain. This problem might be addressed in the future by the development of objective measures of pain.
Spatial variability of force distributions within and between medial and lateral claws may also have complicated lesion-based modeling. It has been demonstrated that the maximal forces are normally borne by the lateral claw rather than the medial claw of the hind limb (Van der Tol et al., 2002, 2003, 2004). Within regions of the claw, forces tended to be maximally applied in regions 3, 5, and 6 (Van der Tol et al., 2002). All things being equal, one could predict that lesions with similar severity scores appearing in regions 3, 5, and 6 of the lateral hind claw might generate greater changes in the 5 LMV than similar score lesions appearing in the medial claw or in regions 1, 2, and 4 of the lateral claws. Thus, modeling lameness by lesion scoring without regard for either pain assessment or regional distribution runs the risk of inappropriate lesion score assignment during modeling. Again, this problem might be addressed by the development of objective measures of pain.
The lack of objectivity in scoring lesions and locomotion presents a fundamental problem for modeling lameness with objective LMV data. One solution would be to develop objective measures of the dependent variables. In this regard, it might be appropriate to develop an objective, graduated measure of claw pain and then relate the 5 LMV to this measure. Another possibility would be to incorporate an exact measurement of lesion size and location, as described by Leach et al. (1998). With independent, continuous measures of lesions and locomotion, it may be possible to model lameness more accurately and relate lameness problems to milk and reproductive losses.
| CONCLUSIONS |
|---|
|
|
|---|
| APPENDIX A: ROC CURVE |
|---|
|
|
|---|
|
![]() |
An ROC is plotted for a binary classification model by applying different classification thresholds, determining the values listed above, and plotting the result as shown in Figure A1
.
|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication February 11, 2005. Accepted for publication June 20, 2006.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. T. Thomsen, L. Munksgaard, and F. A. Togersen Evaluation of a Lameness Scoring System for Dairy Cows J Dairy Sci, January 1, 2008; 91(1): 119 - 126. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Dyer, N. K. Neerchal, U. Tasch, Y. Wu, P. Dyer, and P. G. Rajkondawar Objective Determination of Claw Pain and Its Relationship to Limb Locomotion Score in Dairy Cattle J Dairy Sci, October 1, 2007; 90(10): 4592 - 4602. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. C. Bicalho, S. H. Cheong, G. Cramer, and C. L. Guard Association Between a Visual and an Automated Locomotion Score in Lactating Holstein Cows J Dairy Sci, July 1, 2007; 90(7): 3294 - 3300. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. Pastell and M. Kujala A Probabilistic Neural Network Model for Lameness Detection J Dairy Sci, May 1, 2007; 90(5): 2283 - 2292. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |