Principal component analysis for ataxic gait using a triaxial accelerometer
Journal of NeuroEngineering and Rehabilitation volume 14, Article number: 37 (2017)
It is quite difficult to evaluate ataxic gait quantitatively in clinical practice. The aim of this study was to analyze the characteristics of ataxic gait using a triaxial accelerometer and to develop a novel biomarker of integrated gate parameters for ataxic gait.
Sixty-one patients with spinocerebellar ataxia (SCA) or multiple system atrophy with predominant cerebellar ataxia (MSA-C) and 57 healthy control subjects were enrolled. The subjects were instructed to walk 10 m for a total of 12 times on a flat floor at their usual walking speed with a triaxial accelerometer attached to their back. Gait velocity, cadence, step length, step regularity, step symmetry, and degree of body sway were evaluated. Principal component analysis (PCA) was used to analyze the multivariate gait parameters. The Scale for the Assessment and Rating of Ataxia (SARA) was evaluated on the same day of the 10-m walk trial.
PCA divided the gait parameters into four principal components in the controls and into two principal components in the patients. The four principal components in the controls were similar to those found in earlier studies. The second principal component in the patients had relevant factor loading values for gait velocity, step length, regularity, and symmetry in addition to the degree of body sway in the medio-lateral direction. The second principal component score (PCS) in the patients was significantly correlated with disease duration and the SARA score of gait (ρ = −0.363, p = 0.004; ρ = −0.574, p < 0.001, respectively).
PCA revealed the main component of ataxic gait. The PCS of the main component was significantly different between the patients and controls, and it was well correlated with disease duration and the SARA score of gait in the patients. We propose that this score provides a novel method to assess the severity of ataxic gait quantitatively using a triaxial accelerometer.
Spinocerebellar ataxia (SCA) and multiple system atrophy (MSA) are major neurodegenerative diseases mainly affecting the cerebellum and brainstem. SCA and MSA patients with predominant cerebellar ataxia (MSA-C) usually show ataxic gait as the initial and cardinal symptom [1, 2]. Clinical scales such as the Scale for the Assessment and Rating of Ataxia (SARA) or Unified Multiple System Atrophy Rating Scale are easy to administer and can be used to assess the severity of gait disturbance [3, 4]; however, these scales are mostly qualitative and ordinal. A triaxial accelerometer is a well-known device for measuring human movement . The measured values obtained using a triaxial accelerometer are objective, and test-retest reliability has been validated . Triaxial accelerometers have been applied to several movement disorders to evaluate the characteristics of motor symptoms and the effectiveness of treatment [7,8,9]. Human gait is a highly complicated and integrated bipedal locomotion activity that can be characterized by multivariate data; therefore, simpler and more comprehensive methods are needed for the interpretation of such data. Several of the gait parameters obtained using a triaxial accelerometer have been proven to be useful for the clinical assessment of ataxia in SCA and MSA-C patients [10, 11]; however, a method for the integrated analysis of various gait parameters directly reflecting ataxic gait has not yet been developed. Principal component analysis (PCA) is a statistical approach used to process complicated multivariate data. PCA can reasonably integrate multivariate data covering some principal components (PCs) and can provide continuous variables by calculating the principal component score (PCS), which reflects the aspect of the principal component. PCA has been used recently in patients with dementia and aged patients to reveal the main factors accounting for gait disturbance [12, 13].
The purpose of this study was to reveal the characteristics of ataxic gait and to develop a novel method for the quantitative assessment of ataxic gait in SCA and MSA-C patients. For this purpose, we applied PCA to analyze the multivariate gait parameters obtained using a triaxial accelerometer.
Patients whose SARA score of gait was greater than six points or whose SARA score of stance was greater than three points were excluded because they could not complete the walking task using a triaxial accelerometer. Patients who had comorbid conditions that affect motor function, such as cerebrovascular or orthopedic disorders, were also excluded. After the application of the exclusion criteria, 61 patients clinically diagnosed with SCA or MSA-C and 57 control subjects without gait impairment were enrolled in this study. All patients could stand and walk by themselves, but some used a cane or walker to avoid falling. Fifteen patients (2 with SCA type 6 [SCA6], 8 with SCA31, and 5 with MSA-C) and 18 control subjects were measured twice within an approximately 6-month interval to evaluate chronological changes. When the subjects were recruited, they were provided with all necessary information about the study and informed consent was obtained from all subjects. This study was approved by the Ethics Committee of Shinshu University School of Medicine (No. 2667).
Instrumentation and measurements
A triaxial accelerometer (Jukudai Mate; Kissei Comtec Co., Ltd., Matsumoto, Japan) was used. The device was small (size, 55 mm × 80 mm; thickness, 10 mm) and light (weight, ~90 g). It had a sampling rate of 20 Hz. The range of detection was between −10 G and +10 G (G: acceleration of gravity, 1 G = 9.80665 m/s2) and the resolution power was 0.02 G. The data acquired by the device were analyzed by BIMUTAS II (Kissei Comtec Co., Ltd., Matsumoto, Japan), which was developed for biological processing. The device was attached to the back (median of L3) of the subject by an elastic belt. The subjects were asked to walk on a flat floor at a speed they were comfortable with. The walking distance was 10 m; however, the subjects were instructed to stop walking at 3 m beyond the end of the walkway. The walking test was repeated 12 times (6-times reciprocating walk) consecutively. For subjects who felt it difficult to walk 12 times, the test was aborted at the end of the third reciprocating walk (i.e., they completed the 10-m walk 6 times).
The axes of the direction of acceleration were defined in three dimensions: anterior-posterior direction (AP), medio-lateral direction (ML), and vertical direction (VT). Gait velocity, cadence, step length, auto-correlation coefficient (AC) in each direction, and the root mean square (RMS) in each direction were taken as gait parameters. Velocity was calculated by using the time taken to complete each 10-m walking test, which was measured with a manual stopwatch. Cadence, defined as steps per minute, was calculated by dividing the steps counted from the acceleration data in the 10-m walking by walking time. Step length was calculated using velocity and cadence.
We used an auto-correlation method to calculate AC [14, 15]. To compare the similarity of certain two different acceleration data, a cross-correlation function was used. When the compared two acceleration data were identical, the function was specially called as auto-correlation function [14, 15]. Each value of the auto-correlation function is a correlation coefficient between the raw acceleration data and the data to be shifted at some sample point from the same raw acceleration data. When the shifted sample point was zero, the position was defined as the reference point. The first AC peak next to the reference point represented step regularity [14, 15], and the ratio of the first AC peak next to the reference point to the second peak represented step symmetry . RMS represented the degree of sway during gait . By eliminating the beginning and end of the test, we used the acceleration data from the middle 6.4 s to calculate AC and RMS with as many data points as possible for each subject. As a result, the same number (128) of data points was analyzed in all subjects. The mean value according to the number of measurements was used for gait analysis in all parameters.
Statistical analysis was performed using SPSS for Windows (IBM Statistics 22). The assumption of a normal distribution was assessed by the Kolmogorov-Smirnov test. Between-group differences in age and gender were assessed using an independent t-test (continuous data) and a chi-square test (categorical data), respectively. Differences in gait parameters between both groups were examined with an independent t-test. The correlation between disease duration and gait parameters was assessed by Spearman’s rank correlation coefficient. To validate the reliability of the gait parameters, an intra-class correlation coefficient (ICC) was used. Specifically, to focus on the reliability of a single time measurement, ICC (1, 1) was used.
In PCA, all gait parameters were included and standardized to zero mean and unit variance before PCA was performed. Varimax rotation was used to derive orthogonal factor loading. PCs with eigenvalues greater than one were considered to be relevant. Factor loading values greater than 0.4 as the absolute value were considered relevant . The first and second PCSs were calculated for each subject as the linear combination of the factor loading values in the patients and standardized gait parameters in each subject. The scores were standardized to zero mean and unit variance. For multiple comparisons of PCS divided according to the SARA score of gait, one way analysis of variance (ANOVA) was used. The Bonferroni correction was applied as a post hoc test. The chronological change of PCSs was assessed by a dependent t-test. The level of significance was set at p < 0.05 in all tests.
The characteristics of the subjects are shown in Table 1. There were no significant differences in gender and age distribution between the patients and controls (gender: χ 2 = 0.131, p = 0.71; age: t = −1.881, p = 0.08).
The values of each gait parameter are shown in Table 2. All gait parameters, except for step symmetry, were significantly different between the patients and controls. RMSs for AP and VT were significantly lower, whereas that for ML was significantly higher in the patients than in the controls. The ICC (1, 1) values were approximately or greater than 0.6 for all parameters, except step symmetry, in both groups.
Factor loading values and the proportion of the variance in the patients and controls are shown in Table 3. In the controls, four PCs (PC1, PC2, PC3, and PC4) were found to explain 79% of the variance. In PC1, velocity, step length, and RMS in all directions were relevant. In PC2, 3, and 4, cadence, step symmetry, and step regularity were relevant, respectively.
In the patients, PC1 was similar to that in the controls, except for the involvement of cadence. However, the factor loading values in PC2 were very different. In particular, RMS in ML was relevant and had a negative value.
The distributions of the first and second PCSs are shown in Fig. 1a. Both the first and second PCSs were significantly different between the patients and controls (t = 5.189, p < 0.001; t = 7.527, p < 0.001, respectively), and well correlated with disease duration and the SARA score of gait in the patients (ρ = −0.363, p = 0.004; ρ = −0.574, p < 0.001, respectively). The distributions of the second PCS in the patients and controls are shown in Fig. 1b. The second PCS was significantly different between the patients and controls. There was a significant difference between the subgroups (ANOVA: F (5, 53) = 16.866, p < 0.001; post hoc test: between SARA score of gait 0 and that of 5 and 6; between that of gait 1 and that of 5 and 6; between that of 2 and that of 5). The ICC (1, 1) values of the first PCS were 0.92 (95% confidence interval [CI]: 0.89–0.95) in the patients and 0.85 (95% CI: 0.75–0.89) in the controls; those of the second PCS were 0.63 (95% CI: 0.54–0.72) in the patients and 0.33 (95% CI: 0.24–0.44) in the controls.
Using the first and second PCSs, the chronological change with an interval of approximately 6 months was measured (Fig. 2). However, despite the fact that there were too few MSA-C patients to perform reliable statistical analysis, only the second PCS tended to decline in MSA-C patients (first PCS: t = 1.849, p = 0.138 and second PCS: t = 2.492, p = 0.067), whereas the changes of both the first and second PCSs were not apparent, not only in the controls (first PCS: t = 1.654, p = 0.116 and second PCS: t = 0.086, p = 0.933) but also in the SCA6/SCA31 patients (first PCS: t = 0.603, p = 0.562 and second PCS: t = −0.646, p = 0.535).
As PCA can yield orthogonal factors that account for the variance of data, it is used to divide multivariable data into several components. PCA has also been proven useful to resolve gait into its component elements [12, 13, 17], which were named Pace, Rhythm, Asymmetry, and Variability [12, 17]. Not all of the parameters in this study were the same as those used in previous reports [12, 17], but they were basically interchangeable considering the physiological meaning of each parameter. As for the results of PCA in the controls, it was not difficult to interpret the meanings of PC2, 3, and 4. Cadence had very high factor loading in PC2. In the same manner, only step symmetry was relevant in PC3, and step regularity was relevant in PC4. Thus, it was reasonable that PC2, 3, and 4 were named Rhythm, Asymmetry, and Variability, respectively. In PC1, velocity, step length, and RMSs in each direction were relevant. RMS is closely correlated with gait velocity in normal gait [5, 15]; therefore, it is reasonable that velocity and RMSs in each direction were contained in the same component. As a higher RMS value is caused by a faster gait velocity , gait velocity was more predominant in PC1 than RMS, supporting our contention that PC1 could be named Pace.
The results of PCA in the patients were considerably different from those of the controls. In the patients, the factor loading values of cadence, velocity, step length, and RMS were relevant to PC1. This means that PC1 was not able to differentiate the Pace factor from the Rhythm factor in the patients. In PC2, the difference between both groups was much more evident. To understand the meaning of PC2 in the patients more clearly, we calculated the second PCS for each subject. The second PCS was significantly different between the patients and controls, and was significantly correlated with disease duration and the SARA score of gait in the patients. The second PCS was also different between the subgroups divided according to the SARA score of gait. Hence, it is reasonable to suppose that PC2 in the patients represented the main component of ataxic gait.
RMS increases with the square of gait velocity in normal gait [5, 15]; thus, there should be a positive correlation between RMS and gait velocity in all directions. In fact, the factor loading values of PC1 in the controls showed a positive correlation between RMS and gait velocity. However, the factor loading values of PC2 in the patients demonstrated a negative correlation between RMS in ML and gait velocity. In addition, the factor loading values of PC2 indicated that RMS in ML was more relevant than that in the other directions when calculating the second PCS. As the second PCS had a significant negative correlation with the SARA score and disease duration, a smaller second PCS represented a higher severity of ataxic gait. From the results, PCA showed that a short step, low step regularity and symmetry, and a high degree of body sway in ML were characteristics of ataxic gait.
We considered that the second PCS would be effective for the assessment of chronological changes in the subjects. It is noteworthy that the range of the second PCS was very narrow, and there was no significant chronological change in the controls. Due to the small sample size, the chronological change of the second PCS could not be examined sufficiently in the patients. However, it was likely that the degree of the chronological change of the PCS would be different among the disease subtypes. In particular, the score change in the MSA-C patients, but not the SCA6/SCA31 patients, was relatively large.
A triaxial accelerometer is an objective and reliable tool to study standing and gait [6, 18]. Although the data obtained using a triaxial accelerometer are considered to be highly reliable, the reliability of new parameters calculated from the acceleration data needed to be validated. For this purpose, the intra-rater reliability of each gait parameter was assessed using ICC (1, 1). The ICC (1, 1) values of velocity, step length, and RMS in each direction were greater than 0.8 in the patients and controls. The ICC (1, 1) value of RMS in this study was approximately the same as that in a previous report . We considered that the reliability of the instrument we used was demonstrated by this result. The ICC (1, 1) value of cadence in the controls was also similar to that in a previous report , but it was relatively low in the patients. This difference is due to gait disturbance in the patients. As for step symmetry, there was no significant difference between both groups. In addition, there was no significant correlation between step symmetry and disease duration (AP direction: ρ = −0.147: p = 0.260 and VT direction: ρ = −0.229: p = 0.076). These results indicate that step asymmetry was not very conspicuous in the patients with SCA and MSA-C at the stage at which they could stand and walk unaided.
SARA is used widely to assess the severity of ataxia in clinical neurology settings. However, SARA is an ordinal scale and the change of the total SARA score was approximately 1–2 per year in several subtypes of SCA [20,21,22]. These are the main reasons why SARA does not have sufficient sensitivity for the assessment of ataxic gait in SCA or MSA-C patients. Disease-modifying therapies are still under investigation for SCA or MSA-C, but short-term intensive coordinated rehabilitation has been shown to be effective for ataxia patients [23, 24]. Considering clinical trials of upcoming disease-modifying drugs or of robotics-assisted rehabilitation, additional quantitative methods are needed to investigate short-term (~1 year) changes in the severity of ataxic gait. Our findings indicated that the second PCS is a continuous value reflecting the characteristics of ataxic gait in SCA and MSA-C patients. PCA has been used previously to reveal some of the characteristics of pathological gait not only in neurodegenerative disorders but also in patients with hip fracture or spinal cord injury [25, 26]. Simplified gait parameters such as walking speed or walking distance are not sufficiently sensitive to assess the severity of pathological gait in such conditions, and PCA provides new walking paradigms beyond a simple gait parameter . To the best of our knowledge, this study is the first to apply PCA to gait analysis in SCA or MSA-C patients, and we found that PCA clearly detected the characteristics of ataxic gait. We propose that the second PCS can help us assess ataxic gait more objectively and quantitatively than the available methods.
PCA of gait parameters revealed the main component of ataxic gait. This component was characterized as low gait velocity, short step, low step regularity and symmetry, and a high degree of body sway in the medio-lateral direction. The PCS of the main component was significantly different between the patients and controls, and it was well correlated with disease duration and the SARA score of gait in the patients. Although more detailed validation is needed in future studies regarding how applicable it is to assess the severity of ataxic gait, the PCS will provide a novel quantitative biomarker to assess the severity of ataxia.
Analysis of variance
Intra-class correlation coefficient
Multiple system atrophy
Multiple system atrophy with predominant cerebellar ataxia
Principal component analysis
Principal component score
Root mean square
Scale for the Assessment and Rating of Ataxia
Spinocerebellar ataxia type 31
Spinocerebellar ataxia type 6
Rossi M, Perez-Lloret S, Doldan L, Cerquetti D, Balej J, Millar Vernetti P, et al. Autosomal dominant cerebellar ataxias: a systematic review of clinical features. Eur J Neurol. 2014;21:607–15.
Lin DJ, Hermann KL, Schmahmann JD. Multiple system atrophy of the cerebellar type: clinical state of the art. Mov Disord. 2014;29:294–304.
Schmitz-Hübsch T, du Montcel ST, Baliko L, Berciano J, Boesch S, Depondt C, et al. Scale for the assessment and rating of ataxia: development of a new clinical scale. Neurology. 2006;66:1717–20.
Wenning GK, Tison F, Seppi K, Sampaio C, Diem A, Yekhlef F, et al. Development and validation of the Unified Multiple System Atrophy Rating Scale (UMSARS). Mov Disord. 2004;19:1391–402.
Moe-Nilssen R. A new method for evaluating motor control in gait under real-life environmental conditions. Part 1: the instrument. Clin Biomech (Bristol, Avon). 1998;13:320–7.
Lugade V, Fortune E, Morrow M, Kaufman K. Validity of using tri-axial accelerometers to measure human movement - Part I: posture and movement detection. Med Eng Phys. 2014;36:169–76.
Shull PB, Jirattigalachote W, Hunt MA, Cutkosky MR, Delp SL. Quantified self and human movement: a review on the clinical impact of wearable sensing and feedback for gait analysis and intervention. Gait Posture. 2014;40:11–9.
Godinho C, Domingos J, Cunha G, Santos AT, Fernandes RM, Abreu D, et al. A systematic review of the characteristics and validity of monitoring technologies to assess Parkinson's disease. J Neuroeng Rehabil. 2016;13:24.
Yang CC, Hsu YL. A review of accelerometry-based wearable motion detectors for physical activity monitoring. Sensors (Basel). 2010;10:7772–88.
Matsushima A, Yoshida K, Genno H, Murata A, Matsuzawa S, Nakamura K, et al. Clinical assessment of standing and gait in ataxic patients using a triaxial accelerometer. Cerebellum Ataxias. 2015;2:9.
Shirai S, Yabe I, Matsushima M, Ito YM, Yoneyama M, Sasaki H. Quantitative evaluation of gait ataxia by accelerometers. J Neurol Sci. 2015;358:253–8.
Verghese J, Wang C, Lipton RB, Holtzer R, Xue X. Quantitative gait dysfunction and risk of cognitive decline and dementia. J Neurol Neurosurg Psychiatry. 2007;78:929–35.
Lord S, Galna B, Verghese J, Coleman S, Burn D, Rochester L. Independent domains of gait in older adults and associated motor and nonmotor attributes: validation of a factor analysis approach. J Gerontol A Biol Sci Med Sci. 2013;68:820–7.
Yang CC, Hsu YL, Shih KS, Lu JM. Real-time gait cycle parameter recognition using a wearable accelerometry system. Sensors (Basel). 2011;11:7314–26.
Moe-Nilssen R, Helbostad JL. Estimation of gait cycle characteristics by trunk accelerometry. J Biomech. 2004;37:121–6.
Stevens JP. Applied multivariate statistics for the social sciences. New Jersey: Lawrence Erlbaum Associates; 2002.
Hollman JH, McDade EM, Petersen RC. Normative spatiotemporal gait parameters in older adults. Gait Posture. 2011;34:111–8.
Moe-Nilssen R. Test-retest reliability of trunk accelerometry during standing and walking. Arch Phys Med Rehabil. 1998;79:1377–85.
Fortune E, Lugade V, Morrow M, Kaufman K. Validity of using tri-axial accelerometers to measure human movement - Part II: Step counts at a wide range of gait velocities. Med Eng Phys. 2014;36:659–69.
Jacobi H, du Montcel ST, Bauer P, Giunti P, Cook A, Labrum R, et al. Long-term disease progression in spinocerebellar ataxia types 1, 2, 3, and 6: a longitudinal cohort study. Lancet Neurol. 2015;14:1101–8.
Ashizawa T, Figueroa KP, Perlman SL, Gomez CM, Wilmot GR, Schmahmann JD, et al. Clinical characteristics of patients with spinocerebellar ataxias 1, 2, 3 and 6 in the US; a prospective observational study. Orphanet J Rare Dis. 2013;8:177.
Yasui K, Yabe I, Yoshida K, Kanai K, Arai K, Ito M, et al. A 3-year cohort study of the natural history of spinocerebellar ataxia type 6 in Japan. Orphanet J Rare Dis. 2014;9:118.
Ilg W, Synofzik M, Brötz D, Burkard S, Giese MA, Schöls L. Intensive coordinative training improves motor performance in degenerative cerebellar disease. Neurology. 2009;73:1823–30.
Miyai I, Ito M, Hattori N, Mihara M, Hatakenaka M, Yagura H, et al. Cerebellar Ataxia Rehabilitation Trialists Collaboration. Cerebellar ataxia rehabilitation trial in degenerative cerebellar diseases. Neurorehabil Neural Repair. 2012;26:515–22.
Thingstad P, Egerton T, Ihlen EF, Taraldsen K, Moe-Nilssen R, Helbostad JL. Identification of gait domains and key gait variables following hip fracture. BMC Geriatr. 2015;15:150.
Awai L, Bolliger M, Ferguson AR, Courtine G, Curt A. Influence of spinal cord integrity on gait control in human spinal cord injury. Neurorehabil Neural Repair. 2016;30:562–72.
The authors thank Ms. Asuka Murata, Mrs. Setsuko Matsuzawa, and Mrs. Sonomi Nagasaki for their technical support. We are also grateful to all the subjects for their participation in this study.
This work was supported from Research Committee of the Ataxia, Research on Policy Planning and Evaluation for Rare and Intractable Diseases, Health and Labour Sciences Research Grants, The Ministry of Health, Labour, and Welfare, Japan (K. Yoshida). This work was also supported in part by a grant (H27-3-28) from the Nagano Society for The Promotion of Science, Nagano Prefecture, Japan (A. Matsushima), and by a Grant-in-Aid for Scientific Research (C) (15H03031) from the Japan Society for the Promotion of Science, the Ministry of Education, Culture, Sports, Science, and Technology of Japan (K. Yoshida).
Availability of data and materials
The data that support the findings of this study are available from the corresponding author upon reasonable request.
AM, KY, and HG designed the experiments. AM and KY recruited the subjects. AM, KY, and HG obtained the acceleration data. AM analyzed and interpreted the data. AM and KY wrote the manuscript. All authors provided critical feedback on the manuscript. This study was conducted under the supervision of SI. All authors read and approved the final manuscript.
This study was approved by the COI Committee of Shinshu University, and author HG participated in data collection, but not in data analysis. The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
When the subjects were recruited, they were provided with all necessary information about the study and informed consent was obtained from all subjects. This study was approved by the Ethics Committee of Shinshu University School of Medicine (No. 2667).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Matsushima, A., Yoshida, K., Genno, H. et al. Principal component analysis for ataxic gait using a triaxial accelerometer. J NeuroEngineering Rehabil 14, 37 (2017). https://doi.org/10.1186/s12984-017-0249-7