Next Article in Journal
Entropic Stabilization of Cas4 Protein SSO0001 Predicted with Popcoen
Next Article in Special Issue
Analog Circuit Fault Diagnosis via Joint Cross-Wavelet Singular Entropy and Parametric t-SNE
Previous Article in Journal
Enhancing of Self-Referenced Continuous-Variable Quantum Key Distribution with Virtual Photon Subtraction
Previous Article in Special Issue
Assessing Information Transmission in Data Transformations with the Channel Multivariate Entropy Triangle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis

1
Department of Mechanical Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
2
Department of Physical Therapy, College of Rehabilitation Sciences, University of Manitoba, Winnipeg, MB R3E 0T6, Canada
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(8), 579; https://doi.org/10.3390/e20080579
Submission received: 18 June 2018 / Revised: 27 July 2018 / Accepted: 2 August 2018 / Published: 6 August 2018

Abstract

:
Sample entropy (SampEn) has been used to quantify the regularity or predictability of human gait signals. There are studies on the appropriate use of this measure for inter-stride spatio-temporal gait variables. However, the sensitivity of this measure to preprocessing of the signal and to variant values of template size (m), tolerance size (r), and sampling rate has not been studied when applied to “whole” gait signals. Whole gait signals are the entire time series data obtained from force or inertial sensors. This study systematically investigates the sensitivity of SampEn of the center of pressure displacement in the mediolateral direction (ML COP-D) to variant parameter values and two pre-processing methods. These two methods are filtering the high-frequency components and resampling the signals to have the same average number of data points per stride. The discriminatory ability of SampEn is studied by comparing treadmill walk only (WO) to dual-task (DT) condition. The results suggest that SampEn maintains the directional difference between two walking conditions across variant parameter values, showing a significant increase from WO to DT condition, especially when signals are low-pass filtered. Moreover, when gait speed is different between test conditions, signals should be low-pass filtered and resampled to have the same average number of data points per stride.

1. Introduction

Entropy measures quantify the regularity or predictability of a time series [1,2,3,4,5]. Larger entropy values indicate less regularity or predictability in a time series. These measures have been used in gait analysis, and have been shown to discriminate between fallers and non-fallers [6], older and younger adults [7,8] and walk only (WO) and dual-task (DT) walking condition [7]. Various entropy measures have been proposed based on Shannon’s entropy [9] and its successor method, Approximate Entropy [10]. One of the most commonly methods used to study gait function is Sample Entropy (SampEn) [3].
Earlier application of SampEn was used to examine various inter-stride spatio-temporal gait variables, derived from endpoints of the gait cycle (heel strikes), for example, stride time and step length signals [11]. It has also been shown that there may be temporal scales in changes that occur in spatio-temporal gait variables [12]. However, these signals lack intra-stride information, which represents important passive and active gait control process [7]. Recent gait studies have examined “whole” gait signals (entire time series as opposed to stride-to-stride gait variables), such as trunk linear acceleration [1,6,13,14], joint angular positions [15] and center of pressure of feet displacement [7,16].
Center of pressure displacement in the mediolateral direction (ML COP-D) has been extensively used to examine the balance performance during standing conditions (base of support stationary) [17,18,19]. However, its usage in gait analysis is limited to a few studies [1,7,16,20]. This might be due to the limitation of the facilities that collect data on a treadmill or during overground walking. ML COP-D is capable of representing both single and double support phases of the gait cycle. It has also been shown that SampEn of ML COP-D better distinguishes between different treadmill walking conditions as compared to the SampEn of other commonly used signals (e.g., trunk linear acceleration) [7].
A few studies have investigated how variant template size (m), tolerance size (r), and data length (N) would affect SampEn of short [11] and long [21] time series when using inter-stride spatio-temporal gait variables. It was shown that SampEn values are dependent on the combination of m and r, and not on N [11,21]. However, no study has investigated the effect of parameter selection on SampEn of the human gait whole signals, such as ML COP-D or trunk linear acceleration, over an appropriate amount of continuous strides. During continuous, steady-state gait, these signals are similar in nature with a few dominant frequencies and have consistent fluctuations from stride to stride. Considering the increasing use of SampEn in analyzing human gait whole signals, it is essential to investigate how parameter selection would affect the outcomes. The importance of this investigation stems from the fact that parameter selection for calculating SampEn of whole gait signals, in many studies, is based on those that have analyzed inter-stride gait variables [8,13].
Most studies which have examined the effect of aging or dual-tasking on gait function, use self-paced walking and do not control for gait speed. Self-paced walking results in different walking speeds and, therefore, each walking condition will have a different average number of data points per stride. It has been shown that gait speed is significantly reduced during dual-task walking compared to walk only trials [22] and it has been reported that speed has a significant effect on measures of dynamical systems, such as the largest Lyapunov exponent [23,24,25]. Moreover, researchers have used different sampling frequencies when collecting target whole signals, which in turn have resulted in a different average number of data points per stride. It is unknown whether a different average number of data points per stride caused by varying walking speed or sampling rate would affect SampEn. Furthermore, many researchers have opted to apply SampEn, or other entropy measures, to raw unfiltered signals [8,16,26] to avoid losing or altering information due to filtering. While others have filtered the high-frequency components of trunk linear acceleration signal using a cut-off frequency of 20 Hz [14,27]. Therefore, investigating the effect of filtering would also be beneficial.
The first objective of this study is to systematically examine the sensitivity of SampEn of ML COP-D signals, obtained during treadmill walking, to variant m, r, and sampling rate values. The second objective is to determine the effect of the choice of low-pass filtering and data resampling, to have the same average number of data points per stride, on the SampEn of ML COP-D signals. Discriminatory ability of SampEn will be examined through comparing walk only condition to dual-task walking, which has been shown to adversely affect gait performance [28,29].

2. Materials and Methods

2.1. Experimental Procedure

