banner

Blog

May 18, 2023

Machine learning based estimation of dynamic balance and gait adaptability in persons with neurological diseases using inertial sensors

Scientific Reports volume 13, Article number: 8640 (2023) Cite this article

576 Accesses

2 Altmetric

Metrics details

Poor dynamic balance and impaired gait adaptation to different contexts are hallmarks of people with neurological disorders (PwND), leading to difficulties in daily life and increased fall risk. Frequent assessment of dynamic balance and gait adaptability is therefore essential for monitoring the evolution of these impairments and/or the long-term effects of rehabilitation. The modified dynamic gait index (mDGI) is a validated clinical test specifically devoted to evaluating gait facets in clinical settings under a physiotherapist's supervision. The need of a clinical environment, consequently, limits the number of assessments. Wearable sensors are increasingly used to measure balance and locomotion in real-world contexts and may permit an increase in monitoring frequency. This study aims to provide a preliminary test of this opportunity by using nested cross-validated machine learning regressors to predict the mDGI scores of 95 PwND via inertial signals collected from short steady-state walking bouts derived from the 6-minute walk test. Four different models were compared, one for each pathology (multiple sclerosis, Parkinson's disease, and stroke) and one for the pooled multipathological cohort. Model explanations were computed on the best-performing solution; the model trained on the multipathological cohort yielded a median (interquartile range) absolute test error of 3.58 (5.38) points. In total, 76% of the predictions were within the mDGI's minimal detectable change of 5 points. These results confirm that steady-state walking measurements provide information about dynamic balance and gait adaptability and can help clinicians identify important features to improve upon during rehabilitation. Future developments will include training of the method using short steady-state walking bouts in real-world settings, analysing the feasibility of this solution to intensify performance monitoring, providing prompt detection of worsening/improvements, and complementing clinical assessments.

Healthy people easily adapt to environmental perturbations: they recover from external impulses and learn walking dynamics in different contexts. Locomotion requires a continuous modulation of the coordination within and between limbs to maintain the equilibrium during progression and to accommodate demands coming from the real-world environment (e.g., passing from straight-line walking to walking over an obstacle)1. To achieve this goal, the central nervous system exploits the sensorimotor control that conveys and integrates visual, proprioceptive, and vestibular sensory inputs to detect deviations from the upright posture and generate the appropriate muscle response to correct these deviations1, 2. Because sensorimotor feedback is altered in persons with neurological disorders (PwND), dynamic balance (i.e. the ability to maintain balance while moving the body) and its adaptability to environmental changes are impaired3, leading to a high risk of falling. In fact, the most commonly found predictors of falls are gait impairments and balance disorders4. Thus, many dynamic balance assessment tools have been developed, including clinical scales5, quantitative gait markers6, posturography7, fall prevention protocols8 and single- and dual-task tests9. Regarding clinical scales, the Berg balance scale10, the MiniBESTest10, the timed up and go (TUG)11, and the modified dynamic gait index (mDGI)12, 13 are among the most commonly used in clinical practice to measure balance in PwND. Compared to the Berg scale, which does not assess dynamic balance during locomotion, and the TUG and MiniBESTest, which evaluate this aspect in a few functional tasks (i.e., one task in TUG and five tasks in MiniBESTest), the mDGI is specifically devoted to measuring an individual's capacity to maintain their balance and adapt their gait in the presence of various environmental demands, essential to perform daily-life locomotor activities without falling. The mDGI assesses eight facets of gait; it evaluates distance, temporal, ambient, terrain, physical load, attention, and postural transitions, which are representative of the environmental demands of a walking human. The highest possible score on the mDGI is 64 points. The mDGI has been extensively validated in subjects with mobility impairments14, 15 and in different pathological cohorts such as stroke (ST), vestibular dysfunction, multiple sclerosis (MS), traumatic brain injuries, and Parkinson's disease (PD)12, 13. Recently, Torchio et al. provided cut-off values to identify PwND with no or minimal risk of falls (mDGI score \(> 49\)) and PwND with high risk of falls (mDGI score \(\le 29\))16. The minimal detectable change in mDGI is 5 points in PwND14, 15.

A previous study of the eight items of the mDGI found a strong correlation between the mDGI score and instrumentally-derived indexes17, suggesting that kinematic determinants can predict dynamic balance. Despite these encouraging results, the assessment of equilibrium and gait adaptation to different environments is still performed only by specialised personnel in clinical settings and requires specific tools, such as obstacles or stairs. Such requirements, in turn, limits the number of possible assessment sessions, while a frequent monitoring would better track the possible change caused by pathology course or rehabilitation/pharmacological interventions. In this respect, the opportunity of predicting dynamic balance and locomotor adaptability by using wearable sensors during repeated short (i.e. 10 strides) steady-state walking bouts, easily performed during daily living also by PwND (e.g. during a stroll alone or with a caregiver), would represent a first step in increasing the monitoring frequency and complement the periodic in-clinic evaluations18, 19.

As a preliminary test of this opportunity, this study implemented an interpretable machine learning (ML) model targeting the mDGI score by using inertial measurements units (IMUs) to collect data during short steady-state walking bouts of a 6-minute walk test (6MWT). After statistically confirming the association between instrumental variables and the mDGI score, we hypothesised that ML models based on such variables can predict dynamic balance and gait adaptability (i.e., mDGI scores) in PwND, including multiple sclerosis (MS), Parkinson's disease (PD), or stroke (ST). If this hypothesis is confirmed, the results of the present study could offer a clear starting for assessing the feasibility of this approach in daily life. Determinants of gait temporal aspects, intensity, smoothness, stability, symmetry, and regularity were extracted from a 3-IMU set-up. Then, a regularized Elastic-Net (EN) regression was developed using a nested cross-validation approach. This pipeline was repeated for the multi-pathological cohort (MP) and the single pathology cohorts (SP\(_{MS}\), SP\(_{PD}\), and SP\(_{ST}\)). Furthermore, we integrated the best performing solution model with a Shapley-values based explainability technique (SHapley Additive exPlanations, SHAP20, 21).

The pooled cohort—95 PwND, F = 43, median age = 60 years [IQR = 19]—resulted in a median 6MWT score of 346 m [IQR = 21] and a median mDGI of 46 points [IQR = 21] (Table 1). For the MS, PD and ST groups, the median disease duration corresponded to 19, 4, and 7 years, respectively. The median 6MWT was 316 [IQR = 182], 332 [IQR = 194] and 372 [IQR = 152] m, and the median mDGI score was 40 [IQR = 21], 46 [IQR = 29] and 50 [IQR = 14], respectively. Following the cut-off values defined by Torchio et al.16, 39 of the 95 participants (41%) had a low/minimal fall risk (mDGI score > 49), while 17 (18%) had a high fall risk (mDGI score \(\le\) 29). The number (percentage) of individuals with small/minimal fall risk was 18 (35%) for MS, 17 (59%) for PD, and 4 (27%) for the ST group. The number (percentage) of individuals at high fall risk was 7 (14%) for MS, 3 (10%) for PD, and 4 (27%) for the ST group.

Preliminary univariate Spearman's correlation showed that the 6MWT was significantly associated with the mDGI (\(p < 0.001\)) for the pooled and the individual cohorts (Table 2). In the pooled cohort, longer durations of stride and double support times resulted in a decrease in the mDGI. A positive association with dynamic balance was found for swing and single support time. The same trends were found in the MS and PD groups. Conversely, the ST group exhibited this behavior only for \(T_{d,support}\) and \(T_{s,support}\). Gait regularity (in the form of step and stride regularity) was positively associated with the mDGI score (\(p < 0.01\)) for all axes in the pooled and MS cohorts, and on the vertical axes for the PD group. Moreover, mediolateral (ML) stride regularity was positively correlated with the mDGI in PD patients. Smoothness of walking, expressed as an improved harmonic ratio (iHR), resulted in a positive relationship with the mDGI for all three axes in the MS, PD, and pooled cohorts and in the vertical (VT) direction for patients with ST. Root mean squared acceleration values were significantly associated with the mDGI for all groups and all axes. Lyapunov exponents in the AP direction of the pooled and MS cohorts exhibited a negative association with the outcome (\(p < 0.01\)). The single-pathology group, specifically the MS and PD groups, demonstrated an inversely significant relationship between normalised vertical jerk and the mDGI values. Furthermore, in people with ST and MS, the mDGI was inversely associated with age and disease duration. The presence of either bilateral or monolateral support significantly reduced the mDGI performance, particularly in the MS and PD groups.

In the MP model, the optimized EN model resulted in a regressor with a median absolute validation error of 4.07 points [IQR = 0.07] across the 95 leave-one-subject-out (LOSO) outer folds. The aggregated test absolute error was equal to 3.58 points [IQR = 5.41] with a \(R^{2} = 0.81\) (Fig. 1, panel A). Thus, 76% of the predictions were within the mDGI's minimal detectable change of 5 points. By subdividing MP predictions into the three sub-cohorts (Fig. 1, panel C) and computing correlation coefficients within the groups, an \(R^{2}\) of 0.79, 0.85, and 0.78 was obtained respectively for the MS, PD, and ST sub-cohorts, respectively. Similarly, individual group test errors resulted in 2.92 [IQR = 5.21], 4.11 [IQR = 6.04], and 3.22 [IQR = 4.41] points for the the MS, PD and ST sub-cohorts, respectively. Single-pathology (SP) models resulted in an overall (aggregated) test error of 4.91 points [IQR = 5.09] and an \(R^{2}\) = 0.76 (Fig. 1, panel B). Namely, individual SP models resulted in a median absolute test error of 4.64 points [IQR = 3.86] for the MS group, 4.76 points [IQR = 9.47] for the PD group, and 5.64 points [IQR = 2.92] for the ST group (Fig. 1, panel D). No significant differences (Wilcoxon signed-rank tests) were found between the combined SP models and the MP model or between the individual SP models and their respective groupings of the MP predictions. Models trained using data not derived from inertial sensors resulted in reduced accuracy in both the MP and SP cases. In particular, for the MP and the combined SP models, the \(R^{2}\) decreased from 0.81 to 0.56 and from \(R^{2}\) = 0.76 to 0.46, respectively (Supplementary Fig. 1). Applying a threshold of 5 to the MP model derived from clinical data only, a classification accuracy of 49% was obtained. Similarly, after selecting individual pathologies from the MP model, the \(R^{2}\) decreased from 0.79, 0.85, and 0.78–0.55, 0.75 and 0.22, respectively for the MS, PD, and ST models (Supplementary Fig. 1). In order to assess the model's predictive performance without the 6MWT score (estimate of walking speed), the pipeline was repeated, removing this information from the training data (Supplementary Fig. 2). This approach resulted in a median absolute error of 5.01 [IQR = 5.6] for the MP model (\(R^{2}\) = 0.78), with no significant differences from the MP model with the 6MWT score included (Mann–Whitney, \(p > 0.05\)).

Test predictions plot for the MP model (panel A) and the ST models (panel B). In both panels, the marker color indicates the type of assistive support while the marker type indicates the pathology. The dashed grey line \(y = x\) indicates ideal predictions, while the shaded grey rectangle around it represents a boundary of \(\pm 5\) points (MDC). Below, with violin-plots (and superimposed swarm-plots) the test absolute error distributions are represented. The MP model (MP) test errors were grouped according to the pathology (\(MP_{MS}, \; MP_{PD}, \; MP_{ST}\)) and provided in panel (C), whilst the SP models (\(SP_{combined}\)) combined test errors and the respective individual models are reported in panel (D).

In the MP model (Fig. 2, panel A), the presence of a bilateral or mono-lateral assistive device was the factor that showed the strongest negative association with the outcome. Lyapunov coefficients, calculated on all three axes, were negatively correlated with the mDGI prediction. Additionally, longer \(T_{stride}\) and \(T_{d,supp}\) was related to a reduced dynamic balance (i.e., lower mDGI score). Walking endurance (i.e., the distance covered in 6 min, as measured with the 6MWT), movement intensity (i.e., Acc.RMS), stride regularity computed on the AP/VT axes, and duration of the single support phase \(T_{s,supp}\) were positively associated with the mDGI scores and, consequently, with a strong ability to maintain dynamic balance. In the \(SP_{MS}\) model (Fig. 2, panel B), the presence of an assistive device, longer \(T_{stride}\) and longer \(T_{d,supp}\) were associated with lower mDGI values. Furthermore, higher values of normalized jerk values and higher Lyapunov exponents in the VT/AP directions resulted in reduced dynamic balance. Increased stride regularity on the three axes, longer \(T_{s,supp}\) and higher acceleration values on the AP axes were all related to higher mDGI in the MS group.

Regression \(\beta\) coefficients of the Elastic-Net models for the MP (panel A), \(SP_{MS}\) (panel B), \(SP_{PD}\) (panel C) and \(SP_{ST}\) (panel D). Regression coefficients obtained from the models of the outer leave-one-subject-out split are aggregated. Therefore each variable weight is indicated via a box-plot instead than a bar-plot to account for this variability. Features are sorted in ascending order of median \(\beta\) values across the folds.

In the PD group, the presence of an assistive device was the variable that showed the most negative association with the outcome (Fig. 2, panel C). Lower mDGI scores in people with PD were also negatively impacted by higher values of \(CV_{T,step}\), \(CV_{T,stride}\) (i.e., step and stride variability) and \(T_{stride}\). Conversely, the Lyapunov of the ML/AP axes maintained a strong negative effect on the mDGI in the PD model. People with PD were almost unaffected by changes in Lyapunov coefficients in the vertical direction. Faster acceleration RMS values and higher stride regularity in all three directions were associated with higher mDGI values in the PD model and higher gait symmetry quantified by \(iHR_{AP}\) and \(iHR_{ML}\).

Notably, in the ST model, the disease duration and the presence of assistive devices were the variables most negatively associated variables with the mDGI. Lyapunov coefficients in all directions affected the balance of participants post-stroke (Fig. 2, panel D). As in the case of participants with PD, the ST cohort exhibited a strong negative association between \(CV_{T,step}\) and \(CV_{T,stride}\) and the mDGI score. Stride regularity and RMS acceleration values on the AP axes were the features most strongly correlated with the mDGIs of participants post-stroke.

The SHAP values of the MP model (Fig. 3) indicated that the strongest mDGI predictor is the 6MWT score. With respect to assistive devices, three distinctive effects are observed: no assistive device (blue), monolateral assistance (purple), and bilateral assistance (pink). Compared to monolateral assistance, the use of bilateral assistive devices resulted in stronger limitations of dynamical capabilities. In conjunction with related regression coefficients, high values of stride regularity on the VT and AP axes positively impacted the mDGI. Moreover, higher values of anteroposterior RMS acceleration contributed to increasing predicted values. The SHAP values confirmed that Lyapunov coefficients on all three axes provided a strong contribution to the mDGI prediction with a negative effect (i.e., higher Lyapunov coefficients are predictive of lower predicted mDGI). Furthermore, faster strides and longer single support periods positively influenced the model predictions. The feature contributions to the prediction of patients with an absolute prediction error larger than 10 points were compared to the contributions of the whole cohort. No systematic differences in feature importance were found between the misclassified and the correctly predicted patients (Fig. 3, panel A and B).

SHAP values of the MP model were computed for all features coalitions repeatedly for each of the outer test folds and aggregated together. In Panel (A and B) they are presented ordered according to the \(mean|Shapley_{i}|\) across all test patients. Respectively, in panel (A) all instances are reported whilst on panel (B) only instances mis-classified by more than 10 mDGI points are shown. On the x-axes, the individual features normalize contribution on the mDGI values is reported, with the model prediction for each patient being the sum across all features contribution.

In this study, we trained, optimised, and cross-validated an elastic net regression model capable of predicting the mDGI score from a 3-IMU setup used during short steady-state walking bouts (10 strides each) extracted from the 6MWT. The median [IQR] error was 3.58 [5.38] points, which is lower than the mDGI MDC (i.e., 5 points). Elastic net was chosen as the adopted model because it is a regularised multivariate linear regression and one of the simplest ML models. In particular, regularisation was necessary to address the high dimensionality of the predictors compared to the moderate sample size and to prevent overfitting. Regularisation was chosen instead of statistical feature screening methods to avoid train-test contamination. Lastly, linear regressions such as elastic net also allow for the evaluation of regression coefficients and hence a model-based feature importance assessment. Comparing the developed models with algorithms trained only on features not derived from IMUs (6MWT score and presence of assistive devices) showed that IMU-based features are crucial in the assessment of dynamic balance and gait adaptability. In particular, comparing the classification performances obtained after applying an MDC-based threshold showed how using IMU-related features improved accuracy by 27% (from 49 to 76%), resulting improved classification for 26 out of 95 patients. These encouraging results represent a strong starting point for future studies focusing on the feasibility of this approach in real-world contexts, with the final goal of increasing monitoring frequency and complementing in-clinic assessment.

The multipathology model (MP), trained with a cohort of people with MS, PD, and ST, resulted in the best-performing solution. Unexpectedly, single-pathology (SP) models did not outperform the MP models. The reason may be twofold. First, the PD and ST cohorts were considerably smaller than the MS cohort. Therefore, in people with ST and PD, predictions of models trained on the MP cohort may have been influenced by the presence, in the training set, of the MS patients, who comprised almost 50% of the participants. Second, different neurological impairments could lead to different gait pattern alterations, resulting in comparable effects in terms of dynamic balance and gait adaptability, which may still be interpretable by a MP model. The latter was further confirmed by the presence of significant differences between groups in gait timing (double and single support time, step and swing time) and gait regularity (both step and stride) on all three axes (Supplementary Table 1). This could overcome the need to train and deploy individual models for specific pathologies. In this study, we included several independent variables to predict the mDGI score, including information about participant demographics, disease duration, and type of assistive device, which are known to be associated with dynamic balance17, 22 and fall-risk assessment23, 24. In addition, we included a set of clinical and wearable sensor-based features descriptive of gait spatio-temporal aspects (i.e. 6MWT score, duration of stride, step, single support and double support phases) and gait quality (i.e., movement intensity and smoothness, stride/step regularity and variability, gait symmetry, and local dynamic instability)25,26,27, which are known to be impaired in PwND3, 19, 28,29,30,31,32. Alterations in gait spatio-temporal parameters have been well-documented in MS33, PD34 and ST35. In contrast, the quantification of factors related to gait quality has only recently attracted interest; compared to the spatio-temporal factors, the aspects are (1) more robust to differences in test settings27, (2) more sensitive to mild impairments36, (3) more responsive to rehabilitation effects37, and (4) more strongly associated with patient-reported walking ability38.

Linear models can associate importance levels with the features, retaining transparency in model building (e.g., \(\beta\) coefficient in regressions and test statistics in group comparisons). However, such techniques do not yet provide a patient-wise estimate of the features’ contribution to the predictions. For this reason, embedding the elastic net regularised regression with SHAP results in translatable models, fostering the trust of clinical operators and the interpretability of mistakes. Regarding the most important features according to SHAP (Fig. 3), the 6MWT score had the highest impact on the prediction of the mDGI score and confirmed the strong positive correlation between the 6MWT and clinical balance measures (i.e., DGI, Berg Balance Scale, MiniBESTest, and TUG test) previously found in persons with MS39, PD40 and ST41. Although the 6MWT assesses walking endurance, its score can be considered as an estimate of sustainable gait speed over long periods of time. Thus, these result enforce the importance of gait speed as "sixth vital sign"42 that is associated with several health-related considerations, including dynamic balance and falls43. The use of an assistive device was the second most important feature in mDGI score prediction. This result was expected, since the level of assistance represents one of the mDGI clinical sub-scores. However, it should be highlighted that a previous study on people with neurological disorders17, showed that the IMU-based measures describing locomotion during the mDGI tasks (1) indicated greater impairment in persons using assistive devices than in those who do not, and (2) significantly correlated with mDGI and TUG scores. Combined, these results strengthened the association between walking aids and impaired gait patterns and, as a consequence of this correlation, poorer dynamic balance. Furthermore, since assistive devices are often adopted whenever a dynamic balance disorder occurs, it is reasonable to speculate that assistive devices are associated with a higher fall risk.

Six of the 10 most important contributors were descriptive of gait quality (i.e., regularity, intensity, and dynamic instability), further highlighting the importance of these variables in the characterization of locomotion and dynamic balance. The role of these parameters in the mDGI prediction persisted despite the strong contribution of walking speed (6MWT), suggesting that, despite their dependency on gait velocity44, these metrics provide different additional information. The latter demonstrates that in order to retain sufficient dynamic balance, a fine control of body links’ position or velocity is a crucial necessity, complementary to having a consistent walking speed. Similarly, Carpinella et al. found statistically significant correlations between gait quality metrics and the mDGI—Item 8 subscore, even after correcting for gait speed45.

Among the gait quality metrics, high stride regularity computed from VT and AP trunk accelerations provided positive contributions to the mDGI prediction. This result is consistent with findings in frail older adults, whose balance and locomotor impairments were characterized by higher inter-stride trunk acceleration variability (i.e., lower stride regularity) in the AP and VT directions46. The authors speculated that this finding could be related to the inability of frail subjects to generate a consistent and well-balanced propulsion of the body in the sagittal plane. This hypothesis could also apply for PwND. Significant correlations between stride regularity metrics and balance measures have already been found in previous studies of the elderly47 and PwND45, 48, 49, confirming that gait regularity is a potential rehabilitation target for improving dynamic balance and gait adaptability. Movement intensity in the AP direction (Acc. RMS\(_{AP}\)) was also a strong positive contributor to the mDGI predicted value—i.e., the higher the trunk AP acceleration—the better the dynamic balance and the gait adaptability. This result is supported by previous studies that find a reduction of movement intensity in persons with PD, ST, and MS50,51,52 (characterized by balance deficits) compared to healthy controls. Considering the strong correlation between trunk acceleration and walking speed44, a preliminary explanation for the above result could be that the trunk acceleration is positively associated with dynamic balance because gait velocity is positively associated with dynamic balance. However, a reduction in movement intensity has been previously found in people with MS compared to healthy subjects walking at a comparable speed53. Moreover, in the present model, Acc. RMS\(_{AP}\) represented the fourth most important factor, despite the strong contribution of gait speed (6MWT). A second hypothesis can thus be formulated: people with neurological diseases may minimize upper body motion (i.e., trunk acceleration) during walking in an effort to compensate for lower limb impairments (by damping perturbations of the impaired limb) and maintain dynamic stability50. The local dynamic instability of gait also provided a strong contribution to the mDGI prediction. This result indicated that poorer dynamic balance is associated with higher local dynamic instability, i.e., the locomotor system has greater difficulty coping with small spontaneous perturbations caused by environmental factors (e.g., uneven surfaces) or internal factors (e.g., neurocontrol errors)54. A previous study of early-stage, nondisabled people with MS found that the AP Lyapunov exponent significantly correlated with clinical balance scales such as the Fullerton Advanced Balance Scale and the TUG28. The present findings reinforce this result and demonstrate that local dynamic instability of gait is associated with dynamic balance; this finding is also true for severely impaired persons with different neurological conditions and can thus be considered a valid quantitative measure of balance control during walking. Regarding temporal features of gait, only stride and single support duration have been found among the 10 most important factors for mDGI prediction. In particular, prolonged stride (i.e. reduced cadence) and shorter single support phase (i.e. longer double support duration) are associated to poorer dynamic balance, in line with the typical protective cautious strategy adopted to compensate for balance impairment and maintain a stable gait45. As shown in Fig. 3, gait symmetry, as measured by improved harmonic ratios (iHRs), provided only a minor contribution to the mDGI score estimation, suggesting that this aspect has a stronger association with the energy efficiency of gait than with dynamic balance55.

In summary, this work demonstrates that measures descriptive of the short steady-state walking bouts composing the 6MWT include information sufficient for predicting dynamic balance and gait adaptability to external demands and indicate that the two assessments are not mutually independent. Furthermore, since the mDGI measures the walk adaptability to different environments, it is reasonable that better quality of gait results the focal point of a prompt adaptation, and thus of a healthier walk. These results may inform rehabilitation by indicating the most important features to be addressed during balance training. The deployed models were trained on data recorded during the central 10 strides of a 30-m hospital hallway repeatedly travelled by the subject for 6 min under the supervision of a physiotherapist, as required by the 6MWT. Although this type of walking is not equivalent to that used in everyday life, these results provide a promising basis for future studies that test the possibility of extending the current model validity to more general free walks, including short steady-state walking bouts extracted from a typical stroll. Recently developed technologies based on portable devices already permit the acquisition and processing of data during unconstrained walking and provide results immediately following the recording56 or even in real-time57, 58. Regarding the specific topic of this study, the implementation of a dedicated smartphone app (one that embeds the proposed model with all data-processing procedures) able to automatically send the results to clinicians could allow a repeated assessment (self-administered or minimally supervised by a caregiver) of dynamic balance and gait adaptability in ecological settings during walking tasks that are easily sustained by PwND. Use of such an app would permit increased assessment frequency, allow tracking of performance evolution caused by the disease or by the effects of rehabilitation/pharmacological treatment, and complement in-clinic examinations. One of the study's limitations is the retrospective and monocentric nature of the available data. Thus, given the retrospective nature of the data, the predictive power of IMU-related features on individual subitems of the mDGI scale could not be assessed. The latter would provide insights on the relationship between specific gait determinants/impairments and the different domains of the dynamic equilibrium concept. Moreover, in the proposed model, the strongest contributor to the prediction is gait speed, here estimated with the 6MWT score, which is commonly recorded by the physiotherapist in clinical practice but which is not available in real-life contexts. Future developments of the present approach should include the estimation of gait velocity from inertial sensors on the trunk59 or on the lower limbs60, 61. Another limitation is that the instrumented measures computed from the IMU were related to the central 10 strides of straight-line corridors only, thus excluding indexes descriptive of turns. This choice was made based on previous studies26, 50, 62 and on published results on healthy young adults63, which demonstrated that some of the considered gait quality measures (e.g., harmonic ratio) are affected by directional changes; to the best of our knowledge, however, the effect of turns on other indexes, such as temporal, smoothness, and intensity parameters has not been previously studied so far. Moreover, no data are available regarding the effect of turning while walking on gait parameters in people with neurological disorders. Nevertheless, considering the high impact of turning on balance and falls64, 65, the inclusion of measures descriptive of turns (e.g., turning speed and duration) could further improve the mDGI prediction. Future studies should address this issue and analyse the reliability of the present IMU-based assessment using a different group of 10 strides from the patient's 6MWT. Furthermore, even if the classifiers work well on average, the predictions include outliers. Nevertheless, from a methodological point of view, the SHAP analysis confirms that these outliers did not result from any consistent bias in the prediction and thus were not caused by a systematic error observable in our feature set. The latter, together with the nested cross-validation approach that simulates the testing of new patients entering the model, confirms the reliability of the results. From a clinical point of view, the presence of outliers with an mDGI estimation error greater than 5 points could lead to misinterpretation of the results. For example, if the difference between the mDGI scores collected at two distinct times is greater than 5 points, an actual significant change in performance or a prediction error could be the cause. To mitigate this issue, it is essential that the portable device embedding the prediction model automatically and immediately send the results to the clinician so that the clinician can analyse them and, if necessary, organise an in-clinic examination. However, further analysis concerning the use of a single IMU sensor and its optimal placement is needed before the algorithms are embedded within a smartphone app. In conclusion, we deployed a model using ML techniques to predict the mDGI score in a cohort of patients with neurological impairments and obtained similar accuracies across pathology groups. The nested cross-validation approach ensured that, within the outer leave-one-subject-out test loop, every patient was assigned to the test set once. Thus, the inner k-fold cross-validation loop, used for hyperparameter optimization, avoids train-dev-test contamination for the choice of model parameters. Furthermore, our approach targeted an estimation of the continuous mDGI value instead of classifying groups according to their risk of fall66, 67. This approach allows researcher and clinicians to post-process regression predictions and derive a classification by imposing error boundaries; the inverse is not possible. This valuable tool helps bridge the gap between traditional balance and gait-supervised assessments in clinical settings and automated, self-administered evaluations in real-life contexts, with the aim of reducing the time and costs necessary to track disease evolution or treatment effects.

Ninety-five people suffering from neurological diseases were recruited from the IRCCS Fondazione Don Carlo Gnocchi (Milan, Italy). The cohort was composed of 51 people with multiple sclerosis (MS), 25 people with Parkinson's disease (PD), and 19 people post-stroke (ST). To be included in the study, participants had to be between 20 and 85 years old; be able to walk for 20 m, even with an assistive device; and have a mini-mental state examination (MMSE) score of \(\ge\) 21. People with MS were included only if they had a certain diagnosis and had not relapsed in the previous two months. People with PD were enrolled only if their Hoehn and Yahr score was < 4. Post-stroke participants were included only if the time post-onset was > 2 months. Exclusion criteria included an inability to understand and sign the informed consent, the presence of a psychiatric complication, or major cardiovascular or visual disorders. All participants signed a written informed consent to participate in this study (conformed to the Declaration of Helsinki). All methods were approved by the Ethical Committee of IRCCS Fondazione Don Carlo Gnocchi, Milan (ref 29-03-2017 and 13-02-2019). All procedures were performed in accordance with the relevant guidelines and regulations.

The mDGI was administered by experienced physical therapists (Fig. 4, panel A)14. It consists of 8 items (e.g., walking with head turns, walking around or over obstacles, stairway walking). Each item is evaluated on the basis of three aspects: gait pattern (subscore: 0–3), level of assistance (sub-score: 0–2), and time (subscore: 0–3). The mDGI total score (i.e. the sum of the sub-scores of all item), ranges from 0 to 64, with increasing values indicating better performances. In PwND, the amount of previous activity could have an impact on the following performances; thus, the participants were required to rest after mDGI execution. The resting time was determined by the patient and supervised by the physiotherapist. The participants performed the 6-minute walk test (6MWT68), that measures walking endurance. The test required walking back and forth along a 30-m corridor for 6 min at a fast but safe speed. If needed, the participant could use an assistive device. The presence of an assistive device was here coded as monolateral, bilateral, or no support. Participants with PD were tested while they were on-phase during antiparkinsonian therapy, approximately two hours after medication intake. The distance covered over 6 min was recorded by the examiner and represented the clinical score of the test. The participants executed the 6MWT wearing three IMUs (MTw, XSens, NL) secured to the lower trunk (L5 level) and shanks, about 20 mm above the lateral malleoli. The trunk sensor was placed on the lower back since this position is the most commonly used, according to the literature69. The position of the shank sensors was chosen as it was the one associated with less sensor instability, due to soft tissue artifacts. In particular, report of a decrease of 4–51% artifacts were found compared to placing IMUs in other parts of the shank or of the feet70. The sensors were fixed to the body by elastic bands with Velcro strips, which can also be easily applied by the subject autonomously or with the help of a caregiver. Three-dimensional accelerations and angular velocities were recorded from the three IMUs at a sampling frequency of 75 Hz; this frequency was considered adequate for the purposes of the present study since it was within the range of sampling rates (25–1000 Hz) used in previous studies69, 71. Then trunk accelerations were reoriented to a horizontal–vertical coordinate system72. Only the short steady-state walking bouts, represented by 10 consecutive strides in the middle of each corridor, were considered for the subsequent analysis, after discarding the portions of signals pertaining to the 180\(^\circ\) turns at the end of each 30-m hallway73. Foot-strike and foot-off events were computed from the angular velocity around the medio-lateral axis of each shank74. Next, temporal gait determinants such as mean stride, step, and swing time, and single and double support duration were calculated. These metrics were chosen because they represent traditional measures of gait and show well-documented impairments in PwND33,34,35. Furthermore, for the step and stride times, the respective coefficients of variation were added to the dataset, as measures of step and stride variability, which are usually higher in PwND compared to healthy subjects75. Then, a set of 18 metrics was computed, from all trunk acceleration components (antero-posterior, medio-lateral and vertical), to provide information about gait quality domains (i.e., intensity, regularity, symmetry, stability, and smoothness) proposed by previous literature27,28,29,30,31,32. Gait intensity was quantified through the root mean squared value of the acceleration27. The gait regularity domain was represented by stride and step regularity indexes computed, respectively, as the second and the first peaks of the unbiased autocorrelation function calculated from each acceleration component76. Gait symmetry was quantified by the improved Harmonic Ratio (iHR) computed following Pasciuto et al.62. The stability domain was represented by the short-term Lyapunov exponent. This metric quantifies the local dynamic (in)stability of gait, which reflects the capability of the locomotor system to cope with small perturbations naturally present during walking, such as external disturbances or internal control errors54. The short-term Lyapunov exponent was calculated over the duration of one step, as fully detailed elsewhere28. Briefly, trunk accelerations related to ten consecutive strides in the central part of each walking bout were re-sampled to 1000 frames (10 strides \(\times\) 100 frames)54, 77 to maintain equal data length across walking bouts and participants. This procedure was applied for the computation of the short-term Lyapunov exponent only, since this parameter is highly influenced by signal length54, 77. Hence, the short-term Lyapunov exponent was computed following the Rosenstein method78, with m = 5 and T = 10 samples (m and T were estimated using published algorithms79). Increasing values of the Lyapunov exponent reflect the decreasing ability of the locomotor system to manage small perturbations, thus indicating greater dynamic instability. Finally, the gait smoothness domain was quantified via the logarithm of jerk (first time-derivative of the acceleration), normalized with respect to stride duration and mean acceleration80. All parameters were computed for each short steady-state walking bout (10 strides) derived from the 6MWT; the median values over the whole test were then calculated to reduce the effect of possible outliers. In the present cohort, the number of walking bouts was always greater than or equal to 3. The above gait quality metrics were chosen because they showed a statistically significant correlation with the mDGI score (see Table 2) and because previous literature has demonstrated their robustness to different test settings27, their sensitivity to subtle impairment28, 36, 38 and to rehabilitation effects37, and their ability to discriminate between different levels of disease severity50, 81.

Model pipeline. In panel (A), steps of the data collections protocol are reported. In panel (B and C), respectively the pre-processing steps and model deployment loops are presented.

No missing values were present in the dataset, and no data imputation techniques were adopted. To initially evaluate univariate correlations between IMU-based features and the mDGI, all extracted features were subject to Spearman's correlation analysis, with the dependent variable set to the mDGI. To evaluate whether categorical variables had an effect on the mDGI, binary (e.g., sex) and multi-class (e.g., presence of an assistive device) categorical variables were subject to a Mann–Whitney test and a Kruskal–Wallis test, respectively. This procedure was performed for the pooled cohort and the three single-pathology cohorts. Lastly, to assess whether different pathology groups were associated with different biomechanical determinants, a group analysis (Kruskal–Wallis test) was performed, with independent variables set to the IMU-derived features and grouping variables set to the pathology (MS, PD, ST). The significance threshold was set to 0.05.

To avoid train-test contamination, all extracted features entered subsequent ML models. Namely, a regularised type of regression, the elastic net (EN), was implemented. The EN combines the penalties of the LASSO and Ridge regressions82, overcoming their respective implementation problems. Ridge adds quadratic regularisation via L2 penalties, assigning a non-zero coefficient to all features in the model and keeping the coefficients even if the corresponding independent variable is irrelevant to the prediction. Conversely, LASSO regression is known to deteriorate with multicollinear independent variables83 but neglects specific features. Elastic net combines feature elimination from LASSO and coefficient reduction from Ridge and improves on both, yielding regression parameter estimates as follows:

The special cases \(\lambda _{2} = 0 \; \lambda _{1} \ne 0\) and \(\lambda _{1} = 0 \; \lambda _{2} \ne 0\) correspond to the LASSO and Ridge regressions, respectively, therefore including both LASSO and Ridge in the EN model hypothesis space. In the Scikit-Learn implementation

the \(l1_{ratio}\) describes the tendency toward a LASSO regularization (\(l1_{ratio} \sim 1\)) or a Ridge regularization (\(l1_{ratio} \sim 0\)) and the \(\alpha\) acts a scaling parameter of the regularization process.

To maximize robustness of the model, we implemented a nested cross-validation approach. In brief, such an approach consists of two k-fold cross-validation loops: an outer loop identifies the test set for each of its folds, while the inner loop implements a further split of the dataset for training and validation84. In the outer layer (testing layer), a LOSO testing procedure was used (Fig. 4, panel C). Specifically, one patient at a time was withheld for testing and the remaining N − 1 were used for training and cross-validation. The N − 1 patient were resampled using the Synthetic Minority Oversampling Technique for Regression (SMOTER85), generating training and validation samples equally distributed across the mDGI range. Then, the resampled training set was used to cross-validate and optimize the model's hyper-parameters by minimizing the cross-validation median absolute error \(E_{val} = \frac{\sum _{k=1}^{3} Accuracy_{k}}{K_{inner\, folds}}\) and averaging the accuracies across the K inner folds (Eq. 3). For the EN regression, \(\alpha\) and \(l1_{ratio}\) were optimized together with the majority class down sampling rate (\(\%d\)) and the minority class synthetic samples generation rate (\(\%o\))85. Specifically, both percentages were allowed to vary between the 200 and 800%. The number of inner folds (K) used to optimize all of the aforementioned parameters was set to 3. Then, with optimal hyper-parameters, EN regressions were retrained on all the \(N-1\) training patients, including the synthetic samples. Lastly, the model was tested on the patient withheld from the outer split. This procedure was repeated for the multi-pathological cohort (MP) and the single pathology cohorts (\(SP_{MS}\), \(SP_{PD}\) and \(SP_{ST}\)) and the predictions on the outer test samples were stored and aggregated (Fig. 4). The described pipeline was repeated using three different datasets: (1) IMU-derived features and clinical variables (6MWT score, presence of assistive devices), (2) clinical variables only (Supplementary Fig. 1), and (3) IMU-derived data only (Supplementary Fig. 2). All machine learning pipelines were implemented with the Optuna and Scikit-Learn libraries.

Generalised linear models already enable interpretability and explainability measures by assigning to each independent variable k the magnitude of its related regression coefficient k and therefore calculating the effect of all features on the prediction via the dot product \(\cdot\)x. Nevertheless, given the nested cross-validation implementation, in the outer split, each patient is included N 1 times in a training set and only once in a test set. Consequently, N models trained on the N permutation of N 1 patients will result in N parameter estimates (K, N). Accordingly, evaluating feature importance by averaging the N-coefficients is possible but has two major drawbacks. First, the resulting variability in the parameter estimates can be relevant. Second, estimates are derived from mean trends in the training subset and are not patient-specific. SHAP overcomes these limitations by determining feature contributions to the prediction specifically for the individual subjects, resulting in one value per subject per feature20, 21.

Wilcoxon signed-rank tests were applied between the MP model and the aggregated predictions of the SP models. Furthermore, individual SP models were also compared to the predictions of the same instances made by the MP model with Wilcoxon signed-rank tests.

Data has been provided with the manuscript and codes can be made available upon request to Piergiuseppe Liuzzi and Ilaria Carpinella for replication purposes only.

Modified dynamic gait index

Persons with neurological diseases

Timed up and go

Multiple sclerosis

Parkinson's disease

STroke

6-Minutes walk test

Inertial measurement units

Multi-pathology

Antero-posterior

VerTical

Medio-lateral

SHapley additive explanations

Minimally detectable change

Patla, A. Strategies for dynamic stability during adaptive human locomotion. IEEE Eng. Med. Biol. 22, 48–52. https://doi.org/10.1109/memb.2003.1195695 (2003).

Article Google Scholar

Zampogna, A. et al. Fifteen years of wireless sensors for balance assessment in neurological disorders. Sensors 20, 3247. https://doi.org/10.3390/s20113247 (2020).

Article ADS PubMed PubMed Central Google Scholar

Lencioni, T. et al. Strategies for maintaining dynamic balance in persons with neurological disorders during overground walking. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 235, 1079–1087. https://doi.org/10.1177/09544119211023624 (2021).

Article Google Scholar

Tinetti, M. E. & Kumar, C. The patient who falls: "It's always a trade-off’’. J. Am. Med. Assoc. 2303, 58–266. https://doi.org/10.1001/jama.2009.2024 (2010).

Article Google Scholar

Mancini, M. & Horak, F. B. The relevance of clinical balance assessment tools to differentiate balance deficits. Eur. J. Phys. Rehabil. Med. 46, 239–48 (2010).

CAS PubMed Google Scholar

Verghese, J., Holtzer, R., Lipton, R. B. & Wang, C. Quantitative gait markers and incident fall risk in older adults. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 64A, 896–901. https://doi.org/10.1093/gerona/glp033 (2009).

Article Google Scholar

Edginton Bigelow, K. & Berme, N. Development of a protocol for improving the clinical utility of posturography as a fall-risk screening tool. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 66A, 228–233. https://doi.org/10.1093/gerona/glq202 (2011).

Article Google Scholar

Stevens, J. A. & Phelan, E. A. Development of STEADI: A fall prevention resource for health care providers. Health Promot. Pract. 14, 706–714. https://doi.org/10.1177/1524839912463576 (2013).

Article PubMed Google Scholar

Muir-Hunter, S. & Wittwer, J. Dual-task testing to predict falls in community-dwelling older adults: A systematic review. Physiotherapy 102, 29–40. https://doi.org/10.1016/j.physio.2015.04.011 (2016).

Article CAS PubMed Google Scholar

Godi, M. et al. Comparison of reliability, validity, and responsiveness of the mini-bestest and berg balance scale in patients with balance disorders. Phys. Ther. 93, 158–167. https://doi.org/10.2522/ptj.20120171 (2012).

Article PubMed Google Scholar

Yelnik, A. & Bonan, I. Clinical tool for assessing balance disorders. Neurophysiol. Clin. Clin. Neurophysiol. 38, 439–45. https://doi.org/10.1016/j.neucli.2008.09.008 (2009).

Article Google Scholar

Matsuda, P. N., Taylor, C. S. & Shumway-Cook, A. Evidence for the validity of the modified dynamic gait index across diagnostic groups. Phys. Ther. 9, 996–1004. https://doi.org/10.2522/ptj.20130294 (2014).

Article Google Scholar

Corrini, C. et al. Minimal clinically important difference of modified dynamic gait index in people with neurological disorders. Gait Posture 90, 210–214. https://doi.org/10.1016/j.gaitpost.2021.08.024 (2021).

Article PubMed Google Scholar

Shumway-Cook, A., Taylor, C. S., Matsuda, P. N., Studer, M. T. & Whetten, B. K. Expanding the scoring system for the dynamic gait index. Phys. Ther. 93, 1493–1506. https://doi.org/10.2522/ptj.20130035 (2013).

Article PubMed Google Scholar

Shumway-Cook, A., Matsuda, P. N. & Taylor, C. Investigating the validity of the environmental framework underlying the original and modified dynamic gait index. Phys. Ther. 95, 864–870. https://doi.org/10.2522/ptj.20140047 (2015).

Article PubMed Google Scholar

Torchio, A. et al. Identification of modified dynamic gait index cutoff scores for assessing fall risk in people with Parkinson disease, stroke and multiple sclerosis. Gait Posture 91, 1–6. https://doi.org/10.1016/j.gaitpost.2021.09.201 (2021).

Article PubMed Google Scholar

Anastasi, D. et al. Instrumented version of the modified dynamic gait index in patients with neurologic disorders. PM R 11, 1312–1319. https://doi.org/10.1002/pmrj.12137 (2019).

Article PubMed Google Scholar

Ancona, S. et al. Wearables in the home-based assessment of abnormal movements in Parkinson's disease: A systematic review of the literature. J. Neurol.https://doi.org/10.1007/s00415-020-10350-3 (2021).

Article PubMed Google Scholar

Carpinella, I. et al. Instrumental assessment of stair ascent in people with multiple sclerosis, stroke, and Parkinson's disease: A wearable-sensor-based approach. IEEE Trans. Neural Syst. Rehabil. Eng. 26, 2324–2332. https://doi.org/10.1109/TNSRE.2018.2881324 (2018).

Article PubMed Google Scholar

Lundberg, S. M. et al. Explainable machine-learning predictions for the prevention of hypoxaemia during surgery. Nat. Biomed. Eng. 2, 749. https://doi.org/10.1038/s41551-018-0304-0 (2018).

Article PubMed PubMed Central Google Scholar

Lundberg, S. M. et al. From local explanations to global understanding with explainable AI for trees. Nat. Mach. Intell. 2, 2522–5839. https://doi.org/10.1038/s42256-019-0138-9 (2020).

Article Google Scholar

Demura, S., Yamaji, S. & Kitabayashi, T. Gender and age-related differences of dynamic balancing ability based on various stepping motions in the healthy elderly. J. Hum. Ergol. 34, 1–11 (2005).

Google Scholar

Marschollek, M. et al. Sensor-based fall risk assessment—An expert ‘to go’. Methods Inf. Med. 50, 420–6. https://doi.org/10.3414/ME10-01-0040 (2011).

Article CAS PubMed Google Scholar

van Schooten, K. et al. Ambulatory fall-risk assessment: Amount and quality of daily-life gait predict falls in older adults. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 70, 608–615. https://doi.org/10.1093/gerona/glu225 (2015).

Article CAS Google Scholar

Ihlen, E., Weiss, A., Beck, Y., Helbostad, J. & Hausdorff, J. A comparison study of local dynamic stability measures of daily life walking in older adult community-dwelling fallers and non-fallers. J. Biomech. 49, 1498–1503. https://doi.org/10.1016/j.jbiomech.2016.03.019 (2016).

Article PubMed Google Scholar

Riva, F., Bisi, M. & Stagni, R. Gait variability and stability measures: Minimum number of strides and within-session reliability. Comput. Biol. Med. 50C, 9–13. https://doi.org/10.1016/j.compbiomed.2014.04.001 (2014).

Article Google Scholar

Angelini, L. et al. Is a wearable sensor-based characterisation of gait robust enough to overcome differences between measurement protocols? a multi-centric pragmatic study in patients with multiple sclerosis. Sensors 20, 79. https://doi.org/10.3390/s20010079 (2019).

Article ADS PubMed PubMed Central Google Scholar

Caronni, A. et al. Local dynamic stability of gait in people with early multiple sclerosis and no-to-mild neurological impairment. IEEE Trans. Neural Syst. Rehabil. Eng. 28, 1389–1396. https://doi.org/10.1109/TNSRE.2020.2991636 (2020).

Article PubMed Google Scholar

Plotnik, M., Wagner, J., Adusumilli, G., Gottlieb, A. & Naismith, R. Gait asymmetry, and bilateral coordination of gait during a six-minute walk test in persons with multiple sclerosis. Sci. Rep. 10, 12382. https://doi.org/10.1038/s41598-020-68263-0 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Socie, M. & Sosnoff, J. Gait variability and multiple sclerosis. Multiple Scler. Int. 2013, 645197. https://doi.org/10.1155/2013/645197 (2013).

Article Google Scholar

Patterson, K., Gage, W., Brooks, D., Black, S. & McILROY, W. Evaluation of gait symmetry after stroke: A comparison of current methods and recommendations for standardization. Gait Posture 31, 241–6. https://doi.org/10.1016/j.gaitpost.2009.10.014 (2009).

Article PubMed Google Scholar

Godi, M., Arcolin, I., Giardini, M., Corna, S. & Schieppati, M. A pathophysiological model of gait captures the details of the impairment of pace/rhythm, variability and asymmetry in parkinsonian patients at distinct stages of the disease. Sci. Rep. 11, 21143. https://doi.org/10.1038/s41598-021-00543-9 (2021).

Article ADS PubMed PubMed Central Google Scholar

Chee, J. et al. Influence of multiple sclerosis on spatiotemporal gait parameters: A systematic review and meta-regression. Arch. Phys. Med. Rehabil. 102, 1801–1815. https://doi.org/10.1016/j.apmr.2020.12.013 (2021).

Article PubMed Google Scholar

Morris, M., Iansek, R., Matyas, T. & Summers, J. The pathogenesis of gait hypokinesia in Parkinson's disease. Brain 117, 1169–1181. https://doi.org/10.1093/brain/117.5.1169 (1994).

Article PubMed Google Scholar

Wonsetler, E. & Bowden, M. A systematic review of mechanisms of gait speed change post-stroke. Part 1. Spatiotemporal parameters and asymmetry ratios. Top. Stroke Rehabil. 24, 1–12. https://doi.org/10.1080/10749357.2017.1285746 (2017).

Article Google Scholar

Carpinella, I. et al. Walking with horizontal head turns is impaired in persons with early-stage multiple sclerosis showing normal locomotion. Front. Neurol. 12, 821640. https://doi.org/10.3389/fneur.2021.821640 (2022).

Article PubMed PubMed Central Google Scholar

Hilfiker, R. et al. Local dynamic stability as a responsive index for the evaluation of rehabilitation effect on fall risk in patients with multiple sclerosis: A longitudinal study. BMC Res. Notes 6, 260. https://doi.org/10.1186/1756-0500-6-260 (2013).

Article PubMed PubMed Central Google Scholar

Carpinella, I. et al. Instrumentally assessed gait quality is more relevant than gait endurance and velocity to explain patient-reported walking ability in early stage multiple sclerosis. Eur. J. Neurol. 28, 2259–2268. https://doi.org/10.1111/ene.14866 (2021).

Article PubMed Google Scholar

Bennett, S., Bromley, L., Fisher, N., Tomita, M. & Niewczyk, P. Validity and reliability of four clinical gait measures in patients with multiple sclerosis. Int. J. MS Care 19, 247–252. https://doi.org/10.7224/1537-2073.2015-006 (2017).

Article PubMed PubMed Central Google Scholar

Combs, S. A., Diehl, D. M., Filip, J. & Long, E. Short-distance walking speed tests in people with Parkinson disease: Reliability, responsiveness, and validity. Gait Posture 39, 784–788. https://doi.org/10.1016/j.gaitpost.2013.10.019 (2014).

Article PubMed Google Scholar

Pohl, P. et al. Influence of stroke-related impaiments on performance in 6-minute walk test. J. Rehabil. Res. Dev. 39, 439–44 (2002).

PubMed Google Scholar

Fritz, S. & Lusardi, M. White paper: "Walking speed: the sixth vital sign’’. J. Geriatr. Phys. Ther. 2001(32), 46–9. https://doi.org/10.1519/00139143-200932020-00002 (2009).

Article Google Scholar

Callisaya, M. L. & Verghese, J. The association of clinic-based mobility tasks and measures of community performance and risk. PM R 10, 704-711.e1. https://doi.org/10.1016/j.pmrj.2017.12.008 (2018).

Article PubMed PubMed Central Google Scholar

Menz, H., Lord, S. & Fitzpatrick, R. Age-related differences in walking stability. Age Ageing 32, 137–42 (2003).

Article PubMed Google Scholar

Lencioni, T., Carpinella, I., Rabuffetti, M., Cattaneo, D. & Ferrarin, M. Measures of dynamic balance during level walking in healthy adult subjects: Relationship with age, anthropometry and spatio-temporal gait parameters. Porc. Inst. Mech. Eng. H 65, 131–140. https://doi.org/10.1177/0954411919889237 (2020).

Article Google Scholar

Moe-Nilssen, R. & Helbostad, J. Interstride trunk acceleration variability but not step width variability can differentiate between fit and frail older adults. Gait Posture 21, 164–70. https://doi.org/10.1016/j.gaitpost.2004.01.013 (2005).

Article PubMed Google Scholar

Callisaya, M., Blizzard, L., Mcginley, J., Schmidt, M. & Srikanth, V. Sensorimotor factors affecting gait variability in older people-a population-based study. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 65, 386–92. https://doi.org/10.1093/gerona/glp184 (2009).

Article Google Scholar

Handelzalts, S. et al. Insights into motor performance deficits after stroke: An automated and refined analysis of the lower-extremity motor coordination test (lemocot). J. NeuroEng. Rehabil. 18, 1–10. https://doi.org/10.1186/s12984-021-00950-z (2021).

Article Google Scholar

Kalron, A., Dolev, M. & Givon, U. Further construct validity of the timed up and go test as a measure of ambulation in persons with multiple sclerosis. Eur. J. Phys. Med. Rehabil. 53, 841–847. https://doi.org/10.23736/S1973-9087.17.04599-3 (2017).

Article Google Scholar

Angelini, L. et al. Wearable sensors can reliably quantify gait alterations associated with disability in people with progressive multiple sclerosis in a clinical setting. J. Neurol. 267, 2897–2909. https://doi.org/10.1007/s00415-020-09928-8 (2020).

Article PubMed PubMed Central Google Scholar

Mizuike, C., Ohgi, S. & Morita, S. Analysis of stroke patient walking dynamics using a tri-axial accelerometer. Gait Posture 30, 60–4. https://doi.org/10.1016/j.gaitpost.2009.02.017 (2009).

Article PubMed Google Scholar

Buckley, C., Galna, B., Rochester, L. & Mazzà, C. Upper body accelerations as a biomarker of gait impairment in the early stages of Parkinson's disease. Gait Posture 71, 289–295. https://doi.org/10.1016/j.gaitpost.2018.06.166 (2018).

Article PubMed Google Scholar

Huisinga, J., Mancini, M., St George, R. & Horak, F. Accelerometry reveals differences in gait variability between patients with multiple sclerosis and healthy controls. Ann. Biomed. Eng. 41, 1670–1679. https://doi.org/10.1007/s10439-012-0697-y (2012).

Article PubMed PubMed Central Google Scholar

Bruijn, S., Meijer, O., Beek, P. & Van Dieen, J. Assessing the stability of human locomotion: A review of current measures. J. R. Soc. Interface R. Soc. 10, 20120999. https://doi.org/10.1098/rsif.2012.0999 (2013).

Article CAS Google Scholar

Ellis, R., Howard, K. & Kram, R. The metabolic and mechanical costs of step time asymmetry in walking. Proc. Biol. Sci. R. Soc. 280, 20122784. https://doi.org/10.1098/rspb.2012.2784 (2013).

Article Google Scholar

Alexander, S., Peryer, G., Gray, E., Barkof, F. & Chataway, J. Wearable technologies to measure clinical outcomes in multiple sclerosis: A scoping review. Multiple Scler. 27, 1643–1656. https://doi.org/10.1177/1352458520946005 (2020).

Article Google Scholar

Casamassima, F. et al. A wearable system for gait training in subjects with Parkinson's disease. Sensors (Basel, Switzerland) 14, 6229–46. https://doi.org/10.3390/s140406229 (2014).

Article ADS PubMed Google Scholar

Ginis, P. et al. Feasibility and effects of home-based smartphone-delivered automated feedback training for gait in people with Parkinson's disease: A pilot randomized controlled trial. Parkinsonism Relat. Disord. 22, 28–34. https://doi.org/10.1016/j.parkreldis.2015.11.004 (2015).

Article Google Scholar

Soltani, A. et al. Algorithms for walking speed estimation using a lower-back-worn inertial sensor: A cross-validation on speed ranges. IEEE Trans. Neural Syst. Rehabil. Eng. 29, 1955–1964. https://doi.org/10.1109/TNSRE.2021.3111681 (2021).

Article CAS PubMed Google Scholar

Ferrari, A. et al. A mobile Kalman-filter based solution for the real-time estimation of spatio-temporal gait parameters. IEEE Trans. Neural Syst. Rehabil. Eng. Publ. IEEE Eng. Med. Biol. Soc. 24, 764–773. https://doi.org/10.1109/TNSRE.2015.2457511 (2015).

Article Google Scholar

Trojaniello, D. et al. Estimation of step-by-step spatio-temporal parameters of normal and impaired gait using shank-mounted magneto-inertial sensors: Application to elderly, hemiparetic, parkinsonian and choreic gait. J. Neuroeng. Rehabil. 11, 152. https://doi.org/10.1186/1743-0003-11-152 (2014).

Article PubMed PubMed Central Google Scholar

Pasciuto, I., Bergamini, E., Iosa, M., Vannozzi, G. & Cappozzo, A. Overcoming the limitations of the harmonic ratio for the reliable assessment of gait symmetry. J. Biomech. 53, 84–89. https://doi.org/10.1016/j.jbiomech.2017.01.005 (2017).

Article PubMed Google Scholar

Riva, F., Grimpampi, E., Mazzà, C. & Stagni, R. Are gait variability and stability measures influenced by directional changes?. Biomed. Eng. Online 13, 56. https://doi.org/10.1186/1475-925X-13-56 (2014).

Article PubMed PubMed Central Google Scholar

Caronni, A. et al. Criterion validity of the instrumented timed up and go test: A partial least square regression study. Gait Posture 61, 287–293. https://doi.org/10.1016/j.gaitpost.2018.01.015 (2018).

Article PubMed Google Scholar

Leach, J., Mellone, S., Palumbo, P., Bandinelli, S. & Chiari, L. Natural turn measures predict recurrent falls in community-dwelling older adults: A longitudinal cohort study. Sci. Rep. 8, 1–11. https://doi.org/10.1038/s41598-018-22492-6 (2018).

Article CAS Google Scholar

Lindberg, D. et al. Identification of important factors in an inpatient fall risk prediction model to improve the quality of care using EHR and electronic administrative data: A machine-learning approach. Int. J. Med. Inform. 143, 104272. https://doi.org/10.1016/j.ijmedinf.2020.104272 (2020).

Article PubMed PubMed Central Google Scholar

Noh, B. et al. Xgboost based machine learning approach to predict the risk of fall in older adults using gait outcomes. Sci. Rep. 11, 12183. https://doi.org/10.1038/s41598-021-91797-w (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Brooks, D., Solway, S. & Gibbons, W. ATS statement on six-minute walk test. Am. J. Respir. Crit. Care Med. 167, 1287. https://doi.org/10.1164/ajrccm.167.9.950 (2003).

Article PubMed Google Scholar

Storm, F., Cesareo, A., Reni, G. & Biffi, E. Wearable inertial sensors to assess gait during the 6-minute walk test: A systematic review. Sensors 20, 2660. https://doi.org/10.3390/s20092660 (2020).

Article ADS PubMed PubMed Central Google Scholar

Scalera, G., Ferrarin, M., Marzegan, A. & Rabuffetti, M. Assessment of stability of MIMU probes to skin-marker-based anatomical reference frames during locomotion tasks: Effect of different locations on the lower limb. Front. Bioeng. Biotechnol. 9, 721900. https://doi.org/10.3389/fbioe.2021.721900 (2021).

Article PubMed PubMed Central Google Scholar

Pacini Panebianco, G., Bisi, M., Stagni, R. & Fantozzi, S. Analysis of the performance of 17 algorithms from a systematic review: Influence of sensor position, analysed variable and computational approach in gait timing estimation from imu measurements. Gait Posture 66, 76–82. https://doi.org/10.1016/j.gaitpost.2018.08.025 (2018).

Article PubMed Google Scholar

Moe-Nilssen, R. A new method for evaluating motor control in gait under real-life environmental conditions. Clin. Biomech. 13, 320–327. https://doi.org/10.1016/s0268-0033(98)00089-8 (1998).

Article Google Scholar

Nguyen, H. et al. Auto detection and segmentation of daily living activities during a timed up and go task in people with Parkinson's disease using multiple inertial sensors. J. NeuroEng. Rehabil. 14, 1–13. https://doi.org/10.1186/s12984-017-0241-2 (2017).

Article MathSciNet Google Scholar

Salarian, A. Gait assessment in Parkinson's disease: Toward an ambulatory system for long-term monitoring. IEEE Trans. Biomed. Eng. 51, 1434–43. https://doi.org/10.1109/TBME.2004.827933 (2004).

Article PubMed Google Scholar

Moon, Y., Sung, J., An, R., Hernandez, M. & Sosnoff, J. Gait variability in people with neurological disorders: A systematic review and meta-analysis. Hum. Mov. Sci. 47, 197–208. https://doi.org/10.1016/j.humov.2016.03.010 (2016).

Article PubMed Google Scholar

Moe-Nilssen, R. & Helbostad, J. L. Estimation of gait cycle characteristics by trunk accelerometry. J. Biomech. 37, 121–126. https://doi.org/10.1016/S0021-9290(03)00233-1 (2004).

Article PubMed Google Scholar

England, S. & Granata, K. The influence of gait speed on local dynamic stability of walking. Gait Posture 25, 172–8. https://doi.org/10.1016/j.gaitpost.2006.03.003 (2007).

Article PubMed Google Scholar

Rosenstein, M. T., Collins, J. J. & De Luca, C. J. A practical method for calculating largest Lyapunov exponents from small data sets. Physica D Nonlinear Phenom. 65, 117–134. https://doi.org/10.1016/0167-2789(93)90009-P (1993).

Article ADS MathSciNet MATH Google Scholar

Wallot, S. & Mønster, D. Calculation of average mutual information (AMI) and false-nearest neighbors (FNN) for the estimation of embedding parameters of multidimensional time series in Matlab. Front. Psychol. 9, 1679. https://doi.org/10.3389/fpsyg.2018.01679 (2018).

Article PubMed PubMed Central Google Scholar

Caby, B., Kieffer, S., Hubert, M., Cremer, G. & Macq, B. Feature extraction and selection for objective gait analysis and fall risk assessment by accelerometry. Biomed. Eng. Online 10, 1. https://doi.org/10.1186/1475-925X-10-1 (2011).

Article PubMed PubMed Central Google Scholar

Castiglia, S. et al. Ability of a set of trunk inertial indexes of gait to identify gait instability and recurrent fallers in Parkinson's disease. Sensors 21, 3449. https://doi.org/10.3390/s21103449 (2021).

Article ADS PubMed PubMed Central Google Scholar

Zou, H. & Hastie, T. Regularization and variable selection via the elastic nets. J. R. Stat. Soc. B 67, 301–320 (2015).

Article MathSciNet MATH Google Scholar

Fan, J. & Li, R. Statistical challenges with high dimensionality: Feature selection in knowledge discovery. In Proceedings of the Madrid International Congress of Mathematicians, vol. 3 (2006).

Parvandeh, S., Yeh, H.-W., Paulus, M. P. & McKinney, B. A. Consensus features nested cross-validation. Bioinformatics 36, 3093–3098. https://doi.org/10.1093/bioinformatics/btaa046 (2020).

Article CAS PubMed PubMed Central Google Scholar

Torgo, L., Branco, P., Ribeiro, R. P. & Pfahringer, B. Resampling strategies for regression. Expert Syst. 32, 465–476. https://doi.org/10.1111/exsy.12081 (2015).

Article Google Scholar

Download references

This work was supported in part by the "Ricerca corrente RC2020-RC2021 Program" and in part by the 5 × 1000 funds AF2018: "Data Science in Rehabilitation Medicine" AF2019: "Study and development of biomedical data science and machine learning methods to support the appropriateness and the decision-making process in rehabilitation medicine" by the Italian Ministry of Health. This study was also supported by the Tuscany Region through the Tuscany Network for BioElectronic Approaches in Medicine: AI based predictive algorithms for fine-tuning of eletroceutics treatments in neurological, cardiovascular and endocrinological diseases (TUNE-BEAM, n. H14I20000300002).

These authors contributed equally: Maurizio Ferrarin and Andrea Mannini.

AIRLab, IRCCS Fondazione Don Carlo Gnocchi ONLUS, 50143, Florence, Italy

Piergiuseppe Liuzzi & Andrea Mannini

Scuola Superiore Sant’Anna, Istituto di BioRobotica, 56025, Pontedera, Italy

Piergiuseppe Liuzzi & Maria Chiara Carrozza

LAMoBIR and LaRiCE, IRCCS Fondazione Don Carlo Gnocchi ONLUS, 20148, Milan, Italy

Ilaria Carpinella, Denise Anastasi, Elisa Gervasoni, Tiziana Lencioni, Rita Bertoni, Davide Cattaneo & Maurizio Ferrarin

Dipartimento di Fisiopatologia Medico-Chirurgica e dei Trapianti, Università di Milano, 20122, Milan, Italy

Davide Cattaneo

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

Study conception and design: D.C., M.F., A.M. Data acquisition: E.G., D.A., R.B. Patients’ clinical assessment: E.G., D.A., R.B. Signal elaboration: I.C., T.L. Data analysis: P.L., A.M. Data interpretation: all authors. Manuscript drafting: P.L., I.C. Manuscript critical revision: all authors.

Correspondence to Ilaria Carpinella.

The authors declare no competing interests.

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

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

Reprints and Permissions

Liuzzi, P., Carpinella, I., Anastasi, D. et al. Machine learning based estimation of dynamic balance and gait adaptability in persons with neurological diseases using inertial sensors. Sci Rep 13, 8640 (2023). https://doi.org/10.1038/s41598-023-35744-x

Download citation

Received: 13 February 2023

Accepted: 23 May 2023

Published: 27 May 2023

DOI: https://doi.org/10.1038/s41598-023-35744-x

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

SHARE