Skip to main content

Prediction of gait recovery in spinal cord injured individuals trained with robotic gait orthosis



Motor impairment is a major consequence of spinal cord injury (SCI). Earlier studies have shown that robotic gait orthosis (e.g., Lokomat) can improve an SCI individual’s walking capacity. However, little is known about the differential responses among different individuals with SCI. The present longitudinal study sought to characterize the distinct recovery patterns of gait impairment for SCI subjects receiving Lokomat training, and to identify significant predictors for these patterns.


Forty SCI subjects with spastic hypertonia at their ankles were randomly allocated to either control or intervention groups. Subjects in the intervention group participated in twelve 1-hour Lokomat trainings over one month, while control subjects received no interventions. Walking capacity was evaluated in terms of walking speed, functional mobility, and endurance four times, i.e. baseline, 1, 2, and 4 weeks after training, using the 10-Meter-Walking, Timed-Up-and-Go, and 6-Minute-Walking tests. Growth Mixture Modeling, an analytical framework for stratifying subjects based on longitudinal changes, was used to classify subjects, based on their gait impairment recovery patterns, and to identify the effects of Lokomat training on these improvements.


Two recovery classes (low and high walking capacity) were identified for each clinical evaluation from both the control and intervention groups. Subjects with initial high walking capacity (i.e. shorter Timed-Up-and-Go time, higher 10-Meter-Walking speed and longer 6-Minute-Walking distance) displayed significant improvements in speed and functional mobility (0.033 m/s/week and–0.41 s/week respectively); however no significant change in endurance was observed. Subjects with low walking capacity exhibited no significant improvement. The membership in these two classes—and thus prediction of the subject’s gait improvement trajectory over time—could be determined by the subject’s maximum voluntary torque at the ankle under both plantar-and dorsi-flexion contractions determined prior to any training.


Our findings demonstrate that subjects responded to Lokomat training non-uniformly, and should potentially be grouped based on their likely recovery patterns using objective criteria. Further, we found that the subject’s ankle torque can predict whether he/she would benefit most from Lokomat training prior to the therapy. These findings are clinically significant as they can help individualize therapeutic programs that maximize patient recovery while minimizing unnecessary efforts and costs.


A common focus during rehabilitation after spinal cord injury (SCI) is on promoting improvements in functional walking capacity. Fast and effective enhancement of gait function can enhance a patient's independence, life satisfaction, and subsequent reintegration into society as a fully-participating member. Body weight support treadmill training (BWSTT) has been considered a primary rehabilitative strategy for improving locomotion in patients with incomplete SCI [1, 2]. BWSTT provides intensive, task-specific repetitive walking training, promotes supraspinal neuroplasticity involved in locomotion, and consequently helps patients to regain their walking skills [35]. While conventionally BWSTT is supplemented by manual assistance from therapists, robotic-assisted BWSTT has become widely used over the past decade to provide symmetric walking training to support lower extremity rehabilitation in patients with neurological movement impairment. The most widely-used automated locomotor training system is the Lokomat® (Hocoma AG, Volketswil, Switzerland) [6, 7]. Recent studies have shown that patients can receive positive physical and psychological benefits from Lokomat training [1, 2, 817], such as improved walking capacity [1, 810], improved metabolic performance [2, 12], and increased activity in the cerebellum [15]. However, a recent review from Swinnen [18] has shown that the effectiveness of Lokomat training remains controversial, with some studies showing significant improvements in motor and walking performance after Lokomat trainings, while others showing minimal changes [18]. The authors and previous researchers have pointed out that difficulty arises when translating these research findings to the clinical setting because of the lack of a control group in those studies [19]. The authors also observe that the advantages of Lokomat training over other types of locomotor trainings have yet to be fully demonstrated.

In addition, the high variability in recovery pattern among patients should be explicitly considered in the analysis. Locomotor training, like any other intervention, is not expected to affect all patients equally; instead, different treatment responses are expected for different patient subgroups. Many research studies, however, oversimplify by assuming that all individuals are drawn from a single population with a similar recovery pattern; pre-training and post-training measures are then compared. This approach loses the heterogeneous response of patients to the therapy program as well as the longitudinal information describing the timeline of improvement—information that should be considered when evaluating the efficacy of a particular intervention. In fact, many studies involved with locomotor training have reported results with a standard deviation even larger than the sample mean, indicating that the high inter-individual variance within their sample population results in a less meaningful average value. Thus, it is necessary to address these deficits by first identifying distinct latent classes within each sample and then modeling the recovery patterns for each class separately. In addition, pre-vs. post-training comparisons usually have a lowered effect size and statistical power than longitudinal tracking studies, because the error term has a larger standard deviation than in other statistical models [20].

Recent statistical techniques such as Growth Mixture Modeling (GMM) allow researchers to effectively recognize this high inter-individual variance, and to model the growth pattern of a longitudinal recovery procress. The inclusion of additional time points between the pre-and post-training evaluations can help elucidate the training effect on gait impairment, and further can dramatically improve the statistical power in situations where increasing sample size is not feasible due to costs and patient availability [20].

In the present study, rather than performing a pre-vs. post-training comparison within one experimental treatment group, we examined in detail the growth trajectory of three aspects of gait impairment (speed, functional mobility, and endurance, as evaluated by typical clinical measures) for subjects in an intervention group receiving Lokomat training and those in a control group who received no intervention. Instead of assuming that all subjects are drawn from a single population, we objectively classified the subjects into distinct subgroups based on their growth trajectory of walking capacity, and then modeled the effect of Lokomat training for each subgroup. Finally, we evaluated whether baseline measurements of neuromuscular and clinical performance could serve as statistical predictors for these recovery patterns. In particular, we hypothesized that ankle strength—measured in terms of maximum voluntary dorsiflexion and plantarflexion torque—could predict the recovery trend for these clinical evaluations.


Study participants

Forty incomplete SCI subjects (twenty-seven males, thirteen females) with spastic hypertonia at their lower extremities participated in this single-center, unblinded, and randomized study (Table 1). Each subject was assigned either to the control or intervention (Lokomat training) group according to a permuted block randomization design. All subjects were injured within their cervical or upper thoracic (superior to T10) vertebrae.

Table 1 Baseline characteristics of SCI patients (n = 40) in the Lokomat and control group (mean ± standard deviation), including age, post-injury time, ankle joint tested, gender ratio, walking index for spinal cord injury (WISCI II), plantarflexion MVC torque (Tp) and dorsiflexion MVC torque (Td)

All subjects gave informed consent according to the policies of the Institutional Review Board of Northwestern University.

Lokomat training

The Lokomat was used to train subjects in pre-programmed physiological gait pattern. The system uses a pair of robotic orthoses to move the subject’s legs over a synchronized treadmill while a harness system dynamically adjusts the amount of body-weight support. Gait parameters, such as the range of motion of the hip and knee joints, walking speed, body weight support, and guidance force, were adjusted in real-time by the physical therapist.

Each subject received a one-hour training session three times per week for four consecutive weeks; as it took 15–20 minutes to set up the subject, the gait training lasted up to 45 minutes per session. The goal of the Lokomat training was to increase the walking speed and decrease guidance force and body-weight support over the twelve training sessions based on the subject's improvement. In each session, the training speed was increased gradually from the minimum allowed by the machine to a peak magnitude (up to 3.4 km/hr) where the subject could maintain gait quality without feeling fatigued. The guidance force was then decreased from full assistance of the robot, until the subject was no longer able to maintain the set speed (low to 20% assistance) and the body-weight support, customized to the individual subject’s ability. The same therapist provided all trainings for a given subject.

Outcome measures: walking capacity

Walking capacity was assessed during the training regimen using three standard clinical evaluations: 10-Meter-Walking Test (10MWT) [21], Timed-Up-and-Go (TUG) [22], and 6-Minute-Walking Test (6MWT) [23]. Subjects were evaluated four times for their overground walking capacity: at the baseline (i.e., prior to any Lokomat training), and after 1, 2, and 4 weeks of training. Subjects were permitted to use assistive devices and braces during the tests. Harness systems, parallel bars, and other support systems were not permitted. All tests for an individual subject were measured by the same examiner.

The 10MWT was used to evaluate walking speed over a short distance by measuring the time spent walking over a distance of 10 meters. The 10MWT has been rated as the test that provides the most valid measure of walking speed improvement [24].

The TUG was used to evaluate the functional recovery by measuring time taken by the subject to stand up from an armed chair, walk three meters, turn around, and sit back in the chair. It has been used to evaluate functional mobility and balance for various movement disorders [2527].

The 6MWT test was used to assess walking endurance by measuring the distance walked in six minutes. It has been widely used to evaluate the response to therapeutic treatment for stroke and spinal cord injury, and is believed to provide a more accurate measure of daily activity than other walking tests [28].

These timed walking tests have been confirmed for their high concurrent and discriminant validity for SCI individuals [23, 29], and are widely used as guidelines to both evaluate functional walking capacity and aid in the interpretation of clinical research results.

Muscle strength: walking capacity predictor

The isometric torque resulting from the maximum voluntary contraction (MVC) at the ankle joint was measured after each subject finished all required clinical evaluations. Prior to any testing, the Modified Ashworth Score (MAS) was evaluated at each ankle joint. The ankle joint with the highest modified Ashworth scale and lowest maximum voluntary contraction (both conditions often occurred on the same side) was chosen for the MVC evaluation. Subjects were then seated in an adjustable chair with the ankle attached to a customized shoe mounted securely to a rigid plate. The trunk and thigh were strapped to the chair to minimize their contributions to ankle and foot movement. Subjects were then instructed to generate (MVC) in the dorsiflexion (Td) and plantarflexion (Tp) directions, using only their ankle joint, while the ankle joint was held fixed at 90° and the knee joint was held at 120°. During the experiment, the subject exerted a maximum force by pressing down (plantarflexion) or pulling up (dorsiflexion) against the stationary foot plate, which was in turn attached to a six-axis torque transducer (JR3 Inc., Woodland, CA). The subjects were instructed to practice exerting the MVC torque without using their knee and hip joints by minimizing the EMG activity in quadriceps and hamstring muscles. Electromyography (EMG) activity from the tibialis anterior, gastrocnemius, quadriceps and hamstring muscles was monitored using bi-polar surface electrodes (Ag-Ag/Cl type, Delsys, Boston MA); trials with substantial quadriceps and hamstring activity were repeated. Two trials were performed for each evaluation; the trial with the higher peak torque was chosen for the subsequent data analysis.

Data analysis: recovery pattern identification and prediction

Growth Mixture Modeling (GMM) was used to subdivide participants into multiple latent classes based on the recovery patterns (i.e., the change over time) of their walking measures, and subsequently inspect gait improvement within each class. Logistical regression was then applied to identify potential predictors for subject class membership.

Recovery pattern identification using GMM

GMM is one type of longitudinal mixture model, and has been widely used to analyze data that exhibits heterogeneity in developmental pathways [3037]. GMM assumes that a given sample of subjects can be divided into a finite number of groups, each of which is comprised of subjects with similar recovery patterns. That is, all subjects within a given class were assumed to exhibit similar changes over time. In effect, GMM is analogous to a clustering technique, e.g., k-means, but with consideration of multiple observations over time. The intervention effect was then be modeled within individual subgroups and compared with control or other treatment groups. A two-step GMM procedure was employed in this study. Step 1 classified the control and intervention groups separately, based on the measured growth trajectories of their clinical evaluations. A mixed-effects model (random coefficient regression) was then used to determine if the trajectories showed significant change over time. Step 2 involved a joint analysis to examine the training effect on walking capacity. For Class 1, the difference in the rates of change between the control and intervention groups was compared using the regression slopes; the same procedure was then used for Class 2 [37]. The quality of the resulting classification was evaluated by the posterior probabilities of subject class membership and by the Bayesian information criteria.

Prediction of recovery patterns at baseline using logistic regression analysis

The class membership of subjects in each treatment group was correlated with the latent variables (potential predictors) using logistic regression [3840], in order to determine whether injury level, age, gender, walking index for spinal cord injury (WISCI II) and muscle strength (MVC torque at the ankle joint) evaluated at the baseline could predict clinical assessments of functional improvement over the subsequent four weeks. Such an analysis can help to discover potentially-causal relationships between the predictors and functional improvement.

The group allocation, i.e., control and Lokomat training, was concealed (blinded) to the statistician. The GMM analysis was performed in the R statistical package V 2.12 (open source software,, while the logistic regression was performed using SAS v 9.0 (SAS Institute, Cary NC). Statistical tests with a P < 0.05 were considered significant.


Recovery patterns of gait impairment measures

GMM was used to identify the growth trajectories for both control and intervention groups, and to identify the effect of training on these classes. Figures 1 and 2 show the mean trajectory for each class. Table 2 summarizes the statistical results from the GMM analysis.

Figure 1
figure 1

10MWT speed as a function of training week. Mean ± SE for each latent class shown. The linear regression line for Class 2 in the Lokomat group was obtained from GMM. No other changes in 10MWT speed with time were significant.

Figure 2
figure 2

TUG time as a function of training week. Mean ± SE in each latent class shown. The linear regression line for Class 2 in Lokomat was obtained from GMM. Note that the ordinate scales are different between Class 1 and Class 2, due to their large difference in standard errors.

Table 2 Illustration of classification and joint group analysis from GMM for control and intervention groups

Latent class identification

Control group

The growth trajectory for each of the three clinical measures was examined by GMM over the 4-week training period. Two distinct growth patterns were identified for each walking measure (Figures 1 and 2, Table 2). The latent classes that included subjects with longer 10MWT and TUG times and shorter 6MWT distances at baseline (i.e., the initial test) were defined as the low walking capacity classes (Class 1). Similarly, the latent classes for subjects with shorter 10MWT and TUG times and a longer 6MWT distance were defined as the high walking capacity classes (Class 2).

Intervention group

Similar to the control group, two classes of growth pattern were identified for each clinical evaluation, corresponding to low and high walking capacity at baseline (Figures 1 and 2, Table 2).

Subject class membership agreed across the different clinical evaluations for both groups. The average posterior probabilities of class membership were 0.99 for each walking test, suggesting low classification error.

Training effect identification

For the 10MWT test, the GMM analysis considering joint groups found a significant beneficial training effect for subjects with a high walking capacity (P = 0.02). In the Lokomat group, subjects in Class 2 increased their walking speed at a rate of 0.033 m/s/week (P = 0.01). Subjects in Class 1 did not show significant benefit from the Lokomat training (P = NS).

For TUG, the GMM analysis of joint groups found a beneficial training effect on subjects with high walking capacity (P = 0.05). Subjects in Class 2 of the Lokomat group decreased their time at a rate of 0.41 s/week (P = 0.02). Subjects in Class 1 did not show significant benefit from the Lokomat training (P = NS).

For 6MWT, the training effect was not found to be significant for subjects from either the high or low walking capacity groups, based on the GMM analysis (P = NS).

The change over time for each clinical measure was not significant for any latent class in the control group (P = NS).

Prediction of gait impairment recovery patterns by mvc torque

Figure 3 shows the relationship between baseline MVC torque and subject class membership for the 10MWT evaluation in the Lokomat group. Table 3 presents the results of the logistic regression to determine whether MVC torque is able to predict subject class membership for the Lokomat group.

Figure 3
figure 3

Relationship between plantarflexion (Tp, panel A) torque and dorsiflexion (Td, panel B) MVC torque and subject class membership, for the 10MWT evaluation in the Lokomat group. Note: in each boxplot, the bounds of the box denote the first and third quartiles, while the two whiskers denote the minimum and maximum values.

Table 3 Logistic regression results used to assess whether Tp or Td could predict patient class membership in the Lokomat group

As the plantarflexion MVC torque (Tp) had overlapping ranges for the two latent classes for each clinical evaluation, logistic regression was used to inspect the effects of Tp on the subject class membership, using Class 1 as the reference. For the 10MWT evaluation, subjects with Tp > 18.3 Nm were more likely to be members of Class 2 than subjects with Tp ≤ 18.3 Nm. For each unit increase in Tp, the odds that a subject belonged to Class 2 increased by 18% (P < 0.01). For the TUG test, it was shown that subjects with Tp > 17.5 Nm had a higher likelihood of membership in Class 2 (P < 0.01). For 6MWT, subjects with Tp > 16.1 Nm were more likely to belong to Class 2 than subjects with Tp ≤ 16.1 (P < 0.01).

A similar logistic regression was performed for the dorsiflexion MVC torque (Td). It was shown that subjects with Td > 9.5 Nm had a higher likelihood of Class 2 membership than subjects with Td ≤ 9.5 Nm for both the 10MWT and TUG evaluations (P < 0.01). For the 6MWT, subjects with Td > 8.9 Nm were more likely to belong to Class 2 than subjects with Td ≤ 8.9 Nm (P < 0.01).

Figure 4 shows the predicted probability of a subject to belong to the high walking capacity class (Class 2), based on their measured Tp or Td, assuming pooled subjects of both the control and intervention groups. It further confirmed that Tp and Td were significant predictors for walking-capacity classification (P < 0.01).

Figure 4
figure 4

Logistic regression results to determine ability of Tp (panel A) and Td (panel B) to predict class membership. Regression indicates whether MVC torque predicts whether a subject is in the high walking capacity class (i.e., Class 2) for pooled subjects from both groups. The predicted probability curve (solid line) and 95% confidence limits (dashed line) are shown.

Other possible factors, including injury level, Modified Ashworth Scale score, gender, age and post-injury time were found not to have any significant effect on the class membership for any clinical evaluation (P = NS for all).


The longitudinal growth pattern of gait impairment for two groups of SCI subjects—control and intervention—over four weeks of ongoing Lokomat training, was investigated and analyzed with the GMM. Two classes of growth trajectory in each clinical evaluation were identified for both the control and intervention groups, representing subgroups with low and high walking capacity. The two latent classes responded to the Lokomat training differently; subjects with high walking capacity presented significant and consistent improvements in walking speed and balance with time, while subjects with low walking capacity exhibited no such significant improvement. Interestingly, the results indicated that the baseline (i.e. pre-training) measures of MVC torque (Td and Tp) could predict the differential treatment response, i.e., subjects with high Tp and Td were more likely to have both high walking capacity and receive significant benefit from Lokomat training. These findings can help clinicians to identify the subgroups of individuals with incomplete SCI in terms of their current and potential walking capacity and to estimate the differential treatment effects of the Lokomat intervention on their recovery progress.

The effect of lokomat training on walking improvement

The twelve sessions of Lokomat training for SCI subjects improved their functional ambulation in a selective manner. Subjects with low walking capacity did not show significant improvements. By contrast, subjects with a high walking capacity at baseline presented a consistent linear trend in time for both speed and functional balance over the 4-week training period. Our average improvement in 10MWT speed (0.132 m/s) indicates that the 12-session Lokomat training program can improve the walking speed for SCI subjects with high walking capacity; this improvement exceeds the threshold established by Lam [41] (viz. 0.13 m/s) for significant clinical change in 10MWT. Our results agree with previous reports that Lokomat training improves the overground walking capacity [1, 9, 15, 42], and that subjects with less control impairment achieve greater improvement [19].

While the 10MWT, TUG and 6MWT evaluate different aspects of gait impairment, they have been shown to be highly-correlated in the walking validity test [23, 29, 41, 43, 44]. Despite this, surprisingly we observed no significant change in the 6MWT, whereas the 10MWT and TUG improved significantly. Such a discrepancy was also reported in a similar training for the stroke population [45]. By contrast, other studies reported improvements in endurance for SCI subjects who were provided with two months or more of Lokomat training as well as supplementary pharmaceutical intervention [1, 10]. Thus, lack of significant improvement in 6MWT in our study could be due to the lower number of training sessions in our study than others. In particular, it is likely that the 4-weeks of training was insufficient for the individuals in the low-capacity class, as 4-weeks may be too little to effect the necessary neuroplasticity to improve walking function. Further, our walking endurance test was performed at a self-selected speed that the subject felt was most comfortable. Therefore, trial-to-trial variations during the clinical evaluations might be much higher than the improvements received from the training, i.e., an individual subject did not necessarily perform the clinical evaluations at his/her maximum effort level, complicating the detection of the training effects.

Need for classification of recovery patterns

A few other studies in the literature have classified their participants prior to other analyses, based on a priori or empirical considerations. A walking validity test by van Hedel et. al. grouped subjects into low and high capacity based on their walking index score prior to performing the correlation analysis [23]. Correlation between walking tests and WISCI II score was found to be significant for subgroup with high walking function (WISCI II > 10), while it was little for subgroup with WISCI II ≤ 10. A study by Field-Fote et. al. also classified subjects into two subgroups by their walking speeds [9]. They found that subjects with slower initial walking speed achieved higher improvement in walking speed than those with faster initial walking speeds, in contrast with our finding. The discrepancy could be due to the fact that they used a limited number of subjects who were randomly assigned to 4 different body-weight-support assisted-stepping groups, including treadmill training with manual assistance, robotic assistance, stimulation, and overground training with stimulation. Each slower- and faster- initial walking group included a combination of subjects who received different types of training with different mechanisms of actions, making the interpretation of the effects of Lokomat alone difficult.

By contrast, the GMM analysis in this study classified subjects with the same training in terms of their growth trajectories using a rigorous statistical technique, which provides objective evidence for the validity of classifying subjects into subgroups prior to other statistical analyses.

Further, it can be expected that the recovery manner of SCI subjects during neurorehabilitation might present more complex growth patterns than those found in the current study, e.g., different numbers of latent classes for different clinical tests, and inconsistencies in the class membership between the various clinical evaluations. Classification with objective criteria, such as that provided by GMM, can provide a robust and repeatable analytical framework to reveal distinct subgroups which exhibit substantial differences between subgroups but maintain high within-group similarities.

Prediction of effect of lokomat training on improving walking capacity

To determine the extent to which the recovery pattern of walking capacity can be predicted based on baseline measures, the ability of MVC torque (Td and Tp) to determine subject classification was assessed. Our results showed that subjects with a larger Tp or Td had a higher likelihood of membership in the high walking capacity group for each clinical test, and thus to receive significant training benefit. Such findings indicate that this quantitative measurement of ankle MVC dorsiflexion and plantarflexion torque can potentially be a fast and reliable clinical assessment that can be used to identify which subjects are most likely to make significant progress during a one-month Lokomat training regimen prior to the start of training. The ability of other kinetic and kinematic variables to predict gait recovery, such as knee and hip strength, as well as range of motion of lower limb joints, will be investigated in future studies.

Although clinical assessments of gait impairment, such as the 10MWT, TUG, and 6MWT assessments, have been well-accepted by physicians, their applications are also limited by the requirement of necessary time and space to perform such tests, inter-subject differences due to the choice of route for the 6MWT assessment, and the motivation for subjects to perform the tests [28]. By contrast, our study implies that ankle MVC provides an objective measurement that can be part of a routine office visit with a physician.

Our study also can lend insight into the role of muscle strength on gait impairment. In particular, we expect that weakness of the anterior tibialis (TA) and the gastrocnemius (GS) are likely related to gait impairment. This is expected due to the role played by the neuromuscular properties of the spastic ankle during ambulation. For example, walking capacity after SCI is commonly impaired by equinus foot, partially due to spasticity and weakness of the ankle dorsiflexors. Spasticity of ankle plantarflexors has reciprocal inhibitory effects on TA activity, and additionally causes spastic hyperactivity of GS and hypoactivity of TA. Both of these can reduce the voluntary contraction capacity of these muscles and contribute to equinus foot and walking impairment. It has been shown that, in addition to improving walking capacity, Lokomat training can significantly reduce the neuromuscular abnormalities associated with spasticity and may result in increased ankle MVC torque [46]. Such an explanation is further supported by the finding that the improvements in muscle strength for SCI subjects with spastic hypertonia at the ankle joint are mainly due to reductions in reflex stiffness [47, 48]. Our results extend those findings by showing that MVC torque measured at baseline can be a significant predictor for the recovery pattern of SCI gait impairment.

Study limitations

It was not feasible to blind the researchers who performed the tests to the group allocation since they needed to stay with the subjects during the Lokomat trainings and experimental protocols. All tests and trainings for each individual subject were performed by the same researcher; however, since this longitudinal study occurred over a span of multiple years, it was not feasible for one researcher to test all subjects. The effects of these limitations on the results are unclear.

Directions for future research

Although several studies have compared the therapeutic effects of Lokomat to other physical interventions, it is our interest to use the advanced techniques that we used in this study to compare the therapeutic effects of this relatively new physical intervention to more well-established pharmacological interventions that are developed to improve function and reduce clinical impairments. Determining whether the Lokomat training is equally or more effective than pharmacological treatments is clinically significant, as some patients are able to tolerate physical interventions but not medications, while for others the opposite case is true. Thus, it is important to quantify outcomes for different types of interventions in order to provide optimal treatments. Furthermore, since the mechanisms and actions of these interventions are different, applying a combination of these interventions may provide benefits superior to the individual treatments.


Overall, our study indicates that Lokomat training can be an effective physical intervention for SCI subjects with high walking capacity, but not for those with low walking capacity. Further, our findings indicate that the MVC torque measured prior to training may be an objective and reliable tool to predict which subject groups can benefit the most from Lokomat training, prior to being assigned to that therapy. Such an approach can allow clinicians to prescribe a therapeutic and rehabilitation plan that best optimizes an individual’s recovery outcomes while minimizing unnecessary effort and costs.


  1. Wirz M, Zemon DH, Rupp R, Scheel A, Colombo G, Dietz V, Hornby TG: Effectiveness of automated locomotor training in patients with chronic incomplete spinal cord injury: a multicenter trial. Arch Phys Med Rehabil 2005, 86: 672-680. 10.1016/j.apmr.2004.08.004

    Article  PubMed  Google Scholar 

  2. Hidler J, Hamm LF, Lichy A, Groah SL: Automating activity-based interventions: the role of robotics. J Rehabil Res Dev 2008, 45: 337-344. 10.1682/JRRD.2007.01.0020

    Article  PubMed  Google Scholar 

  3. Barbeau H, Wainberg M, Finch L: Description and application of a system for locomotor rehabilitation. Med Biol Eng Comput 1987, 25: 341-344. 10.1007/BF02447435

    Article  CAS  PubMed  Google Scholar 

  4. Behrman AL, Harkema SJ: Locomotor training after human spinal cord injury: a series of case studies. Phys Ther 2000, 80: 688-700.

    CAS  PubMed  Google Scholar 

  5. Harkema SJ: Plasticity of interneuronal networks of the functionally isolated human spinal cord. Brain Res Rev 2008, 57: 255-264. 10.1016/j.brainresrev.2007.07.012

    Article  PubMed  PubMed Central  Google Scholar 

  6. Hidler J, Wisman W, Neckel N: Kinematic trajectories while walking within the Lokomat robotic gait-orthosis. Clin Biomech 2008, 23: 1251-1259. 10.1016/j.clinbiomech.2008.08.004

    Article  Google Scholar 

  7. Wirz M, Bastiaenen C, de Bie R, Dietz V: Effectiveness of automated locomotor training in patients with acute incomplete spinal cord injury: a randomized controlled multicenter trial. BMC Neurol 2011, 11: 60. 10.1186/1471-2377-11-60

    Article  PubMed  PubMed Central  Google Scholar 

  8. Lam T, Wirz M, Lunenburger L, Dietz V: Swing phase resistance enhances flexor muscle activity during treadmill locomotion in incomplete spinal cord injury. Neurorehabil Neural Repair 2008, 22: 438-446. 10.1177/1545968308315595

    Article  PubMed  Google Scholar 

  9. Field-Fote EC, Lindley SD, Sherman AL: Locomotor training approaches for individuals with spinal cord injury: a preliminary report of walking-related outcomes. J Neurol Phys Ther 2005, 29: 127-137. 10.1097/01.NPT.0000282245.31158.09

    Article  PubMed  Google Scholar 

  10. Hornby TG, Zemon DH, Campbell D: Robotic-assisted, body-weight-supported treadmill training in individuals following motor incomplete spinal cord injury. Phys Ther 2005, 85: 52-66.

    PubMed  Google Scholar 

  11. Mirbagheri MM, Patel C, Quiney K: Robotic-assisted locomotor training impact on neuromuscular properties and muscle strength in spinal cord Injury. Conf Proc IEEE Eng Med Biol Soc 2011, 2011: 4132-4135.

    CAS  PubMed  Google Scholar 

  12. Israel JF, Campbell DD, Kahn JH, Hornby TG: Metabolic costs and muscle activity patterns during robotic- and therapist-assisted treadmill walking in individuals with incomplete spinal cord injury. Phys Ther 2006, 86: 1466-1478. 10.2522/ptj.20050266

    Article  PubMed  Google Scholar 

  13. Hidler JM, Wall AE: Alterations in muscle activation patterns during robotic-assisted walking. Clin Biomech (Bristol, Avon) 2005, 20: 184-193. 10.1016/j.clinbiomech.2004.09.016

    Article  Google Scholar 

  14. Lunenburger L, Bolliger M, Czell D, Muller R, Dietz V: Modulation of locomotor activity in complete spinal cord injury. Exp Brain Res 2006, 174: 638-646. 10.1007/s00221-006-0509-4

    Article  CAS  PubMed  Google Scholar 

  15. Winchester P, McColl R, Querry R, Foreman N, Mosby J, Tansey K, Williamson J: Changes in supraspinal activation patterns following robotic locomotor therapy in motor-incomplete spinal cord injury. Neurorehabil Neural Repair 2005, 19: 313-324. 10.1177/1545968305281515

    Article  PubMed  Google Scholar 

  16. Koenig A, Omlin X, Zimmerli L, Sapa M, Krewer C, Bolliger M, Muller F, Riener R: Psychological state estimation from physiological recordings during robot-assisted gait rehabilitation. J Rehabil Res Dev 2011, 48: 367-385. 10.1682/JRRD.2010.03.0044

    Article  PubMed  Google Scholar 

  17. Nooijen CF, Ter Hoeve N, Field-Fote EC: Gait quality is improved by locomotor training in individuals with SCI regardless of training approach. J Neuroeng Rehabil 2009, 6: 36. 10.1186/1743-0003-6-36

    Article  PubMed  PubMed Central  Google Scholar 

  18. Swinnen E, Duerinck S, Baeyens J-P, Meeusen R, Kerckhofs E: Effectiveness of robot-assisted gait training in persons with spinal cord injury: a systematic review. J Rehabil Med 2010, 42: 520-526. 10.2340/16501977-0538

    Article  PubMed  Google Scholar 

  19. Dobkin BH, Duncan PW: Should body weight-supported treadmill training and robotic-assistive steppers for locomotor training trot back to the starting gate? Neurorehabil Neural Repair 2012, 26: 308-317. 10.1177/1545968312439687

    Article  PubMed  PubMed Central  Google Scholar 

  20. Venter A, Maxwell SE, Bolig E: Power in randomized group comparisons: the value of adding a single intermediate time point to a traditional pretest-posttest design. Psychol Methods 2002, 7: 194-209.

    Article  PubMed  Google Scholar 

  21. Ditunno JF Jr, Ditunno PL, Graziani V, Scivoletto G, Bernardi M, Castellano V, Marchetti M, Barbeau H, Frankel HL, D'Andrea Greve JM, Ko HY, Marshall R, Nance P: Walking index for spinal cord injury (WISCI): an international multicenter validity and reliability study. Spinal Cord 2000, 38: 234-243. 10.1038/

    Article  PubMed  Google Scholar 

  22. Bohannon RW: Reference values for the timed up and go test: a descriptive meta-analysis. J Geriatr Phys Ther 2006, 29: 64-68. 10.1519/00139143-200608000-00004

    Article  PubMed  Google Scholar 

  23. van Hedel HJ, Wirz M, Dietz V: Assessing walking ability in subjects with spinal cord injury: validity and reliability of 3 walking tests. Arch Phys Med Rehabil 2005, 86: 190-196. 10.1016/j.apmr.2004.02.010

    Article  PubMed  Google Scholar 

  24. Jackson AB, Carnel CT, Ditunno JF, Read MS, Boninger ML, Schmeler MR, Williams SR, Donovan WH: Outcome measures for gait and ambulation in the spinal cord injury population. J Spinal Cord Med 2008, 31: 487-499.

    PubMed  PubMed Central  Google Scholar 

  25. Andersson AG, Kamwendo K, Seiger A, Appelros P: How to identify potential fallers in a stroke unit: validity indexes of 4 test methods. J Rehabil Med 2006, 38: 186-191. 10.1080/16501970500478023

    Article  PubMed  Google Scholar 

  26. Zampieri C, Salarian A, Carlson-Kuhta P, Aminian K, Nutt JG, Horak FB: The instrumented timed up and go test: potential outcome measure for disease modifying therapies in Parkinson's disease. J Neurol Neurosurg Psychiatry 2010, 81: 171-176. 10.1136/jnnp.2009.173740

    Article  PubMed  PubMed Central  Google Scholar 

  27. Podsiadlo D, Richardson S: The timed "Up & Go": a test of basic functional mobility for frail elderly persons. J Am Geriatr Soc 1991, 39: 142-148.

    Article  CAS  PubMed  Google Scholar 

  28. ATS statement: Guidelines for the six-minute walk test. Am J Respir Crit Care Med 2002, 166: 111-117.

    Article  Google Scholar 

  29. Scivoletto G, Tamburella F, Laurenza L, Foti C, Ditunno JF, Molinari M: Validity and reliability of the 10-m walk test and the 6-min walk test in spinal cord injury patients. Spinal Cord 2011, 49: 736-740. 10.1038/sc.2010.180

    Article  CAS  PubMed  Google Scholar 

  30. Farina A, Sekizawa A, Ralston SJ, D'Alton ME, Bianchi DW: Latent class analysis applied to patterns of fetal sonographic abnormalities: definition of phenotypes associated with aneuploidy. Prenat Diagn 1999, 19: 840-845. 10.1002/(SICI)1097-0223(199909)19:9<840::AID-PD654>3.0.CO;2-E

    Article  CAS  PubMed  Google Scholar 

  31. Muthen B, Muthen LK: Integrating person-centered and variable-centered analysis: growth mixtue modeling with latent trajectory classes. Alcohol Clin Exp Res 2000, 24: 882-891. 10.1111/j.1530-0277.2000.tb02070.x

    Article  CAS  PubMed  Google Scholar 

  32. Hagenaars JA, McCutcheon AL: Applied latent class analysis. In Book Applied latent class analysis. City: Kluwer; 2002:476.

    Chapter  Google Scholar 

  33. Muthen B: Latent variable analysis: growth mixture modeling and related techniques for longitudinal data. In Handbook of Quantitative Methodology for the Social Sciences. Edited by: Kaplan D. Newbury Park: Sage Publications; 2004:345-368.

    Google Scholar 

  34. Proust-Lima C, Leternneur L, Jacqmin-Gadda H: A nonlinear latent class model for joint analysis of multivariate longitudinal data and a binary outcome. Stat Med 2006, 26: 2229-2245.

    Article  Google Scholar 

  35. Jung T, Wickrama KAS: An introduction to latent class growth analysis and growth mixture modeling. Soc Personal Psychol Compass 2008, 2: 302-317. 10.1111/j.1751-9004.2007.00054.x

    Article  Google Scholar 

  36. Hardy SA, Thiels C: Using latent growth curve modeling in clinical treatment research: an example comparing guided self-change and cognitive behavioral therapy treatments for bulimia nervosa. Int J Clin Health Psychol 2009, 9: 51-71.

    Google Scholar 

  37. Meghani SH, Lee CS, Hanlon AL, Bruner DW: Latent class cluster analysis to understand heterogeneity in prostate cancer treatment utilities. BMC Med Inform Decis Mak 2009, 9: 47. 10.1186/1472-6947-9-47

    Article  PubMed  PubMed Central  Google Scholar 

  38. Krishnapuram B, Carin L, Figueiredo MA, Hartemink AJ: Sparse multinomial logistic regression: fast algorithms and generalization bounds. IEEE Trans Pattern Anal Mach Intell 2005, 27: 957-968.

    Article  PubMed  Google Scholar 

  39. Littell RC: SAS for Mixed Models. 2nd edition. Cary, N.C.: SAS Institute, Inc.; 2006.

    Google Scholar 

  40. Tabachnick BG, Fidell LS: Using Multivariate Statistics. 5th edition. Boston: Pearson/Allyn & Bacon; 2007.

    Google Scholar 

  41. Lam T, Noonan VK, Eng JJ: A systematic review of functional ambulation outcome measures in spinal cord injury. Spinal Cord 2008, 46: 246-254. 10.1038/

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Mehrholz J, Kugler J, Pohl M: Locomotor training for walking after spinal cord injury. Cochrane Database Syst Rev 2008, 2: CD006676.

    PubMed  Google Scholar 

  43. van Hedel HJ, Wirz M, Dietz V: Standardized assessment of walking capacity after spinal cord injury: the European network approach. Neurol Res 2008, 30: 61-73. 10.1179/016164107X230775

    Article  CAS  PubMed  Google Scholar 

  44. Furlan JC, Noonan V, Singh A, Fehlings MG: Assessment of impairment in patients with acute traumatic spinal cord injury: a systematic review of the literature. J Neurotrauma 2011, 28: 1445-1477. 10.1089/neu.2009.1152

    Article  PubMed  PubMed Central  Google Scholar 

  45. Westlake KP, Patten C: Pilot study of Lokomat versus manual-assisted treadmill training for locomotor recovery post-stroke. J Neuroengineering Rehabilitation 2009, 6: 18. 10.1186/1743-0003-6-18

    Article  Google Scholar 

  46. Mirbagheri M, Niu X, Varoqui D: Prediction of Gait Recovery as a Tool to Rationalize Locomotor Training in Spinal Cord Injury. Rome, Italy: Fourth IEEE RAS/EMBS International Conference on Biomedical Robotics and Biomechatronics; 2012.

    Book  Google Scholar 

  47. Mirbagheri MM, Ladouceur M, Barbeau H, Kearney RE: The effects of long-term FES-assisted walking on intrinsic and reflex dynamic stiffness in spastic spinal-cord-injured subjects. IEEE Trans Neural Syst Rehabil Eng 2002, 10: 280-289.

    Article  PubMed  Google Scholar 

  48. Barbeau H, Ladouceur M, Mirbagheri MM, Kearney RE: The effect of locomotor training combined with functional electrical stimulation in chronic spinal cord injured subjects: walking and reflex studies. Brain Res Rev 2002, 40: 274-291. 10.1016/S0165-0173(02)00210-2

    Article  PubMed  Google Scholar 

Download references


This research was supported by the National Institutes of Health [R01HD059895] and the Craig H. Neilsen Foundation awarded to MMM.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Xun Niu or Mehdi M Mirbagheri.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XN and DV conducted the experimental procedure and conducted most of the data processing and statistical analysis. XN primarily drafted the manuscript, with substantial editorial assistance from MMM and MK. MMM conceived and supervised the project, and oversaw all experimental and analytical efforts. All authors read and approved the final manuscript.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Reprints and permissions

About this article

Cite this article

Niu, X., Varoqui, D., Kindig, M. et al. Prediction of gait recovery in spinal cord injured individuals trained with robotic gait orthosis. J NeuroEngineering Rehabil 11, 42 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: