# Sample Entropy Analysis of Horses with and without Lameness at Trot before and after Sedation

Authors: Jing Zhao; Dan B. Marghitu; Amey Rane; John Schumacher
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The horses were trotted in a straight line, and the acceleration data were recorded from the movement of the horse’s head. Inertial sensors were used to collect the information. The sample entropy was used to characterize the dynamics of horse gait. The gait of sound and lame horses was compared, and the sound horses had smaller values of the sample entropy. The effects of sedative detomidine on the gait of lame horses were also examined. The values of sample entropy for sedated horses were smaller than the values of the sample entropy for untreated horses. The sample entropy can provide important information about the dynamics of the equine gait.


---

## Body

## 1. Introduction

The concept of entropy originated from information technology and was introduced to quantify the amount of information obtained in communication signal [1]. The approximate entropy (ApEn) has been developed to classify complex systems that include both deterministic and stochastic processes [2]. The approximate entropy was not suitable for the analysis of the noisy data that are prevalent in biological studies, and inconsistent results were produced. As a modification of the ApEn, the sample entropy (SampEn) was designed to increase the accuracy of ApEn over a broad range of conditions [3]. The sample entropy has been useful in the study of the postural control and the center of pressure displacements [4, 5]. Stergiou et al. adopted entropy to measure the change in postural control for an infant with cerebral palsy and found that entropy values decreased significantly from prop sitting to beginning independent sitting [6]. Jia et al. reported that sample entropy of the functional connections of brain regions in healthy subjects decreased with advancing age [7]. A smaller value of the SampEn also indicates more self-similarity in dataset or less noise. The sample entropy has higher value for more complicated time series. For cardiovascular time series, one fundamental hypothesis is that healthy subjects have more complicated time series than cardiovascular patients [8].Lameness in horses causes serious economic losses to the equine industry [9]. Detection of lameness is not difficult when lameness is moderate to severe, but confusion appears when lameness is subtle. There are disagreements between equine clinicians during subjectively evaluating lameness due to diversity of skills developed for lameness recognition. Clinicians are not able to capture small changes in locomotion pattern because of the relatively low image capturing frequency of human eye [10]. Detomidine, a commonly used sedative in horses, is often administered to unruly horses to facilitate lameness examination [11, 12]. However, there is a concern among veterinarians about the effect of sedation on movement patterns of horses [13, 14]. Many clinicians are reluctant to use sedative during lameness examination as they believe it might affect the identification of subtle lameness because of the inherent analgesic property of this drug [15, 16]. The effect of sedation on the gait of lame horses has been studied using vector sums of maximum and minimum head height and the sum of hip drop and hip hike [17].Multiple methods have been initiated for studying the nonlinear dynamics of gait and human motion which may include model development, dynamic formulation, and simulation solver [18]. The nonlinear dynamics equations are directly integrated or differentiated. The simulation solver demands experimental motion data and optimization-based methods. For the repeated motion of the hip during walking, an analytical approach has been proposed to obtain solution for nonlinear oscillators [19]. Optimal auxiliary functions are employed to approximate the analytical solution of the system. The deep learning technique has been used for the gait patterns [20]. With the help of deterministic learning algorithms, detailed information about gait dynamics can be produced. The authors applied the knowledge fusion strategy for gait aspects.Equine locomotion involves rhythmic and coordinated patterns of body and limb segments [21]. An equine body system can be considered as a dynamic system, and the changes occurring in the body due to lameness can be analyzed according to dynamic laws. For the last two decades, using inertial sensors and high-speed camera capturing devices, researchers have studied the dynamic locomotion of the horse.The concept of entropy was initiated from the energy loss during combustion and not converted into functional work. In thermodynamics, any event in an isolated system will make the entropy larger. The topological entropy estimates, in time, the progression of detectable orbits for an iterated function describing the dynamical system. For a stochastic system, the entropy is the average rate at which the information is produced. The proposed definition of entropy (SampEn) is a chosen option in our study.The objective of this paper was to implement the SampEn for the equine gait. First, we studied the trot of sound and lame horses. Next, we compared the gait of unsedated lame horses with their gait when sedated with detomidine. The vertical movement data of the horse’s head were used for this study. To our best knowledge, the use of SampEn on the analysis of horse gait was not studied. The SampEn can be utilized for the dynamic analysis of the horse gait.


## 2. Materials and Methods

Ten horses that were determined by subjective evaluation from the investigator to be lame and eight to be nonlame were enrolled in this study. The mean weight of the horses was 460 ± 50 kg. In the first study, the gait from sound horses and lame horses was compared by calculating SampEn. Horses were trotted for at least 25 consecutive strides during each of 10 random trials while wearing sensors.In the second study, the effect of sedative detomidine on gait of ten lame horses was studied using two groups: Group 1 (no treatment, control group) and Group 2 (horses administered 10 mg of detomidine HCL). Horses were randomly assigned to two groups [17]. Lame horses were trotted about 25 consecutive strides every 5 minutes for 45 minutes. This study was approved by Auburn University’s Animal Care and Use Committee.Two inertial sensors were attached noninvasively to the horse as shown in Figure1. A single-axis accelerometer was mounted on the head halter to measure the vertical acceleration of the head. A single-axis gyroscope was attached to dorsal surface of the region between the metacarpophalangeal joint and coronary band of the right forelimb to measure the rotation of the foot. The signal from sensors was transmitted to a tablet PC via Bluetooth. The raw data collected from the sensors were processed in MATLAB.Figure 1
Inertial sensors attached noninvasively to the horse: (a) uniaxis vertical accelerometer head sensor and (b) uniaxial gyroscope right forelimb sensor.The data were in tsv format with a unique 32-digit hex address. The raw output value in this file was digitized (8-bit) from the voltage value (+/−5 volts). A conversion factor was used to get the actual acceleration data. Figure2 shows a simultaneous plot of vertical head acceleration and right forefoot angular velocity. Stride split was achieved by using the information from the gyroscope located at the right forelimb. Figure 3 shows the head acceleration data gathered from ten strides.Figure 2
Vertical head acceleration and right forefoot angular velocity.Figure 3
Head acceleration for ten strides.


### 2.1. Sample Entropy

Acceleration data for each stride were time normalized to 50 data points. The whole motion cycle with 20 strides was time normalized with 1000 data points. This data length is considered sufficient for the entropy analysis [22]. After signal is segmented into strides and normalized, the entropy calculation was carried out. The SampEn gives us the probability of two data that are close (predetermined tolerance level) in m dimensional space and remain close in m+1 dimensional space. Larger is the probability less is the SampEn value and vice versa.The sample entropy of a time seriesx1,x2,…xN of length N with zero mean and unit standard deviation is defined by(1)SampEnN,m,r=−logArBr,where m is the length of the subsequences of the time series; r is the maximum distance (or tolerance) between two subsequences of the same length m; Ar is the total number of subsequences in a m+1 dimensional phase space that are within the distance r; and Br is the total number of subsequences in a m dimensional phase space that are within the tolerance r. The complexity level of two signals could be determined from this method. For the computation of the SampEn, the guideline detailed in Ramdani et al. was followed [23].Sample entropy is found to be closely related to selection ofm and r for short datasets. The embedding dimension m and the tolerance r are selected according to Lake et al. [24]. The procedure is based on the convergence criterion which estimates m leading to the selection of r. The median SampEn is calculated for all datasets as a function of r for m=1,2,3,4,5.Figure4 shows that the entropy values start converging when m is equal to 3 which lead us to select m=3. Next step is to compute the median of maximum relative error as a function of r for m=3. The value of r is then estimated to give maximum relative error no higher than 0.05.Figure 4
Median SampEn as a function ofm and r.


### 2.2. Statistical Analysis

The difference of SampEn between sound and lame horses was compared witht test. The SampEn values were subjected to a log transformation to ensure the homogeneity of variance. To test the effect of detomidine on SampEn, linear mixed effects model Y=Xβ+Zμ+e was performed, where Y was the SampEn values, X was the treatment design matrix (treatment with detomidine or no treatment; data collection at 5, 10, 15, 20, 25, 30, 35, 40, and 45 minutes), βwas the vector of fixed treatment effects, Z was the subject design matrix (horse), μwas a vector of random subject-specific effects, and e was a vector of random experimental error. The normality of the data was tested using the Shapiro–Wilk test. First-order autoregression has been used as the correlation structure. If a significant influence was found, paired t tests were performed to determine where significant differences were present. The significance level was set at p<0.05. Statistical analyses were performed using SAS (SAS Studio 3.8, SAS Institute Inc., Cary, NC, USA.).


## 2.1. Sample Entropy

Acceleration data for each stride were time normalized to 50 data points. The whole motion cycle with 20 strides was time normalized with 1000 data points. This data length is considered sufficient for the entropy analysis [22]. After signal is segmented into strides and normalized, the entropy calculation was carried out. The SampEn gives us the probability of two data that are close (predetermined tolerance level) in m dimensional space and remain close in m+1 dimensional space. Larger is the probability less is the SampEn value and vice versa.The sample entropy of a time seriesx1,x2,…xN of length N with zero mean and unit standard deviation is defined by(1)SampEnN,m,r=−logArBr,where m is the length of the subsequences of the time series; r is the maximum distance (or tolerance) between two subsequences of the same length m; Ar is the total number of subsequences in a m+1 dimensional phase space that are within the distance r; and Br is the total number of subsequences in a m dimensional phase space that are within the tolerance r. The complexity level of two signals could be determined from this method. For the computation of the SampEn, the guideline detailed in Ramdani et al. was followed [23].Sample entropy is found to be closely related to selection ofm and r for short datasets. The embedding dimension m and the tolerance r are selected according to Lake et al. [24]. The procedure is based on the convergence criterion which estimates m leading to the selection of r. The median SampEn is calculated for all datasets as a function of r for m=1,2,3,4,5.Figure4 shows that the entropy values start converging when m is equal to 3 which lead us to select m=3. Next step is to compute the median of maximum relative error as a function of r for m=3. The value of r is then estimated to give maximum relative error no higher than 0.05.Figure 4
Median SampEn as a function ofm and r.


