Next Article in Journal
Towards Machine Learning-Based Detection of Running-Induced Fatigue in Real-World Scenarios: Evaluation of IMU Sensor Configurations to Reduce Intrusiveness
Next Article in Special Issue
Smart Devices and Wearable Technologies to Detect and Monitor Mental Health Conditions and Stress: A Systematic Review
Previous Article in Journal
Dimensionally Stable Anode Based Sensor for Urea Determination via Linear Sweep Voltammetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson’s Disease

by
Stefano Filippo Castiglia
1,
Antonella Tatarelli
2,3,
Dante Trabassi
1,
Roberto De Icco
4,5,
Valentina Grillo
4,
Alberto Ranavolo
3,
Tiwana Varrecchia
3,
Fabrizio Magnifica
2,
Davide Di Lenola
1,
Gianluca Coppola
1,
Donatella Ferrari
6,
Alessandro Denaro
6,
Cristina Tassorelli
4,5 and
Mariano Serrao
1,7,*
1
Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza, 04100 Latina, Italy
2
Department of Human Neurosciences, University of Rome Sapienza, 00185 Rome, Italy
3
Department of Occupational and Environmental Medicine, Epidemiology and Hygiene, INAIL, Monte Porzio Catone, 00078 Rome, Italy
4
Movement Analysis Research Unit, IRCCS Mondino Foundation, 27100 Pavia, Italy
5
Department of Brain and Behavioral Sciences, University of Pavia, 27100 Pavia, Italy
6
U.O. Riabilitazione Neurologica, Istituto Chirurgico Ortopedico Traumatologico, Latina 04100, Italy
7
Movement Analysis Laboratory, Policlinico Italia, 00162 Rome, Italy
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(10), 3449; https://doi.org/10.3390/s21103449
Submission received: 22 April 2021 / Revised: 8 May 2021 / Accepted: 13 May 2021 / Published: 15 May 2021
(This article belongs to the Special Issue Feature Papers in Wearables Section 2021)

Abstract

:
The aims of this study were to assess the ability of 16 gait indices to identify gait instability and recurrent fallers in persons with Parkinson’s disease (pwPD), regardless of age and gait speed, and to investigate their correlation with clinical and kinematic variables. The trunk acceleration patterns were acquired during the gait of 55 pwPD and 55 age-and-speed matched healthy subjects using an inertial measurement unit. We calculated the harmonic ratios (HR), percent recurrence, and percent determinism (RQAdet), coefficient of variation, normalized jerk score, and the largest Lyapunov exponent for each participant. A value of ≤1.50 for the HR in the antero-posterior direction discriminated between pwPD at Hoehn and Yahr (HY) stage 3 and healthy subjects with a 67% probability, between pwPD at HY 3 and pwPD at lower HY stages with a 73% probability, and it characterized recurrent fallers with a 77% probability. Additionally, HR in the antero-posterior direction was correlated with pelvic obliquity and rotation. RQAdet in the antero-posterior direction discriminated between pwPD and healthy subjects with 67% probability, regardless of the HY stage, and was correlated with stride duration and cadence. Therefore, HR and RQAdet in the antero-posterior direction can both be used as age- and-speed-independent markers of gait instability.

1. Introduction

Gait abnormality is one of the most invalidating features of persons with Parkinson’s disease (pwPD) [1] and closely parallels disease progression [2], leading to a high risk of falls [2], impaired patient autonomy, and reduced quality of life [3,4]. Based on the significant physical and economic burdens associated with falls [5,6], it is critical to identify sensitive and specific indexes that can represent various aspects of gait imbalance [7]. This would allow us to improve the ability to assess gait stability in pwPD, providing a multifold positive impact in terms of patient care ranging from the prevention of falls to the possibility of monitoring and optimizing both pharmacological and rehabilitation treatments [8,9]. Although many clinical tools have been developed to identify unbalance in pwPD [10,11,12], technology-based objective measures may improve the ability to capture motor behaviors to optimize treatment strategies [8,13].
From a biomechanical perspective, gait stability refers to the ability to maintain functional locomotion in the presence of small kinematic disturbances or control errors [14,15,16,17]. A series of quantitative biomechanical measures have been proposed to assess dynamic gait imbalance in pwPD [18]. One of the most widely explored measures is the coefficient of variation (CV) of spatio-temporal parameters (i.e., stride time and step length) [19,20,21,22,23,24]. Other reported measures have been derived from the accelerative patterns of the trunk and include the harmonic ratio (HR) [25,26], largest Lyapunov exponent (LLE) [27] and normalized jerk score (NJS) [28,29]. Another trunk-acceleration-derived measure, the recurrence quantification analysis (RQA), has been reported to be impaired in patients with vestibulopathy [30] and provides a promising approach to the investigation of center of pressure fluctuations in older fallers [31].
In almost all previous studies, gait stability indexes have been investigated individually to identify preventive measures for falls. However, in general terms, it would be ideal to obtain a set of gait stability indexes to understand which are the most sensitive in terms of distinguishing between pathological and neurotypical gaits, which are the most closely related to the severity of the disease, which aspects of stability are expressed by various indexes, which are the most useful for predicting falls, and which are the most responsive to rehabilitation treatments [7].
Additionally, because the correlations between gait speed, age, and trunk accelerations have been well described and because walking stability is strongly influenced by gait speed in pwPD [32], which is a marker of gait alteration common to most neurological diseases [33], the value of the gait stability indexes as independent markers of gait alteration should be investigated, regardless of these confounding factors [33].
Based on their ease of use and cost effectiveness, wearable sensors are widely used for clinical assessment to obtain more objective measures of walking performance [34,35,36,37,38,39]. Based on the inverted pendulum model of gait [40], mobile inertial measurement units (IMU), which embed accelerometers and gyroscopes, are able to objectively capture the ability to control the body center of mass while moving the base of support, resulting in an effective tool to monitor dynamic balance during gait [36,41]. Ideally, by using several combined IMUs to analyze a subject’s gait, the overall accuracy of the analysis would improve, but this gain would be offset by the wearability burden. Conversely, a single lumbar-mounted IMU provides sensitive information on the gait of pwPD and allows clinicians to monitor the gait of pwPD also in free-living conditions [42]. IMUs directly provide trunk acceleration measures and facilitate the recording of patient gaits for many steps during follow-up clinical assessments in outpatient facilities, making them ideal tools for investigating gait stability in pwPD.
In this study, we measured a set of gait stability indexes in data samples from pwPD that were collected using a wearable device with the goals of (i) determining the accuracy of each index in terms of detecting gait instability in pwPD compared to healthy subjects (HS) by excluding the effects of age and gait speed between groups, (ii) assessing the ability of each identified index to characterize the gait of pwPD who fall recurrently, compared with pwPD who do not fall recurrently, and (iii) exploring the correlations of each index with clinical and biomechanical variables. Our main hypotheses can be summarized as follows: (i) Some gait stability indexes may have sufficient ability to discriminate between pwPD and HS and may be correlated with the spatio-temporal and kinematic gait features of pwPD, regardless of gait speed; (ii) the identified discriminant indexes may identify frequent fallers among pwPD as an expression of walking instability.

2. Materials and Methods

2.1. Subjects

We collected data samples from 62 subjects diagnosed with idiopathic PD (20 females and 42 males, aged 73 ± 6.60 years) who enrolled for our study at “ICOT” in Latina, Italy, and the “IRCCS Mondino Foundation” in Pavia, Italy. The inclusion criteria were defined as follows: (i) diagnosis of idiopathic PD according to the UK bank criteria [43], (ii) Hoehn and Yahr (HY) stages 1 to 3 [44], (iii) ability to walk independently for at least 8 m along a laboratory pathway without exhibiting gait freezing and the ability to perform repeated walking trials with at least 10 consecutive strides [45,46,47], and (iv) on a stable drug program and acclimated to their current medication use for at least two weeks. The exclusion criteria were defined as follows: (i) cognitive deficits (defined as scores less than 26 on the Mini-Mental State Examination) [48,49], (ii) moderate or severe depression (defined as scores greater than 17 on the Beck Depression Inventory) [50,51], and (iii) orthopedic and/or other gait-influencing diseases, such as other neurological diseases, clinically defined osteoarthritis [52,53,54], or total hip joint replacement. Particularly, subjects referring pain in hip or knee joints and reduced range of motion on internal hip rotation, hip flexion and visible anatomic alterations of the joints were excluded [54].
For group comparisons, a group of 55 HS (aged 61.40 ± 11.98 years) were recruited. Each HS repeated the tests twice: once walking at a self-selected speed and once walking at a slower directed speed.
Informed consent was obtained from all participants in compliance with the Helsinki Declaration and local ethics committee approval was obtained (CE Lazio2, protocol number: 0053667/2021).

2.2. Instrumentation

An inertial sensor (BTS GWALK, BTS, Milan, Italy) attached to an ergonomic belt placed around the pelvis at the level of the L5 vertebra and connected to a laptop via Bluetooth was used to acquire data. This sensor is equipped with a tri-axial accelerometer (16 bit/axes), tri-axial magnetometer (13 bits), and tri-axial gyroscope (16 bit/axes). Linear trunk accelerations and trunk angular velocities in the anterior-posterior (AP), mediolateral (ML), and vertical (V) directions were measured at a sampling rate of 100 Hz (Figure 1).