A convenience sample of 29 healthy young participants (eight females, 28.3 ± 2.7 years, 173.4 ± 8.8 cm, 69.7 ± 14.2 kg, mean ± standard deviation (SD)) was recruited. They were screened to ensure that no participant had any illnesses, neuromuscular injuries or previous surgeries that might affect their balance and gait. The University of Manitoba Human Research Ethics Committee approved the study and all participants signed the informed consent form prior to the tests.
Participants were asked to walk on an instrumented Bertec treadmill (Bertec Corporation, Columbus, OH, USA) under three different walking conditions:
(a)
Walk only (WO) trial of 1 min at a speed of 1.0 m/s, and
(b)
Dual-task (DT) walking trial of 1 min at a speed of 1.0 m/s, which is described below, and
(c)
Walk only trial of 1 min at a speed of 1.3 m/s (WO-1.3).
Center of pressure displacements in the mediolateral (ML COP-D) and anteroposterior (AP COP-D) directions (Figure 1) were calculated from the force and moment components, which were sampled at 1000 Hz. Forty seconds of each signal, which contained at least 30 strides [16], were used after discarding approximately the first 4 strides.
During all walking trials, participants viewed an 80 cm computer monitor positioned 1 meter away at eye level. During the WO trials, participants watched a scenery video to maintain gaze and head position relative to the monitor. For the purpose of hands-free interaction with game activities, a commercial inertial-based wireless mouse (Elite mouse, SMK Electronics, Chula Vista, CA, USA) was mounted on a plastic headband worn by each participant. Therefore, during walking, the head rotation was used to control the motion of the computer cursor. The goal of the game was to move a game paddle horizontally to interact with moving game objects. The game objects were categorized as designated targets or designated distractors, with the shape of a soccer ball and dotted sphere, respectively. They appeared at random locations at the top of the display every 2 s and moved diagonally toward the bottom of the display. In response to each “game event” (target appearance), the participant produced a head rotation (i.e., rotation of the motion-sense mouse) to move the game paddle (left/right) to catch the target objects and avoid the distractors. For a full description of the interactive cognitive computer game, see Szturm et al. [29].

2.2. Sample Entropy

SampEn   ( m ,   r ,   N ) [3] of a dataset of length N is the negative natural logarithm of the conditional probability of two successive counts of similar pairs (Chebyshev distance less than a tolerance size of r) of template size m and m + 1 without allowing self-matches. SampEn is calculated as follows [7]; consider a time series of length N given below:
u = u ( 1 ) , u ( 2 ) , , u ( N )   or   u = { u ( j ) : 1 j N }
First, m value is chosen to construct series of pairs, size m as:
X m ( i ) = { u ( i + k ) : 0 k m 1 } ,   1 i N m + 1
Next, matching templates are found by comparing their chebyshev distance to a pre-determined r value while excluding self-comparison. Next, a variable called B i is built which is the number of pairs satisfying the aforementioned criteria:
B i m ( r ) = 1 N m 1   ( #   of   d | X m ( i ) X m ( j ) |   r ,   where   j = 1 : N m   &   i j )
d | X m ( i ) X m ( j ) | = max { | u ( i + k ) u ( j + k ) | : 0 k m 1 }
Next, B m ( r ) is defined as:
B m ( r ) = 1 N m i = 1 N m B i m ( r )
This process is repeated for m + 1 and r to form A m ( r ) :
A i m ( r ) = 1 N m 1   ( #   of   d | X m + 1 ( i ) X m + 1 ( j ) |   r ,   where   j = 1 : N m   &   i j )
A m ( r ) = 1 N m i = 1 N m A i m ( r )
Lastly, SampEn is calculated based on B m ( r ) and A m ( r ) as
SampEn   ( m , r , N ) =   ln A m ( r ) B m ( r )

2.3. Data Analysis

This study consists of two parts. In the first part, the sensitivity of SampEn to changing m, r, and sampling rate was investigated when comparing WO to DT. Two methods were used to downsample signals from 1000 Hz to lower sampling rates (Table 1). The goal was to downsample signals by factors of 1, 2, 4, 8, 16 and 32. The first method, decimation (D) by a factor of f, used an eighth-order low-pass Chebyshev Type I filter, which filtered the signal in forward and reverse directions to remove phase distortions and then selected every fth point (MATLAB command decimate). The filter had a normalized cut-off frequency of 0.8/f. This method was chosen to avoid aliasing distortion that might occur by simply downsampling a signal.
The second method, filtering-and-downsampling (FD) by a factor of f, used a second-order Butterworth low-pass filter with a cut-off frequency of 30 Hz, and then downsampled the signal by a factor of f (MATLAB command downsample). Butterworth low-pass filter is the most common filter used in the literature to reduce the effect of noise [30] along with maintaining the variability in the lower range frequencies where the musculoskeletal motion occurs [31]. A nonparametric PSD estimator, Welch’s algorithm, was used to obtain the cut-off frequency. The dominant peak was at 0.89 ± 0.06 Hz (mean ± SD) for WO, 0.91 ± 0.06 for DT, and 0.99 ± 0.06 for WO-1.3. The last peak before noise floor occurred in the 8–15 Hz frequency range. Therefore, 15 Hz was considered as the highest frequency component and 30 Hz was used as the cut-off frequency.
The two methods yielded approximately the same results with respect to the low-pass filtering for f = 32 . Therefore, the first five f values could shed light on the effect of low-pass filtering prior to the calculation of SampEn.
SampEn was calculated using all combinations of parameter values, m = 2 ,   4 ,   6 ,   8 ,   10 and r = 0.2 and 0.3 × standard deviation (SD) of all the time series, and for all downsampling factors f = 1 ,   2 ,   4 ,   8 ,   16 ,   32 , and for both decimated and filtered-and-downsampled signals of WO and DT condition. The present investigation was based on more m and f values in the selected ranges. However, the necessity for statistical analysis with the purpose of studying the discriminatory ability of SampEn, led to choosing fewer parameter values (levels within a factor); e.g., five levels versus nine levels for template size (m = 2~10). In a previous study [11], m = 2 ,   3 ,   4 were tested when SampEn was applied to the inter-stride spatio-temporal gait variables. The present work included more m values to study the SampEn of the entire gait signals and not just times at heel strike or step distances. It was hypothesized that larger m values could better discern changes when there is a much greater number of data points per gait cycle or stride. Additionally, unlike ApEn, SampEn decreases almost monotonically with increasing r value [3,11] and 0.1–0.3 times the standard deviation has been suggested for inter-stride spatio-temporal gait variables [11]. The current analysis was based on r = 0.1 × SD , r = 0.2 × SD and r = 0.3 × SD . However, when the parameter value r = 0.1 × SD was used, many SampEn values converged to infinity. Therefore this level was not included in the results. Large r values were not included because they result in much smaller SampEn values for each condition, i.e., more matched templates, which diminish the discriminatory ability of SampEn.
In the second part, the effect of low-pass filtering and resampling, to have the same average number of data points per stride, was investigated. SampEn of ML COP-D signal of WO, DT, and WO-1.3 was calculated using m = 4 ,   r = 0.2 × SD , and f = 8 (based on the results of the first part). Four methods of preprocessing were used for each condition;
  • decimation (D),
  • decimation-and-resampling (D-R),
  • filtering-and-downsampling (FD) and,
  • filtering-and-downsampling-and-resampling (FD-R).
The average number of data points per stride for WO, DT and WO-1.3 were 142, 140 and 128, respectively. Therefore, 30 strides of each time series were resampled (MATLAB command resample) so that all of the signals would have an average of 142 data points per stride.

2.4. Statistical Analysis

In the first part of this study, there were 4 factors of within-subject repeated measures, which are 2 levels of walking condition (WO and DT), 2 levels of r, 6 levels of f, and 5 levels of m. The following steps were taken to perform the statistical analysis separately for both decimated and filtered-and-downsampled signals.
  • A two-factor repeated measure ANOVA (walking condition*m) was performed at each f level while considering the first tolerance level.
  • A two-factor repeated measure ANOVA (walking condition*m) was performed at each f level while considering the second tolerance level.
  • A two-factor repeated measure ANOVA (walking condition*f) was performed at each m level while considering the first tolerance level.
  • A two-factor repeated measure ANOVA (walking condition*f) was performed at each m level while considering the second tolerance level.
  • Post hoc pairwise comparisons with Bonferroni correction were performed to examine the effect of dual-tasking at each level.
  • Finally, a two-factor (walking condition*r) repeated measure ANOVA was performed at fixed m = 4 and f = 8 values, which were chosen based on the previous step’s statistical results.
In the second part of this study, two two-factor within-subject ANOVA were used to examine the main and interaction effects of the following factors on SampEn;
  • walking condition (WO versus DT) and preprocessing method (D, D-R, FD, FD-R)
  • gait speed (1.0 m/s versus 1.3 m/s) and preprocessing method (D, D-R, FD, FD-R)
Normality of all dependent variables was checked using the Shapiro-Wilk normality test. Results confirmed that the data was normally distributed. Statistical analyses were carried out using SPSS software version 24. In all the tests, a p-value less than 0.05 was considered significant. A Bonferroni correction was used in the software for multiple comparisons.

3. Results

For the vast majority of combinations of parameter values, SampEn of ML COP-D during dual-task walking was significantly larger than that of walk only. In general, SampEn decreased as m increased, as r increased, and as f factor decreased, i.e., as sampling rate increased or as the number of points per stride increased. However, there were a few exceptions, which will be discussed further. The results of the main and interaction effects of walking condition (WO and DT), m, and f at each r value are presented in Table 2 and Table 3. In addition, the results of pairwise comparisons of the significant main effects of walking condition are presented in Table A1 and Table A2. The detailed results for each downsampling method are presented in the following subsections followed by the results of the effects of the preprocessing methods.

3.1. Sensitivity of SampEn to Variant Parameter Values When Using Filtering-and-Downsampling

Figure 2 and Figure 3 show the effect of changing f and m values on SampEn of the filtered-and-downsampled signals at r = 0.2 × SD . The figures are virtually the same as those of   r = 0.3 × SD . The statistical results are also reported in Table 2, Table 3, Table A1 and Table A2. A statistically significant interaction was found between walking condition and f at m = 2 and m = 4 for both tolerance values. Since the direction of the changes of SampEn with increasing f was increasing for both walking conditions, the interaction effect would signify a difference in the rate of the changes between levels. A statistically significant interaction was found between walking condition and m at all f values except for; (a) f = 4 and f = 32 for r = 0.2 × SD and (b) f = 4 and f = 8 for r = 0.3 × SD . At each f value, the direction of the changes of SampEn of WO, with respect to m, was similar to those of DT. The two-factor repeated measures of ANOVA of walking condition *r showed that there was no significant interaction between the walking condition and the r value (p = 0.813). However, there was a statistically significant decrease of SampEn with increasing the r value (p < 0.001). At each r value, SampEn of DT was statistically significantly larger than that of WO (p < 0.001).
For all m values, there was a statistically significant main effect of walking condition and f value on SampEn. Similarly, for all f values, there was a statistically significant main effect of walking condition and m value on SampEn. SampEn significantly increased from WO to DT for most combinations of r, m, and f values except for; (a) 4 out of 30 combinations of f and m values for r = 0.2 × SD and (b) 3 out of 30 combinations for r = 0.3 × SD . Nevertheless, for these exceptions, there was a trend of increased SampEn from WO to DT. Based on the statistical analysis, the increasing effect of dual-tasking on SampEn of ML COP-D signal can be captured for most combinations of r, m, and f values. The only exceptions are m = 10 at f = 16 and m = 6 ,   8 ,   10 at f = 32 .

3.2. Sensitivity of SampEn to Variant Parameter Values When Using Decimation

Figure 4 and Figure 5 show the effect of changing f and m values on SampEn of the decimated signal at r = 0.2 × SD . The figures are virtually the same as those of r = 0.3 × SD . The statistical results are also reported in Table 2, Table 3, Table A1 and Table A2. A statistically significant interaction was found between walking condition and f at m = 2 , m = 4 , and m = 6 for both tolerance values. The direction of the changes of SampEn with increasing f was increasing for both tasks. A statistically significant interaction was found between walking condition and m at all f values except for; (a) f = 4 , f = 8 and f = 32 for r = 0.2 × SD and (b) f = 8 for r = 0.3 × SD . At each f value, the direction of changes of SampEn of WO, with respect to m, was the same to those of DT. The two-factor repeated measure of ANOVA of walking condition-r showed that there was no significant interaction between the walking condition and the r value (p = 0.980). However, there was a statistically significant decrease of SampEn with increasing the r value (p = 0.001). At each r value, SampEn of DT was statistically significantly larger than that of WO (p < 0.001).
For all m values, there was a statistically significant main effect of walking condition (except for m = 10 ) and f values on SampEn. And for all f values, there was a statistically significant main effect of walking condition (except for f = 1 and f = 2 ) and m values on SampEn. SampEn significantly increased from WO to DT for most combinations of m, r, and f (4, 8, 16, and 32) values with a few exceptions; smaller m values at f = 4 and larger m values at f = 16 and f = 32 . Nevertheless, SampEn was seen to increase from WO to DT for those exceptions.

3.3. Effects of Preprocessing Methods

The descriptive and statistical results of the walking condition (WO versus DT), walking speed (1.0 m/s versus 1.3 m/s) and preprocessing method (D, D-R, FD, FD-R) on SampEn for the combination of f = 8 , m = 4 and r = 0.2 × SD are presented in Figure 6 and Table 4.
There was a significant interaction effect of walking speed and method, while no significant interaction of walking condition and method was found. In addition, all three walking condition, speed, and method had a significant main effect on SampEn. Results revealed that resampling signals to have a larger average number of data points per stride had a decreasing effect on SampEn of WO-1.3 signals. However, there was no significant effect of resampling on SampEn of DT and WO signals. The only exception was the significant decrease of SampEn of WO from FD to FD-R. In addition, SampEn significantly decreased by filtering the high-frequency components. Furthermore, SampEn increased significantly from WO to DT when using all four methods. Finally, there was a significant increase in SampEn with increasing walking speed (WO to WO-1.3) only when using decimation or decimation-and-resampling.

4. Discussion

The goal of this investigation was to identify the sensitivity of SampEn to variant values of parameters (m and r, and sampling rate) and two preprocessing methods when applied to ML COP-D signal obtained during treadmill walking under WO and DT conditions. There were three main observations in this study. First, SampEn of ML COP-D consistently increased significantly from WO to DT for most of the combinations of the parameter values and methods of pre-processing. This finding is in agreement with the previous studies [7,16,32]. The only exceptions were m = 10 at f = 16 and m = 6 ,   8 ,   10 at f = 32 for filtered-and-downsampled ML COP-D. And for decimated ML COP-D, all m values at f = 1 and f = 2 , and smaller m values at f = 4 and larger m values at f = 16 and f = 32 . Second, the results demonstrated that walking speed should be controlled when studying the effect of another factor, such as adding a cognitive task. Finally, using a low-pass filter to eliminate high-frequency signal components and smoothing a time series improved the consistency (a significant increase from WO to DT) of SampEn analysis.
The results of this study showed that SampEn of ML COP-D signal increases as participants perform a concurrent cognitive task while walking. Steady-state gait signals are intrinsically periodic with consistent fluctuation from stride to stride. It can be argued that when human gait is disrupted by performing a secondary task or negatively affected by aging or diseases, these fluctuations would change. The unplanned fluctuations might increase during walking under challenging conditions known to cause gait disturbances and stumbles [33]. Another study has reported that SampEn of trunk linear acceleration signal was smaller in older adults who had reported a fall versus non-fallers when tested during overground walking [6]. However, in that study, speed was not controlled among participants or between groups. Reducing one’s walking speed is a consistent strategy to manage threats to balance and when attending to concurrent cognitive tasks. In this regard, a main finding of the present study was that walking speed had a significant effect on SampEn of ML COP-D signal. One reason is an increase in vibrational noise of the treadmill-mounted force plate at higher walking speeds. Increased noise will cause an increase in SampEn. This problem can be solved by preprocessing the signal using a low-pass filter to eliminate high-frequency vibrational noise. Another issue when comparing signals collected at different speeds is that stride time will be reduced as speed increases, and therefore the number of data points per stride will be reduced. The present results of the effect of different f values on SampEn demonstrated that SampEn increased as the number of data points per stride decreased. One method to deal with this issue is to resample the signals so that the average number of data points per stride would be the same across different walking speeds. Based on these findings, it is important to control the walking speed when comparing SampEn of gait signal between WO and DT conditions; since people, especially older adults, slow down when they perform a secondary task [34,35].
The results of this study suggest that SampEn benefits from a relative consistency (a significant increase from WO to DT) across different combinations of the variant values of m, r, and f. For chaotic signals like Mackey-Glass system, which resemble periodic time series, entropy values decrease with increasing the m values [36]. In the present study, there was a decreasing trend of SampEn of ML COP-D signal with increasing m. However, there were exceptions and SampEn values plateaued only at m = 4 for f = 16 . Nevertheless, for most combinations, there was a statistically significant increase of SampEn from WO to DT. With respect to changes in template size at higher sampling rates, SampEn showed an increasing trend with increasing m. A possible reason for this behavior may be related to the strong periodicity of the signal. To overcome this, larger template sizes might be chosen for signals collected at higher sampling rates. A similar issue exists for lower sampling rates where SampEn values increased with increasing m, but only after a specific m value. This suggests that smaller m values should be chosen for lower sampling rates. The two tolerance values performed similarly since the smaller one was already larger than the noise level of the signal.
There are three limitations to this study. First, the sample size of the current study (29 participants) which led to performing several 2-factor repeated measure ANOVA instead of for example two 3-factor repeated measure ANOVA. A much larger sample size would be more appropriate to study the effect of four factors each with many levels. Second, this study examined the discriminatory ability of SampEn by only comparing WO to DT conditions. Several other factors, such as aging, should be considered to generalize the results of this study. Finally, only SampEn, which is a single-scale entropy measure, was studied. Multi-scale SampEn analysis [12] or modified SampEn analysis (incorporating a time delay greater than one) [37] would likely yield important findings.
For future studies using the ML COP-D signal, it is recommended to use a low-pass filter prior to the calculation of SampEn. In addition, a sampling rate of 125 Hz or 62 Hz with m = 2 ~ 6 and r = 0.2 × SD would be the preferred combinations. For studies testing overground walking where speed is difficult to control, an integer number of strides should be resampled so that the average number of data points per stride remains the same. For research investigations using other gait signals, such as trunk/pelvis linear acceleration, a similar approach should be performed to select the best combinations of m, r, and f values. This study was not designed to investigate the effect of data length on SampEn. Nevertheless, SampEn of whole signals plateaus after a few strides [7] and it has been shown that 30 strides are sufficient to calculate SampEn of whole gait signals [16].

Author Contributions

S.A., T.S., and C.W. conceived and designed this work; S.A. collected and analyzed the data; S.A. and T.S. interpreted the results; S.A. drafted the manuscript; S.A., N.S., T.S. and C.W. revised and finalized the manuscript. All authors read and approved the final manuscript.

Funding

This research was funded by the Natural Science and Engineering Research Council of Canada grant number RGPIN-314032, the Canadian Institute of Health Research grant number CEH-134155, and the University of Manitoba Graduate Fellowship UMGF-7765611.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Descriptive results (mean ± SD) and pairwise comparisons of “walking condition*m” at r = 0.2 × SD and each f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
Table A1. Descriptive results (mean ± SD) and pairwise comparisons of “walking condition*m” at r = 0.2 × SD and each f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
r = 0.2 × SD FDD
WODTp-ValueWODTp-Value
f = 1m = 20.014 ± 0.0020.015 ± 0.0020.0080.096 ± 0.0230.093 ± 0.020-
m = 40.014 ± 0.0020.015 ± 0.0020.0080.058 ± 0.0150.056 ± 0.012
m = 60.015 ± 0.0030.016 ± 0.0030.0080.034 ± 0.0070.033 ± 0.005
m = 80.015 ± 0.0030.016 ± 0.0030.0080.033 ± 0.0070.032 ± 0.005
m = 100.015 ± 0.0030.017 ± 0.0030.0040.032 ± 0.0070.032 ± 0.004
f = 2m = 20.028 ± 0.0050.030 ± 0.0050.0080.117 ± 0.0290.112 ± 0.024-
m = 40.030 ± 0.0060.032 ± 0.0050.0080.059 ± 0.0110.058 ± 0.008
m = 60.031 ± 0.0050.034 ± 0.0060.0020.056 ± 0.0110.057 ± 0.008
m = 80.030 ± 0.0050.033 ± 0.0040.0010.051 ± 0.0090.053 ± 0.006
m = 100.029 ± 0.0040.032 ± 0.004<0.0010.047 ± 0.0080.049 ± 0.005
f = 4m = 20.058 ± 0.0110.063 ± 0.0110.0080.073 ± 0.0140.074 ± 0.0120.568
m = 40.060 ± 0.0100.067 ± 0.0090.0010.068 ± 0.0110.072 ± 0.0100.042
m = 60.055 ± 0.0070.060 ± 0.006<0.0010.059 ± 0.0080.063 ± 0.0070.011
m = 80.045 ± 0.0060.049 ± 0.0060.0020.050 ± 0.0060.053 ± 0.0060.035
m = 100.037 ± 0.0040.040 ± 0.0050.0020.042 ± 0.0040.044 ± 0.0050.058
f = 8m = 20.121 ± 0.0200.134 ± 0.0210.0020.133 ± 0.0220.141 ± 0.0210.042
m = 40.095 ± 0.0120.103 ± 0.0110.0010.100 ± 0.0130.107 ± 0.0110.013
m = 60.068 ± 0.0070.073 ± 0.0080.0020.074 ± 0.0070.078 ± 0.0090.041
m = 80.061 ± 0.0060.064 ± 0.0070.0010.065 ± 0.0060.068 ± 0.0070.008
m = 100.059 ± 0.0060.063 ± 0.0070.0010.063 ± 0.0070.066 ± 0.0070.015
f = 16m = 20.208 ± 0.0250.228 ± 0.022<0.0010.209 ± 0.0260.228 ± 0.0220.001
m = 40.123 ± 0.0120.132 ± 0.0140.0010.128 ± 0.0120.136 ± 0.0150.002
m = 60.119 ± 0.0120.126 ± 0.0130.0010.123 ± 0.0120.128 ± 0.0130.007
m = 80.122 ± 0.0110.128 ± 0.0120.0100.126 ± 0.0120.130 ± 0.0120.049
m = 100.127 ± 0.0110.131 ± 0.0120.0520.130 ± 0.0120.133 ± 0.0120.181
f = 32m = 20.267 ± 0.0270.288 ± 0.030<0.0010.265 ± 0.0260.285 ± 0.0300.001
m = 40.243 ± 0.0230.254 ± 0.0240.0040.241 ± 0.0230.253 ± 0.0230.004
m = 60.262 ± 0.0230.270 ± 0.0250.0510.262 ± 0.0230.270 ± 0.0250.037
m = 80.295 ± 0.0370.302 ± 0.0360.1380.296 ± 0.0370.303 ± 0.0360.105
m = 100.332 ± 0.0520.340 ± 0.0450.2240.335 ± 0.0530.344 ± 0.0470.154
Table A2. Descriptive results (mean ± SD) and pairwise comparisons of “walking condition*m” at r   =   0.3 × SD and each f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
Table A2. Descriptive results (mean ± SD) and pairwise comparisons of “walking condition*m” at r   =   0.3 × SD and each f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
r   =   0.3 × SD FDD
WODTp-ValueWODTp-Value
f = 1m = 20.009 ± 0.0020.010 ± 0.0020.0120.058 ± 0.0140.057 ± 0.011-
m = 40.009 ± 0.0020.010 ± 0.0020.0120.034 ± 0.0080.032 ± 0.006
m = 60.010 ± 0.0020.010 ± 0.0020.0110.019 ± 0.0030.019 ± 0.002
m = 80.010 ± 0.0020.011 ± 0.0020.0110.018 ± 0.0030.018 ± 0.002
m = 100.010 ± 0.0020.011 ± 0.0020.0110.018 ± 0.0030.018 ± 0.002
f = 2m = 20.019 ± 0.0030.020 ± 0.0030.0120.069 ± 0.0160.067 ± 0.013-
m = 40.019 ± 0.0030.021 ± 0.0030.0110.034 ± 0.0060.034 ± 0.004
m = 60.020 ± 0.0040.022 ± 0.0040.0110.032 ± 0.0060.032 ± 0.004
m = 80.021 ± 0.0040.023 ± 0.0040.0050.031 ± 0.0060.032 ± 0.004
m = 100.021 ± 0.0030.023 ± 0.0030.0020.030 ± 0.0050.031 ± 0.004
f = 4m = 20.038 ± 0.0070.041 ± 0.0070.0110.047 ± 0.0090.048 ± 0.0070.578
m = 40.041 ± 0.0080.045 ± 0.0080.0060.046 ± 0.0090.048 ± 0.0080.130
m = 60.041 ± 0.0060.045 ± 0.0060.0010.044 ± 0.0070.047 ± 0.0070.018
m = 80.039 ± 0.0050.043 ± 0.0050.0010.041 ± 0.0060.044 ± 0.0050.007
m = 100.035 ± 0.0040.038 ± 0.0040.0010.037 ± 0.0050.040 ± 0.0040.008
f = 8m = 20.081 ± 0.0150.088 ± 0.0150.0080.089 ± 0.0170.093 ± 0.0160.127
m = 40.079 ± 0.0110.087 ± 0.0100.0010.081 ± 0.0120.088 ± 0.0110.005
m = 60.065 ± 0.0070.070 ± 0.0070.0010.068 ± 0.0080.072 ± 0.0070.010
m = 80.054 ± 0.0040.058 ± 0.005<0.0010.057 ± 0.0050.060 ± 0.0060.003
m = 100.050 ± 0.0040.053 ± 0.005<0.0010.053 ± 0.0040.055 ± 0.0050.002
f = 16m = 20.161 ± 0.0240.178 ± 0.0240.0010.163 ± 0.0250.178 ± 0.0240.002
m = 40.112 ± 0.0100.120 ± 0.011<0.0010.115 ± 0.0100.122 ± 0.0120.001
m = 60.100 ± 0.0080.105 ± 0.010<0.0010.102 ± 0.0090.107 ± 0.0100.002
m = 80.102 ± 0.0090.106 ± 0.0100.0040.104 ± 0.0100.108 ± 0.0110.014
m = 100.106 ± 0.0100.109 ± 0.0100.0420.108 ± 0.0100.111 ± 0.0110.102
f = 32m = 20.247 ± 0.0240.268 ± 0.024<0.0010.246 ± 0.0240.265 ± 0.024<0.001
m = 40.202 ± 0.0180.212 ± 0.0210.0010.201 ± 0.0180.211 ± 0.0200.001
m = 60.219 ± 0.0210.225 ± 0.0220.0690.218 ± 0.0210.224 ± 0.0220.066
m = 80.240 ± 0.0290.245 ± 0.0300.1430.241 ± 0.0300.246 ± 0.0300.156
m = 100.255 ± 0.0350.260 ± 0.0340.2200.257 ± 0.0360.263 ± 0.0350.211

References

  1. Kaptein, R.G.; Wezenberg, D.; IJmker, T.; Houdijk, H.; Beek, P.J.; Lamoth, C.J.; Daffertshofer, A. Shotgun approaches to gait analysis: Insights & limitations. J. Neuroeng. Rehabil. 2014, 11, 120. [Google Scholar] [PubMed]
  2. Pincus, S.M.; Goldberger, A.L. Physiological time-series analysis: What does regularity quantify? Am. J. Physiol. 1994, 266, H1643–H1656. [Google Scholar] [CrossRef] [PubMed]
  3. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [PubMed]
  4. Gow, B.; Peng, C.-K.; Wayne, P.; Ahn, A. Multiscale Entropy Analysis of Center-of-Pressure Dynamics in Human Postural Control: Methodological Considerations. Entropy 2015, 17, 7926–7947. [Google Scholar] [CrossRef] [Green Version]
  5. Lake, D.E.; Richman, J.S.; Griffin, M.P.; Moorman, J.R. Sample entropy analysis of neonatal heart rate variability. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2002, 283, R789–R797. [Google Scholar] [CrossRef] [PubMed]
  6. Ihlen, E.A.F.; Weiss, A.; Bourke, A.; Helbostad, J.L.; Hausdorff, J.M. The complexity of daily life walking in older adult community-dwelling fallers and non-fallers. J. Biomech. 2016, 1–9. [Google Scholar] [CrossRef] [PubMed]
  7. Ahmadi, S.; Wu, C.; Sepehri, N.; Kantikar, A.; Nankar, M.; Szturm, T. The Effects of Aging and Dual Tasking on Human Gait Complexity During Treadmill Walking: A Comparative Study Using Quantized Dynamical Entropy and Sample Entropy. J. Biomech. Eng. 2018, 140, 1–10. [Google Scholar] [CrossRef] [PubMed]
  8. Bisi, M.; Riva, F.; Stagni, R. Measures of gait stability: Performance on adults and toddlers at the beginning of independent walking. J. Neuroeng. Rehabil. 2014, 11, 131. [Google Scholar] [CrossRef] [PubMed]
  9. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef]
  10. Pincus, S. Approximate entropy (ApEn) as a complexity measure. Chaos 1995, 5, 110–117. [Google Scholar] [CrossRef] [PubMed]
  11. Yentes, J.M.; Hunt, N.; Schmid, K.K.; Kaipust, J.P.; McGrath, D.; Stergiou, N. The appropriate use of approximate entropy and sample entropy with short data sets. Ann. Biomed. Eng. 2013, 41, 349–365. [Google Scholar] [CrossRef] [PubMed]
  12. Costa, M.; Peng, C.K.; Goldberger, A.L.; Hausdorff, J.M. Multiscale entropy analysis of human gait dynamics. Phys. A Stat. Mech. Appl. 2003, 330, 53–60. [Google Scholar] [CrossRef]
  13. Van Schooten, K.S.; Pijnappels, M.; Rispens, S.M.; Elders, P.J.M.; Lips, P.; Daffertshofer, A.; Beek, P.J.; Van Dieën, J.H. Daily-life gait quality as predictor of falls in older people: A 1-year prospective cohort study. PLoS ONE 2016, 11, 1–13. [Google Scholar] [CrossRef] [PubMed]
  14. Lamoth, C.J.; van Deudekom, F.J.; van Campen, J.P.; Appels, B.A.; de Vries, O.J.; Pijnappels, M. Gait stability and variability measures show effects of impaired cognition and dual tasking in frail people. J. Neuroeng. Rehabil. 2011, 8, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kurz, M.J.; Hou, J.G. Levodopa influences the regularity of the ankle joint kinematics in individuals with Parkinson’s disease. J. Comput. Neurosci. 2010, 28, 131–136. [Google Scholar] [CrossRef] [PubMed]
  16. Leverick, G.; Szturm, T.; Wu, C.Q. Using Entropy Measures to Characterize Human Locomotion. J. Biomech. Eng. 2014, 136, 121002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Fino, P.C.; Mojdehi, A.R.; Adjerid, K.; Habibi, M.; Lockhart, T.E.; Ross, S.D. Comparing Postural Stability Entropy Analyses to Differentiate Fallers and Non-fallers. Ann. Biomed. Eng. 2016, 44, 1636–1645. [Google Scholar] [CrossRef] [PubMed]
  18. Rhea, C.K.; Kiefer, A.W.; Wright, W.G.; Raisbeck, L.D.; Haran, F.J. Interpretation of postural control may change due to data processing techniques. Gait Posture 2015, 41, 731–735. [Google Scholar] [CrossRef] [PubMed]
  19. Ramdani, S.; Seigle, B.; Lagarde, J.; Bouchara, F.; Bernard, P.L. On the use of sample entropy to analyze human postural sway data. Med. Eng. Phys. 2009, 31, 1023–1031. [Google Scholar] [CrossRef] [PubMed]
  20. Terrier, P.; Dériaz, O. Non-linear dynamics of human locomotion: Effects of rhythmic auditory cueing on local dynamic stability. Front. Physiol. 2013, 4, 1–13. [Google Scholar] [CrossRef] [PubMed]
  21. Yentes, J.M.; Denton, W.; Mccamley, J.; Ra, P.C.; Schmid, K.K. Effect of parameter selection on entropy calculation for long walking trials. Gait Posture 2018, 60, 128–134. [Google Scholar] [CrossRef] [PubMed]
  22. Howcroft, J.; Kofman, J.; Lemaire, E.D.; McIlroy, W.E. Analysis of Dual-Task Elderly Gait in Fallers and Non-Fallers using Wearable Sensors. J. Biomech. 2016, 49, 992–1001. [Google Scholar] [CrossRef] [PubMed]
  23. Bruijn, S.M.; van Dieën, J.H.; Meijer, O.G.; Beek, P.J. Is slow walking more stable? J. Biomech. 2009, 42, 1506–1512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kang, H.G.; Dingwell, J.B. Effects of walking speed, strength and range of motion on gait stability in healthy older adults. J. Biomech. 2008, 41, 2899–2905. [Google Scholar] [CrossRef] [PubMed]
  25. Rispens, S.M.; Dieën, J.H. Van; Schooten, K.S. Van; Lizama, L.E.C.; Daffertshofer, A.; Beek, P.J.; Pijnappels, M. Fall-related gait characteristics on the treadmill and in daily life. J. Neuroeng. Rehabil. 2016, 13, 12. [Google Scholar] [CrossRef] [PubMed]
  26. Katsavelis, D.; Mukherjee, M.; Decker, L.; Stergiou, N. Variability of lower extremity joint kinematics during backward walking in a virtual environment. Nonlinear Dyn. Psychol. Life Sci. 2010, 14, 165–178. [Google Scholar]
  27. Kavanagh, J.; Barrett, R.; Morrison, S. The role of the neck and trunk in facilitating head stability during walking. Exp. Brain Res. 2006, 172, 454–463. [Google Scholar] [CrossRef] [PubMed]
  28. Dingwell, J.B.; Robb, R.T.; Troy, K.L.; Grabiner, M.D. Effects of an attention demanding task on dynamic stability during treadmill walking. J. Neuroeng. Rehabil. 2008, 5, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Szturm, T.; Maharjan, P.; Marotta, J.J.; Shay, B.; Shrestha, S.; Sakhalkar, V. The interacting effect of cognitive and motor task demands on performance of gait, balance and cognition in young adults. Gait Posture 2013, 38, 596–602. [Google Scholar] [CrossRef] [PubMed]
  30. Sloot, L.H.; Houdijk, H.; Harlaar, J. Technical note A comprehensive protocol to test instrumented treadmills. Med. Eng. Phys. 2015, 37, 610–616. [Google Scholar] [CrossRef] [PubMed]
  31. England, S.A.; Granata, K.P. The influence of gait speed on local dynamic stability of walking. Gait Posture 2007, 25, 172–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Riva, F.; Toebes, M.J.P.; Pijnappels, M.; Stagni, R.; van Dieën, J.H. Estimating fall risk with inertial sensors using gait stability measures that do not require step detection. Gait Posture 2013, 38, 170–174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Vaillancourt, D.E.; Newell, K.M. Changing complexity in human behavior and physiology through aging and disease. Neurobiol. Aging 2002, 23, 1–11. [Google Scholar] [CrossRef]
  34. Dorfman, M.; Herman, T.; Brozgol, M.; Shema, S.; Weiss, A.; Hausdorff, J.; Mirelman, A. Dual-Task Training on a Treadmill to Improve Gait and Cognitive Function in Elderly Idiopathic Fallers. J. Neurol. Phys. Ther. 2014, 38, 246–253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lamoth, C.J.C.; Ainsworth, E.; Polomski, W.; Houdijk, H. Variability and stability analysis of walking of transfemoral amputees. Med. Eng. Phys. 2010, 32, 1009–1014. [Google Scholar] [CrossRef] [PubMed]
  36. Restrepo, J.F.; Schlotthauer, G.; Torres, M.E. Maximum approximate entropy and threshold: A new approach for regularity changes detection. Phys. A Stat. Mech. Appl. 2014, 409, 97–109. [Google Scholar] [CrossRef] [Green Version]
  37. Govindan, R.B.; Wilson, J.D.; Eswaran, H.; Lowery, C.L.; Preissl, H. Revisiting sample entropy analysis. Phys. A Stat. Mech. Appl. 2007, 376, 158–164. [Google Scholar] [CrossRef]
Figure 1. Trajectory of the center of pressure under WO, WO-1.3, and DT conditions: (A) The low-pass filtered trajectory of center of pressure displayed as AP COP-D vs. ML COP-D, (B) Filtered ML COP-D, (C) Several strides of unfiltered ML COP-D.
Figure 1. Trajectory of the center of pressure under WO, WO-1.3, and DT conditions: (A) The low-pass filtered trajectory of center of pressure displayed as AP COP-D vs. ML COP-D, (B) Filtered ML COP-D, (C) Several strides of unfiltered ML COP-D.
Entropy 20 00579 g001
Figure 2. Effect of changing f on SampEn of WO and DT at each m value at r = 0.2 × SD . for filtered-and-downsampled ML COP-D. Bottom-right: Effect of r on SampEn at f = 8 and m = 4 . The error bars reflect the standard error of the means.
Figure 2. Effect of changing f on SampEn of WO and DT at each m value at r = 0.2 × SD . for filtered-and-downsampled ML COP-D. Bottom-right: Effect of r on SampEn at f = 8 and m = 4 . The error bars reflect the standard error of the means.
Entropy 20 00579 g002
Figure 3. Effect of changing m on SampEn of WO and DT at each f value at r = 0.2 × SD for filtered-and-downsampled ML COP-D. The error bars reflect the standard error of the means.
Figure 3. Effect of changing m on SampEn of WO and DT at each f value at r = 0.2 × SD for filtered-and-downsampled ML COP-D. The error bars reflect the standard error of the means.
Entropy 20 00579 g003
Figure 4. Effect of changing f on SampEn of WO and DT at each m value at r = 0.2 × SD for decimated ML COP-D. Bottom-right: Effect of r on SampEn at f = 8 and m = 4 . The error bars reflect the standard error of the means.
Figure 4. Effect of changing f on SampEn of WO and DT at each m value at r = 0.2 × SD for decimated ML COP-D. Bottom-right: Effect of r on SampEn at f = 8 and m = 4 . The error bars reflect the standard error of the means.
Entropy 20 00579 g004
Figure 5. Effect of changing m on SampEn of WO and DT at each f value at r = 0.2 × SD for decimated ML COP-D. The error bars reflect the standard error of the means.
Figure 5. Effect of changing m on SampEn of WO and DT at each f value at r = 0.2 × SD for decimated ML COP-D. The error bars reflect the standard error of the means.
Entropy 20 00579 g005
Figure 6. Effects of different preprocessing methods on SampEn of the ML COP-D signal of WO, DT, and WO-1.3 conditions for the combination of f = 8 , m = 4 , and r = 0.2 × SD ; The error bars reflect the standard error of the means. D: Decimation; D-R: Decimation-and-Resampling; FD: Filtering-and-Downsampling; and FD-R: Filtering-and-Downsampling-and-Resampling.
Figure 6. Effects of different preprocessing methods on SampEn of the ML COP-D signal of WO, DT, and WO-1.3 conditions for the combination of f = 8 , m = 4 , and r = 0.2 × SD ; The error bars reflect the standard error of the means. D: Decimation; D-R: Decimation-and-Resampling; FD: Filtering-and-Downsampling; and FD-R: Filtering-and-Downsampling-and-Resampling.
Entropy 20 00579 g006
Table 1. Summary of downsampling factors (f), sampling rates and cut-off frequency for decimation and filtering-and-downsampling methods.
Table 1. Summary of downsampling factors (f), sampling rates and cut-off frequency for decimation and filtering-and-downsampling methods.
fSampling Rate (Hz)Cut-Off Frequency (Hz)
DecimationFiltering-and-Downsampling
1100080030
250040030
425020030
812510030
16625030
32312530
Table 2. Main and interaction effects (p-values) of “walking condition (W-C)*f” at each r and m value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
Table 2. Main and interaction effects (p-values) of “walking condition (W-C)*f” at each r and m value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
W-CfW-C*f W-CfW-C*f
FD
r   =   0.2 × SD
m = 2<0.001<0.001<0.001D
r   =   0.2 × SD
m = 20.009<0.001<0.001
m = 4<0.001<0.0010.006m = 40.009<0.001<0.001
m = 60.001<0.0010.123m = 60.021<0.0010.023
m = 80.004<0.0010.334m = 80.037<0.0010.070
m = 100.020<0.0010.421m =100.075<0.0010.211
W-CfW-C*f W-CfW-C*f
FD
r   =   0.3 × SD
m = 2<0.001<0.001<0.001D
r   =   0.3 × SD
m = 20.008<0.001<0.001
m = 4<0.001<0.001<0.001m = 40.001<0.001<0.001
m = 60.002<0.0010.080m = 60.015<0.0010.034
m = 80.006<0.0010.232m = 80.032<0.0010.154
m = 100.020<0.0010.387m = 100.068<0.0010.279
Table 3. Main and interaction effects (p-values) of “walking condition (W-C)*m” at each r and f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
Table 3. Main and interaction effects (p-values) of “walking condition (W-C)*m” at each r and f value for FD (Filtered-and-Downsampled) and D (Decimated) ML COP-D. The two conditions are WO and DT. p-values in bold indicate a significant difference.
W-CmW-C*m W-CmW-C*m
FD
r   =   0.2 × SD
f = 10.007<0.0010.001D
r   =   0.2 × SD
f = 10.094<0.0010.004
f = 20.002<0.0010.007f = 20.806<0.0010.000
f = 4<0.001<0.0010.128f = 40.049<0.0010.176
f = 8<0.001<0.0010.015f = 80.008<0.0010.145
f = 16<0.001<0.0010.002f = 160.001<0.0010.003
f = 320.007<0.0010.070f = 320.005<0.0010.114
W-CmW-C*m W-CmW-C*m
FD
r   =   0.3 × SD
f = 10.011<0.0010.010D
r   =   0.3 × SD
f = 10.329<0.0010.044
f = 20.006<0.0010.000f = 20.938<0.001<0.001
f = 40.001<0.0010.170f = 40.049<0.0010.020
f = 8<0.001<0.0010.062f = 80.006<0.0010.178
f = 16<0.001<0.0010.002f = 160.001<0.0010.004
f = 320.005<0.0010.006f = 320.005<0.0010.010
Table 4. Statistical results of the walking condition (WO and DT), walking speed (1.0 m/s and 1.3 m/s) and preprocessing method (D, D-R, FD, and FD-R) on SampEn of ML COP-D for the combination of f = 8 , m = 4 , and r = 0.2 × SD ; The top section presents main and interaction effects. The middle section presents the pairwise comparisons between preprocessing methods for each walking condition (WO, WO-1.3, and DT). The bottom section presents the pairwise comparisons between walking conditions for each preprocessing method. p-values in bold indicate a significant difference. D: Decimation; D-R: Decimation-and-Resampling; FD: Filtering-and-Downsampling; and FD-R: Filtering-and-Downsampling-and-Resampling.
Table 4. Statistical results of the walking condition (WO and DT), walking speed (1.0 m/s and 1.3 m/s) and preprocessing method (D, D-R, FD, and FD-R) on SampEn of ML COP-D for the combination of f = 8 , m = 4 , and r = 0.2 × SD ; The top section presents main and interaction effects. The middle section presents the pairwise comparisons between preprocessing methods for each walking condition (WO, WO-1.3, and DT). The bottom section presents the pairwise comparisons between walking conditions for each preprocessing method. p-values in bold indicate a significant difference. D: Decimation; D-R: Decimation-and-Resampling; FD: Filtering-and-Downsampling; and FD-R: Filtering-and-Downsampling-and-Resampling.
Main and Interaction Effects (p-Value)
Condition/SpeedMethodInteraction
WO vs. WO-1.30.017<0.001<0.001
WO vs. DT0.002<0.0010.057
Pairwise Comparisons (p-Value)
D vs. D-RD vs. FDFD vs. FD-R
WO0.104<0.0010.042
WO-1.3<0.001<0.001<0.001
DT0.981<0.0011.000
MethodWO vs. WO-1.3WO vs. DT
D<0.0010.013
D-R0.0010.006
FD0.7010.001
FD-R0.225<0.001

Share and Cite

MDPI and ACS Style

Ahmadi, S.; Sepehri, N.; Wu, C.; Szturm, T. Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis. Entropy 2018, 20, 579. https://doi.org/10.3390/e20080579

AMA Style

Ahmadi S, Sepehri N, Wu C, Szturm T. Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis. Entropy. 2018; 20(8):579. https://doi.org/10.3390/e20080579

Chicago/Turabian Style

Ahmadi, Samira, Nariman Sepehri, Christine Wu, and Tony Szturm. 2018. "Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis" Entropy 20, no. 8: 579. https://doi.org/10.3390/e20080579

APA Style

Ahmadi, S., Sepehri, N., Wu, C., & Szturm, T. (2018). Sample Entropy of Human Gait Center of Pressure Displacement: A Systematic Methodological Analysis. Entropy, 20(8), 579. https://doi.org/10.3390/e20080579

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop