Assessment of Local Dynamic Stability in Gait Based on Univariate and Multivariate Time Series

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 short-term 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 three-dimensional, which enables to construct a state space based on a uni- or multivariate time series. Furthermore, two variants of the short-term LLE are present in the literature which differ in length of a time span, over which the short-term LLE is computed. Both a state space structure and the consistency of the observations based on values of both short-term LLE variants were analyzed using time series representing the joint angles at ankle, knee, and hip joints. The short-term 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 short-term 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.


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]. e 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. e method is based on a trajectory in a state space which is reconstructed from time series generated by a dynamical system. e 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. e authors intended to investigate how a state space structure affects the LDS. e input data for constructing a state space were time series describing the movement at a single joint. e authors created state spaces on the basis of one-or three-dimensional 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.

eoretical 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. is rate, which is called the largest Lyapunov exponent (LLE), is defined as follows: where t 0 is the initial time instant and d(t) represents a distance between corresponding points on initially nearby trajectories at any time instant t. 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). e higher the LLE, the greater the system's sensitivity to extremely small perturbations during gait and thus the lower the LDS. e LLE, which estimates the local stability immediately after a potential perturbation, is called the short-term largest Lyapunov exponent. e short-term LLE is computed over a time span of a length either corresponding to one step [2][3][4] or one stride [1,5,6] using the Rosenstein algorithm [7]. e 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. e reconstruction procedure is based on two parameters: time delay τ (reconstruction delay, lag) and embedding dimension m. For a time series, which is composed of N points x 1 , x 2 , . . . , x N , an m-element vector of delay coordinates of the point x i on the reconstructed trajectory is given by For a multivariate time series, which is composed of K univariate time series of equal length N, the reconstruction parameters are defined by a time delay vector [τ 1 , τ 2 , . . . , τ K ] and an embedding dimension vector [m 1 , m 2 , . . . , m K ]. erefore, the state space dimensionality m is a sum of all m k , k � 1, 2, . . . , K, and the m delay coordinates of the point [10,11].

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 fall-prone adults [16]. LDS turned out to be sensitive to age-related degeneration [2]. e 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]. e 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. e 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. e FNN2 method is an exhaustive algorithm. e 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. e 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 (one-dimensional and threedimensional), and location of an inertial sensor (trunk and forefoot). Piórek et al. [24] used a quaternion-based 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. e same set of experimental data was also analyzed using a new quaternion-based 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. Computational and Mathematical Methods in Medicine 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. e 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. ree 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. e 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. e 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. e 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. e 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 short-term LLE is computed, influence the difference in the local dynamic stability between compared gait variants?

e Research
Procedure. e 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 short-term 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. e UniS state space is reconstructed on the basis of a univariate time series describing a movement at a joint in the sagittal plane. e 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.
e MultiFull space is built on the basis of three pairs of independently determined parameters (m k , e 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. e 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 (Rtol) and the second (Atol) criterion of the FNN for designating a point as a "false" neighbor: Rtol � 15, Atol � 2. 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]. e 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. e coordinates are taken from cyclically changed component time series S, F, and T (the symbols "S," "F," and "T" stand for sagittal, frontal, and transverse planes, respectively), with appropriate time delay (τ S , τ F , τ T ) for each time series. e time series describing a movement in the sagittal plane is the first one used in each cycle. e 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 P in a space of dimensionality d turns out to be "false," when it is no longer a neighbor of P in a space of dimensionality d + 1). As a result, the number of coordinates taken from the S time series (m S ) cannot be lower than the number of coordinates from the F series (m F ) which, in Computational and Mathematical Methods in Medicine turn, cannot be lower than the number of coordinates from the T series (m T ): m S ≥ m F ≥ m T . By that means, the method to some extent takes into account the domination of an anterior-posterior movement in gait.
is research is a part of an extensive project carried out in cooperation with the University of the ird Age (U3A). e 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).
e 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. e 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. e participants performed six scenarios of self-paced 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. e CAREN treadmill's self-paced mode enables the subject to initiate gait and walk at her/his own pace which determines the instant walking speed. e treadmill adjusts then its speed to adapt to the subjects' pace. e self-paced mode was used in Normal, Perturbation, Up, and Down scenarios.
e CAREN output data include the instant walking speed, so the values referring to the Normal scenario were averaged, thereby determining the basis (PWS s -subject's mean preferred walking speed) for the imposed constant walking speed used in Faster and Slower scenarios (1.2 · PWS s and 0.8 · PWS s , 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. e perturbation was induced unexpectedly by the staff member, who supervised the experiment, at time instants which were similar for all the subjects. e 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.
e recorded data were initially filtered and optionally repaired (e.g., in view of occluded markers) using the Vicon software. e beginning of each stride was demarcated based on precisely marked occurrence of the "heel-strike" 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 short-term 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 short-term LLE is then labeled by λ S0.5 ) (ii) of a length equal to 100 which is equivalent to one stride (λ S1 ) e 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. e length of the segments was set to 1000. e 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

Results
Examples of times series representing a movement in the sagittal plane, recorded for a 75-year-old 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. e corresponding 3D projections of the reconstructed trajectories are presented in Figure 1 e 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 Q 1 -1.5 · (Q 3 − Q 1 ) and the largest value that is less than or equal to Q 3 + 1.5 · (Q 3 − Q 1 ), where Q 1 and Q 3 denote the 25 th and 75 th percentiles, respectively, and the outliers are individually marked by circles.

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 short-term 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 (λ S0.5 > λ S1 > 0) 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.
e 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 short-term LLE values were computed in two variants using different state space structures.
e null hypothesis assumes that with respect to joint J, measure M, and state space S, the probability of an observation randomly selected from the group related to the Sc 1 scenario exceeding an observation randomly selected from the group representing the Sc 2 scenario equals the probability of an observation randomly selected from the group related to the Sc 2 exceeding an observation randomly selected from the group representing the Sc 1 .
e 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., H 1 (knee, λ S0.5 , UniS): P(Faster > Slower) > P(Slower > Faster) (see Figures 2(c) and 2(d)). e selection of an appropriate statistical test for verification of hypotheses should be preceded by the analysis of distribution for both measures λ S0. 5 and λ S1 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). e 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%). us, 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 M, state space S) and includes "�" if for a given pair of scenarios Sc 1 , Sc 2 and a given joint J the null hypothesis H 0 (J, M, S): P(Sc 1 > Sc 2 ) � P(Sc 2 > Sc 1 ) was not rejected. Otherwise, one of the symbols "<," ">" determines a one-sided alternative hypothesis (H 1 (J, M, S): P(Sc 1 > Sc 2 ) < P(Sc 2 > Sc 1 ) or H 1 (J, M, S): P(Sc 1 > Sc 2 ) > P(Sc 2 > Sc 1 )), in favor of which the corresponding null hypothesis was rejected. e total number of rejected null hypotheses for individual pairs of scenarios is as follows: Normal-Perturbation: 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. e lower values of the short-term 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 Up-Down pair where the hip joint is responsible for bending and straightening the torso. e lower values of the short-term LLE for the Down scenario are supposedly caused by a slightly rigid way of walking downwards and involuntary straightening.
e lower values for hip joints for MultiFull and Mul-tiFNN spaces in comparison to UniS space result from using a multivariate time series and thus adding information about movement in frontal and transverse planes. e 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. e 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, λ S0.5 , Faster-Slower}, {ankle, λ S1 , Faster-Slower}), 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 λ S0. 5 and λ S1 separately). e 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,  Computational and Mathematical Methods in Medicine  (f ) Down. All values are positive. e platform jerks seem to have no effect. e walking speed seems to affect the knee joint primarily, while the platform inclination seems to have the greatest impact on the hip joint.    All values are positive. e platform jerks seem to have no effect. e walking speed seems to affect the knee joint primarily, while the platform inclination seems to have the greatest impact on the hip joint.   (f ) Down. All values are positive but smaller than the corresponding λ S0.5 values: locally unstable behavior is gradually attenuated by the locomotor system. e platform jerks seem to have no effect. e walking speed seems to affect the knee joint primarily, while the platform inclination seems to have the greatest impact on the hip joint.   Up; (f ) Down. All values are positive but smaller than the corresponding λ S0.5 values: locally unstable behavior is gradually attenuated by the locomotor system. e platform jerks seem to have no effect. e walking speed seems to affect the knee joint primarily, while the platform inclination seems to have the greatest impact on the hip joint. (e) Up; (f ) Down. All values are positive but smaller than the corresponding λ S0.5 values: locally unstable behavior is gradually attenuated by the locomotor system. e platform jerks seem to have no effect. e previously indicated dependencies for individual joints are no longer visible.  A similar remark refers to another test in which the accepted hypotheses based on λ S0.5 were juxtaposed with their counterparts based on λ S1 . 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 short-term 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 short-term LLE. e mean computation time values for individual state spaces are presented in Table 10.
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). e 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 two-dimensional 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 data-based 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.

Conclusions
e method of computation of the local dynamic stability requires a proper reconstruction of a state space trajectory. e authors intended to investigate how a state space structure affects the short-term 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 short-term LLE was estimated. e improved method of LDS assessment will be used in experiments focused on finding easy-to-measure, 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 e motion capture data (joint angles) used to support the findings of this study are available from the corresponding author upon request. Disclosure e data used in this project were obtained from the Centre for Research and Development of the Polish-Japanese