# Sample entropy characteristics of movement for four foot types based on plantar centre of pressure during stance phase

- Zhanyong Mei
^{1, 2}, - Guoru Zhao
^{1}, - Kamen Ivanov
^{1, 2}, - Yanwei Guo
^{1, 2}, - Qingsong Zhu
^{1}, - Yongjin Zhou
^{1}and - Lei Wang
^{1}Email author

**12**:101

https://doi.org/10.1186/1475-925X-12-101

© Mei et al.; licensee BioMed Central Ltd. 2013

**Received: **7 May 2013

**Accepted: **27 September 2013

**Published: **10 October 2013

## Abstract

### Background

Motion characteristics of CoP (Centre of Pressure, the point of application of the resultant ground reaction force acting on the plate) are useful for foot type characteristics detection. To date, only few studies have investigated the nonlinear characteristics of CoP velocity and acceleration during the stance phase. The aim of this study is to investigate whether CoP regularity is different among four foot types (normal foot, pes valgus, hallux valgus and pes cavus); this might be useful for classification and diagnosis of foot injuries and diseases. To meet this goal, sample entropy, a measure of time-series regularity, was used to quantify the CoP regularity of four foot types.

### Methods

One hundred and sixty five subjects that had the same foot type bilaterally (48 subjects with healthy feet, 22 with pes valgus, 47 with hallux valgus, and 48 with pes cavus) were recruited for this study. A Footscan® system was used to collect CoP data when each subject walked at normal and steady speed. The velocity and acceleration in medial-lateral (ML) and anterior-posterior (AP) directions, and resultant velocity and acceleration were derived from CoP. The sample entropy is the negative natural logarithm of the conditional probability that a subseries of length m that matches pointwise within a tolerance r also matches at the next point. This was used to quantify variables of CoP velocity and acceleration of four foot types. The parameters *r* (the tolerance) and *m* (the matching length) for sample entropy calculation have been determined by an optimal method.

### Results

It has been found that in order to analyze all CoP parameters of velocity and acceleration during the stance phase of walking gait, for each variable there is a different optimal r value. On the contrary, the value m=4 is optimal for all variables.

Sample entropies of both velocity and acceleration in AP direction were highly correlated with their corresponding resultant variables for r>0.91. The sample entropy of the velocity in AP direction was moderately correlated with the one of the acceleration in the same direction (r≥0.673), as well as with the resultant acceleration (r≥0.660). The sample entropy of resultant velocity was moderately correlated with the one of the acceleration in AP direction, as well as with the resultant acceleration (for the both r≥0.689). Moderate correlations were found between variables for the left foot and their corresponding variables for the right foot.

Sample entropies of AP velocity, resultant velocity, AP acceleration, and resultant acceleration of the right foot as well as AP velocity and resultant velocity of the left foot were, respectively, significantly different among the four foot types.

### Conclusions

It can be concluded that the sample entropy of AP velocity (or the resultant velocity) of the left foot, ML velocity, resultant velocity, ML acceleration and resultant acceleration could serve for evaluation of foot types or selection of appropriate footwear.

## Keywords

## Background

Foot problems are often seen among people of all ages over the world. The situation worsens in case of society aging since the incidence of foot diseases grows as age increases [1]. The common foot types are normal foot, pes valgus, pes cavus and hallux valgus. Pes valgus (also known as pes planus, flat foot etc.) is characterized with a collapse of the medial longitudinal arch. Hallux valgus exhibits medial deviation of the first metatarsal. Pes cavus (also known as high arch, cavoid foot etc.) has a higher medial longitudinal arch. Each of these abnormalities leads to a higher incidence of specific complications. For example, people suffering from high arched feet are more likely to develop tibial and femoral stress fractures, and low-arched people are more likely to develop metatarsal injuries [2]; hallux valgus may lead to metatarsalgia, plantar callosities, hammer toe deformities; all these conditions worsen life quality [3, 4].

However, in the most countries the available podiatric resources are not enough to solve all foot problems [5]. This leads to a lack of early foot problem recognition. Furthermore, most of the people are not aware of their own foot types [6], and if unsuitable footwear is selected the foot condition will worsen the situation. For these reasons, it is necessary to investigate the characteristics of different foot types to automatically detect foot types.

In gait and posture research, plantar pressure measurement is important to determine the condition of foot and ankle. Such measurement is used in footwear design [7, 8], gait identification [9], gait alteration recognition [10], etc. It is also used to investigate the loading characteristics and classification of foot types.

Bavornrit et al. [11] found that, compared with the normal foot, the low arch foot presents an increased contact area beneath the medial midfoot and a decreased peak pressure beneath the lateral forefoot. Pes planus presents significantly more force at the subhallucal area with no difference seen under the other areas [12]. Significant difference between normal foot and pes cavus exists for pressure–time integrals [13]. The peak pressure under first metatarsal and second metatarsal region and the pressure–time integral at the first metatarsal region in hallux valgus group were higher than those in the control group. People with lesser toes deformity exhibited higher pressure–time integral and peak pressure under the 2nd and 3rd metatarsal joints and 2nd toe and 3rd-5th toes compared with people who had normal feet [14]. There are also differences between hallux valgus and normal foot in duration of the forefoot loading and medial-to-lateral transition and in balance pattern [15].

Features obtained from CoP were also used to analyze motion characteristics of different foot types. Ledoux et al. and Hillstrom et al. [12, 16] found that the pes planus has lower center of pressure excursion index (CPEI, a measure of CoP concavity in the metatarsal area divided by foot width during stance phase) in comparison with the neutrally aligned foot and cavus foot. The CoP displacement and velocity were investigated for high arch, normal and low arch feet, and the low arch foot exhibited a more lateral COP course [17].

Spatio-temporal variables of plantar pressure were also used to classify the foot types. For example, Bertani et al. [18] used a heuristic optimization (based on Discard-Insert-Exchange optimization method) to choose the most discriminative samples of GRF (ground reaction force) as features for classifying foot types into normal foot and flat foot. Xu et al. [19] proposed a neural network combined with fuzzy logic for classifying feet into normal foot, pes cavus and pes valgus based on four static characteristic indexes: Staheli index (the ratio of the minimum width of midfoot to the maximum width of rearfoot), Chippaux-Smirak index (the ratio of the minimum width of midfoot to maximum width of forefoot), arch index (the ratio of the area of the middle third of the toeless footprint to the entire footprint area) and modified arch index (it is similar to the Arch Index, but it refers the ratio of the sum of the pressures). De Cock et al. [20] selected as features peak pressures, regional impulses and relative regional impulses under hallux, medial and lateral heel and five metatarsal joints; then they used K-mean cluster algorithm to classify pressure patterns into four categories. The K-mean cluster algorithm is a partition-based cluster method; observations are partitioned into k clusters according to the distance between each observation and each centroid.

Until now, linear features are usually extracted from CoP of different foot types. To some degree, they can indicate the biomechanical characteristics of different foot types. However, the linear features which were always extracted from one or several important data points (e.g. maximal force, CPEI) are easily contaminated with noise. Averaging is always used to deal with multiple measurements of each subject; this leads to suppress the time-evolving structure of CoP variability that is regarded as deterministic property [21]. In fact, human gait is a rhythmical oscillation and the movements do not repeat precisely in every gait cycle, thus in every human there is a stride-to-stride variability. Examining this variability allows identifying pathologies of neuromusculoskeletal system.

It is possible for CoP progression patterns of the four foot types to be different from one to another because each of the four foot types has its typical anatomical structure. To investigate the CoP differences of the four foot types, sample entropy has been used in this study. Sample entropy is a complexity measure of time series without counter of self-matches. This was proposed by Richman and Moorman [22]. It is not sensitive to data corrupted by noise [23] and can be used to quantify the stride-to-stride fluctuation of CoP. Sample entropy has been used to analyze CoP data during standing condition [24, 25]. To the best of our knowledge, it has not been used to analyze CoP data during the stance phase, especially for the CoP velocity and acceleration of the four foot types. It has been hypothesized that there are differences in sample entropy values for CoP variables among the four foot types. The aim of this study is 1) to investigate whether different foot types could present difference in sample entropy for every CoP variable; 2) to study the relationship between CoP components; 3) to study whether sample entropy of left feet and the one of the right feet have the same statistical characteristics.

The rest of the paper is organized as follows: in Method Section the experimental platform and the way to acquire the data are described; then, sample entropy algorithm is introduced of which two important parameters m and r are determined. Statistical analysis methods are given in Statistical analysis Section. Results from comparison of entropy values of the four foot types are reported in Result Section. Discussion and conclusion are presented in Discussion Section.

## Method

### Experimental data acquisition

In total, 165 subjects (88 females, 77 males, foot size: 24.4±1.4 cm, weight: 58.7±15.4 kg, height: 162.9±9.5 cm, age: 32.6±14.2 years) were recruited in the study. Among them there were 22 subjects with pes valgus, 48 with pes cavus, 47 with hallux valgus and 48 with normal feet. Subjects who had pes valgus and hallux valgus present together on the same foot were excluded from the study. The both feet of each recruited subject were of the same type.

**Subjects’ demographic characteristics: allocation of the foot types among subjects, mean age and standard deviation of age, height, weight, and foot size**

Foot Type | Age, years Mean (std.dev.) | Females/males | Height, cm Mean (std.dev.) | Weight, kg Mean (std.dev.) | Foot size, cm Mean (std.dev.) |
---|---|---|---|---|---|

Normal foot | 36.5 (12.1) | 23/25 | 163.9 (8.3) | 60.4 (11.8) | 24.6 (1.4) |

Pes cavus | 32.1 (12.8) | 20/28 | 166.6 (8.7) | 62.4 (16.0) | 24.7 (1.5) |

Hallux valgus | 35.8 (15.3) | 37/10 | 160.5 (6.4) | 54.6 (10.3) | 23.9 (1.1) |

Pes valgus | 26.3 (16.3) | 8/14 | 157.8 (14.6) | 55.8 (25.0) | 24.4 (1.9) |

Whole group | 32.6 (14.2) | 88/77 | 162.9 (9.5) | 58.7 (15.4) | 24.4 (1.4) |

People tend to spend most of the time in their daily life walking at their normal or comfortable speed. Thus the biomechanical characteristics and those of the plantar pressure for normal walking speed are representative. For this reason plantar pressure data at normal walking speed and a mid-step protocol was used (i.e. the subject was required to make at least three steps before and after contacting the Footscan [26, 27]). A trial was considered valid when the following criteria were met: (1) Subjects were instructed to look forward and not to look at the pressure plate and the walkway. (2) Each subject has not suddenly changed his/her gait before and after accessing the plate i.e. there were no changes in the step length, the cadence, etc. (3) Plantar contact area was confined into the measurement area. (4) For each subject qualified plantar data was acquired six times for each foot.

### Sample entropy

*X*

_{ i },

*Y*

_{ i }and

*F*

_{ i }are the coordinates of CoP in ML and AP directions, and the force at a given time

*i*, respectively.

*Δt*represents the sampling interval. ${V}_{M{L}_{i}}$, ${V}_{A{P}_{i}}$, ${V}_{{R}_{i}}$ and ${V}_{{F}_{i}}$ are the velocities in ML and AP directions, the resultant velocity, and the force change rate, respectively. ${A}_{M{L}_{i}}$, ${A}_{A{P}_{i}}$ and ${A}_{{R}_{i}}$ represent the accelerations in ML and AP directions and the resultant acceleration in both directions. CoP data of all trials for the left foot and right foot of each subject were concatenated as whole time series. In Figure 2 and Figure 3 the time series for all variables related to velocity and acceleration are illustrated. The variables CoP velocity and acceleration for left foot and right foot are listed in Table 2. Before calculating the sample entropy, each time series was normalized to unit variance by dividing it by its respective standard deviation.

**The abbreviated variables**

The abbreviated variables | The meaning of CoP variables | |
---|---|---|

Left foot | L_ML_V | Medial-lateral velocity of CoP for left foot |

L_AP_V | Anterior-posterior velocity of CoP for left foot | |

L_Res_V | Resultant velocity of CoP for left foot | |

L_F_V | Force change rate of vertical ground reaction force for left foot | |

L_ML_A | Medial-lateral acceleration of CoP for left foot | |

L_AP_A | Anterior-posterior acceleration of CoP for left foot | |

L_Res_A | Resultant acceleration of CoP for left foot | |

Right foot | R_ML_V | Medial-lateral velocity of CoP for right foot |

R_AP_V | Anterior-posterior velocity of CoP for right foot | |

R_Res_V | Resultant velocity of CoP for right foot | |

R_F_V | Force change rate of ground vertical reaction force for right foot | |

R_ML_A | Medial-lateral acceleration of CoP for right foot | |

R_AP_A | Anterior-posterior acceleration of CoP for right foot | |

R_Res_A | Resultant acceleration of CoP for right foot |

Sample entropy is the negative natural logarithm of the conditional probability that a subseries of length m that matches pointwise within a tolerance r will also match at the next point [28]. Moreover, sample entropy derived from approximate entropy without counting self-matches. Lower sample entropy value indicates more self-similarity in the time series. In addition, for eliminating self-matches, sample entropy algorithm is simpler than approximate entropy algorithm, and its execution time is approximately the half of the one that is required by the approximate entropy algorithm. Sample entropy is largely independent of record length and demonstrates relative consistency under circumstances where approximate entropy does not [22, 29].

The algorithm to calculate sample entropy (SamEn) is as follows [28]:

Given a time sequence data *x*_{
1
}, *x*_{
2
}, *…*, *x*_{
n
}, n is the total number of data points. First, it is necessary to construct a subseries (template vector) of length m: *X*_{
i
}=[*x*_{
i
}*, x*_{
i+1
}*,…, x*_{
i+m-1
}], where i=1, 2, …, n-m+1; m represents embedding dimension.

*Xi*is calculated:

*n*

_{ i }

*(m, r)*is the number of vectors

*X*

_{ j }that are similar to

*X*

_{ i }with the following constrain:

where *d(X*_{
i
}*, X*_{
j
}*)* is defined as the maximal absolute difference between vectors *X*_{i} and *Y*_{j} in their respective scalar components; r specifies the filter level (tolerance). If the distance between *X*_{
i
} and *X*_{
j
} is less than r, the counter of vectors which are similar to *X*_{
i
} will increase by one. In the above counter of similar matches, *j=i* cases are not included to avoid self-match counting.

Similar process is repeated for an embedded dimension m+1 to calculate *Φ*(*m* + 1, *r*).

To calculate the sample entropy for each of the CoP measurement variables, the parameters *m* and *r* should be determined. As CoP data were acquired during the stance phase, *m* and *r* may be different from the recommended values (*m* can take value 1 or 2, *r* can take values between 0.1 and 0.25 [30]). The method proposed in [31] and [32] was used to determine the values of *m* and *r*.

If the number of matches of length m and m+1 increase, the accuracy and confidence of sample entropy estimate will increase. Small m values and large r values result into an increased number of matches. However, as r increases, the probability of match will tend to 1, thus the quantified sample estimate will lose the discriminative ability; as m decreases, underlying physical process may be obscured [32]. Therefore, m values and r values should be rationally chosen for reliable estimation of sample entropy.

*m*values, first the sample entropy of each time-series with combination of m from 1 to 6 with step 1, and r from 0.1 to 1 with step 0.05 has been calculated; then median sample entropy has been obtained for given m and r. The median sample entropy for each CoP variable is illustrated in Figure 4 (a)-(g). It can be seen that the median sample entropy of all variables converges when m≥4 for almost all r values.