2.3. Task Description

To familiarize themselves with the experimental procedure, the participants were asked to walk on the ground along a predetermined pathway before the experimental session. Both the pwPD and HS were asked to walk barefoot at a self-selected speed along a corridor (approximately 3 m wide and 10 m long). To avoid the influence of gait speed on other gait parameters and collect the largest possible sample size for speed-matched comparisons, the HS were also asked to walk at a slower speed [1,55,56,57,58,59]. Because natural locomotion was the main focus of this study, only general and qualitative instructions were provided, leaving the subjects free to choose their speed level without interfering with their pacing and rhythm based on external sensory cues.

2.4. Inertial Sensor Data Processing

The “Walk+” protocol of the G-STUDIO software (G-STUDIO, BTS, Milan, Italy) was used to detect trunk acceleration patterns, right and left heel strikes, toe-off, spatiotemporal parameters, and pelvis kinematics. The HR, RQA, CV, NJS, and largest Lyapunov exponent (LLE) were calculated (Figure 1) using MATLAB (MATLAB 7.4.0, MAthWorks, Natick, MA, USA).

2.4.1. HR Calculation

The HR [60] provides an indication of the acceleration patterns, smoothness, and rhythm of the gait. The greater the HR value, the smoother the walking pattern.
A stable and rhythmic gait pattern should be characterized by multiples of two acceleration patterns within any given stride. Irregular acceleration signals during a walking trial reflect unstable gait patterns and should be characterized by acceleration patterns that do not repeat in multiples of two, resulting in out-of-phase accelerations. The harmonic content of the trunk acceleration signals was analyzed in the AP (HRAP), ML (HRML), and V (HRV) directions by using stride frequency as the principal frequency component. Because most of the power in a normal cadence occurs below 10 Hz, 20 harmonics were calculated for each subject according to their stride times. A discrete Fourier transform was used to break down the trunk accelerations of each stride into individual sinusoidal waveforms. In a stable and smooth gait pattern, the acceleration signals in the AP and V directions repeat in multiples of two during a single stride. Therefore, HRAP and HRV were calculated as the ratio of the sum of the amplitudes of the first 10 even harmonics divided by the sum of the amplitudes of the first 10 odd harmonics. HRML was calculated as the sum of the amplitudes of the odd harmonics divided by the sum of the amplitudes of the even harmonics. A high-pass filter with a 20 Hz cutoff was used to eliminate noise signals and the HR values were calculated as follows [61]:
H R A P ,   V   = i A i * 2 i A i * 2 1 ,
H R M L   = i A i * 2 1 i A i 2 ,
where A i represents the amplitudes of the first 20 even harmonics and A2i−1 represents the amplitudes of the first 20 odd harmonics.

2.4.2. RQA Calculation

RQA is a nonlinear technique that can provide useful information regarding the patterns and structures of system dynamics. It provides a characterization of a variety of features of a given time series, including the quantification of deterministic structures and non-stationarity, based on the construction of recurrence plots. A detailed description of RQA calculation was provided by Webber et al. (1994) [62]. Acceleration and angular velocity data are embedded in “m” dimensions using “m” copies of the original time series, where each copy is shifted in time by integer multiples of “τ” samples. The embedding dimension “m” is the first recurrence parameter and is estimated using the nearest-neighbor method [63], which compares the distances between neighboring trajectories at successively higher dimensions. ‘‘False neighbors’’ occur when trajectories that overlap in dimension m i are distinguished in dimension m i + 1 . As i increases, the total percentage of false neighbors decrease, and m is chosen where this percentage approaches 0. False neighbors analysis was performed using values of Rtol = 17 and Atol = 2 [14,30,46,64,65]. We considered a maximum embedding dimension of 10 (range: 2–10), and m = 5 was considered as the optimal embedding dimension. Here, “τ” represents the second recurrence parameter and is selected to minimize the interaction between the points in the measured time series. Two common methods for selecting a proper delay are finding the first minimum in either the (linear) autocorrelation function or (nonlinear) mutual information function of the continuous time series. τ was calculated from the first minimum of the average mutual information (AMI) function [66], which evaluates the shared amount of information in bits between two data sets over a range of time delays. By choosing the first minimum of the AMI function, adjacent delay coordinates with a minimum of redundancy are provided. A range of τ between 7 and 18 was evaluated in this study. The time delay computed by the first minimum of the AMI considered as optimal was 10 samples [46,62,64,65].
A distance matrix is then computed by calculating the Euclidean distances between all embedded vectors. A recurrence matrix is computed by selecting a threshold (radius) of 10% of the maximum distance, where all cells with values below this threshold are identified as recurrent points. RQA variables are used to quantify the structure of the recurrence matrix. Percent recurrence (RQArec) can be calculated to understand how often a trajectory visits similar locations in the state space (time independent) and is computed as the percentage of recurrent points in the recurrence matrix. Percent determinism (RQAdet) can be calculated to understand how often a trajectory repeatedly revisits similar state space locations (time dependent) and is quantified as the percentage of recurrent points in the diagonal line structures (at least four consecutive points in length) parallel to the main diagonal. Therefore, RQArec quantifies the number of potentially recurrent points, where only a portion of these points recur periodically and are related to the predictability of the target dynamical system. The higher the RQArec and RQAdet values, the higher the predictability of the system [30].

2.4.3. CV

To compute the step length CV, the step length was estimated from the upward and downward movements of the trunk, as proposed by Zijlstra and Hof [67]. Assuming a compass gait type, the body’s center of mass (CoM) movements in the sagittal plane follow a circular trajectory during each single-support phase. In this inverted pendulum model, changes in the height of the CoM depend on the step length [67]. Therefore, the step length is calculated as follows:
s t e p   l e n g t h = 2 2 l h h 2 ,
where h denotes the change in the height of the CoM and l represents the pendulum length. A double integration of the vertical acceleration was implemented to calculate changes in the vertical position. A fourth-order zero-lag Butterworth high-pass filter (0.1 Hz) was used to avoid integration drift. The amplitude of the changes in the vertical position (h) was calculated as the difference between the highest and lowest positions during a step cycle. The leg length was considered to be the pendulum length (l). The step length was calculated as the mean of the observed step lengths during seven consecutive steps for each subject. The stride length CV was computed as follows:
CV = 100SD/mean,
where mean is the mean step length and SD is the standard deviation over all step lengths for each subject [68]. The higher the CV, the higher the variability in step length.

2.4.4. NJS

The NJS measures the time-normalized rate of change in the acceleration signal during stepping [28,29]. The acceleration data were first bandpass filtered using a fourth-order zero-lag Butterworth filter with a cutoff frequency of 20 Hz. We then calculated the NJS from the time duration between each foot contact by applying the following formula [29]:
G a i t   N J S = 1 N i = 1 N ( h s i + 1 h s i ) 5 2 h s i h s i + 1 ( a ) 2 d t ,
where hsi is the time of the ith heel strike and a is the acceleration. The next step consists of bandpass filtering of the NJS using a fourth-order zero-lag Butterworth filter with a cutoff frequency of 5 Hz [28]. The higher the NJS, the greater the gait smoothness.

2.4.5. LLE

The LLE characterizes the behavior of a dynamical system. It quantifies gait stability as the average logarithmic rate of divergence following infinitesimal perturbations. An LLE less than zero represents the rate of convergence of the system’s trajectory to its nearest neighboring trajectory, whereas an LLE greater than zero represents the rate of divergence. When trajectories converge, the observed system is considered to have local dynamic stability, whereas divergence indicates local dynamic instability. The procedure described by Van Schooten et al. (2014) [69] was used to estimate the LLE in our study. To avoid the loss of spatiotemporal fluctuations and nonlinearities, no filtering was applied to the triaxial accelerations and the accelerations were time-normalized to obtain 100 data points per stride, thereby excluding the effects of the time duration of the data series on dynamic stability measures. Local dynamic stability was computed for each triaxial trunk acceleration over the strides considered in each trial. We assessed the short-term maximum finite-time LE ( λ m a x ) for each stride based on the AP, ML, and V trunk accelerations. The value of λ m a x was determined according to Rosenstein’s algorithm for short time series by using the Lyaprosen MATLAB toolbox for nonlinear time series analysis. A multidimensional state space, whose dimensions were determined using the classical global false-nearest-neighbor method, was reconstructed from the recorded one-dimensional time series data by juxtaposing the original data and delayed copies (the time delay was determined according to the first minimum of the average mutual information function) to evaluate dynamic perturbations [70]. Low values of this index indicate more stable trunk dynamics and high values indicate less stable trunk dynamics.

2.5. Clinical Assessment

The severity of PD was assessed using the HY disease staging system and the motor examination section of the Unified Parkinson’s Disease Rating Scale (UPDRS-III) [71]. Clinical scales were administered by an assessor unaware of the stability index reports. To identify fallers among the pwPD, participants were also asked to report the number of falls they experience over the past year.

2.6. Statistical Analysis

A needed sample of 62 subjects (31 pwPD and 31 HS) was calculated to identify gait stability indexes with good ability (area under the receiver operating characteristic (ROC) curve (AUC ≥ 0.70) for discriminating between pwPD and HS at a 95% significance level with 80% power under the null hypothesis of AUC ≤ 0.50. Based on the assumption of a high dropout rate as a result of the number of recorded strides required to be included in the analysis and matching procedures, we initially recruited 117 participants.
A one-to-one optimal data matching procedure using the propensity score difference method was implemented to match pwPD with HS [72]. We excluded 21 out of the 110 recorded HS trials based on an insufficient number of recorded strides. Therefore, 89 HS trials were used for the matching procedure. Propensity scores were calculated for each HY subgroup through logistic regression analysis using age [73,74] and speed [32,75,76,77] as covariates. Univariate analysis of variance (ANOVA) with Dunnett’s post-hoc analysis [78] was performed to verify the effectiveness of the age and speed matching procedures between pwPD at each HY stage and matched HS (HSmatched).
After checking the normality of the distributions and homoscedasticity of the variances using the Shapiro-Wilk test for all variables and Levene’s test, respectively, univariate ANOVA with Bonferroni’s post-hoc analysis and the Kruskal-Wallis H test with Dunn’s post-hoc analysis was performed to identify differences in the gait stability indexes between pwPD and HSmatched, as well as across the HY stages within the pwPD.
To assess the ability of the identified indexes to discriminate between pwPD and HSmatched, characterize the gait of pwPD according to the HY stages, and characterize the gait of recurrent fallers, ROC curves were plotted and AUCs were calculated with the presence of PD, HY stages, and a reported number of falls ≥5 as anchors [79,80,81], respectively. AUC values greater than 0.60 with a confidence interval lower bound greater than 0.50 were considered for sufficient overall discriminative ability [82]. The optimal cutoff points (OCPs) for each discriminative index were calculated as the points on the ROC curves that maximized the sum of sensitivity (Se) and specificity (Sp) values. Positive and negative likelihood ratios (LR+ and LR−) were also computed and transformed into post-test probabilities by using Fagan’s nomogram to analyze the probability of being correctly classified by a given index at the OCP [83]. Because the prevalence of recurrent fallers in our sample was low, to improve the generalizability of our results, the previously reported 39% prevalence of recurrent fallers in the general PD population [80,84] was used to calculate post-test probabilities for the identification of recurrent fallers [85].
Partial correlation coefficients (r) after removing the effects of age and gait speed were calculated to identify correlations between the identified indexes and UPDRS-III scores, history of falls, and spatiotemporal and kinematic parameters. Statistical analyses were performed using the IBM SPSS ver. 27 and NCSS 2019 software.

3. Results

3.1. Subjects and Matching Procedure Results

After the matching procedure, 55 walking trials from pwPD and 55 age-and-speed-matched walking trials from 30 HS (HSmatched) were included in our analysis. The final pwPD group consisted of 16 females and 39 males diagnosed with PD 8.20 ± 5.37 years ago and aged 71.17 ± 4.74 years. Twelve subjects were assessed at the HY = 1 disability stage, 19 subjects at the HY = 2 stage, and 24 at the HY = 3 stage. The mean UPDRS-III score was 39.48 ± 16.91. Nineteen subjects (34.4%) had experienced at least one fall in the previous year with 1.63 ± 4.54 falls an average. Seven subjects (12.7%) reported more than five falls during the past six months and were classified as recurrent fallers [79]. The mean gait speed of the included pwPD was 0.76 ± 0.23 m/s. The final HSmatched group consisted of 55 walking trials from a sample of 30 subjects (15 females and 15 males aged 70.11 ± 5.62 years). Eighteen of the 55 HSmatched trials were recorded at self-selected speeds (0.93 ± 0.38 m/s) and 37 were recorded at a reduced speed (0.65 ± 0.17 m/s). The mean speed of the final HSmatched group was 0.77 ± 0.25 m/s.
No significant overall differences between group means were identified for age (F(3, 105) = 2.81; p = 0.29) or gait speed (F(3, 106) = 2.49; p = 0.06) after the matching procedure. Post-hoc analysis revealed no significant differences for age (HY=1 vs. HS: p = 1.00; HY = 2 vs. HS: p = 1.00; HY = 3 vs. HS: p = 0.09; HY = 1 vs. HY = 2: p = 1.00; HY = 1 vs. HY = 3: p = 0.11; HY = 2 vs. HY = 3: p = 0.25) and gait speed (HY=1 vs. HS: p = 1.00; HY = 2 vs. HS: p = 1.00; HY = 3 vs. HS: p = 0.29; HY = 1 vs. HY = 2: p = 1.00; HY = 1 vs. HY = 3: p = 0.12; HY = 2 vs. HY = 3: p = 0.21) between the HY stage subgroups and HS, neither between the HY subgroups.

3.2. Discriminative Ability Results

Significant main effects of the group on HRAP (H = 9.48, p = 0.024), HRML (H = 10.24, p = 0.017), HRV (H = 14.30, p = 0.003), RQAdet in the AP direction (RQAdetAP) (H = 8.27, p = 0.041), and CV (H = 10.41, p = 0.015) were identified (Figure 2, Table 1).
Post-hoc analysis revealed significant differences in HRAP between pwPD at HY stage = 3 and HSmatched (p = 0.003), between pwPD at HY = 3 and HY = 1 (p = 0.043), and between pwPD at HY = 3 and HY = 2 (p = 0.032). A significant difference was identified in HRML between pwPD at HY = 3 and HSmatched (p = 0.001). Significant differences were identified in HRV between pwPD at HY = 3 and HSmatched, (p = 0.000), and between pwPD at HY = 3 and HY = 2 (p = 0.037). Significant differences were identified in RQAdetAP between pwPD at all HY stages and HSmatched (HY = 1, p = 0.015; HY = 2, p = 0.041, and HY = 3, p = 0.039). A significant difference in CV was identified between pwPD at HY = 3 and HSmatched (p = 0.004) (Table 1).
A good ability (AUC > 0.70) to discriminate between pwPD at HY = 3 and HSmatched was identified for HRAP, HRML, HRV, CV (Table 2, Figure 3 and Figure 4). A moderate ability (AUC = 0.65) to discriminate between pwPD at HY = (1, 2, 3) and HSmatched was identified for RQAdetAP (Table 2, Figure 4). HRAP values ≤ 1.50 identified pwPD at HY = 3 with 67% probability, HRML values ≤ 1.58 identified pwPD at HY = 3 with 54% probability, HRV values ≤ 1.74 identified pwPD at HY = 3 with 57% probability, and CV values ≥ 38.06 identified pwPD at HY = 3 with 58% probability. RQAdetAP values ≤ 38.85 identified pwPD with 67% probability, regardless of the HY stage.
HRAP exhibited a good ability to discriminate between pwPD at HY = 3 and HY = (1, 2) (AUC = 0.70) (Table 2, Figure 3). HRAP values ≤ 1.50 discriminated pwPD at HY = 3 from subjects at lower HY stages with 73% probability. Furthermore, HRAP exhibited a strong ability to identify recurrent fallers (AUC = 0.80) (Table 2, Figure 3). HRAP values ≤ 1.50 identified frequent fallers with 77% probability.

3.3. Partial Correlation Analysis Results

After removing the effects of gait speed and age, HRAP was negatively correlated with the history of falls (r = −0.45, p = 0.004) and positively correlated with pelvic obliquity (r = 0.37, p = 0.024) and pelvic rotation (r = 0.31, p = 0.040). RQAdetAP was negatively correlated with gait cadence (r = −0.35, p = 0.031) and positively correlated with stride time (r = 0.36, p = 0.030) and UPDRS III score (r = 0.385, p = 0.004). HRML, HRV, and CV exhibited no correlations with clinical features or gait kinematics.

4. Discussion

The main goal of this study was to determine the ability of 16 gait stability indexes derived from trunk acceleration signals to differentiate between pwPD and HS. By excluding the effects of age and gait speed, we found that the HRs calculated for the three spatial planes (AP, ML, and V) and the step length CV were able to discriminate between pwPD at moderate stages of disease progression (HY = 3) and HS with 67%, 54%, 57%, and 58% probabilities, respectively. HRAP was also able to discriminate between pwPD at the moderate disability stage (HY = 3) and pwPD at lower disability stages with 73% probability. Furthermore, HRAP exhibited good ability to characterize the gait of recurrent fallers with 77% probability. RQAdetAP was able to discriminate between pwPD and HSmatched, regardless of the HY score, with 67% probability. Overall, these results are in line with those from previous studies reporting that pwPD exhibit disruptions in the rhythmicity of pelvic acceleration [25,26] and a divergence from the typical quasi-periodic time-dependent recurrence of gait [86,87], particularly in the AP direction.
HR is a gait index that quantifies the smoothness of the trunk acceleration patterns in the three spatial directions during the gait and is defined as the ability to synchronize changes in trunk mechanics with those in the lower extremities [18]. In this study, the ability to control the trunk smoothly during the gait was found to be altered in pwPD at HY = 3. In particular, the HR in the AP direction was found to be correlated with a decrease in pelvic mobility. Additionally, we found that HRAP differentiated pwPD according to disease stages and was correlated with the history of falls, exhibiting a strong ability to identify recurrent fallers. Considering the progression of PD, HY stage 3 classifies subjects with PD with mild bilateral and axial involvement and initial postural instability. Therefore, HRAP may reflect the appearance of postural instability leading to falls and is correlated with axial rigidity. In this study, we could not directly measure trunk kinematics based on the localization of the inertial measurement unit. Therefore, information regarding lower trunk kinematics was inferred from pelvic kinematics [35]. In particular, the correlation coefficients revealed that HRAP was lower in subjects with lower pelvic obliquity and pelvic rotation. Alterations in trunk rotation and pelvic kinematics have been consistently described as characteristics of pwPD [55,88,89]. Based on our results, we can argue that the trunk rotation rigidity reflected by pelvic rotation and pelvic obliquity leads to alterations in AP trunk smoothness during the gait. This deficit in AP trunk control may be considered as a marker of postural instability that is independent of age and gait speed, as confirmed by our partial correlation analysis with a history of falls.
Although we determined that the HR in the ML and V directions can differentiate between pwPD at HY = 3 and HSmatched, the low discriminative power (<60%) of these variables for correctly identifying pwPD and their lack of correlations with clinical features disqualify them as accurate markers of gait instability in PD. These results are in line with previous research [25,26] evaluating the ability of HRML [26] and HRV [25] to differentiate between pwPD and HS, yet revealing a loss in their informative value when adjusting for gait speed [25,26]. Additionally, we did not find any correlations between HRML or HRV and a history of falls. In contrast, Cole et al. (2017) [32] identified significant differences in HR values in all planes between pwPD and age-matched healthy subjects, as well as a strong association between the HR values of recurrent fallers and gait speed. However, in their study, the experimental groups were not matched by gait speed. Therefore, direct comparisons to our results are difficult. The objective of this study was to determine the discriminative ability of HR values irrespective of gait speed. Therefore, our results indicate that HRML and HRV cannot be considered as independent markers of falling risk.
RQAdet is an index that expresses the predictability of acceleration trajectories during the gait. It represents how often a subject’s trunk accelerations revisit similar locations in the three spatial planes during their gait. Previous studies [86,87,90] have shown that RQAdet is the most convenient RQA index for objectively separating pwPD from HS. Accordingly, we identified lower RQAdet values in the AP direction in pwPD compared to HSmatched, regardless of the disability stage. Therefore, we can argue that pwPD exhibit disruptions in the quasi-periodic recurrence of their gaits in the early stages of PD, which reflects the early temporal gait alterations experienced by pwPD [91]. This consideration is further reinforced by the correlation we identified between RQAdetAP and temporal gait parameters, such as cadence and stride duration. Since the earliest stages of PD, pwPD exhibit shorter step lengths and stride times with an increased cadence [55,92] as a part of the abnormal gait pattern that progressively deteriorates into festination [7]. The reduced predictability of AP trunk accelerations may reflect systemic instability caused by altered temporal gait patterns. Therefore, RQAdetAP can be considered as a temporal marker of gait stability that can identify pwPD independently of gait speed. Furthermore, RQAdetAP was the only index correlated with the motor aspects of PD measured using UPDRS-III scores. Considering the variety of motor symptoms other than gait impairment assessed by UPDRS-III scores, a lack of correlation with gait-oriented measures is reasonable [93]. Therefore, the positive correlation between RQAdetAP and UPDRS-III may reflect its ability to capture gait impairment besides other symptoms such as reduced hand dexterity, facial mimicry, or altered posture.
The step length CV is a measure of gait spatial variability. Significant differences in the step length CV between pwPD and HS have been observed by other researchers [18,74,94]. These alterations have been attributed to compensatory behaviors based on the reduction of arm swing or to expressions of problems in terms of regulating the amplitude of step length during walking [95]. Regarding HRs, we found specific alterations in the step length CV only in pwPD with HY = 3. This result is in line with the results of previous studies that did not observe step length CV alterations in subjects at milder disease stages [94]. We speculate that in terms of HRs, CV alterations occur concomitantly with the progression of bilateral and axial involvement. However, the low probability of correctly discriminating pwPD (58%) at the cutoff value and a lack of correlations with clinical variables and spatiotemporal gait parameters do not allow us to consider this index as a marker of gait stability that is independent of speed.
Previous studies have reported [29,96] significantly lower NJS values in the ML and V directions in pwPD compared to age-matched HS. These differences have been interpreted as a reflection of PD-related bradykinesia, which parallels lower gait speed [96] and reduced arm swing movements [97]. Furthermore, the NJS in the ML and V directions has been reported to be responsive to dopaminergic medication in subjects with improved gait speed. In this study, we identified no differences in NJS between pwPD and age-and speed-matched HS, indicating that NJS is a marker of a general loss in complexity of the motor control system [29] that depends on gait speed.
We also identified no differences between groups in the LLE. In contrast, other studies have reported significant group differences between pwPD and HS [59] and between walking at self-selected speeds and dual-task walking [98] in terms of LLE. These contrasting results can be attributed to varying matching procedures and testing conditions. In one of these studies [98], no speed matching was performed between pwPD and HS because they were focused on observing changes in local dynamic stability induced by dual-task conditions, which does not correspond to our aim. Another study [59] observed differences when walking at comfortable speeds between pwPD and non-matched HS, as well as an increase in local dynamic instability when speeding up the gait of pwPD to match the HS gait speed. Given that the goals of the aforementioned studies were different, combining their interpretations with our results may suggest that the gait of pwPD is characterized by dynamic instability, but that LLE is speed dependent.
One of the main strengths of this study is that we investigated the discriminatory ability of the target indexes as independent diagnostic tools by excluding the influence of one of the gait features that differentiates pwPD from HS at the same age, namely gait speed [7]. Therefore, combined with the post-test probability results, the relatively low AUC values we observed may suggest that the identified indexes cannot be used as diagnostic markers for PD. However, these indexes may represent useful tools for assessing subtle gait alterations that clinical observations cannot detect. Furthermore, our results demonstrated that HRAP is a marker of gait alteration that can potentially identify recurrent fallers. Falls at a rate of more than one per year are frequently reported by elderly people [80]. The ability to detect recurrent fallers represents a key feature that a gait stability index should address because it can identify subtle alterations leading to falls in pwPD when compared to HS [80]. However, these results should be interpreted with caution. In this study, we only investigated gait patterns during the steady state. In pwPD, most falls occur during transitions such as turning or during the execution of concomitant cognitive tasks [84]. We attempted to overcome this limitation by using a high number of reported falls (greater than five) [80] as an anchor to discriminate recurrent fallers. Therefore, we can argue that HRAP is able to characterize the gait of recurrent fallers among the most disabled subjects, even during the analysis of linear walking without additional tasks. Regardless, the prevalence of fallers in our sample was lower than that in the general PD population, which can be attributed to our inclusion criteria. Because many predictors of falls, such as cognitive deficits, freezing of gait, depression, and orthopedic and peripheral conditions, were considered as exclusion criteria, this may have led to a lower prevalence of potential fallers, which could influence the calculation of post-test probabilities [85]. We compensated for this limitation by using the reported 39% prevalence of recurrent fallers in the general PD population in our calculations of post-test probabilities, thereby increasing the generalizability of our results to larger samples. Another limitation of this study is the retrospective self-reported history of falls, which could have led to recall bias. Therefore, additional longitudinal studies investigating the ability of HRAP to predict falls during the gait and more complex tasks should be conducted. Noteworthy, other acceleration-derived gait stability indexes, such as entropy measures or the maximal longest diagonal in recurrence quantification analysis, which we did not consider in this study, can be calculated to assess human gait instability. However, as these indexes are highly dependent on the number of recorded consecutive strides, other studies involving long-time gait series recordings are needed to analyze their ability to characterize the gait instability of pwPD.

5. Conclusions

In conclusion, the results of this study demonstrated that HRAP and RQAdetAP are age- and speed-independent markers of gait instability in subjects with PD. Both HRAP and RQAdetAP are able to discriminate pwPD from HS matched for age and gait speed. HRAP is altered in subjects at HY = 3, correlated with alterations in pelvic movements, and are able to characterize the gait instability of recurrent fallers. RQAdetAP captures early temporal gait alterations in pwPD but is not able to characterize the gait of fallers. Given the ease of use and relatively low cost of inertial motion sensors, these indices could provide useful clinical information regarding gait instability or alterations in the pace of subjects with PD, which are not directly observable through routine clinical assessments. Additional studies should be conducted to assess their ability to predict future falls, identify alterations during more demanding gait tasks, and their responsiveness to pharmacological and rehabilitative treatments.

Author Contributions

Conceptualization. S.F.C., M.S. and C.T.; methodology. S.F.C., M.S., A.R.; formal analysis. S.F.C., T.V. and M.S.; investigation. A.T., T.V., V.G., R.D.I. and F.M.; resources. A.D., D.D.L., D.F., C.T., M.S.; data curation. A.T., D.T., V.G. and R.D.I.; writing—original draft preparation. S.F.C., D.T. and M.S.; writing—review and editing. C.T., M.S., G.C., and A.R.; supervision. C.T. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee Lazio/2 (protocol N. 0053667/2021, date of approval: 03/10/2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author and stored in a password-protected PC located in the Department of Medico-Surgical Sciences and Biotechnologies, University of Rome Sapienza.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Morris, M.E.; Iansek, R.; Matyas, T.A.; Summers, J.J. The pathogenesis of gait hypokinesia in parkinson’s disease. Brain 1994, 117, 1169–1181. [Google Scholar] [CrossRef] [PubMed]
  2. Amboni, M.; Iuppariello, L.; Iavarone, A.; Fasano, A.; Palladino, R.; Rucco, R.; Picillo, M.; Lista, I.; Varriale, P.; Vitale, C.; et al. Step length predicts executive dysfunction in Parkinson’s disease: A 3-year prospective study. J. Neurol. 2018, 265, 2211–2220. [Google Scholar] [CrossRef]
  3. De Boer, A.G.E.M.; Wijker, W.; Speelman, J.D.; De Haes, J.C.J.M. Quality of life in patients with Parkinson’s disease: Development of a questionnaire. J. Neurol. Neurosurg. Psychiatry 1996, 61, 70–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Martínez-Martín, P. An introduction to the concept of “quality of life in Parkinson’s disease”. J. Neurol. 1998, 245 (Suppl. S1). [Google Scholar] [CrossRef]
  5. McCrone, P.; Allcock, L.M.; Burn, D.J. Predicting the cost of Parkinson’s disease. Mov. Disord. 2007, 22, 804–812. [Google Scholar] [CrossRef] [PubMed]
  6. Martinez-Martin, P.; Macaulay, D.; Jalundhwala, Y.J.; Mu, F.; Ohashi, E.; Marshall, T.; Sail, K. The long-term direct and indirect economic burden among Parkinson’s disease caregivers in the United States. Mov. Disord. 2019, 34, 236–245. [Google Scholar] [CrossRef] [PubMed]
  7. Mirelman, A.; Bonato, P.; Camicioli, R.; Ellis, T.D.; Giladi, N.; Hamilton, J.L.; Hass, C.J.; Hausdorff, J.M.; Pelosin, E.; Almeida, Q.J. Gait impairments in Parkinson’s disease. Lancet Neurol. 2019, 18, 697–708. [Google Scholar] [CrossRef]
  8. Espay, A.J.; Bonato, P.; Nahab, F.B.; Maetzler, W.; Dean, J.M.; Klucken, J.; Eskofier, B.M.; Merola, A.; Horak, F.; Lang, A.E.; et al. Technology in Parkinson’s disease: Challenges and opportunities. Mov. Disord. 2016, 31, 1272–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Espay, A.J.; Hausdorff, J.M.; Sánchez-Ferro, Á.; Klucken, J.; Merola, A.; Bonato, P.; Paul, S.S.; Horak, F.B.; Vizcarra, J.A.; Mestre, T.A.; et al. A roadmap for implementation of patient-centered digital outcome measures in Parkinson’s disease obtained using mobile health technologies. Mov. Disord. 2019, 34, 657–663. [Google Scholar] [CrossRef]
  10. Winser, S.J.; Kannan, P.; Bello, U.M.; Whitney, S.L. Measures of balance and falls risk prediction in people with Parkinson’s disease: A systematic review of psychometric properties. Clin. Rehabil. 2019, 33, 1949–1962. [Google Scholar] [CrossRef]
  11. Jacobs, J.V.; Horak, F.B.; Tran, V.K.; Nutt, J.G. Multiple balance tests improve the assessment of postural stability in subjects with Parkinson’s disease. J. Neurol. Neurosurg. Psychiatry 2006, 77, 322–326. [Google Scholar] [CrossRef] [Green Version]
  12. Bloem, B.R.; Marinus, J.; Almeida, Q.; Dibble, L.; Nieuwboer, A.; Post, B.; Ruzicka, E.; Goetz, C.; Stebbins, G.; Martinez-Martin, P.; et al. Measurement instruments to assess posture, gait, and balance in Parkinson’s disease: Critique and recommendations. Mov. Disord. 2016, 31, 1342–1355. [Google Scholar] [CrossRef] [PubMed]
  13. Sangarapillai, K.; Norman, B.M.; Almeida, Q.J. Rehabilitation of falls in parkinson’s disease: Self-perception vs. objective measures of fall risk. Brain Sci. 2021, 11, 320. [Google Scholar] [CrossRef] [PubMed]
  14. Dingwell, J.B.; Cusumano, J.P. Nonlinear time series analysis of normal and pathological human walking. Chaos 2000, 10, 848–863. [Google Scholar] [CrossRef] [PubMed]
  15. England, S.A.; Granata, K.P. The influence of gait speed on local dynamic stability of walking. Gait Posture 2007, 25, 172–178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hamacher, D.; Singh, N.B.; Van Dieën, J.H.; Heller, M.O.; Taylor, W.R. Kinematic measures for assessing gait stability in elderly individuals: A systematic review. J. R. Soc. Interface 2011, 8, 1682–1698. [Google Scholar] [CrossRef] [PubMed]
  17. Bruijn, S.M.; Meijer, O.G.; Beek, P.J.; Van Dieen, J.H. Assessing the stability of human locomotion: A review of current measures. J. R. Soc. Interface 2013, 10, 20120999. [Google Scholar] [CrossRef]
  18. Siragy, T.; Nantel, J. Quantifying Dynamic Balance in Young, Elderly and Parkinson’s Individuals: A Systematic Review. Front. Aging Neurosci. 2018, 10, 387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Frenkel-Toledo, S.; Giladi, N.; Peretz, C.; Herman, T.; Gruendlinger, L.; Hausdorff, J.M. Effect of gait speed on gait rhythmicity in Parkinson’s disease: Variability of stride time and swing time respond differently. J. Neuroeng. Rehabil. 2005, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Baltadjieva, R.; Giladi, N.; Gruendlinger, L.; Peretz, C.; Hausdorff, J.M. Marked alterations in the gait timing and rhythmicity of patients with de novo Parkinson’s disease. Eur. J. Neurosci. 2006, 24, 1815–1820. [Google Scholar] [CrossRef]
  21. Auriel, E.; Hausdorff, J.M.; Herman, T.; Simon, E.S.; Giladi, N. Effects of methylphenidate on cognitive function and gait in patients with Parkinson’s disease: A pilot study. Clin. Neuropharmacol. 2006, 29, 15–17. [Google Scholar] [CrossRef]
  22. Plotnik, M.; Giladi, N.; Hausdorff, J.M. A new measure for quantifying the bilateral coordination of human gait: Effects of aging and Parkinson’s disease. Exp. Brain Res. 2007, 181, 561–570. [Google Scholar] [CrossRef] [PubMed]
  23. Herman, T.; Giladi, N.; Gruendlinger, L.; Hausdorff, J.M. Six Weeks of Intensive Treadmill Training Improves Gait and Quality of Life in Patients with Parkinson’s Disease: A Pilot Study. Arch. Phys. Med. Rehabil. 2007, 88, 1154–1158. [Google Scholar] [CrossRef] [PubMed]
  24. Cole, M.H.; Silburn, P.A.; Wood, J.M.; Worringham, C.J.; Kerr, G.K. Falls in Parkinson’s disease: Kinematic evidence for impaired head and trunk control. Mov. Disord. 2010, 25, 2369–2378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Latt, M.D.; Menz, H.B.; Fung, V.S.; Lord, S.R. Acceleration patterns of the head and pelvis during gait in older people with Parkinson’s disease: A comparison of fallers and nonfallers. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2009, 64, 700–706. [Google Scholar] [CrossRef] [Green Version]
  26. Lowry, K.A.; Smiley-Oyen, A.L.; Carrel, A.J.; Kerr, J.P. Walking stability using harmonic ratios in Parkinson’s disease. Mov. Disord. 2009, 24, 261–267. [Google Scholar] [CrossRef]
  27. Thumm, P.C.; Maidan, I.; Brozgol, M.; Shustak, S.; Gazit, E.; Shema Shiratzki, S.; Bernad-Elazari, H.; Beck, Y.; Giladi, N.; Hausdorff, J.M.; et al. Treadmill walking reduces pre-frontal activation in patients with Parkinson’s disease. Gait Posture 2018, 62, 384–387. [Google Scholar] [CrossRef]
  28. Miller Koop, M.; Ozinga, S.J.; Rosenfeldt, A.B.; Alberts, J.L. Quantifying turning behavior and gait in Parkinson’s disease using mobile technology. IBRO Rep. 2018, 5, 10–16. [Google Scholar] [CrossRef]
  29. Palmerini, L.; Mellone, S.; Avanzolini, G.; Valzania, F.; Chiari, L. Quantification of motor impairment in Parkinson’s disease using an instrumented timed up and go test. IEEE Trans. Neural Syst. Rehabil. Eng. 2013, 21, 664–673. [Google Scholar] [CrossRef]
  30. Sylos Labini, F.; Meli, A.; Ivanenko, Y.P.; Tufarelli, D. Recurrence quantification analysis of gait in normal and hypovestibular subjects. Gait Posture 2012, 35, 48–55. [Google Scholar] [CrossRef]
  31. Ramdani, S.; Tallon, G.; Bernard, P.L.; Blain, H. Recurrence quantification analysis of human postural fluctuations in older fallers and non-fallers. Ann. Biomed. Eng. 2013, 41, 1713–1725. [Google Scholar] [CrossRef]
  32. Cole, M.H.; Sweeney, M.; Conway, Z.J.; Blackmore, T.; Silburn, P.A. Imposed Faster and Slower Walking Speeds Influence Gait Stability Differently in Parkinson Fallers. Arch. Phys. Med. Rehabil. 2017, 98, 639–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fukuchi, C.A.; Fukuchi, R.K.; Duarte, M. Effects of walking speed on gait biomechanics in healthy participants: A systematic review and meta-analysis. Syst. Rev. 2019, 8, 153. [Google Scholar] [CrossRef] [Green Version]
  34. Zampogna, A.; Mileti, I.; Palermo, E.; Celletti, C.; Paoloni, M.; Manoni, A.; Mazzetta, I.; Costa, G.D.; Pérez-López, C.; Camerota, F.; et al. Fifteen years of wireless sensors for balance assessment in neurological disorders. Sensors 2020, 20, 3247. [Google Scholar] [CrossRef] [PubMed]
  35. Brognara, L.; Palumbo, P.; Grimm, B.; Palmerini, L. Assessing Gait in Parkinson’s Disease Using Wearable Motion Sensors: A Systematic Review. Diseases 2019, 7, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Schlachetzki, J.C.M.; Barth, J.; Marxreiter, F.; Gossler, J.; Kohl, Z.; Reinfelder, S.; Gassner, H.; Aminian, K.; Eskofier, B.M.; Winkler, J.; et al. Wearable sensors objectively measure gait parameters in Parkinson’s disease. PLoS ONE 2017, 12, e0183989. [Google Scholar] [CrossRef]
  37. Rovini, E.; Maremmani, C.; Cavallo, F. How wearable sensors can support parkinson’s disease diagnosis and treatment: A systematic review. Front. Neurosci. 2017, 11, 555. [Google Scholar] [CrossRef]
  38. Maetzler, W.; Domingos, J.; Srulijes, K.; Ferreira, J.J.; Bloem, B.R. Quantitative wearable sensors for objective assessment of Parkinson’s disease. Mov. Disord. 2013, 28, 1628–1637. [Google Scholar] [CrossRef]
  39. Ramdhani, R.A.; Khojandi, A.; Shylo, O.; Kopell, B.H. Optimizing clinical assessments in Parkinson’s disease through the use of wearable sensors and data driven modeling. Front. Comput. Neurosci. 2018, 12, 72. [Google Scholar] [CrossRef] [Green Version]
  40. Kuo, A.D.; Donelan, J.M. Dynamic principles of gait and their clinical implications. Phys. Ther. 2010, 90, 157–174. [Google Scholar] [CrossRef] [Green Version]
  41. Horak, F.B.; Mancini, M. Objective biomarkers of balance and gait for Parkinson’s disease using body-worn sensors. Mov. Disord. 2013, 28, 1544–1551. [Google Scholar] [CrossRef] [Green Version]
  42. Czech, M.; Demanuele, C.; Erb, M.K.; Ramos, V.; Zhang, H.; Ho, B.; Patel, S. The impact of reducing the number of wearable devices on measuring gait in parkinson disease: Noninterventional exploratory study. JMIR Rehabil. Assist. Technol. 2020, 7. [Google Scholar] [CrossRef] [PubMed]
  43. Hughes, A.J.; Daniel, S.E.; Kilford, L.; Lees, A.J. Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: A clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatry 1992, 55, 181–184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hoehn, M.M.; Yahr, M.D. Parkinsonism: Onset, progression, and mortality. Neurology 1967, 17, 427–442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kroneberg, D.; Elshehabi, M.; Meyer, A.C.; Otte, K.; Doss, S.; Paul, F.; Nussbaum, S.; Berg, D.; Kühn, A.A.; Maetzler, W.; et al. Less is more—Estimation of the number of strides required to assess gait variability in spatially confined settings. Front. Aging Neurosci. 2019, 11, 435. [Google Scholar] [CrossRef] [Green Version]
  46. Riva, F.; Bisi, M.C.; Stagni, R. Gait variability and stability measures: Minimum number of strides and within-session reliability. Comput. Biol. Med. 2014, 50, 9–13. [Google Scholar] [CrossRef] [PubMed]
  47. 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. 2017, 53, 84–89. [Google Scholar] [CrossRef] [PubMed]
  48. Folstein, M.F.; Folstein, S.E.; McHugh, P.R. “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res. 1975, 12, 189–198. [Google Scholar] [CrossRef]
  49. Perneczky, R.; Wagenpfeil, S.; Komossa, K.; Grimmer, T.; Diehl, J.; Kurz, A. Mapping scores onto stages: Mini-mental state examination and clinical dementia rating. Am. J. Geriatr. Psychiatry 2006, 14, 139–144. [Google Scholar] [CrossRef] [PubMed]
  50. Beck, A.T.; Ward, C.H.; Mendelson, M.; Mock, J.; Erbaugh, J. An Inventory for Measuring Depression. Arch. Gen. Psychiatry 1961, 4, 561–571. [Google Scholar] [CrossRef] [Green Version]
  51. Goodarzi, Z.; Mrklas, K.J.; Roberts, D.J.; Jette, N.; Pringsheim, T.; Holroyd-Leduc, J. Detecting depression in Parkinson disease: A systematic review and meta-analysis. Neurology 2016, 87, 426–437. [Google Scholar] [CrossRef]
  52. Altman, R.; Alarcón, G.; Appelrouth, D.; Bloch, D.; Borenstein, D.; Brandt, K.; Brown, C.; Cooke, T.D.; Daniel, W.; Feldman, D.; et al. The American College of Rheumatology criteria for the classification and reporting of osteoarthritis of the hip. Arthritis Rheum. 1991, 34, 505–514. [Google Scholar] [CrossRef] [PubMed]
  53. Fitzgerald, G.K.; Hinman, R.S.; Zeni, J.; Risberg, M.A.; Snyder-Mackler, L.; Bennell, K.L. OARSI Clinical Trials Recommendations: Design and conduct of clinical trials of rehabilitation interventions for osteoarthritis. Osteoarthr. Cartil. 2015, 23, 803–814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Damen, J.; Van Rijn, R.M.; Emans, P.J.; Hilberdink, W.K.H.A.; Wesseling, J.; Oei, E.H.G.; Bierma-Zeinstra, S.M.A. Prevalence and development of hip and knee osteoarthritis according to American College of Rheumatology criteria in the CHECK cohort. Breast Cancer Res. 2019, 21. [Google Scholar] [CrossRef] [Green Version]
  55. Serrao, M.; Chini, G.; Caramanico, G.; Bartolo, M.; Castiglia, S.F.; Ranavolo, A.; Conte, C.; Venditto, T.; Coppola, G.; Di Lorenzo, C.; et al. Prediction of responsiveness of gait variables to rehabilitation training in Parkinson’s disease. Front. Neurol. 2019, 10, 826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Rinaldi, M.; Ranavolo, A.; Conforto, S.; Martino, G.; Draicchio, F.; Conte, C.; Varrecchia, T.; Bini, F.; Casali, C.; Pierelli, F.; et al. Increased lower limb muscle coactivation reduces gait performance and increases metabolic cost in patients with hereditary spastic paraparesis. Clin. Biomech. 2017, 48, 63–72. [Google Scholar] [CrossRef]
  57. Mari, S.; Serrao, M.; Casali, C.; Conte, C.; Martino, G.; Ranavolo, A.; Coppola, G.; Draicchio, F.; Padua, L.; Sandrini, G.; et al. Lower limb antagonist muscle co-activation and its relationship with gait parameters in cerebellar ataxia. Cerebellum 2014, 13, 226–236. [Google Scholar] [CrossRef] [PubMed]
  58. Cofré, L.E.; Lythgo, N.; Morgan, D.; Galea, M.P. Aging modifies joint power and work when gait speeds are matched. Gait Posture 2011, 33, 484–489. [Google Scholar] [CrossRef]
  59. Peterson, D.S.; Mancini, M.; Fino, P.C.; Horak, F.; Smulders, K. Speeding Up Gait in Parkinson’s Disease. J. Parkinsons Dis. 2020, 10, 245–253. [Google Scholar] [CrossRef]
  60. Smidt, G.L. Methods of studying gait. Phys. Ther. 1974, 54, 13–17. [Google Scholar] [CrossRef]
  61. Iosa, M.; Picerno, P.; Paolucci, S.; Morone, G. Wearable inertial sensors for human movement analysis. Expert Rev. Med. Devices 2016, 13, 641–659. [Google Scholar] [CrossRef]
  62. Webber, C.L.; Zbilut, J.P. Dynamical assessment of physiological systems and states using recurrence plot strategies. J. Appl. Physiol. 1994, 76, 965–973. [Google Scholar] [CrossRef] [PubMed]
  63. Kennel, M.B.; Abarbanel, H.D.I. False neighbors and false strands: A reliable minimum embedding dimension algorithm. Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Top. 2002, 66. [Google Scholar] [CrossRef] [PubMed]
  64. 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. 2018, 9, 1679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Toebes, M.J.P.; Hoozemans, M.J.M.; Furrer, R.; Dekker, J.; Van Dieën, J.H. Local dynamic stability and variability of gait are associated with fall history in elderly subjects. Gait Posture 2012, 36, 527–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Fraser, A.M.; Swinney, H.L. Independent coordinates for strange attractors from mutual information. Phys. Rev. A 1986, 33, 1134–1140. [Google Scholar] [CrossRef] [PubMed]
  67. Zijlstra, W. Assessment of spatio-temporal parameters during unconstrained walking. Eur. J. Appl. Physiol. 2004, 92, 39–44. [Google Scholar] [CrossRef]
  68. Serrao, M.; Pierelli, F.; Ranavolo, A.; Draicchio, F.; Conte, C.; Don, R.; Di Fabio, R.; Lerose, M.; Padua, L.; Sandrini, G.; et al. Gait pattern in inherited cerebellar ataxias. Cerebellum 2012, 11, 194–211. [Google Scholar] [CrossRef]
  69. Van Schooten, K.S.; Rispens, S.M.; Elders, P.J.M.; van Dieën, J.H.; Pijnappels, M. Toward ambulatory balance assessment: Estimating variability and stability from short bouts of gait. Gait Posture 2014, 39, 695–699. [Google Scholar] [CrossRef] [Green Version]
  70. Chini, G.; Ranavolo, A.; Draicchio, F.; Casali, C.; Conte, C.; Martino, G.; Leonardi, L.; Padua, L.; Coppola, G.; Pierelli, F.; et al. Local Stability of the Trunk in Patients with Degenerative Cerebellar Ataxia During Walking. Cerebellum 2017, 16, 26–33. [Google Scholar] [CrossRef]
  71. Goetz, C.G.; Tilley, B.C.; Shaftman, S.R.; Stebbins, G.T.; Fahn, S.; Martinez-Martin, P.; Poewe, W.; Sampaio, C.; Stern, M.B.; Dodel, R.; et al. Movement Disorder Society-Sponsored Revision of the Unified Parkinson’s Disease Rating Scale (MDS-UPDRS): Scale presentation and clinimetric testing results. Mov. Disord. 2008, 23, 2129–2170. [Google Scholar] [CrossRef] [PubMed]
  72. Yao, X.I.; Wang, X.; Speicher, P.J.; Hwang, E.S.; Cheng, P.; Harpole, D.H.; Berry, M.F.; Schrag, D.; Pang, H.H. Reporting and Guidelines in Propensity Score Analysis: A Systematic Review of Cancer and Cancer Surgical Studies. J. Natl. Cancer Inst. 2017, 109, djw323. [Google Scholar] [CrossRef] [PubMed]
  73. Nemanich, S.T.; Duncan, R.P.; Dibble, L.E.; Cavanaugh, J.T.; Ellis, T.D.; Ford, M.P.; Foreman, K.B.; Earhart, G.M. Predictors of gait speeds and the relationship of gait speeds to falls in men and women with parkinson disease. Parkinsons Dis. 2013, 2013. [Google Scholar] [CrossRef] [PubMed]
  74. Noh, B.; Youm, C.; Lee, M.; Cheon, S.M. Gait characteristics in individuals with Parkinson’s disease during 1-minute treadmill walking. PeerJ 2020, 2020. [Google Scholar] [CrossRef]
  75. Lindemann, U. Spatiotemporal gait analysis of older persons in clinical practice and research: Which parameters are relevant? Z. Gerontol. Geriatr. 2020, 53, 171–178. [Google Scholar] [CrossRef] [PubMed]
  76. Huijben, B.; van Schooten, K.S.; van Dieën, J.H.; Pijnappels, M. The effect of walking speed on quality of gait in older adults. Gait Posture 2018, 65, 112–116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Craig, J.J.; Bruetsch, A.P.; Huisinga, J.M. Coordination of trunk and foot acceleration during gait is affected by walking velocity and fall history in elderly adults. Aging Clin. Exp. Res. 2019, 31, 943–950. [Google Scholar] [CrossRef]
  78. Lee, S.; Lee, D.K. What is the proper way to apply the multiple comparison test? Korean J. Anesthesiol. 2018, 71, 353–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Paul, S.S.; Allen, N.E.; Sherrington, C.; Heller, G.; Fung, V.S.C.; Close, J.C.T.; Lord, S.R.; Canning, C.G. Risk factors for frequent falls in people with Parkinson’s disease. J. Parkinsons Dis. 2014, 4, 699–703. [Google Scholar] [CrossRef] [PubMed]
  80. Allen, N.E.; Schwarzel, A.K.; Canning, C.G. Recurrent falls in parkinson’s disease: A systematic review. Parkinsons Dis. 2013, 2013, 906274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Thomas, A.A.; Rogers, J.M.; Amick, M.M.; Friedman, J.H. Falls and the falls efficacy scale in Parkinson’s disease. J. Neurol. 2010, 257, 1124–1128. [Google Scholar] [CrossRef] [PubMed]
  82. Carter, J.V.; Pan, J.; Rai, S.N.; Galandiuk, S. ROC-ing along: Evaluation and interpretation of receiver operating characteristic curves. Surgery 2016, 159, 1638–1645. [Google Scholar] [CrossRef] [PubMed]
  83. Kallner, A. Bayes’ theorem, the roc diagram and reference values: Definition and use in clinical diagnosis. Biochem. Med. 2018, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Fasano, A.; Canning, C.G.; Hausdorff, J.M.; Lord, S.; Rochester, L. Falls in Parkinson’s disease: A complex and evolving picture. Mov. Disord. 2017, 32, 1524–1536. [Google Scholar] [CrossRef]
  85. Eusebi, P. Diagnostic accuracy measures. Cerebrovasc. Dis. 2013, 36, 267–272. [Google Scholar] [CrossRef] [PubMed]
  86. Afsar, O.; Tirnakli, U.; Marwan, N. Recurrence Quantification Analysis at work: Quasi-periodicity based interpretation of gait force profiles for patients with Parkinson disease. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef] [PubMed]
  87. Afşar, Ö. Recurrence Quantification Analysis on Gait Reaction Forces of Elderly Adults for Determination of Pathological States. Celal Bayar Univ. Fen Bilim. Derg. 2018, 14, 309–314. [Google Scholar] [CrossRef]
  88. Varrecchia, T.; Castiglia, S.F.; Ranavolo, A.; Conte, C.; Tatarelli, A.; Coppola, G.; Di Lorenzo, C.; Draicchio, F.; Pierelli, F.; Serrao, M. An artificial neural network approach to detect presence and severity of Parkinson’s disease via gait parameters. PLoS ONE 2021, 16, e0244396. [Google Scholar] [CrossRef]
  89. Creaby, M.W.; Cole, M.H. Gait characteristics and falls in Parkinson’s disease: A systematic review and meta-analysis. Park. Relat. Disord. 2018, 57, 1–8. [Google Scholar] [CrossRef]
  90. Pham, T.D. Pattern analysis of computer keystroke time series in healthy control and early-stage Parkinson’s disease subjects using fuzzy recurrence and scalable recurrence network features. J. Neurosci. Methods 2018, 307, 194–202. [Google Scholar] [CrossRef]
  91. Djurić-Jovičić, M.; Belić, M.; Stanković, I.; Radovanović, S.; Kostić, V.S. Selection of gait parameters for differential diagnostics of patients with de novo Parkinson’s disease. Neurol. Res. 2017, 39, 853–861. [Google Scholar] [CrossRef]
  92. Kwon, K.Y.; Lee, H.M.; Kang, S.H.; Pyo, S.J.; Kim, H.J.; Koh, S.B. Recuperation of slow walking in de novo Parkinson’s disease is more closely associated with increased cadence, rather than with expanded stride length. Gait Posture 2017, 58, 1–6. [Google Scholar] [CrossRef]
  93. Mancini, M.; Carlson-Kuhta, P.; Zampieri, C.; Nutt, J.G.; Chiari, L.; Horak, F.B. Postural sway as a marker of progression in Parkinson’s disease: A pilot longitudinal study. Gait Posture 2012, 36, 471–476. [Google Scholar] [CrossRef] [Green Version]
  94. Bovonsunthonchai, S.; Vachalathiti, R.; Pisarnpong, A.; Khobhun, F.; Hiengkaew, V. Spatiotemporal Gait Parameters for Patients with Parkinson’s Disease Compared with Normal Individuals. Physiother. Res. Int. 2014, 19, 158–165. [Google Scholar] [CrossRef]
  95. Mak, M.K.Y. Reduced step length, not step length variability is central to gait hypokinesia in people with Parkinson’s disease. Clin. Neurol. Neurosurg. 2013, 115, 587–590. [Google Scholar] [CrossRef] [PubMed]
  96. Warabi, T.; Furuyama, H.; Kato, M. Gait bradykinesia: Difficulty in switching posture/gait measured by the anatomical y-axis vector of the sole in Parkinson’s disease. Exp. Brain Res. 2020, 238, 139–151. [Google Scholar] [CrossRef] [PubMed]
  97. Zampier, V.C.; Vitório, R.; Beretta, V.S.; Jaimes, D.A.R.; Orcioli-Silva, D.; Santos, P.C.R.; Gobbi, L.T.B. Gait bradykinesia and hypometria decrease as arm swing frequency and amplitude increase. Neurosci. Lett. 2018, 687, 248–252. [Google Scholar] [CrossRef] [PubMed]
  98. Fino, P.C.; Mancini, M.; Curtze, C.; Nutt, J.G.; Horak, F.B. Gait stability has phase-dependent dual-task costs in Parkinson’s disease. Front. Neurol. 2018, 9, 373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Graphical representation of the accelerations-derived gait indexes of a representative healthy subject: (a) amplitudes of the filtered acceleration signals in the antero-posterior (AP), medio-lateral (ML), and vertical (V) direction as a function of the time range data; (b) Harmonic Ratio values for each of the 20 considered strides; (c) 2D-reconstructed state space of the acceleration and its time-delayed copies (time delay of 12 data samples). The distance (d) of two neighboring trajectories at a one-time sample, which is needed to calculate the Lyapunov exponent, is highlighted; (d) representation of the jerks during the whole gait cycle; (e) recurrence matrix. Based on the percent of the recurrent points in the diagonal line structure parallel to the main diagonal (i.e., the blue circled points), the RQAdet was calculated.
Figure 1. Graphical representation of the accelerations-derived gait indexes of a representative healthy subject: (a) amplitudes of the filtered acceleration signals in the antero-posterior (AP), medio-lateral (ML), and vertical (V) direction as a function of the time range data; (b) Harmonic Ratio values for each of the 20 considered strides; (c) 2D-reconstructed state space of the acceleration and its time-delayed copies (time delay of 12 data samples). The distance (d) of two neighboring trajectories at a one-time sample, which is needed to calculate the Lyapunov exponent, is highlighted; (d) representation of the jerks during the whole gait cycle; (e) recurrence matrix. Based on the percent of the recurrent points in the diagonal line structure parallel to the main diagonal (i.e., the blue circled points), the RQAdet was calculated.
Sensors 21 03449 g001
Figure 2. (a) Graphical representation of the Harmonic Ratios in the antero-posterior, medio-lateral, and vertical directions of a representative age-and-speed-matched healthy subject (blue) and a subject with PD at Hoehn and Yahr stage = 3 (red); (b) recurrence matrices in the antero-posterior direction of the same representative subjects.
Figure 2. (a) Graphical representation of the Harmonic Ratios in the antero-posterior, medio-lateral, and vertical directions of a representative age-and-speed-matched healthy subject (blue) and a subject with PD at Hoehn and Yahr stage = 3 (red); (b) recurrence matrices in the antero-posterior direction of the same representative subjects.
Sensors 21 03449 g002
Figure 3. ROC curves for the HRs in identifying pwPD vs. HSmatched, pwPD at HY = 3 from milder HY and recurrent fallers. The red line represents the HRAP, the blue line the HRML, and the green line the HRV.
Figure 3. ROC curves for the HRs in identifying pwPD vs. HSmatched, pwPD at HY = 3 from milder HY and recurrent fallers. The red line represents the HRAP, the blue line the HRML, and the green line the HRV.
Sensors 21 03449 g003
Figure 4. ROC curves of the CV (a) and RQAdetAP (b) in discriminating pwPD from HSmatched.
Figure 4. ROC curves of the CV (a) and RQAdetAP (b) in discriminating pwPD from HSmatched.
Sensors 21 03449 g004
Table 1. Comparison of the gait stability indexes between subjects with PD and HSmatched.
Table 1. Comparison of the gait stability indexes between subjects with PD and HSmatched.
HY 1HY 2HY 3HSmatched
Mean (SD)Mean (SD)Mean (SD)Mean (SD)
HRAP2.00 (0.55)1.94 (0.51)1.65 (0.34) * #2.08 (0.67)
HRML1.82 (0.49)1.84 (0.48)1.63 (0.37) *1.97 (0.51)
HRV1.93 (0.53)1.86 (0.41)1.64 (0.37) *2.05 (0.50)
RQArecAP3.13 (2.04)4.21 (2.91)4.35 (4.10)5.24 (3.54)
RQArecML4.43 (4.83)3.80 (4.08)2.53 (3.77)3.86 (3.31)
RQArecV4.73 (5.50)3.44 (4.65)5.38 (7.62)2.76 (2.25)
RQAdetAP23.20 (15.94) *30.23 (17.89) *28.37 (22.74) *41.75 (27.07)
RQAdetML31.09 (26.35)33.93 (26.27)22.63 (25.10)34.93 (25.38)
RQAdetV27.17 (28.83)20.61 (20.05)25.40 (23.67)23.78 (16.04)
CV39.36 (17.85)35.56 (16.94)40.92 (18.30) *28.72 (14.08)
NJSAP2261.67 (1387.79)4157.58 (3465.20)4061.89 (2839.03)3512.49 (3047.39)
NJSML1169.38 (959.77)1518.38 (1115.65)1622.29 (1217.89)1558.43 (1221.20)
NJSV1127.09 (808.44)1526.17 (1166.67)1591.34 (1925.34)1656.37 (1169.45)
LLEAP0.53 (0.26)0.49 (0.20)0.60 (0.23)0.53 (0.26)
LLEML0.64 (0.23)0.58 (0.18)0.63 (0.20)0.54 (0.28)
LLEV0.82 (0.28)0.88 (0.21)0.88 (0.21)0.86 (0.29)
* significant differences between subjects with PD and HSmatched; #, significant differences between subjects with PD at HY 1,2 and HY 3; HY, Hoehn and Yahr disease staging classification system; HR, harmonic ratio; RQArec, recurrence quantification analysis, percent recurrence; RQAdet, recurrence quantification analysis, percent determinism; CV, coefficient of variation of the step length; NJS, normalized Jerk score; LLE, largest Lyapunov exponent.
Table 2. ROC curve and cutoff analysis results.
Table 2. ROC curve and cutoff analysis results.
Gait IndexSubjectsAUC (95% CI)OCPSe (95% CI)Sp (95% CI)LR+LR−PTP+PTP−
HRAPHY = 3 vs. HSmatched0.71 (0.56–0.81)≤1.500.42 (0.22–0.63)0.91 (0.80–0.97)4.580.6467%22%
HY = 3 vs. HY = 1.20.70 (0.53–0.80)≤1.500.46 (0.25–0.67)0.87 (0.70–0.96)3.550.6273%32%
recurrent fallers (≥5)0.80 (0.52–0.92)≤1.500.85 (0.42–0.98)0.83 (0.70–0.92)5.140.1777%10%
HRMLHY = 3 vs. HSmatched0.72 (0.56–0.82)≤1.580.58 (0.36–0.78)0.78 (0.65–0.88)2.670.5354%19%
HRVHY = 3 vs. HSmatched0.76 (0.61–0.86)≤1.740.71 (0.49–0.87)0.76 (0.63–0.87)3.000.3857%14%
RQAdetAPAll subjects with PD vs. HSmatched0.65 (0.53–0.75)≤38.850.74 (0.61–0.85)0.63 (0.49–0.76)2.050.4067%29%
CVHY = 3 vs. HSmatched0.72 (0.57–0.82)≥38.060.58 (0.36–0.78)0.81 (0.69–0.91)3.210.5158%18%
AUC, area under the ROC curve; OCP, optimal cutoff point; Se, sensitivity; Sp, specificity; LR+, positive likelihood ratio; LR−, negative likelihood ratio; PTP+, positive post-test probability; PTP−, negative post-test probability; HR, harmonic ratio; RQAdetAP. Recurrence quantification analysis, percent determinism in the antero-posterior direction; CV, coefficient of variation. of the step length.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Castiglia, S.F.; Tatarelli, A.; Trabassi, D.; De Icco, R.; Grillo, V.; Ranavolo, A.; Varrecchia, T.; Magnifica, F.; Di Lenola, D.; Coppola, G.; et al. Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson’s Disease. Sensors 2021, 21, 3449. https://doi.org/10.3390/s21103449

AMA Style

Castiglia SF, Tatarelli A, Trabassi D, De Icco R, Grillo V, Ranavolo A, Varrecchia T, Magnifica F, Di Lenola D, Coppola G, et al. Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson’s Disease. Sensors. 2021; 21(10):3449. https://doi.org/10.3390/s21103449

Chicago/Turabian Style

Castiglia, Stefano Filippo, Antonella Tatarelli, Dante Trabassi, Roberto De Icco, Valentina Grillo, Alberto Ranavolo, Tiwana Varrecchia, Fabrizio Magnifica, Davide Di Lenola, Gianluca Coppola, and et al. 2021. "Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson’s Disease" Sensors 21, no. 10: 3449. https://doi.org/10.3390/s21103449

APA Style

Castiglia, S. F., Tatarelli, A., Trabassi, D., De Icco, R., Grillo, V., Ranavolo, A., Varrecchia, T., Magnifica, F., Di Lenola, D., Coppola, G., Ferrari, D., Denaro, A., Tassorelli, C., & Serrao, M. (2021). Ability of a Set of Trunk Inertial Indexes of Gait to Identify Gait Instability and Recurrent Fallers in Parkinson’s Disease. Sensors, 21(10), 3449. https://doi.org/10.3390/s21103449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop