Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2017; 12(11): e0188711.
Published online 2017 November 28. doi:  10.1371/journal.pone.0188711
PMCID: PMC5705160

Fractal analyses reveal independent complexity and predictability of gait

Frédéric Dierick, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing,1,2 Anne-Laure Nivard, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft,1 Olivier White, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing,3,4 and Fabien Buisseret, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing1,5,*
Steven Allen Gard, Editor


Locomotion is a natural task that has been assessed for decades and used as a proxy to highlight impairments of various origins. So far, most studies adopted classical linear analyses of spatio-temporal gait parameters. Here, we use more advanced, yet not less practical, non-linear techniques to analyse gait time series of healthy subjects. We aimed at finding more sensitive indexes related to spatio-temporal gait parameters than those previously used, with the hope to better identify abnormal locomotion. We analysed large-scale stride interval time series and mean step width in 34 participants while altering walking direction (forward vs. backward walking) and with or without galvanic vestibular stimulation. The Hurst exponent α and the Minkowski fractal dimension D were computed and interpreted as indexes expressing predictability and complexity of stride interval time series, respectively. These holistic indexes can easily be interpreted in the framework of optimal movement complexity. We show that α and D accurately capture stride interval changes in function of the experimental condition. Walking forward exhibited maximal complexity (D) and hence, adaptability. In contrast, walking backward and/or stimulation of the vestibular system decreased D. Furthermore, walking backward increased predictability (α) through a more stereotyped pattern of the stride interval and galvanic vestibular stimulation reduced predictability. The present study demonstrates the complementary power of the Hurst exponent and the fractal dimension to improve walking classification. Our developments may have immediate applications in rehabilitation, diagnosis, and classification procedures.


The stride interval of normal human walking is the time period between consecutive heel strikes of the same foot [1]. For more than two decades, a line of research focused on the understanding of the nature of the subtle variations observed in stride intervals and the origin of typical long-range structures in these variations. Today, these investigations are of paramount importance since they could provide a better understanding of the physiological mechanisms involved in normal human walking and in alterations observed in clinical practice. The nature of these stride interval variations could arise either from noisy neural processes that result in errors in the motor output or from alterations in the motor command that account for balance instabilities [2].

Normal gait is characterized by the presence of autocorrelations in the stride interval when considering walking on a sufficiently long time scale [1, 3]. The origin of these autocorrelations may be attributed to neural central pattern generators (CPGs) [1, 3] or a super CPG coupled to a forced Van der Pol oscillator [4], and/or to the biomechanics of walking [5, 6]. For many years, gait analysis has been studied with classical methods adopting biomechanical models in which variability was not of interest. More recent techniques derived from chaos theory are well adapted to analyse time series that exhibit long-range autocorrelation. Importantly, they treat variability as a meaningful interpretable signal. Since the pioneering works of Hausdorff et al. [1, 3], long-range autocorrelations in time series are estimated by the Hurst or fractal exponent (α). A fractal, introduced in 1975 by the French mathematician Benoît Mandelbrot (1924–2010) [7], is defined as a geometrical structure that has a regular or an uneven shape repeated over all scales of measurement. It is characterized by a fractal dimension (D) greater than the spatial dimension of the structure [8]. A famous example of such object is a snow flake. Objects that are statistically self-similar—parts of it show the same statistical properties at many scales—exhibit strong autocorrelation. The Hurst exponent α is a statistical measure of long-term memory of time series (see e.g. [9] for a review) and is usually associated to fractal-like behaviour. In particular, the peculiar behavior of the stride interval may be referred to as “fractal behavior” [3].

The theoretical model of optimal movement complexity [10] is based on the complementary concepts of predictability and complexity. Nature let us find optimal behavior in terms of skills and variability through evolution. As proposed by Lipsitz and Glodberger in their pioneering work, the optimal state of a biological system may be characterized by chaotic temporal variations in the steady state output that correspond to maximal complexity [11]. Any deviation from healthy state, like senescence and disease, causes a loss in complexity (see also [12]). Too few practice results in high disorder (randomness, no predictability) and excessive practice leads to high order (periodic signal, maximal predictability). Adaptation of a system to external stimuli is maximal only at an intermediate state of predictability. Furthermore, a signal from a dynamical system also holds some inherent complexity. A decrease of complexity of a physiological system results from either a reduction in the number of structural components or an alteration in the coupling function between these components. For instance, a joint can become rigid with senescence, hence decreasing the degree of freedom of the system and consequently, its complexity. A holistic approach to study these mechanisms requires to associate specific measurements to these two concepts. The Hurst exponent α captures part of the story and is well suited to reflect predictability. While the Minkowski fractal dimension D provides good measurability of the “apparent rugosity” of fractals [13] and reflects complexity. Note that the quantification of a concept such as complexity may not be linked to a unique measure; entropy-related measures have also been shown to be relvant in that domain [12]. Here, we use these parameters to complement the usual quantification of autocorrelation α in unusual and perturbed gait conditions in an attempt to probe adaptability in the framework of the model of optimal movement complexity [10].

As of today, the vast majority of studies explored autocorrelation in the stride interval during natural forward walking. In one notable exception however, Bollens et al. [14] also tested backward walking in a small sample of young healthy subjects. The authors did not find significant differences in long-range autocorrelation between both walking directions. However, backward walking measures revealed to be more sensitive than forward walking measures to classify elderly fallers compare to non-fallers [15]. The study of backward walking under the perspective of fractal analyses is therefore promising to provide more reliable predictive index of fallers, as previously proposed for forward walking [16]. Backward walking is also frequently used in sports and in rehabilitation settings, and a better understanding of the variability of stride interval in this condition is needed since it is believed that backward walking is at least partly controlled by specialized neural circuits [17].

The vestibular system provides an essential sensory contribution to the maintenance of balance during human walking [18]. Individuals with vestibular disorders show a decreased walking stability accompanied by an increased risk to fall [19]. Therefore, perturbing the vestibular system of healthy subjects with galvanic vestibular stimulation (GVS) is a well targeted mean to probe gait: it is standardized, well tolerated by subjects, generated by currently affordable electrostimulators, and easy to implement when a large number of stride intervals are recorded with an instrumented treadmill. The use of GVS is also an increasingly common clinical intervention on locomotion [2022].

Previously, autocorrelations in stride interval time series have been identified not only in healthy young adults [3] but also in children [23] and elderly [24], and even—although significantly modified—in several neurodegenerative conditions. In particular, the cases of Huntington’s disease [24], amyotrophic lateral sclerosis [25], and Parkinson’s disease have been studied [26, 27], with a hope of connecting the observed modifications of fractal behavior to some relevant evaluation of the risk of falling [16]. Here, we hypothesize that the combined effects of walking direction and the application of GVS on long-range autocorrelations in the stride interval could enhance the sensitivity of fractal analysis to identify impaired gait. We measured α and D during forward and backward walking, with and without the application of binaural and monaural GVS. We speculate that these two indexes should be able to capture differences between experimental conditions and therefore provide better indexes to classify patients.

Material and methods


Thirty-four undergraduate and graduate healthy students (18 males, 16 females) in physiotherapy took part to this study and were recruited at Haute Ecole Louvain en Hainaut (Charleroi, Belgium). Mean age was 23 years (standard deviation, SD = 2), height was 173 cm (SD = 9), mass was 69 kg (SD = 10), body mass index was 23 kg m−2 (SD = 3), and lower limb length (L), measured in standing position as the distance between the floor and great trochanter, was 88 cm (SD = 5).

Subjects were not medicated and did not exhibit any neuromusculoskeletal, orthopaedic, respiratory, or cardiovascular disorders that could influence their gait. Exclusion criteria included vestibular disorders in addition to specific GVS exclusion criteria: presence of a heart pacemaker, pregnancy, metallic brain implants, epilepsy, and skin damage behind the ears or forehead. Eligible participants were required to be able to respond to verbal questions, comprehend questionnaires, and understand instructions during the procedures of the study. Prior to participating, subjects read and signed an informed consent form. The study was approved by the ethics committee of Grand Hôpital de Charleroi and conducted in accordance with the declaration of Helsinki.

Experimental procedure

Subjects walked on an instrumented treadmill (70 cm wide, 185 cm long) with an integrated force plate and an overhead safety frame (N-Mill, Motekforce Link, The Netherlands). They wore comfortable running shoes, a safety harness, and were asked to keep their eyes fixed straight ahead. Four walking conditions were studied during two measurement sessions on two different days: forward walking without GVS (FWS0) and with GVS (FWS+) and backward walking without GVS (BWS0) and with GVS (BWS+). Session 1 included FWS0 and FWS+ conditions and session 2 included BWS0 and BWS+ conditions. During each condition, subjects walked on the treadmill for 15 minutes. Before each session, subjects were given five minutes to familiarize themselves with the treadmill and the conditions.

Subjects walked at their comfortable speed that was determined during the familiarization procedure by the same experimenter by tuning the speed of the belt while the subject was walking without GVS. The same speed was then imposed when GVS was applied. Vertical ground reaction force (Fv) and centre of pressure (CoP) of each foot was recorded at a sampling rate of 500 Hz using the manufacturer’s software (CueFors 2, Motekforce Link, The Netherlands). Time series stride interval were computed from heel strikes (during forward walking) or toe strikes (during backward walking) of the right foot identified on Fv-time histories and time series step width from maximal medio-lateral displacement of CoP of two consecutive steps. At completion of both sessions, four time series containing the values of the stride intervals in the different conditions were obtained for each subject. Typical plots are displayed in Fig 1.

Fig 1
Typical stride interval time series in the different experimental conditions.

Bipolar GVS was applied with a regulated, direct-current device (Compex 3 Professional, Compex Medical SA, Switzerland) with a maximum output current of 20 mA by steps of 0.125 mA. The carbon electrodes (20 cm2) were covered with a saline-soaked sponge held in place over the mastoids or forehead with a strap. The 34 subjects were randomly exposed to one of the three different transcranial stimulation conditions: binaural (n = 12), unilateral left (n = 11), and unilateral right (n = 11). For the binaural stimulation, electrodes were randomly placed over the mastoids with a cathode-left anode-right or cathode-right anode-left montage. For the monoaural stimulations, the cathode was randomly placed over the right mastoid and the anode on the right part of the forehead (right stimulation) or the cathode was placed over the left mastoid and the anode on the left part of the forehead (left stimulation). Fig 2 shows a schematic representation of the location of the electrodes over the head for the 3 stimulation conditions. The intensity was set at the highest sensory tolerance threshold, that was determined by increasing the current intensity slowly by 0.125-mA steps. That intensity was maintained constant throughout the walking period. Mean current density for subjects was 0.07 mA cm−2 (range: 0.04–0.08 mA cm−2). The intensity and duration of the GVS adhered to the safety criteria for transcranial direct current stimulation [28]. After FWS+ and BWS+ conditions, each subject completed a home-made French translation of the Motion Sickness Assessment Questionnaire (MSAQ) [29], that consists of 16 questions, allowing to differentiate motion sickness symptoms along the gastrointestinal, central, peripheral, and sopite-related dimensions.

Fig 2
The three types of galvanic stimulations used.

Data analysis

The treadmill software directly computed the mean stride interval, T, and the mean step width, w, for each subject in each condition. The stride amplitude (θ0), i.e. the angle between the leg and the vertical at heel strike, has then been computed from the relation vT4=Ltanθ0, displayed in Fig 3, where v is the walking speed and L is the lower limb length of the subject.

Fig 3
Schematic representation of the body (circle) and the lower limb (solid lines) when the heel strikes the ground in forward direction (FW) or the toes in backward direction (BW).

The temporal analysis of our experimental data has to go a step beyond mere descriptive statistics to study the information contained in stride interval variability. Let T = {Ti: i = 1, …, n} be a time series, where Ti is the stride interval of cycle i and where n is the number of cycles recorded during 15 minutes. The first indicator of variability is the coefficient of variation CVT = SD(T)/T. Because CVT provides no information on the dynamics of the stride interval fluctuations, several indexes characterizing the dynamical structure of the time series have been computed such as the Hurst exponent, α, the spectral exponent, γ, and the Minkowski fractal dimension, D.

We assess the presence of long-range autocorrelations with the Hurst exponent α and the spectral exponent γ. These parameters quantify the “predictability” of the time series. As pointed out in e.g. [30, 31], using a single parameter may not be generally enough to assess the presence or not of such autocorrelations. The Hurst exponent, computed by using the Detrended Fluctuation Analysis (DFA) with a linear detrending [32], provides a diagnostic on the long-range trend of the time series. DFA consists in several steps. First, one has to compute the shifted time series Tτ = {Ti+τ: i = 1, …, nτ} and the cumulated time series Z={Zi=j=1i(Tj-T):i=1,,n}. One has then to divide the cumulated time series Z into windows of length l, leading to the samples Z(m)(l), m labelling the window. For each window, a local least squares linear fit is calculated, leading to the fitted values Z¯(m)(l). Second, one computes the fluctuation function F(l)=1ni=1n[Zi(m)(l)Z¯i(m)(l)]2. The Hurst exponent is then defined as the scaling exponent of F(l), i.e. F(l) [proportional, variant] lα. Stationary times series originating from long-range (anti)correlated processes correspond to 0.5 < α [less, similar] 1 (0 ≤ α < 0.5), respectively. When α = 0.5, the process is random. Values larger than 1 correspond to unbounded, unstable, processes [9]. A strongly autocorrelated signal can be denoted “predictable”: Its value at a given step is strongly dependent on the system’s previous state.

The spectral exponent γ can be extracted from the low-frequency behaviour of the power spectral density P(f) of T, f being the frequency: P(f) [proportional, variant] fγ. Actually P(f) is the Fourier transform of the autocorrelation function C(τ) = E(TTτ), where E denotes the average value. The parameter γ is expected to take values between 0 and 1 for long-range autocorrelated processes. For large enough time series, the asymptotic relation


relates α and γ [9].

We finally compute the Minkowski fractal dimension D of the time series, defined through the box-counting method stating that, if N(ϵ) is the number of square boxes of size ϵ needed to fully cover the time series once plotted, then one has N(ϵ) [proportional, variant] ϵD for small ϵ. For time series such as the ones we deal with, D will typically lie between 1 (differentiable curve) and 2 (surface with differentiable boundary). Even if T does not define a fractal in a rigorous mathematical way, D can be thought of as a relevant estimator of the “apparent roughness” of the corresponding curve [33]. The roughness of the stride interval time series—i.e. the variability of fluctuations from one stride to another—may be associated to the “complexity” of the process [34].

It is worth mentioning at this stage that D and α may be seen as independent variables characterizing a time series. The link α = 2 − D is actually valid only for some widely studied random walks, but processes with arbitrary values of α and D can be built and could be more representative of realistic time series [35].

For completeness we show in Fig 4 typical plots of the necessary steps performed to compute the parameters α, γ and D. All these quantities are the slopes of the different linear regressions performed.

Fig 4
Typical logarithmic plots showing the linear regressions performed to find α (upper left panel), γ (upper right panel) and D (lower left panel) in the four conditions.

Statistical analysis

All data were checked for normality (Shapiro-Wilk) and equal variance tests. A two-way (WD × GVS) repeated measures ANOVA with post hoc Holm-Sidak method for pairwise multiple comparisons has been performed and used to examine the effects of walking direction (FW or BW), GVS (S+ or S0), and their interaction on the computed parameters. The significance level has been set at p = .05 for all analyses and post hoc statistical power has been calculated. The correlation between parameters are provided as Spearman’s rank correlation coefficient ρ, Pearson’s correlation coefficient r, and prinicipal component analysis (PCA).

Finally, the median (Mdn) scores and interquartile ranges (IQR) related to the four symptom dimensions assessed by the MSAQ questionnaire have been computed from the 16 items scores according to [29].

All statistical procedures were performed with SigmaPlot software version 11.0 (Systat Software, San Jose, CA). Indexes have been computed using R free software environment (v. 3.2.2) [36].


Subjects adopted mean comfortable forward speed of 4.4 km h−1 (SD = 0.4) and backward speed of 2.2 km h−1 (SD = 0.5). The rather small SDs indicate that our sample was quite homogeneous with respect to that parameter. Results of all parameters analysed according to the four experimental conditions are reported in Table 1. Some parameters particularly relevant for the Discussion are also shown as box plots in Fig 5. In any of the four experimental conditions, stride interval time series are such that 0 < γ < 1 and 0.5 < α < 1. These values confirm that time series characterise long-range autocorrelated processes with memory [30]. Note that our value of α in the FWS0 condition (95% CI of [0.78–0.90]) are fully compatible with the interval reported in the meta-analysis [37] in healthy subjects (95% CI of [0.85–0.97]). Moreover, Pearson’s correlation coefficient between α and γ (computed over the complete data set) is equal to r = 0.287, p < .01. Parameters α and γ can then be seen as linearly correlated as expected from Eq (1). We decided to keep α since it is the most widely used autocorrelation-related parameter in the literature [37]. The parameter γ will no longer be considered since it only confirms the presence of long-range autocorrelations. In contrast, Pearson’s and Spearman’s correlation coefficients between α and D are equal to r = − 0.081, p = .350 and ρ = −0.031, p = .719 respectively. Moreover, a PCA conducted on our full data set, including the parameters listed in Table 1, shows that nearly 66% of the total variance is carried by the first two dimensions, the angle between α and D being equal to 131° while that between alpha and γ being equal to 14°. To sum up, these analyses demonstrate that α and D can be considered as independent parameters and are therefore both retained in our analysis.

Table 1
Mean ± standard deviations values for the parameters computed in the four experimental conditions.
Fig 5
Box plots of the mean stride time (T, upper left panel), the coefficient of variation (CVT, upper right panel), the stride amplitude (θ0, middle left panel), the mean step width (w, middle right panel), the Hurst exponent α (lower left ...

We first report significant differences found in the two-way RM ANOVA related to WD and GVS factors. WD induces significant differences in all variables (see also Fig 5): T (F = 216, p < .001, partial effect size ηp2=.991), CVT (F = 324, p < .001, ηp2=.999), w (F = 13.6, p < .001, ηp2=.851), θ0 (F = 429, p < .001, ηp2=.997), α (F = 5.594, p = .024, ηp2=.222) and D (F = 12.1, p = .001, ηp2=.276). GVS has a less important influence but has nevertheless a significant impact on T (F = 5.58, p = .024, ηp2=.135), CVT (F = 6.24, p = .018, ηp2=.872) and w (F = 11.0, p = .002, ηp2=.190). The two-way RM ANOVA also reported significant interaction between WD and GVS for T (F = 1.2, p = .003, ηp2=.236), w (F = 4.75, p = .036, ηp2=.126), θ0 (F = 11.4, p = .002, ηp2=.257) and D (F = 7.69, p = .009, ηp2=.189).

Post hoc pairwise multiple comparisons results are listed in Table 2. Only w and D are significantly modified when the GVS was active during forward walking. The GVS had a stronger impact in backward walking, with a significant modification of T, CVT, θ0, α. The factor WD had a major impact on all variables. All parameters were significantly different when comparing FWS0 and BWS0 and all parameters but α and D were significantly different when comparing FWS+ and BWS+.

Table 2
Significance level (p–values) of the post hoc (Holm-Sidak) pairwise multiple comparisons between the different experimental conditions for the parameters analyzed.

Table 3 summarizes MSAQ results according to the four dimensions listed in [29]. The highest median scores were reached for the items related to the central nervous system dimension. It is worth noticing that beyond these four scores, the third item of the questionnaire [29]—“I felt annoyed/irritated” scored highest in the FWS+ condition (Mdn = 4, IQR = [2 − 6]) and in the BWS+ condition (Mdn = 3, IQR = [1 − 5]).

Table 3
Median scores (in %) and interquartile range (IQR) results of the MSAQ.

Finally, one can ask the question as to how basic kinematic parameters relate to high level indexes such as α and D? Spearman’s coefficient between α and θ0, calculated from each subject’s data, is equal to ρ = −0.207, p = .017. It is therefore relevant to consider that both parameters are correlated. Fig 6 depicts the evolution of α vs θ0 in the FWS0 and BWS0 conditions and compares it to that of data presented in Ahn and Hogan’s model (full circles). Fig 7 reports the evolution of the fractal dimension versus the Hurst exponent in the four experimental conditions. It shows that changing walking direction with respect to the control condition FWS0 leads to an increase of α and a decrease of D (independently of GVS). However, turning on GVS during forward walking only leads to a decrease of D. There is some decoupling between α and D.

Fig 6
Evolution of the Hurst exponent (α) versus stride amplitude (θ0) in the FWS0 and BWS0 conditions (empty open circles).
Fig 7
Evolution of the fractal dimension (D) versus the Hurst exponent (α) taken as complexity and predictability indices respectively for the four experimental conditions.


Locomotion has been used as a tool to identify and characterise diverse impairments. Here, we set out to use techniques beyond classical linear analyses of spatio-temporal gait parameters in order to define more sensitive indexes. We asked participants to walk under perturbed conditions induced either by reversing the direction of walking or perturbing the vestibular system, and measured proxies of “complexity” and “predictability” through the Hurst exponent and the fractal dimension.

Walking direction and galvanic vestibular stimulation

Walking direction has a major impact on our results. Stride interval increased in backward walking compared to forward walking which is in line with previous work e.g. [38]. This fact is coherent with the lower walking speed spontaneously adopted by the subjects when walking backward. The kinematic parameters CVT and w were larger in the backward walking conditions than in forward walking. Fluctuations in the stride interval during backward walking were indeed larger than in forward walking. The increase of CVT has been reported previously [14, 15] in young and elderly adults while, to our knowledge, the corresponding increase in the step width has not been reported.

While the existence of autocorrelations in backward and forward walking has been acknowledged in [14], that study did not find a significant difference in the Hurst exponent, presumably because that experiment included a small sample (12 subjects). The Hurst exponent appeared to be larger in backward compared to forward walking. While neurological diseases generally decrease α [37] (more random motion), an increase in α has been reported in children up to typically 14 years old [23]. Both backward walking in adults and forward walking in children can be related to learning processes, with stereotyped, more predictable, motion.

In their recent paper, Ahn and Hogan showed that long-range autocorrelations may emerge from the dynamics of a particular pendular model of walking described in detail in [6, 39]. It is, to our knowledge, the only model linking kinematic variables and autocorrelations indexes. Here, we compared experimental data with that model for the first time. One of the key plots of [6] shows the variation of the Hurst exponent computed on a range of 500 strides for realistic values of the different parameters versus stride amplitude θ0, which is fixed in their model. A comparison of this plot to our results is shown in Fig 6. Results from the FWS0 condition are in agreement with the model. The BWS0 condition can be compared too since the dynamical equations presented in [39] are invariant with respect to time reversal. Our results match quite well with the trend of the model. Moreover, GVS had no significant influence on α. This is in favour of a mostly mechanical origin of long-range stride interval autocorrelation. Galvanic vestibular stimulation is a known procedure to electrically stimulate vestibular afferents [4042]. Here, we used continuous GVS with an average current of 1.4 mA, which was well tolerated by all subjects (low MSAQ median scores).

Previously, vestibular inputs were thought to be primarily required for stabilizing the head to ensure stable gaze control during gait and for spatial orientation in navigational tasks [40, 43, 44]. More recently, it has also been suggested that vestibular inputs play a role in maintaining dynamic walking stability since they generate phase-dependent influences on lower body control during walking by fine tuning the timing and magnitude of foot displacement [18, 45, 46]. In agreement with those recent findings, our results show that GVS significantly modifies T, CVT and w. The magnitude of CVT is regularly associated to the risk of fall [47]. We observed that applying GVS significantly decreased CVT in backward walking; hence training in BWS+ may be relevant to decrease the risk of fall. Previous research demonstrated that GVS mostly affects stability in the medio-lateral direction [4042, 48]. This is in accordance with our results that showed an increased w during forward walking.

Interestingly, we also found that GVS induced larger D in FWS0 compared to FWS+ condition. This indicates a less complex stride interval time series in the non-stimulated condition of forward walking. This result is compatible with [45] showing that planning of the foot placement at heel contact is modulated by vestibular information. Here, we provide another evidence of the influence of GVS on walking variability. It is known that the vestibular system is essential to the maintenance of balance throughout the stepping cycle, with phase-velocity/cadence-dependent modulation on the activity of hip, knee and ankle muscles [18]. Vestibular-muscle coupling is specific for each muscle, probably organised according to each muscle’s functional role in whole-body stabilization during walking. Our analyses suggest that the less complex nature of the stride interval time series reflects the disruption of dynamic balance evoked by GVS. It is worth mentioning that treatments are aleady known to have a positive impact on gait complexity: The use of a stochastic resonance-based therapy in the elderly indeed increases complexity in the center of pressure displacement [49].

The optimal complexity model

Extending on the maximal complexity hypothesis of Lipsitz and Goldberger [11], Stergiou and Decker [10] proposed that time series originating from human motion could be classified by using two indices. The first catches signal complexity and the other measures its predictability. In this context, a healthy motion should be chaotic, characterized by a maximal complexity reflecting the adaptability of the subject to external perturbation, and an intermediate predictability. Pathological motion should be characterized by a lower complexity (fewer adaptability) and a predictability that could be either lower (random, “drunken-sailor-like”, motion) or higher (robotic motion) than the healthy motion.

Following on that line, we interpret D as a measure of the complexity of gait time series. Indeed, a large fractal dimension is associated to an apparently rough time series, with abrupt relative changes of values stepwise. A complex time series may be the signature of an adaptable behavior. The more a subject is able to change her/his stride interval from one cycle to another, the more s/he should be able to modify her/his pattern. Therefore, D could be a good indicator of complexity during walking. Moreover, we think that the Hurst exponent—that was independent of D—could be a relevant predictability index. Indeed, α can discriminate between a random motion (α = 0.5) and a far more predictable, strongly autocorrelated, time series (0.5 < α [less, similar] 1). So α may provide an answer to the question as to how much a stride interval depends on history? This is exactly what predictability stands for. Healthy subjects should be characterised by a maximal value of D (high complexity, good adaptation skills) and an optimal value of α (good but not too high predictability). Any significant deviations from these values could indicate pathological motion linked to one or both dimensions (predictability or complexity).

Our results are displayed in a (α, D)-plane in Fig 7. It clearly appears that the FWS0 condition—the healthy motion—has the higher complexity and an intermediate predictability in agreement with [10, 11]. The other conditions, non-standard but not pathological either, have lower complexities. Walking backward without GVS leads to a larger value of α, that is a more stereotyped, more predictable walking. GVS slightly decreased α in backward walking. In that condition, walking gets closer to a random process, presumably because of the perturbation of the vestibular system. It is worth noticing that at least one of the two parameters (α, D) is significantly modified when going from one condition to another. A two-dimensional representation is necessary to classify all the experimental conditions we study. Considering only α may be too restrictive to discriminate pathologies. We suggest that other non-linear indexes are worth to be added.

Previous studies only computed the Hurst exponent and implicitly considered that the fractal dimension and the exponent were related. Here, we computed D beside α and show for the first time that these two parameters are actually decoupled in some conditions. As already pointed out, children walking forward also have a larger α than healthy young adults walking forward. We used open access dataset [23] to calculate the average D for the 50 children having participated to the study. We found values equal to 1.45 ± 0.21, a significantly lower value than in our FWS0 condition (t = −6.77, p < .001), as expected. In Parkinson’s disease, α are smaller than in young healthy adults. It has been shown that α decreases with disease’s severity [27]. We have computed D = 1.69 ± 0.10 from the data of [27] (20 patients with Parkinson’s disease walking for 10 min, with α = 0.70 ± 0.09). Although lower than our maximal, FWS0, the difference between both values is not significant (t = −0.941, p = .351). Similarly, neurodegenerative pathologies have actually been shown to generally decrease α with respect to its optimal value [37]. As can be deduced from the above discussion, the parameters D and α are good candidates to disentangle and characterise the main long-term features of walking. The Hurst exponent α is a widely used indicator of long-term autocorrelations, and adding D opens new classification perspectives. Our findings may have immediate applications in rehabilitation, diagnosis, and classification procedures.

We also think that this field could benefit in a near future from new techniques such as a representation of the stride interval time series in terms of complex networks (visibility graphs) [50, 51]. This technique has already proven its efficiency to distinguish healthy from epileptic EEG signals [52]. Hence it can reasonably be assumed that visibility graphs could provide relevant information on the structure of stride interval time series. New classification schemes resting on machine learning could also shed new light on walking dynamics. Algorithms like random forests could help to find better indices to disentangle the different experimental conditions [53]. There is hope that such new techniques could better classify the stride interval time series, but with less common indices, either less intuitive or less easily compared to the literature in walking analysis. Such a research program is beyond the scope of the present study and we leave it for future investigations.


Our findings show that stride interval dynamics behave as a chaotic system exhibiting long-range autocorrelations independently of walking direction. The Hurst exponent is increased when walking backward, suggesting that the more predictable fluctuations of the stride interval reflect more stereotyped motion adopted by subjects in this unusual condition. The magnitude of these fluctuations are however larger in backward walking, due to the weaker stability of the subjects. The Minkowski fractal dimension complements the characterisation of stride interval variability by considering complexity, or more intuitively, the adaptive capacities of the subject. Any nonstandard condition reduced complexity. The present study thus paves the way toward more accurate classification methods of healthy or pathological walking by considering complexity and predictability of simple stride interval time series.


The authors thank C Detrembleur and T Warlop for stimulating discussions at early stages of this work and for having provided us the data of [27], and M Scohier for useful comments.

Funding Statement

The authors received no specific funding for this work.

Data Availability

Data Availability

Data are available from Open Science Framework (


1. Hausdorff JM, Purdon PL, Peng CK, Ladin ZVI, Wei JY, Golberger AL. Fractal dynamics of human gait: stability of long-range correlations in stride interval fluctuations. J Appl Physiol. 1996;80: 1448–1457. [PubMed]
2. Kurz MJ, Wilson TW, Arpin DJ. Stride-time variability and sensorimotor cortical activation during walking. Neuroimage. 2012;59: 1602–1607. doi: 10.1016/j.neuroimage.2011.08.084 [PubMed]
3. Hausdorff JM, Peng CK, Ladin Z, Wei JY, Goldberger AL. Is walking a random walk? Evidence for long-range correlations in stride interval of human gait. J Appl Physiol. 1995; 78: 349–358. [PubMed]
4. West BJ, Scafetta N. Nonlinear dynamical model of human gait. Phys Rev E Stat Nonlin Soft Matter Phys. 2003;67: 051917 doi: 10.1103/PhysRevE.67.051917 [PubMed]
5. Gates DH, Su JL, Dingwell JB. Possible biomechanical origins of the long-range correlations in stride intervals of walking. Physica A. 2007;380: 259–270. doi: 10.1016/j.physa.2007.02.061 [PMC free article] [PubMed]
6. Ahn J, Hogan N. Long-Range correlations in stride intervals may emerge from non-chaotic walking dynamics. PLoS One. 2013;8: e73239 doi: 10.1371/journal.pone.0073239 [PMC free article] [PubMed]
7. Mandelbrot B. Les objets fractals: forme, hasard et dimension. 1st ed Paris: Flammarion; 1975.
8. Random House Kernerman Webster’s College Dictionary, 2010. Available from:
9. Kantz H, Schreiber T. Nonlinear Time Series Analysis. 2nd edition Cambridge: Cambridge University Press; 2004.
10. Stergiou N, Decker LM. Human movement variability, nonlinear dynamics, and pathology: is there a connection?. Hum Mov Sci. 2011;30: 869–888. doi: 10.1016/j.humov.2011.06.002 [PMC free article] [PubMed]
11. Lipsitz LA, Goldberger AL. Loss of ‘Complexity’ and Aging. JAMA 1992; 67: 1806–1809. doi: 10.1001/jama.1992.03480130122036 [PubMed]
12. Goldberger AL, Peng CL, Lipsitz LA. What is physiologic complexity and how does it change with aging and disease?. Neurobiology of Aging 23 (2002) 23–26. [PubMed]
13. Falconer KJ. On the Minkowski measurability of fractals. Proc Am Math Soc. 1995;123: 1115–1124. doi: 10.1090/S0002-9939-1995-1224615-4
14. Bollens B, Crevecoeur F, Detrembleur C, Warlop T, Lejeune TM. Variability of human gait: effect of backward walking and dual-tasking on the presence of long-range autocorrelations. Ann Biomed Eng. 2013;42: 742–750. doi: 10.1007/s10439-013-0961-9 [PubMed]
15. Fritz NE, Worstell AM, Kloos AD, Siles AB, White SE, Kegelmeyer DA. Backward walking measures are sensitive to age-related changes in mobility and balance. Gait Posture. 2013; 37:593–597. doi: 10.1016/j.gaitpost.2012.09.022 [PubMed]
16. Hausdorff JM. Gait dynamics, fractals and falls: finding meaning in the stride-to-stride fluctuations of human walking. Hum Mov Sci. 2007;26: 555–589. doi: 10.1016/j.humov.2007.05.003 [PMC free article] [PubMed]
17. Hoogkamer W, Meyns P, Duysens J. Steps forward in understanding backward gait: from basic circuits to rehabilitation. Exerc Sport Sci Rev. 2014;42: 23–29. doi: 10.1249/JES.0000000000000000 [PubMed]
18. Dakin CJ, Inglis JT, Chua R, Blouin JS. Muscle-specific modulation of vestibular reflexes with increased locomotor velocity and cadence. J Neurophysiol. 2013;110: 86–94. doi: 10.1152/jn.00843.2012 [PubMed]
19. Whitney SL, Hudak MT, Marchetti GF. The dynamic gait index relates to self-reported fall history in individuals with vestibular dysfunction. J Vestib Res. 2000;10: 99–105. [PubMed]
20. Buzzi UH, Stergiou N, Kurz MJ, Hageman PA, Heidel J. Nonlinear dynamics indicates aging affects variability during gait. Clin Biomech. 2003;18: 435–443. doi: 10.1016/S0268-0033(03)00029-9 [PubMed]
21. Dingwell JB, Cusumano JP. Nonlinear time series analysis of normal and pathological human walking. Chaos. 2000;10: 848–863. doi: 10.1063/1.1324008 [PubMed]
22. Dingwell JB, Marin LC. Kinematic variability and local dynamic stability of upper body motions when walking at different speeds. J Biomech. 2006;39: 444–452. doi: 10.1016/j.jbiomech.2004.12.014 [PubMed]
23. Hausdorff JM, Zemany L, Peng C, Goldberger AL. Maturation of gait dynamics: stride-to-stride variability and its temporal organization in children. J Appl Physiol. 1999; 86: 1040–1047. [PubMed]
24. Hausdorff JM, Mitchell SL, Firtion R, Peng CK, Cudkowicz ME, Wei JY, et al. Altered fractal dynamics of gait: reduced stride-interval correlations whit aging and Huntington’s disease. J Applied Physiol. 1997; 82: 262–269. [PubMed]
25. Hausdorff JM, Lertratanakul A, Cudkowicz ME, Peterson AL, Kaliton D, Goldberger AL. Dynamic markers of altered gait rhythm in amyotrophic lateral sclerosis. J Appl Physiol. 2000; 88: 2045–2053. [PubMed]
26. Hausdorff JM. Gait dynamics in Parkinson’s disease: common and distinct behavior among stride length, gait variability, and fractal-like scaling. Chaos. 2009;19: 026113 doi: 10.1063/1.3147408 [PubMed]
27. Warlop T, Detrembleur C, Bollens B, Stoquart G, Crevecoeur F, Jeanjean A, et al. Temporal organization of stride duration variability as a marker of gait instability in Parkinson’s disease. J Rehabil Med. 2016;48: 865–871. doi: 10.2340/16501977-2158 [PubMed]
28. Utz SK, Dimova V, Oppenlander K, Kerkhoff G. Electrified minds: transcranial direct current stimulation (tDCS) and galvanic vestibular stimulation (GVS) as methods of non-invasive brain stimulation in neuropsychology—a review of current data and future implications. Neuropsychologia. 2010;48: 2789–2810. doi: 10.1016/j.neuropsychologia.2010.06.002 [PubMed]
29. Gianaros PJ, Muth ER, Mordkoff JT, Levine ME, Stern RM. A questionnaire for the assessment of the multiple dimensions of motion sickness. Aviat Space Environ Med. 2001;72: 115–119. [PMC free article] [PubMed]
30. Rangarajan G, Ding M. Integrated approach to the assessment of long range correlation in time series data. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 2000;61: 4991–5001. [PubMed]
31. Crevecoeur F, Bollens B, Detrembleur C, Lejeune TM. Towards a “gold-standard” approach to address the presence of long-range auto-correlation in physiological time series. J Neurosci Methods. 2010;192: 163–172. doi: 10.1016/j.jneumeth.2010.07.017 [PubMed]
32. Peng CK, Buldyrev SV, Havlin S, Simons M, Stanley HE, Goldberger AL. Mosaic organization of DNA nucleotides. Phys Rev E Stat Phys Plasmas Fluids Relat Interdiscip Topics. 1994;49: 1685–1689. [PubMed]
33. Gneiting T, Ševčíková H, Percival DB. Estimators of fractal dimension: assessing the roughness of time series and spatial data. Stat Sci. 2012;27: 247–277. doi: 10.1214/11-STS370
34. Vaillancourt DE, Newell KM. Changing complexity in human behavior and physiology through aging and disease. Neurobiol Aging. 2002;23: 1–11. doi: 10.1016/S0197-4580(01)00247-0 [PubMed]
35. Gneiting T, Schlather M. Stochastic models that separate fractal dimension and the Hurst effect. SIAM Rev. 2004;46: 269–282. doi: 10.1137/S0036144501394387
36. More specifically, the packages nonlinearTseries, fractaldim have been used to compute the DFA Hurst exponent and the Minkowski fractal dimension respectively.
37. Moon Y, Sung J, An R, Hernandez ME, Sosnoff JJ. Gait variability in people with neurological disorders: a systematic review and meta-analysis. Hum Mov Sci. 2016;47: 197–208. doi: 10.1016/j.humov.2016.03.010 [PubMed]
38. Katsavelis D, Mukherjee M, Decker L, Stergiou N. Variability of lower extremity joint kinematics during backward walking in a virtual environment. Nonlinear Dynamics Psychol Life Sci. 2010;14: 165–178. [PubMed]
39. Ahn J, Hogan N. A simple state-determined model reproduces entrainment and phase-locking of human walking. PloS One. 2012;7: e47963 doi: 10.1371/journal.pone.0047963 [PMC free article] [PubMed]
40. Jahn K, Strupp M, Schneider E, Dieterich M, Brandt T. Differential effects of vestibular stimulation on walking and running. Neuroreport. 2000;11: 1745–1748. doi: 10.1097/00001756-200006050-00029 [PubMed]
41. Scinicariello AP, Eaton K, Inglis JT, Collins JJ. Enhancing human balance control with galvanic vestibular stimulation. Biol Cybern. 2001;84: 475–480. doi: 10.1007/PL00007991 [PubMed]
42. Balter SGT, Stokroos RJ, De Jong I, Boumans R, Van De Laar M, Kingma H. Background on methods of stimulation in galvanic-induced body sway in young healthy adults. Acta Otolaryngol. 2004;124: 262–271. doi: 10.1080/00016480310015245 [PubMed]
43. Fitzpatrick RC, Wardman DL, Taylor JL. Effects of galvanic vestibular stimulation during human walking. J Physiol. 1999;517: 931–939. doi: 10.1111/j.1469-7793.1999.0931s.x [PubMed]
44. Pozzo T, Berthoz A, Lefort L. Head stabilization during various locomotor tasks in humans. I. Normal subjects. Exp Brain Res. 1990;82: 97–106. doi: 10.1007/BF00230842 [PubMed]
45. Bent LR, Inglis JT, McFadyen BJ. When is vestibular information important during walking? J Neurophysiol. 2004;92: 1269–1275. doi: 10.1152/jn.01260.2003 [PubMed]
46. Blouin JS, Dakin CJ, van den Doel K, Chua R, McFadyen BJ, Inglis JT. Extracting phase-dependent human vestibular reflexes during locomotion using both time and frequency correlation approaches. J Appl Physiol. 2011;111: 1484–1490. doi: 10.1152/japplphysiol.00621.2011 [PubMed]
47. van Schooten KS, Sloot LH, Bruijn SM, Kingma H, Meijer OG, Pijnappels M, et al. Sensitivity of trunk variability and stability measures to balance impairments induced by galvanic vestibular stimulation during gait. Gait Posture. 2011;33: 656–660. doi: 10.1016/j.gaitpost.2011.02.017 [PubMed]
48. Kennedy PM, Cressman EK, Carlsen AN, Chua R. Assessing vestibular contributions during changes in gait trajectory. Neuroreport. 2005;16: 1097–1100. doi: 10.1097/00001756-200507130-00013 [PubMed]
49. Costa M, Priplata AA, Lipsitz LA, Wu Z, Huang NE, Goldberger AL, et al. Noise and poise: Enhancement of postural complexity in the elderly with a stochastic-resonance-based therapy. Europhys Lett. 2007. March; 77: 68008 doi: 10.1209/0295-5075/77/68008 [PMC free article] [PubMed]
50. Gao ZK, Yang YX, Fang PC, Zou Y, Xia CY, Du M. Multiscale complex network for analyzing experimental multivariate time series. EPL. 2015;109: 30005 doi: 10.1209/0295-5075/109/30005
51. Gao ZK, Small M, Kurths J. Complex network analysis of time series. EPL. 2016;116: 50001 doi: 10.1209/0295-5075/116/50001
52. Gao ZK, Cai O, Yang YX, Dong N, Zhang SS. Visibility graph from adaptive optimal kernel time-frequency representation for classification of epileptiform EEG. Int J Neural Syst. 2017;27: 1750005 doi: 10.1142/S0129065717500058 [PubMed]
53. Ho TK. The random subspace method for constructing decision forests. IEEE Trans Pattern Anal Mach Intell. 1998;20: 832–844. doi: 10.1109/34.709601

Articles from PLoS ONE are provided here courtesy of Public Library of Science