## 2.2. Statistical Analysis

The difference of SampEn between sound and lame horses was compared witht test. The SampEn values were subjected to a log transformation to ensure the homogeneity of variance. To test the effect of detomidine on SampEn, linear mixed effects model Y=Xβ+Zμ+e was performed, where Y was the SampEn values, X was the treatment design matrix (treatment with detomidine or no treatment; data collection at 5, 10, 15, 20, 25, 30, 35, 40, and 45 minutes), βwas the vector of fixed treatment effects, Z was the subject design matrix (horse), μwas a vector of random subject-specific effects, and e was a vector of random experimental error. The normality of the data was tested using the Shapiro–Wilk test. First-order autoregression has been used as the correlation structure. If a significant influence was found, paired t tests were performed to determine where significant differences were present. The significance level was set at p<0.05. Statistical analyses were performed using SAS (SAS Studio 3.8, SAS Institute Inc., Cary, NC, USA.).


## 3. Results

When applied fast Fourier transform, no major differences were found between lame and sound horses (Figure5). The phase planes are affected by the duration of the detomidine treatment but there are no limit cycles or clear attractors (Figure 6). For a periodic solution, the phase plane is a closed curve. For the return maps, there are differences between 5 min and 45 min following detomidine treatment but the systems are not periodic or quasiperiodic (Figure 7). The return map of a periodic system consists of a finite set of points. For a quasiperiodic system, a closed orbit tends to display on the return map. If a system has chaotic motion, an infinite number of points with a well-defined structure usually make up the return map. We cannot find distinct points, segments, or curves on the return maps.Figure 5
Fast Fourier transform for (a) a sound horse and (b) a lame horse.
(a)(b)Figure 6
Phase plane for a lame horse at (a) 5 min and (b) 45 min following treatment with detomidine.
(a)(b)Figure 7
The return map for a lame horse at (a) 5 min and (b) 45 min following treatment with detomidine.
(a)(b)Gait from sound horses was associated with significantly smaller SampEn than lame horsesp<0.0001 (Figure 8).Figure 8
Comparison of SampEn between sound and lame horses.Treatment with detomidine had a significant effect on the SampEn of lame horses compared with no treatmentp=0.0225 (Figure 9). At the 15 min data collection, the SampEn was significantly lower after treatment with detomidine compared with no treatment p=0.0438. No significant differences were found at other times of data collection.Figure 9
Mean SampEn of ten horses at 5, 10, and every 5 min over a total period of 45 min following no treatment and treatment with detomidine.


## 4. Discussion

Based on the results of this study, the sample entropies of lame horses are larger than that of sound horses. The sample entropies were significantly smaller at 15 min after the intravenous administration of 10 mg detomidine compared with no treatment. Our results are consistent with the results of Buchner et al. which observed a decrease of head height and increase of stride duration and stride length after detomidine sedation, with the maximum effect of detomidine at 15 min recording [11].The investigation of sample entropy may reveal important information about the characteristics of locomotion pattern. Higher entropy values are associated with a functional decline of the locomotion control system, indicating a greater need for adaptability to perturbations in the system [5]. In contrast, the signal emanated from a system with low entropy values generates less additional information over time and is more predictable [25]. The less complex gait was observed in the sound horses in this study.It has been reported that gait complexity of human adults walking on a treadmill increased with increasing walking speed, possibly related to the mechanism of movement control [26]. In this study, however, horses were trotted by a single handler for consistency of speed. Evaluation of the effect of speed on sample entropy of horses should be considered for future study to help understanding the control system in horses. Furthermore, future studies of different gait of sound horses would benefit from using sample entropy as a new measurement variable.


## 5. Conclusion

In spite of achievements in the area, there is no common measure to evaluate the dynamics of animal and human locomotion. In the present study, we investigated the equine gait of sound and lame horses. By comparison, smaller SampEn was identified on sound horses. Moreover, we observed the effect of detomidine on gait of lame horses. The sample entropy values for horses with detomidine were smaller than the SampEn of horses that did not receive a treatment. In the present work, we have shown that, using an appropriate methodology, the SampEn algorithm was able to characterize the dynamic features of equine gait. The technique in this study gives the clinical team a useful measure of the gait. It is expected that this research will have straight implementation to the analysis of gait in humans and animals.



---

*Source: 4827587-2020-07-13.xml*