*m*and

*r*values, conditional probability has been calculated:

*A(r)*and

*B(r)*are, respectively, the number of matches of length m+1 and m within tolerance r. The variance of CP can be estimated as:

*K*

_{ A }and

*K*

_{ B }are, respectively, the number of pairs of matching templates of length m+1 and m that overlap within tolerance r. The r values are determined by minimizing the quantity:

where *Q(m, r)* is the maximum relative error of SampEn and the CP estimate.

The metric simultaneously penalizes CP near 0 and 1 and it is a tradeoff between accuracy and discriminative capability. The maximal relative error is set less than 0.05; this value corresponds to 95% confidence interval, which is 10% sample entropy estimation.

The median between the maximum relative error for m≥4 and all r values for each CoP variable has been calculated. In Figure 4 (a’)-(g’), the medians of maximum relative error for m=4, 5, 6 are illustrated. Under the condition that the median of maximum relative error is set to 0.05, in order to reach the best discriminative ability of sample entropy, the values have been set as follows: for m a value equal to 4 has been chosen. On the contrary the r values of all variables for the left foot should be: ML velocity (L_ML_V) 0.16, AP velocity (L_AP_V) 0.19, resultant velocity (L_Res_V) 0.18, force change rate (L_F_V) 0.28, ML acceleration (L_ML_A) 0.19, AP acceleration (L_AP_A) 0.22, resultant acceleration (L_Res_A) 0.22; the corresponding variables for the right foot should have values of 0.17, 0.18, 0.17, 0.29, 0.18, 0.21and 0.20, respectively.

### Statistical analysis

To investigate the relationship between variables for left foot or right foot, Spearman Correlation analysis was performed because some variables did not follow normal distribution. If the coefficient |r|≥0.75, the two variables highly correlate; if 0.25≤|r|<0.75, the two variables moderately correlate; if |r|<0.25, the two variables weakly correlate [33]. The sample entropy of each variable was log-transformed to achieve normal distribution and equality of variance. Because of the inequality of variables R_ML_V, R_AP_V, L_ML_V and L_ML_A, Kruskal-Wallis test was applied to investigate whether sample entropy is statistically different for each of the three variables among the four foot types. Analysis of variance was performed with pairwise comparisons with Bonferroni adjustment for the rest variables (L_AP_V, L_Res_V, L_F_V, L_AP_A, L_Res_A, R_Res_V, R_F_V, R_ML_A, R_ AP_A, and R_Res_A).

With regard to the corresponding variables of left foot and right foot, if there is no high correlation between them and analysis of variance or Kruskal-Wallis test shows that there is a statistical difference between them, then data from the left foot and the right foot should be collected and processed separately. On the other hand, with regard to the CoP variables of the same-side foot, if there is a high correlation between variables and analysis of variance or Kruskal-Wallis test shows that they have similar statistical characteristics, then one of the variables cannot be used for classification of foot types or evaluation of foot function. All statistics were calculated with SPSS 17.0 and P<0.05 was taken as significant level.

## Result

### Relationships between variables

**Spearman correlation coefficients for sample entropy estimate between variables of the left foot**

L_ML_V | L_AP_V | L_Res_V | L_F_V | L_ML_A | L_AP_A | L_Res_A | ||
---|---|---|---|---|---|---|---|---|

L_ML_V | r | 1.000 | .278 | .271 | .191 | .662 | .192 | .196 |

p-value | .000 | .000 | .015 | .000 | .014 | .012 | ||

L_AP_V | r | .278 | 1.000 | .971 | .259 | .335 | .723 | .714 |

p-value | .000 | .000 | .001 | .000 | .000 | .000 | ||

L_Res_V | r | .271 | .971 | 1.000 | .246 | .348 | .742 | .746 |

p-value | .000 | .000 | .002 | .000 | .000 | .000 | ||

L_F_V | r | .191 | .259 | .246 | 1.000 | .243 | .329 | .324 |

p-value | .015 | .001 | .002 | .002 | .000 | .000 | ||

L_ML_A | r | .662 | .335 | .348 | .243 | 1.000 | .394 | .427 |

p-value | .000 | .000 | .000 | .002 | .000 | .000 | ||

L_AP_A | r | .192 | .723 | .742 | .329 | .394 | 1.000 | .985 |

p-value | .014 | .000 | .000 | .000 | .000 | .000 | ||

L_Res_A | r | .196 | .714 | .746 | .324 | .427 | .985 | 1.000 |

p-value | .012 | .000 | .000 | .000 | .000 | .000 |

**Spearman correlation coefficients for sample entropy estimate between variables of the right foot**

R_ML_V | R_AP_V | R_Res_V | R_F_V | R_ML_A | R_AP_A | R_Res_A | ||
---|---|---|---|---|---|---|---|---|

R_ML_V | r | 1.000 | .255 | .302 | .254 | .697 | .187 | .243 |

p-value | .001 | .000 | .001 | .000 | .017 | .002 | ||

R_AP_V | r | .255 | 1.000 | .963 | .314 | .229 | .673 | .660 |

p-value | .001 | .000 | .000 | .003 | .000 | .000 | ||

R_Res_V | r | .302 | .963 | 1.000 | .265 | .279 | .669 | .681 |

p-value | .000 | .000 | .001 | .000 | .000 | .000 | ||

R_F_V | r | .254 | .314 | .265 | 1.000 | .340 | .408 | .427 |

p-value | .001 | .000 | .001 | .000 | .000 | .000 | ||

R_ML_A | r | .697 | .229 | .279 | .340 | 1.000 | .422 | .511 |

p-value | .000 | .003 | .000 | .000 | .000 | .000 | ||

R_AP_A | r | .187 | .673 | .669 | .408 | .422 | 1.000 | .975 |

p-value | .017 | .000 | .000 | .000 | .000 | .000 | ||

R_Res_A | r | .243** | .660** | .681** | .427** | .511** | .975** | 1.000 |

p-value | .002 | .000 | .000 | .000 | .000 | .000 |

**Spearman correlation coefficients for sample entropy estimate between variables of left foot and right foot**

R_ML_V | R_AP_V | R_Res_V | R_F_V | R_ML_A | R_AP_A | R_Res_A | ||
---|---|---|---|---|---|---|---|---|

L_ML_V | r | .637 | .058 | .067 | .204 | .393 | .056 | .089 |

p-value | .000 | .458 | .399 | .009 | .000 | .476 | .256 | |

L_AP_V | r | .214 | .586 | .564 | .329 | .246 | .443 | .464 |

p-value | .006 | .000 | .000 | .000 | .002 | .000 | .000 | |

L_Res_V | r | .192 | .583 | .574 | .328 | .236 | .454 | .481 |

p-value | .014 | .000 | .000 | .000 | .002 | .000 | .000 | |

L_F_V | r | .163 | .227 | .193 | .738 | .320 | .330 | .348 |

p-value | .038 | .004 | .014 | .000 | .000 | .000 | .000 | |

L_ML_A | r | .389 | -.008 | -.007 | .292 | .453 | .182 | .223 |

p-value | .000 | .920 | .934 | .000 | .000 | .020 | .004 | |

L_AP_A | r | .097 | .398 | .405 | .396 | .260 | .601 | .611 |

p-value | .217 | .000 | .000 | .000 | .001 | .000 | .000 | |

L_Res_A | r | .097 | .381 | .393 | .387 | .267 | .587 | .601 |

p-value | .217 | .000 | .000 | .000 | .001 | .000 | .000 |

Moderate correlations were found between variables for left foot and corresponding variables for the right foot. With regard to the variables, AP velocity of the left foot (L_AP_V) was moderately correlated with the resultant velocity of the right foot (R_Res_V) (r=0.564, p=0.000); resultant velocity of the left foot (L_Res_V) was also moderately correlated with AP velocity of the right foot (R_AP_V) (r=0.583, p=0.000); AP acceleration of the left foot (L_AP_A) was moderately correlated with resultant acceleration of the right foot (R_Res_A) (r=0.611, p=0.000); resultant acceleration of the left foot (L_Res_A) was also moderately correlated with AP acceleration of the right foot (R_AP_A) (r=0.587, p=0.000).

### Comparisons of sample entropy for velocity and acceleration variables

**Kruskal-Wallis test for variables of velocity and acceleration for left foot and right foot**

Mean rank |
χ
| p-value | |||||
---|---|---|---|---|---|---|---|

Normal Foot | Pes Valgus | Hallux Valgus | Pes Cavus | ||||

Left foot | L_ML_V | 75.38 | 76.10 | 88.28 | 85.13 | 2.312 | .510 |

L_ML_A | 75.21 | 75.05 | 86.09 | 87.96 | 2.550 | .466 | |

Right foot | R_ML_V | 71.10 | 74.57 | 94.36 | 84.09 | 6.394 | .094 |

R_AP_V | 85.17 | 96.43 | 95.83 | 58.49 | 17.876 | .000 |

**Mean and standard deviations, and pairwise comparisons of sample entropy between normal foot, pes cavus, pes valgus, and hallux valgus**

Mean(STD) | |||||
---|---|---|---|---|---|

Normal Foot | Pes Valgus | Hallux Valgus | Pes Cavus | ||

Left foot | L_AP_V | −2.000 (.358) | −1.890 (.358) | −1.998(.324) | −2.180(.303) |

L_Res_V | −1.970(.363) | −1.866(.324) | −1.959(.330) | −2.151(.310) | |

L_F_V | −1.359(.267) | −1.426(.352) | −1.332(.310) | −1.429(.317) | |

L_AP_A | −1.914 (.790) | −2.050 (.941) | −1.890(.710) | −2.104(.635) | |

L_Res_A | −1.917 (.774) | −2.066 (.886) | −1.873 (.703) | −2.104(.634) | |

Right foot | R_Res_V | −1.968(.290) | −1.960(.418) | −1.878 (.341) | −2.149(.260) |

R_F_V | −1.439 (.280) | −1.452 (.289) | −1.408 (.306) | −1.504(.324) | |

R_ML_A | −1.750(.621) | −1.917(.686) | −1.650(.729) | −1.652(.816) | |

R_AP_A | −1.825(.652) | −2.180(.889) | −1.826(.681) | −2.149(.621) | |

R_Res_A | −1.781(.650) | −2.159(.832) | −1.790(.653) | −2.106(.639) |

## Discussion

This study was conducted to investigate whether there is a difference in sample entropy of each CoP variable among four foot types. Optimal values for m and r (which should be determined before calculation of the sample entropy) for every variable were found and they were different from those that were used in the analysis of the CoP sway data (m=3, r=0.3 [24]). Statistical differences were found for sample entropy of CoP velocity and acceleration.

### Foot structure and function of different foot types result in different CoP patterns

Above-mentioned statistical characteristics are based on the biomechanics of the different foot types during the stance phase. During this phase, the roll-off process of foot and ankle is accompanied with complex movements including dorsiflexion/plantarflexion on frontal plane, enversion/inversion on transverse plane and adduction/abduction on sagittal plane. Abnormity of foot anatomical structure (e.g. pes cavus, pes valgus, hallux valgus) will exhibit different patterns of the complex movement, and this will result into difference in CoP displacement, velocity and acceleration among the four foot types. For example, the forefoot of the pes planus was less adducted at toe-off in the transverse plane compared with normal foot [34]; at the ankle, the high-arched athletes exhibited significantly smaller peak eversion angles in walking compared to the low-arched athletes [35]. Therefore, complexity of CoP velocity and acceleration could indicate the functional difference of the foot and ankle among the four foot types.

In this study it has been found that sample entropy of ML acceleration and resultant acceleration were different between normal foot and pes valgus. Ledoux et al. [36] found that there is no difference in the acceleration at heel strike between normal foot (neutrally aligned foot) and pes valgus; however, in our study it has been found that sample entropy of AP and resultant acceleration were different between those of the two foot types. This may indicate that, compared with some of the linear methods, sample entropy has certain advantages in providing information about condition of foot and ankle.

### Similarities in sample entropy between velocity/acceleration variables

Spearman correlation indicated that two pairs of variables (AP velocity vs. resultant velocity, AP acceleration vs. resultant acceleration) of the left foot were highly correlated (p>0.96 for the both). Sample entropies for AP velocity and resultant velocity of the left foot exhibited similar statistical characteristics: pes cavus was different from any of the other three foot types. There was no statistical difference for sample entropy of AP acceleration (L_AP_A) and resultant acceleration (L_Res_A) among the four foot types. For the right foot, AP velocity (R_AP_V) and resultant velocity (R_Res_V) also demonstrated similar statistical characteristics. Sample entropies of AP acceleration (R_AP_A) and resultant acceleration (R_Res_A) were also similar with the only exception that the sample entropy of resultant acceleration (R_Res_A) for pes valgus and hallux valgus were different.

In conclusion, sample entropy for AP velocity and resultant velocity, AP acceleration and resultant acceleration had the same statistical characteristics. Actually, CoP displacement in AP direction is larger than the one in ML direction during the rollover process, so at the same time, velocity and acceleration in AP direction are larger than those in ML direction. In other words, ML velocity (and ML acceleration) components almost had no effect on the resultant velocity (and resultant acceleration). The above result seems similar with the one observed by De Cock et al. [17] that the curve of the resultant velocity resembles the curve of velocity in AP direction.

### Differences in sample entropy between right and left foot

There was no significant difference for sample entropy of AP acceleration (L_AP_A) and resultant acceleration (L_Res_A) of the left foot among the four foot types. However, for the right foot, statistically significant differences were found for pairwise comparisons of sample entropy of AP acceleration. It seems that CoP of the right foot can provide more discriminative information than the one of the left foot. The possible reason is that there is a functional asymmetry between left foot and right foot (the left foot is responsible for supporting the body while the right foot is responsible for propulsion during walking [37]). Since the right foot is in charge of propulsion, its sample entropy may exhibit more abnormalities of structure and function among the foot types. Considering the difference between statistical characteristics of the left foot and the right foot, data sets of left and right foot should be collected independently for sample entropy analysis of CoP velocity and acceleration.

### Limitations of the study

In this study, each subject walked many times in order to collect enough plantar pressure information to perform analysis. After this, the CoP time series of the different measurements of each subject were concatenated. If in-shoe plantar pressure device was used, it reduces the number of trails for each subject, and it is also of benefit for the subsequent data processing. In this study only four foot types have been explored. However, each foot type could be divided into sub-groups according to aetiology or symptom, and the pattern of plantar pressure of one sub-group may differ from those of the other sub-groups of the same foot type. For example, plantar pressures of idiopathic and neurogenic pes cavus are different [13], as well as hallux valgus with pain is different than the one without pain [15]. Merging all subgroups of one foot type into a common group may not lead to significant difference when comparing normal foot with a foot from another type. In future studies, more sub-groups will be considered. As a limitation of our method appears the dependence of the sample entropy on the sampling rate. It results from the fact that at higher sampling rates the self-match counter value will be higher. The optimal sampling rate should be determined for each particular system.

## Conclusion

- 1.
In order to analyse CoP data during the stance phase, values of m and r different from the recommended should be used.

- 2.
The statistical characteristics of the corresponding variables differ between the left and the right foot. For this reason CoP data for the left foot and the right foot should be collected and processed separately; in this way the study of CoP non-linear characteristics of different foot types is facilitated.

- 3.
Statistical difference was found for the following variables: sample entropies of AP velocity (R_AP_V), resultant velocity (R_Res_V), resultant acceleration (R_Res_A) of the right foot, sample entropies of AP velocity (L_AP_V) and resultant velocity (L_Res_V).

The results can be potentially used for evaluation of foot function, automatic classification of the foot type, or selection of footwear.

## Declarations

### Acknowledgements

This study was financed partially by the Projects of National Natural Science Foundation of China (Grant Nos. 60932001, 51105359 and 61072031), the National 863 Program of China (Grant No. 2012AA02A604), the National 973 Program of China (Grant No. 2010CB732606), the ‘Low-cost Healthcare’ Programs of Chinese Academy of Sciences, the Guangdong Innovation Research Team Fund for Low-cost Healthcare and International Science and Technology Cooperation Program of Guangdong Province (2012B050200004).

## Authors’ Affiliations

## References

- Burzykowski T, Molenberghs G, Abeck D, Haneke E, Hay R, Katsambas A, Roseeuw D, van de Kerkhof P, van Aelst R, Marynissen G: High prevalence of foot diseases in Europe: results of the Achilles Project.
*Mycoses*2003, 46: 496–505. 10.1046/j.0933-7407.2003.00933.xView ArticleGoogle Scholar - Simkin A, Leichter I, Giladi M, Stein M, Milgrom C: Combined effect of foot arch structure and an orthotic device on stress fractures.
*Foot Ankle Int*1989, 10(1):25–29. 10.1177/107110078901000105View ArticleGoogle Scholar - Saro C, Andrén B, Felländer-Tsai L, Lindgren U, Arndt A: Plantar pressure distribution and pain after distal osteotomy for hallux valgus.
*Foot*2007, 17: 84–93. 10.1016/j.foot.2006.11.002View ArticleGoogle Scholar - Lorei TJ, Kinast C, Klärner H, Rosenbaum D: Pedographic, clinical and functional outcome after scarf osteotomy.
*Clin Orthop Relat Res*2006, 451: 161–166.View ArticleGoogle Scholar *Future of podiatric medicine: an international perspective*. 2013. https://podiatry.com/ezines/index.php?action=viewPublication&nopopout=true&confirmOff=true&SearchText=&id=543&keepThis=true&TB_iframe=true&height=700&width=768&full=true- Hohmann E, Reaburn P, Imhoff A: Runner’s knowledge of their foot type: do they really know?
*Foot (Edinb)*2012, 22: 205–210. 10.1016/j.foot.2012.04.008View ArticleGoogle Scholar - Salles AS, Gyi DE: An evaluation of personalised insoles developed using additive manufacturing.
*J Sports Sci*2013, 31: 442–450. 10.1080/02640414.2012.736629View ArticleGoogle Scholar - Praet SF, Louwerens JW: The influence of shoe design on plantar pressure in Neuropathic Feet.
*Diabetes care*2003, 26: 441–445. 10.2337/diacare.26.2.441View ArticleGoogle Scholar - Pataky TC, Mu T, Bosch K, Rosenbaum D, Goulermas JY: Gait recognition: highly unique dynamic plantar pressure patterns among 104 individuals.
*J R Soc Interface*2012, 9: 790–800. 10.1098/rsif.2011.0430View ArticleGoogle Scholar - Giacomozzi C, Martelli F, Nagel A, Schmiegel A, Rosenbaum D: Cluster analysis to classify gait alterations in rheumatoid arthritis using peak pressure curves.
*Gait Posture*2009, 29: 220–224. 10.1016/j.gaitpost.2008.08.004View ArticleGoogle Scholar - Chuckpaiwong B, Nunley JA, Mall NA, Queen RM: The effect of foot type on in-shoe plantar pressure during walking and running.
*Gait Posture*2008, 28: 405–411. 10.1016/j.gaitpost.2008.01.012View ArticleGoogle Scholar - Ledoux WR, Hillstrom HJ: The distributed plantar vertical force of neutrally aligned and pes planus feet.
*Gait Posture*2002, 15: 1–9. 10.1016/S0966-6362(01)00165-5View ArticleGoogle Scholar - Burns J, Crosbie J, Hunt A, Ouvrier R: The effect of pes cavus on foot pain and plantar pressure.
*Clin Biomech (Bristol, Avon)*2005, 20: 877–882. 10.1016/j.clinbiomech.2005.03.006View ArticleGoogle Scholar - Mickle KJ, Munro BJ, Lord SR, Menz HB, Steele JR: Gait, balance and plantar pressures in older people with toe deformities.
*Gait Posture*2011, 34: 347–351. 10.1016/j.gaitpost.2011.05.023View ArticleGoogle Scholar - Wen J, Ding Q, Yu Z, Sun W, Wang Q, Wei K: Adaptive changes of foot pressure in hallux valgus patients.
*Gait Posture*2012, 36: 344–349. 10.1016/j.gaitpost.2012.03.030View ArticleGoogle Scholar - Hillstrom HJ, Song J, Kraszewski AP, Hafer JF, Mootanah R, Dufour AB, Chow BS, Deland JT: Foot type biomechanics part 1: Structure and function of the asymptomatic foot.
*Gait Posture*2013, 37(3):445–451. 10.1016/j.gaitpost.2012.09.007View ArticleGoogle Scholar - De Cock A, Vanrenterghem J, Willems T, Witvrouw E, De Clercq D: The trajectory of the centre of pressure during barefoot running as a potential measure for foot function.
*Gait Posture*2008, 27: 669–675. 10.1016/j.gaitpost.2007.08.013View ArticleGoogle Scholar - Bertani A, Cappello A, Benedetti MG, Simoncini L, Catani F: Flat foot functional evaluation using pattern recognition of ground reaction data.
*Clinical Biomechanics*1999, 14: 484–493. 10.1016/S0268-0033(98)90099-7View ArticleGoogle Scholar - Xu S, Zhou X, Sun YN: A novel gait analysis system based on adaptive neuro-fuzzy inference system.
*Expert Syst Appl*2010, 37: 1265–1269. 10.1016/j.eswa.2009.06.026View ArticleGoogle Scholar - De Cock A, Willems T, Witvrouw E, Vanrenterghem J, De Clercq D: A functional foot type classification with cluster analysis based on plantar pressure distribution during jogging.
*Gait Posture*2006, 23: 339–347. 10.1016/j.gaitpost.2005.04.011View ArticleGoogle Scholar - Cavanaugh JT, Guskiewicz KM, Stergiou N: A nonlinear dynamic approach for evaluating postural control, new directions for the management of sport-related cerebral concussion.
*Sports Med*2005, 35: 935–950. 10.2165/00007256-200535110-00002View ArticleGoogle Scholar - Richman JS, Joshua JR: Physiological time-series analysis using approximate entropy and sample entropy.
*Am J Physiol Heart Circ Physiol*2000, 278: 2039–2048.Google Scholar - Ramdani S, Bouchara F, Lagarde J: Influence of noise on the sample entropy algorithm.
*Chaos*2009, 19(1):031231–031236.MathSciNetView ArticleMATHGoogle Scholar - Ramdani S, Seigle B, Lagarde J, Bouchara F, Bernard PL: On the use of sample entropy to analyze human postural sway data.
*Med Eng Phys*2009, 31: 1023–1031. 10.1016/j.medengphy.2009.06.004View ArticleGoogle Scholar - Roerdink M, De Haart M, Daffertshofer A, Donker SF, Geurts AC, Beek PJ: Dynamical structure of center-of-pressure trajectories in patients recovering from stroke.
*Exp Brain Res*2006, 174: 256–269. 10.1007/s00221-006-0441-7View ArticleGoogle Scholar - Miller CA, Verstraete MC: A mechanical energy analysis of gait initiation.
*Gait Posture*1999, 9: 158–166. 10.1016/S0966-6362(99)00008-9View ArticleGoogle Scholar - Lindemann U, Najafi B, Zijlstra W, Hauer K, Muche R, Becker C, Aminian K: Distance to achieve steady state walking speed in frail elderly persons.
*Gait Posture*2008, 27: 91–96. 10.1016/j.gaitpost.2007.02.005View ArticleGoogle Scholar - Pham TD: GeoEntropy: A measure of complexity and similarity.
*Pattern Recogn*2010, 43: 887–896. 10.1016/j.patcog.2009.08.015View ArticleMATHGoogle Scholar - Yentes JM, Hunt N, Schmid KK, Kaipust JP, McGrath D, Stergiou N: The appropriate use of approximate entropy and sample entropy with short data sets.
*Ann Biomed Eng*2013, 41: 349–365. 10.1007/s10439-012-0668-3View ArticleGoogle Scholar - Aboy M, Cuesta-Frau D, Austin D, Mico-Tormos P:
*Characterization of Sample Entropy in the Context of Biomedical Signal Analysis*. Lyon, France: presented at the IEEE EMBS; 2007.View ArticleGoogle Scholar - Ramdani S, Seigle B, Varoqui D, Bouchara F, Blain H, Bernard PL: Characterizing the dynamics of postural sway in humans using smoothness and regularity measures.
*Ann Biomed Eng*2011, 39: 161–171. 10.1007/s10439-010-0137-9View ArticleGoogle Scholar - Lake DE, Richman JS, Griffin MP, Moorman JR: Sample entropy analysis of neonatal heart rate variability.
*Am J Physiol Regul Integr Comp Physiol*2002, 283: R789-R797.View ArticleGoogle Scholar *Pearson’s Correlation Coefficient*. http://www.statisticssolutions.com/academic-solutions/resources/directory-of-statistical-analyses/pearsons-correlation-coefficient/- Hunt AE, Smith RM: Mechanics and control of the flat versus normal foot during the stance phase of walking.
*Clin Biomech (Bristol, Avon)*2004, 19: 391–397. 10.1016/j.clinbiomech.2003.12.010View ArticleGoogle Scholar - Powell DW, Long B, Milner CE, Zhang S: Frontal plane multi-segment foot kinematics in high- and low-arched females during dynamic loading tasks.
*Hum Mov Sci*2011, 30: 105–114. 10.1016/j.humov.2010.08.015View ArticleGoogle Scholar - Ledoux WR, Willianm HJ: Acceleration of the calcaneus at heel strike in neutrally aligned and pes planus feet.
*Clinical Biomechanics*2001, 16: 608–613. 10.1016/S0268-0033(01)00041-9View ArticleGoogle Scholar - Sadeghi H, Allard P, Prince F, Labelle H: Symmetry and limb dominance in able-bodied gait: a review.
*Gait Posture*2000, 12: 34–45. 10.1016/S0966-6362(00)00070-9View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.