Skip to main content

Regularity of kinematic data between single and dual-task treadmill walking in people with Parkinson’s disease



Regularity, quantified by sample entropy (SampEn), has been extensively used as a gait stability measure. Yet, there is no consensus on the calculation process and variant approaches, e.g. single-scale SampEn with and without incorporating a time delay greater than one, multiscale SampEn, and complexity index, have been used to calculate the regularity of kinematic or kinetic signals. The aim of the present study was to test the discriminatory performance of the abovementioned approaches during single and dual-task walking in people with Parkinson’s disease (PD).


Seventeen individuals with PD were included in this study. Participants completed two walking trials that included single and dual-task conditions. The secondary task was word searching with twelve words randomly appearing in the participants’ visual field. Trunk linear acceleration at sternum level, linear acceleration of the center of gravity, and angular velocity of feet, shanks, and thighs, each in three planes of motion were collected. The regularity of signals was computed using approaches mentioned above for single and dual-task conditions.


Incorporating a time delay greater than one and considering multiple scales helped better distinguish between single and dual-task walking. For all signals, the complexity index, defined as the summary of multiscale SampEn analysis, was the most efficient discriminatory index between single-task walking and dual-tasking in people with Parkinson's disease. Specifically, the complexity index of the trunk linear acceleration of the center of gravity distinguished between the two walking conditions in all three planes of motion.


The significant results observed across the 24 signals studied in this study are illustrative examples of the complexity index’s potential as a gait feature for classifying different walking conditions.


Researchers have proposed and studied a plethora of measures to quantitatively analyze walking [1]. Among them, entropy measures, such as Sample Entropy (SampEn), have shown great potential to successfully discern interrupted gait from uninterrupted [2]. Entropy measures represent regularity or predictability of signals and have been shown to distinguish between fallers and non-fallers [3], single and dual-task walking conditions [2], patients with various neurological and/or physical impairments and control groups [4]. Although first calculated for discrete stride-interval time series [5], entropy measures are increasingly applied to continuous human signals [3, 4, 6,7,8,9]. Trunk linear acceleration [4, 8, 10,11,12,13], lower extremity’s angular velocity [10, 13], and center of foot pressure [6] have been used with various entropy measures. Unlike inter-stride signals, such as step time, kinematic signals offer important information on motor control strategies as they comprise intra-stride details of gait cycles [6]. A previous study suggests that the within-stride phases, such as  the weight transfer phase, are critical for gait stability [14]. These phases are actively engaged with intra-stride patterns of cortical activity [15] which are involved in steady-speed human locomotion. [15, 16].

SampEn, among different types of entropy measures, is the most commonly used measure for human gait analysis. This single-scale regularity measure is representative of difficultness of describing the patterns of the trajectory of a system [17]. Single-scale analysis, however, is not capable of representing the multiscale richness of information embedded in a signal. Therefore, multiscale sample entropy (MSE) [5], which quantifies regularity on multiple scales, was introduced. One limitation of MSE is that it is reported as multiple SampEn values over a range of scales and does not produce a single value. This limitation hinders its application in gait analysis, especially when single quantitative features are rather desirable. Complexity index (CI) [18] calculated as the sum of SampEn values over a range of scales was introduced as a solution. Both single-scale and multiscale SampEn have shown great potential to capture the increased risk of fall despite not having incorporated a time-delay greater than one [2, 8]. For time series generated by nonlinear dynamics that have a long range correlation, it has been shown that using a higher time delay can provide additional information regarding the characteristics of a time series when conducting comparative analyses [19, 20]. In finite noisy time-series data, an appropriate choice of time delay [21] is necessary to properly reconstruct the state space as a too small time delay would produce a compressed attractor along the identity line [22]. In other words, the time delay is chosen such that the components of an embedding window or a template (which are the building blocks of nonlinear measures) are independent [23]. Incorporating a non-unity time delay has been used to calculate other gait stability measures, such as largest Lyapunov exponents, when applied to continuous kinematic gait signals [24, 25]. Nonetheless, the influence of incorporating a non-unity time delay, i.e. a time delay greater than one, when studying continuous signals with entropy measures is yet to be determined. The three approaches mentioned earlier, i.e. single-scale SampEn with and without a time delay greater than one, multiscale SampEn, and complexity index when calculated for continuous signals could somehow yield different results.

As Parkinson’s disease (PD) is the second most common neurodegenerative disease worldwide, after Alzheimer’s Disease, and has an exacerbated fall risk during walking, PD provides an opportunity to examine the efficacy of the aforementioned entropy metrics [26,27,28]. Epidemiological evidence demonstrates that the incidence of PD ranges between 10 and 50 per 100,000 individuals and has a prevalence between 100 and 300 per 100,000 individuals [29, 30]. Additionally, the number of individuals diagnosed with PD is expected to double worldwide by 2030 [31]. PD is a progressive neurodegeneration of dopaminergic neurons within the Basal Ganglia [26, 27]. In addition to causing PD’s cardinal symptoms (bradykinesia, rigidity, and tremor), the Basal Ganglia’s reduced functionality impairs automatic motor processes that contribute to effective locomotion [1, 32,33,34]. To overcome the impaired automaticity of the Basal Ganglia, previous research proposes that individuals with PD attempt to recruit relatively intact higher-level cortical structures (pre-frontal and frontal cortices) to direct conscious attentional control on their locomotion [32, 35]. However, when attention is divided between simultaneous tasks, the attentional resources required for effective compensation of impaired neuromuscular pathways become strained [32, 35]. For instance, dual-tasking evidence demonstrates that during the presence of a concurrent mental task, there is a breakdown in the gait pattern for individuals with PD which increases their fall risk [32, 33, 35, 36]. Hillel et al. discussed that changes in PD gait performance during dual-tasks are more representative of how these individuals ambulate in their everyday lives where multitasking is ever present in constantly changing environments [37]. To-date, however, it remains unexamined whether incorporating a non-unity time delay and multiscale analysis on SampEn results affects the sensitivity of entropy to respond to dual-tasks.

The main objective of the current study was to investigate the effect of incorporating a non-unity time delay and multiscale analysis on SampEn results. It was hypothesized that these approaches would increase the ability of SampEn to discriminate between single and dual-task trials. In addition, another objective of this paper was to determine whether the SampEn of eight different signals, each in three directions, would result in the same outcome. These eight different signals are trunk linear acceleration at sternum level, linear acceleration of the center of gravity, and angular velocity of least-affected and most-affected feet, shanks, and thighs. It was hypothesized that the discriminatory ability of SampEn would be greater when calculated for signals in the mediolateral direction based on previous studies showing an increased demand on mediolateral stability during visual perturbations [38,39,40].



Twenty (13 males; 7 females) individuals with idiopathic PD (confirmed by a neurologist), aged between 48 and 79 years old (age: 63.8 ± 9.0 years, mass: 72.3 kg ± 19.4, height: 172.8 cm ± 8.0) were recruited from the Ottawa-Gatineau region. Due to incomplete data collection, only 17 participants (11 males and 6 females, 64.8 ± 7.4) were used in the analysis. Participants were assessed with the original Unified Parkinson’s Disease Rating Scale Motor examination (11 ± 6) and were between I-III on the Hoehn and Yahr scale. Average disease duration (8.0 ± 5.1 years) and age at onset (56.8 ± 9.6 years) data were collected. Further, seven participants reported freezing of gait based on the Freezing of Gait Questionnaire. Participants were tested on their optimally medicated state. Exclusion criteria encompassed if participants reported any physical discomfort using a virtual reality system, reported any injuries and/or orthopedic surgeries that could interfere with walking, could only ambulate with the assistance of a walking aid and had any additional conditions other than PD. All participants provided written informed consent and the study was approved by the University of Ottawa Research Ethics Boards (REB) and Ottawa Hospital Research Institute (OHRI) in accordance with the Helsinki protocol.


Participants completed two treadmill walking trials, where they walked on a straight path (i.e. the treadmill belt) with their preferred walking speed, that included single and dual-task conditions. These two trials were part of a larger protocol [41, 42]. The single-task trial had a duration of three minutes while the dual-task trial lasted two minutes. Prior to data collection, participants walked for approximately 20 s to achieve steady-state walking. The dual-task consisted of a word searching task with twelve words randomly appearing in the participants’ visual field. Words were projected one at a time on a screen directly in front of the participants. Words appeared randomly on either the left or right side of the screen at angles varying between 20 and 70 degrees. The duration of each word was 3 s with a 2–4 s pause between following words. Participants verbally called out the words as they appeared [41]. The word searching task was chosen as a more ecologically representation of a potential secondary task that individuals would encounter in their daily lives when walking. Specifically, this task reflects scanning the environment for relevant visual cues such as when passing by signs which provide direction.

Three-dimensional motion analysis was completed using the CAREN-Extended System (Motek Medical, Amsterdam NL). This system combines a six degrees-of-freedom motion platform with an embedded dual-belt treadmill, 12 camera Vicon motion capture system, 180-degree projector screen, and a safety harness. A 57-marker set for tracking full body kinematics was used. Kinematic data were collected at 100 Hz.

Regularity measures

We briefly describe SampEn and define its single-scale and multiscale implementation with and without a non-unity time delay. Next, complexity index as a summarizing index of MSE is explained.

Consider a time series of length \(N\) given below:

$$X=\left\{x\left(1\right),x\left(2\right),\dots ,x\left(N\right)\right\} or X=\left\{x\left(j\right):1\le j\le N\right\}$$

The value for template size (\(m\)) is chosen to construct series of pairs, size \(m\) as:

$${X}_{m}\left(i\right)=\left\{x\left(i+k\right):0\le k\le m-1\right\}, 1\le i\le N-m+1$$

Next, matching templates are found by comparing their Chebyshev distance (denoted as \(\mathrm{d}\left|.\right|\)) to a pre-determined threshold size (\(r\)) while excluding self-comparison. Next, the variable \({B}_{i}\), which is the number of pairs satisfying the aforementioned criteria, is built.

$${B}_{i}^{m}\left(\mathrm{r}\right)=\frac{1}{N-m-1} \left(\#\mathrm{ of d}\left|{X}_{m}\left(i\right)-{X}_{m}\left(j\right)\right| \le r, where j=1:N-m \& i\ne j\right)$$
$$\mathrm{d}\left|{X}_{m}\left(i\right)-{X}_{m}\left(j\right)\right|=\mathrm{max}\left\{\left|x\left(i+k\right)-x\left(j+k\right)\right|:0\le k\le m-1\right\}$$

Next, \({B}^{m}\left(r\right)\) is defined as:

$${B}^{m}\left(r\right)=\frac{1}{N-m}{\sum }_{i=1}^{N-m}{B}_{i}^{m}\left(r\right)$$

This process is repeated for \(m+1\) and \(r\) to form \({A}^{m}\left(r\right)\):

$${A}_{i}^{m}\left(\mathrm{r}\right)=\frac{1}{N-m-1} (\#\mathrm{ of d}\left|{X}_{m+1}\left(i\right)-{X}_{m+1}\left(j\right)\right|\le r), \mathrm{where} j=1:N-m \& i\ne j$$
$${A}^{m}\left(r\right)=\frac{1}{N-m}{\sum }_{i=1}^{N-m}{A}_{i}^{m}\left(r\right)$$

Lastly, single-scale SampEn [43] is calculated based on \({B}^{m}\left(r\right)\) and \({A}^{m}\left(r\right)\) as

$$\mathrm{SampEn}\left(m,r,N\right)= -\mathrm{ln}\frac{{A}^{m}\left(r\right)}{{B}^{m}\left(r\right)}$$

where \(m\), \(r\) and \(N\) are the template size (i.e., the length of template vector), tolerance size and the length of time series, respectively [6].

For MSE, consecutive coarse-grained time series are constructed by averaging a successively increasing number of data points in non-overlapping windows. The number of data points in each window is determined by the scale factor which ranges from 1 to 20 in this study. MSE is usually presented as SampEn for each one of the coarse-grained time series plotted as a function of scale factor. And the complexity index is defined as the sum of the SampEn values of different scales [18].

In order to incorporate a non-unity time delay, a time delay \((\tau )\), which has been calculated using the minimum average mutual information method [22], is used when constructing the templates in [2] as below while the rest of the calculation process remain the same.

$${X}_{m}\left(i\right)=\left\{x\left(i+k\tau \right):0\le k\le m-1\right\}, 1\le i\le N-m+1$$

Data processing

Markers data were processed in Vicon Nexus (Nexus 2.6, Oxford, UK) and 3D kinematic calculations were performed in Visual 3D. A 4th order, low-pass Butterworth filter with a 10 Hz cut-off frequency was used to filter marker data. Visual 3D calculates the center of gravity based on the anthropometric data provided by Hanavan [44]. Center of gravity velocity and acceleration were calculated as the first and second derivative of the center of gravity’s position. For this study, only trunk linear acceleration at sternum level (SternumLA), center of gravity linear acceleration (COG-LA), foot angular velocity (FootAV), shank angular velocity (ShankAV), and thigh angular velocity (ThighAV), each in mediolateral (ML), anteroposterior (AP), and vertical (V) directions were exported from Visual 3D (see Fig. 1). Subsequently, the exported data were analyzed in MATLAB R2019b. For angular velocity signals, results are reported for the most affected (MA) and the least affected (LA) foot, shank, and thigh instead of left and right categorization.

Fig. 1
figure 1

A 7-s segment of biological data for continuous foot angular velocity (FootAV), shank angular velocity (ShankAV), thigh angular velocity (ThighAV), COG linear acceleration (COG-LA), and sternum linear acceleration (SternumLA) in ML (left), AP (middle) and V (right) directions

70 strides of each signal were resampled to have the same \(108\times 70\) data points, where 108 was the average number of data points per stride for 17 participants. This data manipulation was recommended to account for variant gait speed among participants and across trials [6]. Next, SampEn was calculated for each signal using a template size of 4, a tolerance size of 0.2 times the standard deviation of all time-series, and a time delay of one. Similarly, SampEn was recalculated for all signals using a template size of 4, a tolerance size of 0.2 times the standard deviation of all time-series, and a specified non-unity time delay. A systematic parameter selection [6] was followed to select template size and tolerance size, and the minimum average mutual information method was used to calculate the time delay [22]. A range of time delay values was obtained for each signal when calculated for different participants [45,46,47], and the median values were selected for the analysis (see Table 1). The dominant frequency of intrinsically periodic human gait signals is different across body segments, and it affects the calculated time delay [47]. Signals with less stride-to-stride changes have larger time delay values; i.e. fewer data points are needed to accurately reconstruct the state-space using a single time-series data [22].

Table 1 Median time delay in mediolateral, anteroposterior, and vertical directions

Subsequently, consecutive coarse-grained signals were constructed by averaging a successively increasing number of data points (up to 20) in non-overlapping windows. Next, SampEn was calculated for each coarse-grained signal using a template size of 4, a tolerance size of 0.2 times the standard deviation of all time-series of the same body segment, and a specified time delay greater than one. Time delay at each scale was the quotient of the division of the time delay at scale 1 by the scale values. For quotient values below 2, a time delay of 2 was chosen. This approach to time delay selection was developed after calculating the time delay for select signals using the minimum average mutual information method. Finally, the complexity index was computed by summing up SampEn values over scales 1 to 20.

Statistical analysis

The normality of all dependent variables was checked using the Shapiro–Wilk normality test. Paired t-test and Wilcoxon signed ranks test were used respectively for variables normally distributed and non-normally distributed to investigate the effect of dual-tasking. A p-value less than 0.05 was considered significant in all tests.


For single-scale SampEn with a time delay of one (see Table 2), there was no significant change from single-task to dual-task walking conditions, except for vertical MA-FootAV (Z = − 2.580, p = 0.010) where SampEn increased significantly. For single-scale SampEn with a larger time delay (see Table 3), dual-tasking elicited a statistically significant change in LA-ThighAV in AP direction (t = − 3.898, p = 0.001), MA-ShankAV in AP direction (t = − 3.031, p = 0.008), COG-LA in ML direction (t = − 2.818, p = 0.012), and SternumLA in ML (t = − 2.574, p = 0.020) and V direction (t = − 2.767, p = 0.014).

Table 2 Descriptive results (mean ± SD) and pairwise comparisons of the SampEn (with a time delay of one) of single-task and dual-task trials
Table 3 Descriptive results (mean ± SD) and pairwise comparisons of the SampEn (with a non-unity time delay) of single-task and dual-task trials

Figure 2, 3, 4, and 5 show that SampEn values vary with each scale. Nonetheless, the directional change from single-task to dual-task walking condition remains consistent for the examined signals. In other words, the SampEn of various signals tends to increase when participants engage in dual-tasking. A closer look reveals that LA-ShankAV and MA-ShankAV in AP direction (see Fig. 3), LA-ThighAV and MA-ThighAV in AP direction, LA-ThighAV in V direction (see Fig. 4), and COG-LA and SternumLA in ML direction (see Fig. 5) better distinguish between single-task and dual-task walking condition. These signals show significant changes across many time scales unlike the rest of the signals which discriminate between single and dual-task conditions only at sparse time scales.

Fig. 2
figure 2

The SampEn of the least-affected (left column) and most-affected (right column) foot angular velocity (FootAV) in the ML (top row), AP (middle row), and V (bottom row) directions versus time scale for single-task (ST) and dual-task (DT) walking condition. Scales highlighted in green elicited significant increase (p < 0.05) in SampEn from ST to DT

Fig. 3
figure 3

SampEn of the least-affected (left column) and most-affected (right column) shank angular velocity (ShankAV) in the ML (top row), AP (middle row), and V (bottom row) directions versus time scale for single-task (ST) and dual-task (DT) walking condition. Scales highlighted in green elicited significant increase (p < 0.05) in SampEn from ST to DT

Fig. 4
figure 4

SampEn of the least-affected (left column) and most-affected (right column) thigh angular velocity (ThighAV) in the ML (top row), AP (middle row), and V (bottom row) directions versus time scale for single-task (ST) and dual-task (DT) walking condition. Scales highlighted in green elicited significant increase (p < 0.05) in SampEn from ST to DT

Fig. 5
figure 5

SampEn of the linear acceleration of COG (COG-LA) (left column) and linear acceleration of sternum (SternumLA) (right column) in the ML (top row), AP (middle row), and V (bottom row) directions versus time scale for single-task (ST) and dual-task (DT) walking condition. Scales highlighted in green elicited significant increase in SampEn from ST to DT

Finally, the complexity index is presented in Table 4. The values for foot, shank, and thigh angular velocity in AP and V directions are significantly greater when dual-tasking. No significant changes were observed in the ML direction. For the COG-LA signal, the complexity index increased significantly in all three directions. SternumLA signal, however, was able to discriminate between the single and dual-task walking conditions in ML and V directions.

Table 4 Descriptive results (mean ± SD) and pairwise comparisons of the complexity index of single-task and dual-task trials 


In this paper, a specific condition was used to investigate the efficacy of SampEn when a non-unity time delay and multiscale analysis were employed. The specific condition was comparing single and dual-task walking trials in people with PD. Both approaches enhanced the efficacy of single-scale SampEn measure regarding discriminating between single and dual-task trials. Complexity index, which is the summary of multiscale SampEn analysis, was the most effective index of increased risk of fall, showing a significant increase from single-tasking to dual-tasking. Furthermore, while the complexity index of angular velocity signals yielded similar results, they were different from the results of trunk linear acceleration signals. Trunk linear acceleration was the only signal that was able to distinguish between the two walking conditions in three directions (ML, AP, and V).

It has previously been recommended that a time-delay greater than one be incorporated for discrete [19] and continuous time series [48]. However, it was not clear whether SampEn results would be any different if a non-unity time delay was used. In the present study, when SampEn was calculated with a time delay of one, only one out of 24 signals was able to distinguish between the two conditions. This was improved when a non-unity time delay was used, and 6 out of 24 signals showed a significant increase from single to dual-task walking conditions. Interestingly, however, the SampEn of MA-FootAV in the vertical direction, which increased significantly when using a time delay of one, did not change significantly when incorporating a non-unity time delay. Nevertheless, there was an increasing trend from single-task to dual-task. This observation might suggest that the significant difference found for the SampEn of MA-FootAV with a unity time delay is not a true finding and is rather a false positive or a type I error.

Kinematic and kinetic human signals resemble periodic signals with fluctuations from stride to stride. The dependency of the points of kinematic or kinetic signals to their neighboring points might jeopardize the discriminatory ability of the entropy measures. By incorporating a non-unity time delay, templates are constructed such that there is the least dependence between elements of each template [22]. As a result, this preprocessing would help reveal small differences between the two conditions (e.g. single vs. dual tasking) that might have been masked because of the strong periodicity of the signals. In other words, periodic signals result in smaller entropy values compared to non-periodic signals. A non-unity time delay addresses this issue, which yields greater entropy values, and, consequently, increases the discriminatory ability of entropy measures.

Similar to incorporating a non-unity time delay, our study confirms that multiscale SampEn [5] rather than single-scale SampEn can better discriminate between different walking conditions. In addition to a consistent increasing trend from walk only to dual-task walking condition, there were some time scales where values calculated for the dual-task walking condition were significantly larger. The accumulation of these pronounced values resulted in the superior efficacy of the complexity index which is the summary of multi-scale SampEn values. This index was able to successfully detect the effects of dual-tasking using 17 out of 24 signals. The results of our study were similar to the results by Ihlen et al. [3] in terms of the sparsity of significant results in multiscale SampEn plots. Furthermore, in both our study and in the study by Ihlen et al. [3], the direction of change in entropy measure was consistent for all scales. In other words, if one condition resulted in smaller entropy values, this decreasing effect held throughout all time scales.

Each signal in this study was presented in three planes of motion, mediolateral, anteroposterior, and vertical plane. It was hypothesized that the discriminatory ability of SampEn values in the mediolateral direction would be greater using either approach. This was confirmed only for linear acceleration signals obtained from COG and sternum. Bauby and Kuo [49] demonstrated that unlike the AP direction, lower extremity movement in the ML is controlled by higher level cortical structures actively processing visual information. Therefore, the differences in the ML direction for the COG and sternum plausibly arose when visual attention was divided between controlling the upper extremity and performing the visual word searching task. Interestingly, however, no concurrent differences occurred in lower extremity ML motion. In healthy young adults, disruptions to the trunk’s trajectory during walking can be offset by appropriate adaptations in the lower extremity to maintain postural control and dynamic balance [42]. This is primarily due to the fact that the neuromuscular system determines the stepping foot’s final placement (in the AP and ML directions) on its ability to predict the trunk’s future position and velocity [50, 51]. Lack of task differences in our lower extremity ML parameters may suggest a lack of appropriate adaptation in individuals with PD to assess trunk movement changes and alter their base of support [41, 42]. In other words, since COG changes between single and dual-task condition, but there is no parallel change in lower extremity ML parameters, this might suggest that participants were unable to adapt their foot placement to modify their base of support in the ML direction. This may have arisen from differences in task prioritization in people with PD compared to both aged matched controls and young adults [35, 52]. Dual-tasking evidence posits that performance during dual-tasking is not only dictated by two simultaneous tasks competing for attentional resources, but also from an individual’s ability to correctly prioritize the tasks [35]. Indeed, healthy young adults, and to some extent elderly adults, prioritize gait stability (“posture first” strategy) over performance of a secondary cognitive tasks, while people with PD utilize a “posture second” strategy where they prioritize performance of the secondary task over their gait stability [35]. The inability to correctly prioritize gait over the secondary tasks has been associated with deterioration in the prefrontal and anterior cingulate cortices, two higher level cortical structures implicated in the compensation of impaired basal ganglia functionality, in people with PD [35]. However, as this study did not examine differences between people with PD and aged matched controls, the extent to which prioritization affected Entropy measures cannot be elaborated on.

This study investigated the efficacy of SampEn in discriminating between different walking conditions in people with PD when a non-unity time delay and multiscale analysis were used. The results of the complexity index were promising as it was able to highlight small changes made to the signals due to performing a simple secondary task while walking. However, there are a few limitations that need to be considered when interpreting the findings. The secondary task, although being representative of similar situations in every-day life, was not very challenging. 12 out of 17 participants were able to successfully complete the secondary task, i.e. calling out 12 out of 12 words. Moreover, the remaining 5 participants managed to verbally call out the majority of words (mean ± SD, 10 ± 1.5). Additionally, the effect of dual-tasking on entropy measures may vary depending on the relative complexity of the secondary mental task as well as the stage of disease progression in PD.

The question of whether SampEn with and without a non-unity time delay would show more significant results if a more challenging task was used should be investigated in a future study. In addition, although treadmill walking is an efficient and commonly used method to collect data from several strides [41, 53], it may reduce the effect of dual-tasking on gait characteristics [54]. However, it has also been reported that the difference between overground walking and treadmill walking are generally small [53]. Therefore, as the objective of this study was to find the most effective entropy metric rather than generalizing the results to overground walking, this limitation would have a minimal effect on the findings. Finally, performing the secondary task while walking on the treadmill would make participants rotate their head left and right. However, the majority of head rotations were less than 20 degrees as the far left/right displayed words were still in the peripheral vision of our participants. Therefore, they only needed to slightly rotate their head when needed. Nevertheless, not studying the possible effects of this left/right head rotations on the results could be considered as a limitation of this study.


In summary, the three approaches used in the current study, i.e. single-scale SampEn with and without incorporating a non-unity time delay, multiscale SampEn, and summarizing feature called complexity index, showed different discriminatory performance across 24 continuous signals when comparing single and dual-task walking in people with Parkinson's disease. Incorporating a non-unity time delay improved the discriminatory efficacy of SampEn of 6 signals and the multiscale approach proved that not all scales could distinguish between the two walking conditions. The complexity index of the majority of signals (17 out of 24), however, successfully distinguished normal single-task walking from dual-task walking. Therefore, the complexity index of continuous signals can be used as a gait feature for classifying different treadmill walking conditions.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.



Sample entropy


Parkinson’s disease


Multiscale sample entropy


Complexity index


Approximate entropy


Trunk linear acceleration at sternum level


Center of gravity linear acceleration


Foot angular velocity


Shank angular velocity


Thigh angular velocity








Most affected


Least affected


  1. Siragy T, Nantel J. Quantifying dynamic balance in young, elderly and Parkinson’s individuals: a systematic review. Front Aging Neurosci. 2018;10(11):1–22.

    Google Scholar 

  2. Ahmadi S, Sepehri N, Wu C, Szturm T. Comparison of selected measures of gait stability derived from center of pressure displacement signal during single and dual-task treadmill walking. Med Eng Phys. 2019;74:49–57.

    Article  PubMed  Google Scholar 

  3. Ihlen EAF, Weiss A, Bourke A, Helbostad JL, Hausdorff JM. The complexity of daily life walking in older adult community-dwelling fallers and non-fallers. J Biomech. 2016;49(9):1420–8.

    Article  PubMed  Google Scholar 

  4. Lamoth CJC, Ainsworth E, Polomski W, Houdijk H. Variability and stability analysis of walking of transfemoral amputees. Med Eng Phys. 2010;32(9):1009–14.

    Article  PubMed  Google Scholar 

  5. Costa M, Peng CK, Goldberger AL, Hausdorff JM. Multiscale entropy analysis of human gait dynamics. Phys A Stat Mech Appl. 2003;330:53–60.

    Article  Google Scholar 

  6. 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.

    Article  PubMed Central  Google Scholar 

  7. Kaptein RG, Wezenberg D, IJmker T, et al. Shotgun approaches to gait analysis: insights and limitations. J Neuroeng Rehabil. 2014;11:120.

  8. Riva F, Toebes MJP, Pijnappels M, Stagni R, van Dieën JH. Estimating fall risk with inertial sensors using gait stability measures that do not require step detection. Gait Posture. 2013;38:170–4.

    Article  CAS  PubMed  Google Scholar 

  9. Tochigi Y, Segal NA, Vaseenon T, Brown TD. Entropy analysis of tri-axial leg acceleration signal waveforms for measurement of decrease of physiological variability in human gait. J Orthop Res. 2012;30(6):897–904.

    Article  PubMed  Google Scholar 

  10. 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):011006.

    Article  Google Scholar 

  11. Rispens SM, Van DJH, Van SKS, Lizama LEC, Daffertshofer A, Beek PJ, et al. Fall-related gait characteristics on the treadmill and in daily life. J Neuroeng Rehabil. 2016.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bisi MC, Stagni R. Complexity of human gait pattern at different ages assessed using multiscale entropy: from development to decline. Gait Posture. 2016;47:37–42.

    Article  CAS  PubMed  Google Scholar 

  13. Kavanagh JJ. Lower trunk motion and speed-dependence during walking. J Neuroeng Rehabil. 2009;6:9.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Ihlen EAF, Goihl T, Wik PB, Sletvold O, Helbostad J, Vereijken B. Phase-dependent changes in local dynamic stability of human gait. J Biomech. 2012;45(13):2208–14.

    Article  PubMed  Google Scholar 

  15. Gwin JT, Gramann K, Makeig SFD. Electrocortical activity is coupled to gait cycle phase during treadmill walking. Neuroimage. 2011;54(2):1289–96.

    Article  PubMed  Google Scholar 

  16. Bruijn SM, Van Dieën JH, Daffertshofer A. Beta activity in the premotor cortex is increased during stabilized as compared to normal walking. Front Hum Neurosci. 2015;9(October):1–13.

    Google Scholar 

  17. Kinsner W. Single-scale measures for randomness and complexity. In: 6th IEEE international conference on cognitive informatics, Lake Tahoo, CA, 2007, p. 554–568.

  18. Costa MD, Goldberger AL. Generalized multiscale entropy analysis: application to quantifying the complex volatility of human heartbeat time series. Entropy. 2015;17(3):1197–203.

    Article  PubMed  Google Scholar 

  19. Kaffashi F, Foglyano R, Wilson CG, Loparo KA. The effect of time delay on approximate and sample entropy calculations. Phys D Nonlinear Phenom. 2008;237(23):3069–74.

    Article  Google Scholar 

  20. Govindan RB, Wilson JD, Eswaran H, Lowery CL, PreiBl H. Revisiting sample entropy analysis. Phys A Stat Mech Appl. 2007;376:158–64.

    Article  Google Scholar 

  21. Cao L. Practical method for determining the minimum embedding dimension of a scalar time series. Phys D Nonlinear Phenom. 1997;110:43–50.

    Article  Google Scholar 

  22. Fraser AM, Swinney HL. Independent coordinates for strange attractors from mutual information. Phys Rev A. 1986;33:1134–40.

    Article  CAS  Google Scholar 

  23. Kim HS, Eykholt R, Salas JD. Nonlinear dynamics, delay times, and embedding windows. Phys D Nonlinear Phenom. 1999;127(1–2):48–60.

    Article  Google Scholar 

  24. Kang HG, Dingwell JB. Dynamic stability of superior vs inferior segments during walking in young and older adults. Gait Posture. 2009;30(2):260–3.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Dingwell JB, Robb RT, Troy KL, Grabiner MD. Effects of an attention-demanding task on dynamic stability during treadmill walking. J Neuroeng Rehabil. 2008;5:12.

  26. Blandini F, Nappi G, Tassorelli C, Martignoni E. Functional changes of the basal ganglia circuitry in Parkinson’s disease. Prog Neurobiol. 2000;62(1):63–88.

    Article  CAS  PubMed  Google Scholar 

  27. Poewe W, Seppi K, Tanner CM, Halliday GM, Brundin P, Volkmann J, et al. Parkinson disease. Nat Rev Dis Prim. 2017;3:17013.

    Article  PubMed  Google Scholar 

  28. Findley LJ. The economic impact of Parkinson’s disease. Parkinsonism Relat Disord. 2007;13(Suppl):S8–12.

    Article  PubMed  Google Scholar 

  29. Von Campenhausen S, Bornschein B, Wick R, Bötzel K, Sampaio C, Poewe W, et al. Prevalence and incidence of Parkinson’s disease in Europe. Eur Neuropsychopharmacol. 2005;15(4):473–90.

    Article  CAS  Google Scholar 

  30. Pringsheim T, Jette N, Frolkis A, Steeves TDL. The prevalence of Parkinson’s disease: a systematic review and meta-analysis. Mov Disord. 2014;29(13):1583–90.

    Article  PubMed  Google Scholar 

  31. Dorsey ER, Constantinescu R, Thompson JP. Projected number of people with Parkinson disease in the most populous nations, 2005 through 2030. Neurology. 2007;69(2):223–4.

    Article  Google Scholar 

  32. Hausdorff JM, Zitser J, Mirelman AGN. Interaction between cognition and gait in patients with Parkinson’s disease. In: Cognitive impairment and dementia in Parkinson’s disease. Oxford: Oxford University Press; 2015. p. 01.

  33. Plotnik M, Dagan Y, Gurevich T, Giladi N, Hausdorff JM. Effects of cognitive function on gait and dual tasking abilities in patients with Parkinson’s disease suffering from motor response fluctuations. Exp Brain Res. 2011;208(2):169–79.

    Article  PubMed  Google Scholar 

  34. Jehu D, Nantel J. Fallers with Parkinson’s disease exhibit restrictive trunk control during walking. Gait Posture. 2018;65:246–50.

    Article  PubMed  Google Scholar 

  35. Yogev-Seligmann G, Hausdorff JM, Giladi N. The role of executive function and attention in gait. Mov Disord. 2008;23(3):329–42.

    Article  PubMed  Google Scholar 

  36. Yogev G, Plotnik M, Peretz C, Giladi N, Hausdorff JM. Gait asymmetry in patients with Parkinson’s disease and elderly fallers: when does the bilateral coordination of gait require attention? Exp Brain Res. 2007;177(3):336–46.

    Article  PubMed  Google Scholar 

  37. Hillel I, Gazit E, Nieuwboer A, Avanzino L, Rochester L, Cereatti A, et al. Is every-day walking in older adults more analogous to dual-task walking or to usual walking? Elucidating the gaps between gait performance in the lab and during 24/7 monitoring. Eur Rev Aging Phys Act. 2019;16:6.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Bruijn SM, van Dieen JH. Control of human gait stability through foot placement. J R Soc Interface. 2018;15(20170816).

  39. Verrel J, Lövdén M, Lindenberger U. Motor-equivalent covariation stabilizes step parameters and center of mass position during treadmill walking. Exp Brain Res. 2010;207:13–26.

    Article  PubMed  Google Scholar 

  40. O’Connor SM, Kuo AD. Direction-dependent control of balance during walking and standing. J Neurophysiol. 2009;102(3):1411–9.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Siragy T, Nantel J. Absent arm swing and dual tasking decreases trunk postural control and dynamic balance in people with Parkinson’s disease. Front Neurol. 2020;11(4):1–10.

    Google Scholar 

  42. Siragy T, Mezher C, Hill A, Nantel J. Active arm swing and asymmetric walking leads to increased variability in trunk kinematics in young adults. J Biomech. 2020;99:109529.

    Article  PubMed  Google Scholar 

  43. Richman JS, Moorman JR. Physiological time-series analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol. 2000;278(6):H2039–49.

    Article  CAS  PubMed  Google Scholar 

  44. Hanavan EP Jr. A Mathematical model of the human body. AMRL-TR-64-102. AMRL TR. 1964 Oct:1-149. PMID: 14243640.

  45. Bruijn SM, van Dieën JH, Meijer OG, Beek PJ. Is slow walking more stable? J Biomech. 2009;42(10):1506–12.

    Article  PubMed  Google Scholar 

  46. McAndrew PM, Wilken JM, Dingwell JB. Dynamic stability of human walking in visually and mechanically destabilizing environments. J Biomech. 2011;44(4):644–9.

    Article  PubMed  Google Scholar 

  47. Dingwell JB, Marin LC. Kinematic variability and local dynamic stability of upper body motions when walking at different speeds. J Biomech. 2006;39(3):444–52.

    Article  PubMed  Google Scholar 

  48. McCamley JD, Denton W, Arnold A, Raffalt PC, Yentes JM. On the calculation of sample entropy using continuous and discrete human gait data. Entropy. 2018;20(10):1–20.

    Article  Google Scholar 

  49. Bauby CE, Kuo AD. Active control of lateral balance in human walking. J Biomech. 2000;33(11):1433–40.

    Article  CAS  PubMed  Google Scholar 

  50. Winter D. The biomechanics and motor control of human gait. Waterloo: University of Waterloo Press; 1987.

    Google Scholar 

  51. Winter D. Human balance and posture control during standing and walking. Gait Posture. 1995;3(4):193–214.

    Article  Google Scholar 

  52. Kelly VE, Eusterbrock AJ, Shumway-Cook A. A review of dual-task walking deficits in people with Parkinson’s disease: motor and cognitive contributions, mechanisms, and clinical implications. Parkinsons Dis. 2012.

    Article  PubMed  PubMed Central  Google Scholar 

  53. van der Krogt MM, Sloot LH, Harlaar J. Overground versus self-paced treadmill walking in a virtual environment in children with cerebral palsy. Gait Posture. 2014;40(4):587–93.

    Article  PubMed  Google Scholar 

  54. Wrightson JG, Schäfer L, Smeeton NJ. Dual-task prioritization during overground and treadmill walking in healthy adults. Gait Posture. 2020;75(2):109–14.

    Article  PubMed  Google Scholar 

Download references


We would like to acknowledge Courtney Bridgewater, Hossein Gholizadeh, Louis Goudreau, Allen Hill, Jonathan Lommen, Mary-Elise MacDonald and Andrew Smith for their assistance during data collection and processing.


Funding for this project was provided by the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery grant (RGPIN-2016-04928), NSREC Accelerator supplement (RGPAS 493045-2016) and by the Ontario Ministry of Research, Innovation and Science Early Researcher Award (ERA 16-12-206).

Author information

Authors and Affiliations



SA designed the study, analyzed the data, and wrote the manuscript; JN and TS designed and carried out the experiments, performed data analysis and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Julie Nantel.

Ethics declarations

Ethics approval and consent to participate

All participants provided written informed consent to participate and local ethics was obtained.

Consent for publication

Written informed consent was obtained from all the participants for the publication of this report.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ahmadi, S., Siragy, T. & Nantel, J. Regularity of kinematic data between single and dual-task treadmill walking in people with Parkinson’s disease. J NeuroEngineering Rehabil 18, 20 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: