Research Article  Open Access
Henryk Josiński, Adam Świtoński, Agnieszka Michalczuk, Piotr Grabiec, Magdalena Pawlyta, Konrad Wojciechowski, "Assessment of Local Dynamic Stability in Gait Based on Univariate and Multivariate Time Series", Computational and Mathematical Methods in Medicine, vol. 2019, Article ID 6917658, 13 pages, 2019. https://doi.org/10.1155/2019/6917658
Assessment of Local Dynamic Stability in Gait Based on Univariate and Multivariate Time Series
Abstract
The ability of the locomotor system to maintain continuous walking despite very small external or internal disturbances is called local dynamic stability (LDS). The importance of the LDS requires constantly working on different aspects of its assessment method which is based on the shortterm largest Lyapunov exponent (LLE). A state space structure is a vital aspect of the LDS assessment because the algorithm of the LLE computation for experimental data requires a reconstruction of a state space trajectory. The gait kinematic data are usually one or threedimensional, which enables to construct a state space based on a uni or multivariate time series. Furthermore, two variants of the shortterm LLE are present in the literature which differ in length of a time span, over which the shortterm LLE is computed. Both a state space structure and the consistency of the observations based on values of both shortterm LLE variants were analyzed using time series representing the joint angles at ankle, knee, and hip joints. The shortterm LLE was computed for individual joints in three state spaces constructed on the basis of either univariate or multivariate time series. Each state space revealed walkers’ locally unstable behavior as well as its attenuation in the current stride. The corresponding conclusions made on the basis of both shortterm LLE variants were consistent in ca. 59% of cases determined by a joint and a state space. Moreover, the authors present an algorithm for estimation of the embedding dimension in the case of a multivariate gait time series.
1. Introduction
Stability means the ability to return to a stable state after having been subjected to some form of perturbation. Focusing on gait, if infinitesimally small perturbations, naturally occurring tiny variations in the walking surface and/or natural noise in the neuromuscular system, are concerned, then the ability of the locomotor system to keep the gait smooth by attenuating them is called local dynamic stability (LDS) [1]. The aforementioned disturbances are the cause of slightly different conditions at the beginning of successive strides. As a consequence, the LDS can be assessed using a measure of the extreme sensitivity to initial conditions.
Gait stability is of great importance for older people who are considered prone to falls. It requires constantly working on different aspects of LDS assessment method, which is derived from the dynamical systems theory. The method is based on a trajectory in a state space which is reconstructed from time series generated by a dynamical system. The dynamical properties of a system in the true state space are preserved under the reconstruction process, which enables to analyze the system’s behavior using the reconstructed trajectory, with particular emphasis on system’s sensitivity to initial conditions.
The authors intended to investigate how a state space structure affects the LDS. The input data for constructing a state space were time series describing the movement at a single joint. The authors created state spaces on the basis of one or threedimensional time series for hip, knee, and ankle joints separately and used the reconstructed trajectory for the LDS assessment according to the approach briefly described in the next section.
2. Materials and Methods
2.1. Theoretical Background
A symptom of extreme sensitivity to initial conditions is the exponential rate of divergence of trajectories from their starting points which are located in a state space very close to each other. This rate, which is called the largest Lyapunov exponent (LLE), is defined as follows:where is the initial time instant and represents a distance between corresponding points on initially nearby trajectories at any time instant . From the perspective of gait analysis, a positive LLE value indicates locally unstable behavior (i.e., trajectories diverge; however, due to the presence of the attractor the distance between them cannot grow without limit). The higher the LLE, the greater the system’s sensitivity to extremely small perturbations during gait and thus the lower the LDS. The LLE, which estimates the local stability immediately after a potential perturbation, is called the shortterm largest Lyapunov exponent. The shortterm LLE is computed over a time span of a length either corresponding to one step [2–4] or one stride [1, 5, 6] using the Rosenstein algorithm [7]. The idea behind this method is that pairs of segments of the state space trajectory reconstructed on the basis of experimental data repeatedly imitate two initially neighboring trajectories, which makes it possible to trace the divergence of them. A comprehensive and precise description of the component methods leading finally to the determination of LLE was included in the work of Perc [8].
A state space structure is an important aspect of the LDS assessment. The reconstruction procedure is based on two parameters: time delay (reconstruction delay, lag) and embedding dimension . For a time series, which is composed of points , an element vector of delay coordinates of the point on the reconstructed trajectory is given by , where (the reconstructed trajectory consists of points) [9].
For a multivariate time series, which is composed of univariate time series of equal length , the reconstruction parameters are defined by a time delay vector and an embedding dimension vector . Therefore, the state space dimensionality is a sum of all , , and the delay coordinates of the point form a vector where [10, 11].
2.2. Review of the Previous Work
Concentrating on the importance of analyzing LDS and utilizing its results, some other pieces of research deserve a mention, e.g., Terrier et al. [12] investigated LDS in patients with chronic impairments after foot and ankle injuries, Bruijn et al. [13] discussed the relationship between gait stability and arm swing, and Dingwell and Marin [5], as well as England and Granata [6], analyzed the influence of gait speed on LDS. Several studies [14, 15] indicate that LDS is associated with the fall risk. Moreover, LDS may be used as a potential fall predictor to differentiate fallprone adults [16]. LDS turned out to be sensitive to agerelated degeneration [2]. The authors of [3] show that dance training might improve LDS of normal walking of the elderly. Another important aspect of LDS is its suitability for monitoring of geriatric or neurological pathologies in their early phases [2]. Comprehensive reviews of measures assessing the stability of human locomotion were prepared by Hamacher et al. [17], Bruijn et al. [18], and Van Emmerik et al. [19].
The discussion of state space structures in the context of motion data was initiated by Gates and Dingwell [20] who focused on univariate time series based on Euler angles describing rotational motion of a shoulder. The authors conclude that the comparison between outcomes for different state spaces should be made with caution; however, the trends identified in the analyzed data remain relevant. Besides, the authors do not recommend the tested PCAbased reduction of a state space dimensionality. In [21], a multivariate time series, which represented movement at hip, knee, and ankle joints in the sagittal plane, was applied to analyze quiet standing balance. As far as the total embedding dimension for a multivariate time series is concerned, Vlachos and Kugiumtzis [10] presented two modified variants (FNN1 and FNN2) of the false nearest neighbors method [22], which were adjusted to a multivariate time series. According to FNN1, the same embedding dimension is applied to all the component time series. The FNN2 method is an exhaustive algorithm. The third method proposed in [10] was based on the criterion of prediction error minimization (PEM), whereas Zhang et al. [23] suggested applying a maximal joint entropy criterion. The consequences of a fixed time delay and/or a fixed embedding dimension were investigated by van Schooten et al. [4]. Hamacher et al. [2] evaluated multiple state space definitions differing in signal type (linear acceleration and angular velocity), signal dimension (onedimensional and threedimensional), and location of an inertial sensor (trunk and forefoot). Piórek et al. [24] used a quaternionbased interpretation of body segments’ rotations and replaced a multivariate time series of Euler angles by a quaternion angle time series. Moreover, they showed a correlation between LLE values computed for time series consisting of (1) quaternion angles and (2) joint angles in a group of young individuals for hip, knee, and ankle joints in different variants of walking speed and ground inclination. The same set of experimental data was also analyzed using a new quaternionbased variant of the approximate entropy measure [25]. A systematic review of methodological approaches of the LLE quantification was prepared by Mehdizadeh [26].
It should also be pointed out that the range of applications of the LLE as a measure of sensitivity to infinitesimal changes in initial conditions goes beyond the gait analysis. For instance, Jagrič et al. [27] analyzed the irregularity in short electrocardiographic (ECG) recordings in a similar manner to predict successful defibrillation in patients with ventricular fibrillation. A higher level of irregularity was interpreted as an indicator of patients who may be subjected to effective defibrillation.
2.3. The Goal of the Research
As mentioned above, the examined state space structures were constructed on the basis of time series built of joint angles at hip, knee, and ankle joints. Three time series, which are related to the given joint, represent specific types of movement in sagittal, frontal, and transverse planes. For instance, movements at hip joint are called flexion/extension, abduction/adduction, and internal/external rotation, respectively. However, analysis of human gait focuses often on the sagittal plane to which the vast majority of the work during gait is assigned (ca. 74%, 85%, and 93% in case of hip, knee, and ankle joints, respectively) [28]. All the planes can be included using a state space based on a multivariate times series. At the previous research stage, which was extensively described in [29], the authors only used one state space that was based on multivariate time series composed of experimental data recorded in the CAREN extended environment (http://www.motekforcelink.com/product/caren/). Various experiments’ scenarios (i.e., variants of gait) were proposed which differed from each other with respect to walking speed, platform slope, and optional external perturbation. The results presented here are based on the same set of experiments. However, this time both the comparison of the LDS in three pairs of the “opposed” scenarios (i.e., gait variants which differ in one of the aforementioned aspects) and the statistical analysis were made for ankle, knee, and hip joints separately. Moreover, two additional structures of a state space were taken into consideration. The state spaces based either on a multivariate times series or on a univariate time series, which represents joint angles in the sagittal plane, will be described in Section 2.4.
Finally, three state space structures, which were constructed for each joint separately, were used to verify if the differences in LLE values between the opposed scenarios are significant for individual joints.
The authors also present a modification of the LDS computation method, i.e., an algorithm for estimation of one of its crucial parameters, embedding dimension, for the case of a multivariate gait time series, in which the parameter is not estimated for each of the component time series separately, but holistically.
The following research questions are addressed in the paper:(i)Are there any significant differences in the local dynamic stability between compared gait variants, which can be revealed using the individual state spaces?(ii)Is the predominant role of sagittal plane preserved in a state space which is based on a multivariate time series?(iii)Does the length of the time span, over which the shortterm LLE is computed, influence the difference in the local dynamic stability between compared gait variants?
2.4. The Research Procedure
The research procedure was composed of the following steps:(1)Data acquisition and preprocessing.(2)Estimation of the reconstruction parameters.(3)Trajectory reconstruction.(4)Estimation of the shortterm LLE, taking into consideration both the aforementioned time span variants.
Inspired by reports in the literature, the authors decided to incorporate three different state space structures into research. The UniS state space is reconstructed on the basis of a univariate time series describing a movement at a joint in the sagittal plane. The next two spaces—MultiFull and MultiFNN—are reconstructed on the basis of a multivariate time series which is composed of three univariate time series. Each of the series is related to movement at a given joint in one of the motion planes: sagittal, frontal, and transverse. The MultiFull space is built on the basis of three pairs of independently determined parameters , , resulting in space dimensionality . The average mutual information (AMI) method [30] was used in each case to determine time delays, whereas the false nearest neighbors (FNN) method was utilized to estimate embedding dimensions. The number of bins required by the AMI was determined according to the Sturges formula [31], and according to Kennel et al.’s example [22], the following values were assigned to the first () and the second () criterion of the FNN for designating a point as a “false” neighbor: , . However, in the case of MultiFNN, a variant of the FNN adjusted to multivariate time series was applied which takes into consideration the quota of work done during gait in individual motion planes. Based on this criterion, the planes are ordered descending as follows: sagittal, frontal, and transverse [28].
The dimensionality of MultiFNN space is determined holistically, i.e., on all three time series treated as a whole. Starting from 1, the embedding dimension is gradually increased by inserting successive elements to the vector of delay coordinates. The coordinates are taken from cyclically changed component time series , , and (the symbols “,” “,” and “” stand for sagittal, frontal, and transverse planes, respectively), with appropriate time delay for each time series. The time series describing a movement in the sagittal plane is the first one used in each cycle. The stop criterion is met when the percentage of the “false” nearest neighbors falls below a given threshold (e.g. 1%) (a neighbor of a given point in a space of dimensionality turns out to be “false,” when it is no longer a neighbor of in a space of dimensionality ). As a result, the number of coordinates taken from the time series cannot be lower than the number of coordinates from the series which, in turn, cannot be lower than the number of coordinates from the series : . By that means, the method to some extent takes into account the domination of an anteriorposterior movement in gait.
This research is a part of an extensive project carried out in cooperation with the University of the Third Age (U3A). The project focuses on elderly people that would like to remain active over the age of 65. Table 1 presents characteristics of 14 U3A students, who agreed to participate in the experiments (12 women, 2 men), including median, mean, and standard deviation (SD) values of age, height, weight, and the body mass index (BMI).

The authors state that the study has been approved by the Ethical Committee and all the subjects gave informed written consent to participate in the research after they were briefly introduced to the research protocol.
The CAREN extended system, which was used as the research environment, guaranteed not only fully immersive virtual scenery and 6 DOF motion platform but also safety and comfort of the participants. Besides, during the experiments, the walkers were under constant medical supervision. It is worth mentioning that the U3A students willingly took part in the experiments, especially when the research environment turned out to be so immersive, attractive, and safe at the same time, as the CAREN extended system is.
The participants performed six scenarios of selfpaced or fixed speed treadmill walking on level ground or on inclined platform, which are briefly presented in Table 2.

In each scenario, the subjects walked through a virtual forest. The CAREN treadmill’s selfpaced mode enables the subject to initiate gait and walk at her/his own pace which determines the instant walking speed. The treadmill adjusts then its speed to adapt to the subjects’ pace. The selfpaced mode was used in Normal, Perturbation, Up, and Down scenarios. The CAREN output data include the instant walking speed, so the values referring to the Normal scenario were averaged, thereby determining the basis (—subject’s mean preferred walking speed) for the imposed constant walking speed used in Faster and Slower scenarios ( and , respectively). A single external disruption used in the Perturbation scenario was a sudden vertical jerk of a platform with a constant amplitude for all the participants. The perturbation was induced unexpectedly by the staff member, who supervised the experiment, at time instants which were similar for all the subjects.
The participants practiced each scenario until they were able to walk comfortably. Next, three trials were recorded using the integrated Vicon motion capture system at the frequency of 100 Hz giving together 18 gait sequences for every subject (several exceptions were caused by fatigue). Every time series, which is analyzed by means of the LLE, should include the equal number of strides as well as the equal number of data points [18]. So, the sequences were long enough to contain 50 strides as the assumed final length.
The recorded data were initially filtered and optionally repaired (e.g., in view of occluded markers) using the Vicon software. The beginning of each stride was demarcated based on precisely marked occurrence of the “heelstrike” event. A stride interval varied not only across subjects but also across experiments’ scenarios. Mean and standard deviation values of the stride interval for different scenarios are included in Table 3.

Where necessary, the time series were cropped to 50 strides. Next, every stride was separately normalized using linear interpolation to contain 100 points. Subsequently, the time series were subject to estimation of reconstruction parameters. It deserves a mention that the most frequently occurring values for the dimensionality of the MultiFNN state space are 7, 8 (UniS: 4, 5; MultiFull: 14, 15). Afterward, the shortterm LLE values were calculated using the reconstructed trajectory and taking into account two variants of a time span(i)of a length equal to 50 which is equivalent to one step, i.e., a half of a stride (the shortterm LLE is then labeled by )(ii)of a length equal to 100 which is equivalent to one stride ()
The mean period parameter as the threshold for temporal separation of the nearest neighbors on two different segments of the reconstructed state space trajectory, which repeatedly imitate two initially neighboring trajectories, was estimated as the reciprocal of the mean frequency of the power spectrum [7]. In each iteration, one segment starts from the next point of the reconstructed trajectory. Taking into account the temporal separation mentioned above, the nearest neighbor of this starting point on the adjacent orbit becomes the first point of the second segment. The length of the segments was set to 1000. The Euclidean distance for each pair of two corresponding points on both segments is computed and stored to finally determine the average logarithmic divergence of the neighboring trajectories which is required to estimate the LLE value.
Finally, the results were aggregated across all the six scenarios, three state space structures, and individual joints. All the computations were performed using MATLAB and MySQL DBMS. The statistical analysis of the results focused on investigating if differences between the shortterm LLE values in three pairs of compared scenarios (NormalPerturbation, FasterSlower, UpDown) are significant for individual joints, taking into account that the shortterm LLE values were computed using different state space structures. The complete analysis was performed in Excel using the Real Statistics Resource Pack (https://www.realstatistics.com/).
3. Results
Examples of times series representing a movement in the sagittal plane, recorded for a 75yearold woman performing the Normal scenario are presented in Figure 1(a) (the red line corresponds to the left ankle joint, the blue line to the left knee joint, and the green one to the left hip joint). Each of these series was independently used for trajectory reconstruction in UniS state space. The corresponding 3D projections of the reconstructed trajectories are presented in Figure 1(b) (left ankle joint), Figure 1(c) (left knee joint), and Figure 1(d) (left hip joint).
(a)
(b)
(c)
(d)
The final results are presented as box plots. On each box, the boundary between the areas of different colors indicates the median, the × symbol denotes the mean, the edges of the box are the 25^{th} and 75^{th} percentiles, the “whiskers” indicate the most extreme values which are not outliers, i.e., the smallest value that is larger than or equal to and the largest value that is less than or equal to , where and denote the 25^{th} and 75^{th} percentiles, respectively, and the outliers are individually marked by circles. Figures 2–4 include box plots for computed in UniS, MultiFull, and MultiFNN spaces, respectively. Analogically aggregated results for are presented in Figures 5–7. Each figure comprises six subfigures corresponding to six trial scenarios: (a) Normal, (b) Perturbation, (c) Faster, (d) Slower, (e) Up, and (f) Down. Each subfigure includes separate box plots for ankle, knee, and hip joints on both sides of a body (“,” left; “,” right).
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
(a)
(b)
(c)
(d)
(e)
(f)
4. Discussion
Some general remarks, which were formulated at the earlier stage of research [29], will be briefly reminded here for clarity of further considerations. First, it is worth mentioning that the positive values of both shortterm LLE variants calculated for a uni or multivariate time series utilizing each of the considered state spaces confirm that the elderly are locally unstable during gait; however, the LLE values are lower than our results for young subjects reported in [24]. Besides, the relationship between the corresponding LLE values indicates that locally unstable behavior is gradually attenuated by the locomotor system. Secondly, comparisons between the Faster and Slower scenarios (e.g., in Figures 2(c) and 2(d), 3(c) and 3(d), 4(c) and 4(d), 5(c) and 5(d)) as well as between Up and Down scenarios (e.g., in Figures 2(e) and 2(f), 4(e) and 4(f)) suggest that both walking slower and walking downwards are more stable than their opposed variants.
The main goal of the statistical analysis is to investigate if the differences between the opposed scenarios for individual joints are significant, taking into consideration that the shortterm LLE values were computed in two variants using different state space structures.
For a given joint (ankle, knee, and hip), measure ( and ), state space (UniS, MultiFull, and MultiFNN), and the pair of opposed scenarios and (NormalPerturbation, FasterSlower, and UpDown), a null hypothesis is formulated as follows: . The null hypothesis assumes that with respect to joint , measure , and state space , the probability of an observation randomly selected from the group related to the scenario exceeding an observation randomly selected from the group representing the scenario equals the probability of an observation randomly selected from the group related to the exceeding an observation randomly selected from the group representing the .
The expected conclusions are placed in the alternative hypotheses, according to which observations in one group tend to be greater than observations in the other group, e.g., (knee, , UniS): () () (see Figures 2(c) and 2(d)).
The selection of an appropriate statistical test for verification of hypotheses should be preceded by the analysis of distribution for both measures and independently, taking into account each dataset related to a pair (scenario, joint) separately. Totally, 216 cases were considered (3 joints 2 body sides 6 scenarios 3 state spaces 2 measures = 216 datasets).
The outcome of the Shapiro–Wilk test of normality at significance level of 5% turned out to be negative in 74 from 216 cases (ca. 34%). Thus, the verification of the hypotheses requires a nonparametric test without any assumptions related to the distribution of scores.
With regard to the outcomes of the Shapiro–Wilk test, the hypotheses were verified for both measures independently using the nonparametric Mann–Whitney–Wilcoxon test at significance level of 5%. Because the number of null hypotheses is large (3 joints 3 pairs of scenarios 3 state spaces 2 measures = 54), the results of their verification are presented symbolically. Each of Tables 4–9 represents a pair (measure , state space ) and includes “=” if for a given pair of scenarios , and a given joint the null hypothesis : was not rejected. Otherwise, one of the symbols “,” “” determines a onesided alternative hypothesis (: or : ), in favor of which the corresponding null hypothesis was rejected.






The total number of rejected null hypotheses for individual pairs of scenarios is as follows: NormalPerturbation: 2, FasterSlower: 11, and UpDown: 10. Thus, there are no significant differences in the local dynamic stability between Normal and Perturbation scenarios. The platform jerks were probably too mild and too infrequent to have an impact on the shortterm LLE.
As far as two other pairs of scenarios are concerned, the high number of rejected null hypotheses indicates that the sensitivity to tiny local perturbations in both compared scenarios is different. The lower values of the shortterm LLE for the Slower scenario suggest that this variant of walking is conducive to greater stability and this observation is consistent with a common impression that elderly people, scared of fall, try to walk more carefully. During slow gait, the movement of limbs is more accurately controlled by the central nervous system.
Independently of the LDS measure and the applied state space, the significant differences between Faster and Slower scenarios are always visible in the case of a knee joint. LDS for the knee joint is lower in faster walking due to more frequent loss of body balance, more frequent bending of the knee, and the largest lateral knee movements.
A similar remark refers to a hip joint as regards the UpDown pair where the hip joint is responsible for bending and straightening the torso. The lower values of the shortterm LLE for the Down scenario are supposedly caused by a slightly rigid way of walking downwards and involuntary straightening.
The lower values for hip joints for MultiFull and MultiFNN spaces in comparison to UniS space result from using a multivariate time series and thus adding information about movement in frontal and transverse planes.
The verification outcomes for UniS and MultiFNN spaces are very consistent (in 8 from 9 cases for each measure). It can be explained by the domination of movement in the sagittal plane during gait. The MultiFNN space is based only on a subset of delay coordinates from a multivariate time series as opposed to the MultiFull space which includes all of them. Consequently, in the case of MultiFNN space, the delay coordinates, which were determined on the basis of a component time series describing a movement at a joint in the sagittal plane, play a predominant role. In both cases, in which the result of verification is different for UniS and MultiFNN ({hip, , FasterSlower}, {ankle, , FasterSlower}), the null hypothesis was not rejected for UniS, while it happened for MultiFNN which was presumably caused by additional information related to other motion planes. Greater discrepancies between UniS and MultiFull (7 different verification results) as well as between MultiFNN and MultiFull (5 different verification results) suggest that MultiFull includes redundant information.
Inspired by the close connection between UniS and MultiFNN, the authors conducted next tests to assess a statistical significance of differences between datasets related to individual state spaces (for and separately). The deviation of individual datasets from normal distributions and the lack of homogeneity of their variances were showed by the Shapiro–Wilk test and the Levene test, respectively. In consequence, a significant difference between all three datasets was identified by means of the nonparametric Kruskal–Wallis test. Afterward, the Nemenyi test and the pairwise Mann–Whitney test indicated that the following pairs of datasets: UniSMultiFull and MultiFullMultiFNN are significantly different and confirmed that the pair UniSMultiFNN is not significantly different; however, only in the case of ( values were 0.99 for the Nemenyi test and 0.50 for the pairwise Mann–Whitney test). The risk of committing a type I error was reduced by using the Bonferroni correction.
Therefore, although each examined state space reveals locally unstable behavior during gait as well as its attenuation in the current stride, a direct comparison of the shortterm Lyapunov exponents computed for different state spaces bears the risk of a wrong conclusion.
A similar remark refers to another test in which the accepted hypotheses based on were juxtaposed with their counterparts based on . In 16 of 27 (ca. 59%) cases (where each case was determined by a joint and a state space), the hypotheses were consistent. It means that the length of a time span, over which the shortterm LLE is computed, must be taken into account when comparing the results of different studies.
It should also be noted that the structure of a state space affects the computation time which comprises the following operations on a uni or multivariate time series: interpolation, estimation of time delay(s) and embedding dimension(s), reconstruction of the state space trajectory, and computation of the shortterm LLE. The mean computation time values for individual state spaces are presented in Table 10. The computer parameters were as follows: Aspire X5950, Intel Core i5, 3.2 GHz, RAM 4 GB.

As expected, the mean computation time is maximal for the MultiFull space (owing to the estimation of three embedding dimensions) and almost three times smaller for the MultiFNN space (embedding dimension was estimated only once). The mean computation time for the MultiFNN space turned out to be even slightly smaller than in the case of the UniS which is based on a univariate time series.
It is worth mentioning that some other studies which require the LLE computation and are based on a univariate time series could benefit from adding one or more time series to input data used for a state space construction. For instance, dynamic features of eye movement [32] were originally analyzed only using the values of the first derivative of horizontal eye displacements during a fixation. A twodimensional time series could be created which would also include velocity in the vertical direction while maintaining the dominant role of horizontal eye displacements. Another example is motion databased quality assessment of completion of rehabilitation exercises maintaining mobility of the hip in case of coxarthrosis [33] where the complex procedure of the LLE computation using a multivariate time series could replace or expand the results of the RQA (recurrence quantification analysis) measures' application for a time series describing a hip joint movement in the sagittal plane.
5. Conclusions
The method of computation of the local dynamic stability requires a proper reconstruction of a state space trajectory. The authors intended to investigate how a state space structure affects the shortterm LLE as the measure related to the local dynamic stability. One of the state spaces was constructed using an algorithm which estimates the embedding dimension holistically in the case of a multivariate gait time series while taking into account the quota of work done during gait in individual motion planes. It should also be mentioned here that a direct comparison of the shortterm LLE computed for different state spaces is burdened with the risk of a wrong conclusion. Due attention should also be paid to the length of the time span, over which the shortterm LLE was estimated.
The improved method of LDS assessment will be used in experiments focused on finding easytomeasure, objective biomarkers that could classify PD (Parkinson’s disease) patients in early (preclinical) stages of the disease. Identification of the first deviation from the norm in patient’s physical movement like walking, which is often unobservable to a neurologist, might help follow disease progression, make more adjusted treatment, and lead to modification of disease course.
Data Availability
The motion capture data (joint angles) used to support the findings of this study are available from the corresponding author upon request.
Disclosure
The data used in this project were obtained from the Centre for Research and Development of the PolishJapanese Academy of Information Technology (PJAIT) (http://bytom.pja.edu.pl/).
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the Statutory Research Funds of Institute of Informatics, Silesian University of Technology, Gliwice, Poland (BK/204/RAU2/2019).
References
 J. B. Dingwell and J. P. Cusumano, “Nonlinear time series analysis of normal and pathological human walking,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 10, no. 4, pp. 848–863, 2000. View at: Publisher Site  Google Scholar
 D. Hamacher, D. Hamacher, N. B. Singh, W. R. Taylor, and L. Schega, “Towards the assessment of local dynamic stability of levelgrounded walking in an older population,” Medical Engineering & Physics, vol. 37, no. 12, pp. 1152–1155, 2015. View at: Publisher Site  Google Scholar
 D. Hamacher, D. Hamacher, K. Rehfeld, and L. Schega, “Motorcognitive dualtask training improves local dynamic stability of normal walking in older individuals,” Clinical Biomechanics, vol. 32, pp. 138–141, 2016. View at: Publisher Site  Google Scholar
 K. S. van Schooten, S. M. Rispens, M. Pijnappels, A. Daffertshofer, and J. H. van Dieën, “Assessing gait stability: the influence of state space reconstruction on inter and intraday reliability of local dynamic stability during overground walking,” Journal of Biomechanics, vol. 46, no. 1, pp. 137–141, 2013. View at: Publisher Site  Google Scholar
 J. B. Dingwell and L. C. Marin, “Kinematic variability and local dynamic stability of upper body motions when walking at different speeds,” Journal of Biomechanics, vol. 39, no. 3, pp. 444–452, 2006. View at: Publisher Site  Google Scholar
 S. A. England and K. P. Granata, “The influence of gait speed on local dynamic stability of walking,” Gait & Posture, vol. 25, no. 2, pp. 172–178, 2007. View at: Publisher Site  Google Scholar
 M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,” Physica D: Nonlinear Phenomena, vol. 65, no. 12, pp. 117–134, 1993. View at: Publisher Site  Google Scholar
 M. Perc, “The dynamics of human gait,” European Journal of Physics, vol. 26, no. 3, pp. 525–534, 2005. View at: Publisher Site  Google Scholar
 M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “Reconstruction expansion as a geometrybased framework for choosing proper delay times,” Physica D: Nonlinear Phenomena, vol. 73, no. 12, pp. 82–98, 1994. View at: Publisher Site  Google Scholar
 I. Vlachos and D. Kugiumtzis, “State space reconstruction for multivariate time series prediction,” Nonlinear Phenomena in Complex Systems, vol. 11, no. 2, pp. 241–249, 2008. View at: Google Scholar
 L. Cao, A. Mees, and K. Judd, “Dynamics from multivariate time series,” Physica D: Nonlinear Phenomena, vol. 121, no. 12, pp. 75–88, 1998. View at: Publisher Site  Google Scholar
 P. Terrier, F. Luthi, and O. Deriaz, “Do orthopaedic shoes improve local stability of gait? An observation study in patients with chronic foot and ankle injuries,” BMC Musculoskeletal Disorders, vol. 14, no. 1, 2013. View at: Publisher Site  Google Scholar
 S. M. Bruijn, O. G. Meijer, P. J. Beek, and J. H. van Dieën, “The effects of arm swing on human gait stability,” Journal of Experimental Biology, vol. 213, no. 23, pp. 3945–3952, 2010. View at: Publisher Site  Google Scholar
 M. J. P. Toebes, M. J. M. Hoozemans, R. Furrer, J. Dekker, and J. H. van Dieën, “Local dynamic stability and variability of gait are associated with fall history in elderly subjects,” Gait & Posture, vol. 36, no. 3, pp. 527–531, 2012. View at: Publisher Site  Google Scholar
 F. Reynard, P. Vuadens, O. Deriaz, and P. Terrier, “Could local dynamic stability serve as an early predictor of falls in patients with moderate neurological gait disorders? a reliability and comparison study in healthy individuals and in patients with paresis of the lower extremities,” PLoS One, vol. 9, no. 6, Article ID e100550, 2014. View at: Publisher Site  Google Scholar
 T. E. Lockhart and J. Liu, “Differentiating fallprone and healthy adults using local dynamic stability,” Ergonomics, vol. 51, no. 12, pp. 1860–1872, 2008. View at: Publisher Site  Google Scholar
 D. Hamacher, N. B. Singh, J. H. van Dieën, M. O. Heller, and W. R. Taylor, “Kinematic measures for assessing gait stability in elderly individuals: a systematic review,” Journal of the Royal Society Interface, vol. 8, no. 65, pp. 1682–1698, 2011. View at: Publisher Site  Google Scholar
 S. M. Bruijn, O. G. Meijer, P. J. Beek, and J. H. van Dieёn, “Assessing the stability of human locomotion: a review of current measures,” Journal of the Royal Society Interface, vol. 10, no. 83, Article ID 20120999, 2013. View at: Publisher Site  Google Scholar
 R. E. A. van Emmerik, S. W. Ducharme, A. C. Amado, and J. Hamill, “Comparing dynamical systems concepts and techniques for biomechanical analysis,” Journal of Sport and Health Science, vol. 5, no. 1, pp. 3–13, 2016. View at: Publisher Site  Google Scholar
 D. H. Gates and J. B. Dingwell, “Comparison of different state space definitions for local dynamic stability analyses,” Journal of Biomechanics, vol. 42, no. 9, pp. 1345–1349, 2009. View at: Publisher Site  Google Scholar
 K. Liu, H. Wang, and J. Xiao, “The multivariate largest Lyapunov exponent as an agerelated metric of quiet standing balance,” Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 309756, 11 pages, 2015. View at: Publisher Site  Google Scholar
 M. B. Kennel, R. Brown, and H. D. I. Abarbanel, “Determining embedding dimension for phasespace reconstruction using a geometrical construction,” Physical Review A, vol. 45, no. 6, pp. 3403–3411, 1992. View at: Publisher Site  Google Scholar
 C.T. Zhang, J. Guo, Q.L. Ma, H. Peng, and X.D. Zhang, “Phase space reconstruction and prediction of multivariate chaotic time series,” in Proceedings of the 9th International Conference on Machine Learning and Cybernetics, pp. 2428–2433, Qingdao, China, July 2010. View at: Google Scholar
 M. Piórek, H. Josiński, A. Michalczuk, A. Świtoński, and A. Szczęsna, “Quaternions and joint angles in an analysis of local stability of gait for different variants of walking speed and treadmill slope,” Information Sciences, vol. 384, pp. 263–280, 2017. View at: Publisher Site  Google Scholar
 A. Szczęsna, “Quaternion entropy for analysis of gait data,” Entropy, vol. 21, no. 1, p. 79, 2019. View at: Publisher Site  Google Scholar
 S. Mehdizadeh, “The largest Lyapunov exponent of gait in young and elderly individuals: a systematic review,” Gait & Posture, vol. 60, pp. 241–250, 2018. View at: Publisher Site  Google Scholar
 T. Jagrič, M. Marhl, D. Štajer et al., “Irregularity test for very short electrocardiogram (ECG) signals as a method for predicting a successful defibrillation in patients with ventricular fibrillation,” Translational Research, vol. 149, no. 3, pp. 145–151, 2007. View at: Publisher Site  Google Scholar
 J. J. Eng and D. A. Winter, “Kinetic analysis of the lower limbs during walking: what information can be gained from a threedimensional model?” Journal of Biomechanics, vol. 28, no. 6, pp. 753–758, 1995. View at: Publisher Site  Google Scholar
 H. Josiński, A. Świtoński, A. Michalczuk, P. Grabiec, M. Pawlyta, and K. Wojciechowski, “Analysis of chaotic behaviors in gait of the elderly using the CAREN extended system,” in Proceedings of the 20th International Conference on eHealth Networking, Applications and Services (2018 IEEE Healthcom), Ostrava, Czech, August 2018. View at: Google Scholar
 B. Henry, N. Lovell, and F. Camacho, “Nonlinear dynamics time series analysis,” in Nonlinear Biomedical Signal Processing: Dynamic Analysis and Modeling, WileyIEEE Press, New York, NY, USA, 2001. View at: Google Scholar
 H. A. Sturges, “The choice of a class interval,” Journal of the American Statistical Association, vol. 21, no. 153, pp. 6566, 1926. View at: Publisher Site  Google Scholar
 K. Harężlak, “Eye movement dynamics during imposed fixations,” Information Sciences, vol. 384, pp. 249–262, 2017. View at: Publisher Site  Google Scholar
 H. Josiński, A. Michalczuk, R. Mucha, A. Świtoński, A. Szczęsna, and K. Wojciechowski, “Analysis of human motion data using recurrence plots and recurrence quantification measures,” in Proceedings of the 8th Asian Conference on Intelligent Information and Database Systems (ACIIDS 2016), LNAI 9622, pp. 397–406, Springer, Da Nang, Vietnam, March 2016. View at: Google Scholar
Copyright
Copyright © 2019 Henryk Josiński et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.