Skip to main content

Patient-cooperative control increases active participation of individuals with SCI during robot-aided gait training



Manual body weight supported treadmill training and robot-aided treadmill training are frequently used techniques for the gait rehabilitation of individuals after stroke and spinal cord injury. Current evidence suggests that robot-aided gait training may be improved by making robotic behavior more patient-cooperative. In this study, we have investigated the immediate effects of patient-cooperative versus non-cooperative robot-aided gait training on individuals with incomplete spinal cord injury (iSCI).


Eleven patients with iSCI participated in a single training session with the gait rehabilitation robot Lokomat. The patients were exposed to four different training modes in random order: During both non-cooperative position control and compliant impedance control, fixed timing of movements was provided. During two variants of the patient-cooperative path control approach, free timing of movements was enabled and the robot provided only spatial guidance. The two variants of the path control approach differed in the amount of additional support, which was either individually adjusted or exaggerated. Joint angles and torques of the robot as well as muscle activity and heart rate of the patients were recorded. Kinematic variability, interaction torques, heart rate and muscle activity were compared between the different conditions.


Patients showed more spatial and temporal kinematic variability, reduced interaction torques, a higher increase of heart rate and more muscle activity in the patient-cooperative path control mode with individually adjusted support than in the non-cooperative position control mode. In the compliant impedance control mode, spatial kinematic variability was increased and interaction torques were reduced, but temporal kinematic variability, heart rate and muscle activity were not significantly higher than in the position control mode.


Patient-cooperative robot-aided gait training with free timing of movements made individuals with iSCI participate more actively and with larger kinematic variability than non-cooperative, position-controlled robot-aided gait training.


Body weight supported treadmill training (BWSTT) has become a widely used rehabilitation technique for individuals with walking disabilities due to neurological disorders such as stroke and spinal cord injury [14].

Robotic devices have been developed to relieve physical therapists from the straineous and unergonomical burden of manual BWSTT [5]. The Lokomat (Hocoma AG, Switzerland) [6], the ReoAmbulator (Motorika, USA), and the Gait Trainer (Reha-Stim, Germany) are used in clinical practice to automate BWSTT by moving patients repetitively along pre-defined walking trajectories.

A growing body of studies shows that both manual BWSTT and robot-aided treadmill training improve gait quality [715]. While some of these studies found advantages of robot-aided treadmill training compared to BWSTT [9, 11, 14], others found conventional treadmill training to be more effective [12, 13].

The studies in favor of robot-aided treadmill training focused more closely on non-ambulatory patients, while the studies reporting better outcome of conventional treadmill training included mainly ambulatory patients. These results suggest that currently, robot-aided treadmill training is most effective for severely affected, non-ambulatory patients, whereas it may not be ideal for more advanced, ambulatory patients. In contrast to these ambulatory patients, who may benefit more from other approaches like over-ground training, patients in the transition phase between being non-ambulatory and ambulatory still require much physical support during training. This situation demonstrates the need to improve current rehabilitation robots in a way that extends their spectrum of effective treatment to functionally more advanced patients. Such an improvement would allow patients to benefit from robot-aided treadmill training up to a point where they can safely and efficiently perform over-ground training. Thus, rehabilitation robots would be able to optimally support patients in their progression through their different stages of recovery.

In most of the studies mentioned above, the rehabilitation robots were controlled in a very simple way. A pre-recorded gait pattern was replayed by the robot as accurately as possible. This position control approach allows the patient to remain passive during the training [16] and reduces kinematic variability to a minimum [17]. However, both active participation and kinematic variability are considered as important promotors of motor learning and rehabilitation [1823]. fMRI studies comparing training tasks with active and passive movements have shown stronger cortical activation and subsequently also more cortical reorganization leading to more effective formation of motor memory when subjects where contributing actively to the trained movements compared to being passively moved [18, 19]. In a review of robotic therapy approaches based on these findings, Dromerick et al. conclude that these approaches are effective, but rigorous comparisons with traditional techniques still need to be performed [20].

Bernstein emphasized the crucial role of kinematic variability during motor learning ("repetition without repetition") based on practical experience and theoretical considerations [21]. Lewek et al. have shown that kinematic variability as introduced by conventional treadmill training improved the coordination of intralimb kinematics in ambulatory stroke patients while position-controlled robot-aided treadmill training with little kinematic variability did not [22]. Huang and Krakauer argue that from a computational motor-learning perspective, robots should ensure the successful completion of movements, allowing the adapting human nervous system to identify combinations of sensory states and their transitions associated with the motor commands required for the movements [23].

Therefore, researchers in the field of rehabilitation robotics believe that robotic control approaches, which increase active participation of the patients and allow more kinematic variability while still guaranteeing successful task execution, have the potential to substantially boost the efficacy of robot-aided rehabilitation, especially in functionally more advanced patients. Numerous research groups have been working on these patient-cooperative control strategies [2434]. While there have been extensive tests of control strategies that increase patient participation during training for upper-extremity robots [35, 36], most of the approaches for lower extremity-robots have only been evaluated in single case studies with patients or in proof-of-concept experiments with healthy volunteers.

In a recent publication, our group has demonstrated a patient-cooperative control strategy ("Path Control") for the Lokomat which allows free timing of leg movements while ensuring that the spatial kinematics of the legs stay within definable desired limits [37]. We could show that healthy volunteers participated more actively and with more--especially temporal--variability than in a classical, position controlled training mode. Moreover, we were able to modulate the level of activity by an additional supportive "flow" that did not reduce the amount of movement variability when providing more support. We assume that the ability to modulate the level of required activity will be an important feature to adapt the controller to the individual capabilities of patients, particularly of patients transitioning from a non-ambulatory to an ambulatory state during their rehabilitation process. Finally, we evaluated the feasibility of the path control strategy with 15 individuals with chronic incomplete spinal cord injury (iSCI). Assuming a minimal level of voluntary motor control, the patients were able to train with the patient-cooperatively controlled Lokomat.

In the present paper, we have investigated if the short-term effects found for healthy volunteers do also translate to spinal cord injured patients. More specifically, we have posed the following research questions: (1) Does patient-cooperative robot-aided treadmill training lead to more active participation of individuals with iSCI than classical, position-controlled training? (2) Can we deliberately modulate the activity required by the iSCI patient during the training? (3) Can we increase the variability of the iSCI patients' leg movements while still maintaining functional gait?


Gait training robot

Experiments were performed with the gait rehabilitation robot Lokomat. The robot automates body weight-supported treadmill training of patients with locomotor dysfunctions in the lower extremities such as spinal cord injury and hemiplegia after stroke [38]. It comprises two actuated leg orthoses that are attached to the patients' legs. Each orthosis has one linear drive in the hip joint and one in the knee joint to induce flexion and extension movements of hip and knee in the sagittal plane. Knee and hip joint torques can be determined from force sensors between actuators and orthosis. Passive foot lifters can be added to induce ankle dorsiflexion during swing phase. A body weight support system with a harness attached to the patients' trunk reduces the effective body weight by a definable amount.

Control algorithms

Position control

The first approach implemented for the Lokomat was position control [38]. In this approach, the control algorithm tries to match the pre-defined reference trajectory qref(t) as closely as possible1.

Impedance control

A first step towards patient-cooperative behavior of the robot was the implementation of an impedance control algorithm [26]. The actual joint positions qact are virtually coupled to the reference positions qref(t) by a simulated spring and damper system with spring stiffness K and damping constant B. If Δq denotes the control deviation,

Δ q = q ref ( t ) q act ,

the desired joint torques τc for the robot drives are

τ c = K Δ q + B d d t Δ q = ( K hip 0 0 K knee ) Δ q + ( B hip 0 0 B knee ) d d t Δ q .

By adjusting the parameters of the virtual impedance, the therapist can make the training more or less demanding for the patient. With a very low virtual stiffness, the patient has to participate more actively to maintain a functional gait pattern. In practice, only K is adjusted by therapist, and B is adapted automatically as a function of K[26]. The classical position control mode is included as a special case with K set to the maximally achievable stiffness (Fig. 1, left side).

Figure 1
figure 1

Control algorithms. Control algorithms. Impedance control (with its special case position control) is illustrated on the left side. Path control is illustrated on the right side: (1) control action to bring the patient's leg back to the inside of the virtual tunnel, (2) "flow" of supportive torques, (3) "moving window" around time-dependent reference.

Path control

A prominent feature of the position and impedance control approaches is the direct coupling of temporal and spatial guidance. The path control strategy [37] and related approaches [35, 39, 40] overcome this limitation by providing a virtual tunnel. Within this tunnel, patients can move their legs with their own desired timing of movements. The boundaries of the virtual tunnel provide spatial guidance to make sure that the movements still follow a physiologically meaningful pattern in space.

Details about the implementation of the path control strategy for the Lokomat are given in [37]. In the context of the impedance control algorithm described above, the time-dependent reference qref(t) is replaced by the nearest neighbor qNN(qact) on the spatial pattern template. The modified control deviation Δ q ƒ is then the difference between qNN and qact, reduced by a dead zone around the tunnel center (Fig. 1, right side, (1)). The spring stiffness rendering the tunnel wall is linearly scaled from zero at the tunnel border to a maximum of Khip = 720 Nm/ rad, Kknee = 540 Nm/rad.

For the supporting "flow", a torque vector is calculated by differentiating the reference trajectory qref with respect to the relative position in the gait cycle S. Thus, the direction of the torque vector is tangential to the movement path in joint space (Fig. 1, right side, (2)).

τ ƒ s ( S ) = d d S q ref ( S ) d d S q ref ( S )

The actual supportive torques are

τ s ( S ) = τ ƒ s ( S ) ( 1 d c ) k s , d c [ 0 , 1 ]

where ks is a scalar factor that determines the amount of support in Nm, and dc is the relative distance of the current position qact to the center of the path. The relative distance dc is normalized to the width of the tunnel and saturated to the upper limit 1 for positions qact outside the tunnel. Thus, supportive torques are only provided within the tunnel.

Finally, a "moving window" can limit free timing to a definable range wwindow around the timed reference qref(t) as it is used by the impedance controller. qNN is then constrained to be maximally a definable percentage of the gait cycle ahead or behind the timed reference qref(t) (Fig. 1, right side, (3)).

Experimental design

Fifteen patients with chronic iSCI (Table 1) participated in a test training session to evaluate if they were able to train successfully with patient-cooperative controllers. Two out of these 15 patients were not able to train with the path control strategy because they had very weak control over their extensor muscles. Hence, they were not able to induce sufficient knee extension at the end of swing phase to move along the desired path. Two other patients dropped out because of personal reasons. The 11 remaining patients volunteered to participate in further experiments.

Table 1 Patient characteristics

All experimental procedures were approved by the Ethics Committee of the Canton of Zurich, Switzerland, and all participants provided informed consent before the experiments.

The 11 chronic iSCI patients trained with the Lokomat at a walking speed of 2 km/h (0.55 m/s) and with 30-50% body weight support under four different conditions:

  1. 1.

    POS: Position control with the stiffness of the Lokomat controller set to K hip = 1200 Nm/rad, K knee = 900 Nm/rad 2.

  2. 2.

    SOFT: Impedance control with the stiffness set to K hip = 192 Nm/rad, K knee = 144 Nm/rad

  3. 3.

    COOP: Path control with w window set to 20% of the gait cycle and the support gain k s adjusted individually for each patient3

  4. 4.

    COOP+: Path control with w window set to 20% of the gait cycle and the support gain k s increased to 130% of the value used in the previous condition

Prior to the experiment, surface EMG electrodes were attached to the patients' gastrocnemius medialis (GM), tibialis anterior (TA), vastus medialis (VM), rectus femorisi (RF), and biceps femoris (BF) muscles of the left leg. The electrodes were placed according to the SENIAM guidelines [42]. Custom-built foot-switches were taped under the heel of the left foot of the patients to determine heel strikes.

Two additional surface electrodes were placed over the electrical dipole axis of the heart, one below the right clavicle and one below the left pectoral muscle to record a simplified ECG for heart rate extraction. Before each condition, the patients were quietly standing in the Lokomat for 60 seconds. During the final 30 seconds of this period, ECG was recorded to determine the heart rate prior to each condition. After these resting period, patients walked for two minutes to get used to the respective controller. Afterwards, data was recorded during one minute of walking. In addition to the EMG and ECG signals, joint angles from the left hip and knee joints were recorded by sensors at the joint axes of the Lokomat.

Data analysis

Spatiotemporal variability

To quantify the amount of temporal and spatial variations in the gait patterns during walking in the different conditions, we computed the spatio-temporal characteristics of the recorded trajectories qact(t) according to the procedure described by Ilg et al. [43].

The recorded joint angles of each condition were cut into single strides triggered by the heel strike signal of the foot switches. The single strides were normalized in time to the interval [0, 1), with S denoting the normalized stride time. The trajectory of the kth normalized stride is referred to as q(k)(S), and the number of recorded strides is denoted N. The average trajectory qavg(S) was determined as a reference for the spatio-temporal analysis:

q avg ( S ) = 1 N k = 1 N q ( k ) ( S )

Each trajectory q(k)was mapped to the reference trajectory qavg by a spatial shift function ξ(k)(S) and a time shift function τ shift ( k ) ( S ) .

q ( k ) ( S ) = q avg ( S + τ shift ( k ) ( S ) ) + ξ ( k ) ( S )

The values of the shift functions ξ(k)(S) and τ shift ( k ) ( S ) were determined by optimization as described in [44].

The weighting factor for the optimization was determined according to the rules suggested in [43].

Finally, the spatial variability var ξ and the temporal variability var τ as defined in [43] were computed using the following equations:

var ξ = 1 N k = 1 N ( 0 1 | ξ ( k ) ( S ) | d S )
var τ = 1 N k = 1 N ( 0 1 | τ shift ( k ) ( S ) | d S )

The resulting spatial and temporal variability were compared by a Friedman test (nonparametric equivalent to a repeated measures ANOVA) at the 5% significance level [45]. Multiple comparisons were accounted for by the Bonferroni adjustment.

Interaction torques

To better understand the interactions between robot and patient, the interaction torques in the joints of the robot have been calculated. The robot's force sensors are located between drives and exoskeleton and not directly at the interaction points with the human, such that a model of the exoskeleton's dynamics has to be used to derive the interaction torques τint from the torques τmot, which are measured at the robot's drives:

τ int = τ mot M exo ( q act ) q ¦ act + n exo ( q act , q ˙ act )

with Mexo being the mass matrix capturing the inertia of the Lokomat exoskeleton and nexo subsuming the gravitational, friction, and Coriolis torques of the exoskeleton. Static friction in the joints has been identified in a separate experiment to be below 0.5 Nm and has thus been neglected in the dynamic model. To allow comparisons of the interaction torques under the different conditions, we computed the root mean square over whole recording time Trec:

τ ª int = 1 T rec 0 T rec ( τ int ( t ) ) 2 d t .

The root mean square values under the different conditions were compared by a Friedman test (nonparametric equivalent to a repeated measures ANOVA) at the 5% significance level with Bonferroni adjustment.

Heart rate

Heart rate was extracted from the simplified ECG recordings by custom Matlab code which determined the length of the RR intervals IRR. The reciprocal of the median of all RR intervals during the 30 seconds prior to each condition constitutes the pre-condition heart rate

HR pre = l median pre ( I RR ) .

Analogously, the heart rate during a condition HRduring was defined as the reciprocal of the median of all RR intervals during the last 30 seconds of each condition. The absolute heart rate increase ΔHR for each condition was then defined as

Δ HR = HR during HR pre .

We defined the maximal heart rate increase ΔHRmax for a specific patient as the maximum of the values for ΔHR under the four different training conditions. Finally, we normalized the absolute heart rate increase for the different conditions with respect to ΔHRmax to account for the variable cardiovascular reactions of the different patients. The normalization results in the relative heart rate increase Δ HRrel

Δ HR rel = Δ HR Δ HR max .

The values for ΔHRrel under the different conditions were compared by a Friedman test (nonparametric equivalent to a repeated measures ANOVA) at the 5% significance level with Tukey-Kramer adjustment.

Muscle activity

EMG signals were band-pass filtered between 15 and 300 Hz, rectified, and cut into single strides triggered by the heel strike signal of the foot switches. The single strides were normalized in time to 1001 samples each. All strides of a patient under a given condition were then averaged. Next, the average strides were broken up into seven phases (initial loading, mid stance, terminal stance, pre-swing, initial swing, mid swing, terminal swing) according to Perry [46]. The root mean square (RMS) of the EMG signals was calculated for each muscle within each of these phases.

The RMS values of the EMG signals showed high inter-subject variability, and the repeated measurements for a single subject were not independent of each other. Linear mixed models [47] are a statistical tool that can account for such circumstances. In these models, random variables can capture the covariance of multiple data values originating from different individual sources. The remaining subject-independent effects can be described as the linear influence of fixed factors.

To investigate the influence of the different conditions on muscle activity, we fitted a separate linear mixed model to the logarithm of the RMS values of the EMG signals of each muscle. For a given muscle, we define the logarithmized RMS for an observation j in a subject i as EMG ij . An observation is a combination of one of the four conditions and one of the seven gait phases. Hence, there were 7 × 4 = 28 observations j (j = 1, 2, ..., 28) per subject. We included the factors "condition" and "gait phase" as fixed effects. Thus, the value of EMG ij for a given observation j on the i-th subject was modeled as

EMG i j = β 0 + β 1 × COND 1 i j + β 2 × COND 2 i j + β 3 × COND 3 i j + β 4 × PHASE 1 i j + β 5 × PHASE 2 i j + ... + β 9 × PHASE 6 i j + u 0 i + ε i j .

The indicator variables COND1 ij to PHASE6 ij were set to one, if the observation j belonged to the respective condition or gait phase, otherwise to zero. To account for the correlation of repeated measurements within a subject i, a random intercept u 0i was assumed for each subject. The residual ε ij captures the difference between the measured value EMG ij and the prediction of the model.

In order to compare the different conditions, we computed the estimated marginal means for each condition by averaging the model predictions across the different gait phases. These estimated marginal means were then compared with post-hoc tests at the 5% significance level. In these tests, multiple comparisons were accounted for by the Bonferroni adjustment. A similar statistical analysis of EMG data has been performed in [37] and in [30].


Kinematics and spatiotemporal variability

Patients changed their gait kinematics notably under the different training conditions (Fig. 2). The virtual tunnel in the path control modes allowed for a less extended knee at initial contact, and consequently, patients reduced their peak knee extension. Patients also increased their maximal hip flexion during swing phase in the path control modes.

Figure 2
figure 2

Kinematic data. Resulting kinematic data. Trajectories in joint space for one exemplary patient (P12) under the different conditions POS (a), SOFT (b), COOP+ (c), COOP (d).

Spatial variability under conditions SOFT (soft impedance control mode), COOP (path control mode), and COOP+ (path control mode with increased supportive flow) was significantly higher than under condition POS (stiff position control mode). There were no significant differences between the conditions SOFT, COOP+, and COOP (Fig. 3, left).

Figure 3
figure 3

Spatiotemporal variability. Spatial variabilty var ξ (°) and temporal variability var τ (% gait cycle).

Temporal variability under the conditions COOP+ and COOP was significantly higher than under condition POS. Condition SOFT was not significantly different from any other condition (Fig. 3, right).

Interaction torques

Interaction torques in the hip joint between patient and robot were significantly smaller under conditions COOP+ and COOP than under condition POS. No significant differences between the conditions could be found for the interaction torques in the knee joint (Fig. 4).

Figure 4
figure 4

Interaction torques. Interaction torques τint (Nm) for hip and knee joint.

Heart rate

The relative heart rate increase ΔHRrel was significantly larger under condition COOP than under condition POS. No other significant differences could be identified (Fig. 5).

Figure 5
figure 5

Relative change of heart rate. Relative change of heart rate ΔHRrel while walking under the different conditions.

Muscle activity

Activity of the Tibialis anterior muscle was significantly increased under the COOP+ and COOP conditions compared to the POS and the SOFT conditions. No significant differences could be found for the activity of the Gastrocnemius medialis muscle. Activity of the Rectus femoris muscle was significantly increased under the COOP+ and COOP conditions compared to the POS condition. For the Vastus medialis muscle, conditions SOFT, COOP+, and COOP caused significantly higher activity than POS. Activity of the Biceps femoris muscle was significantly higher under the COOP condition than under the POS condition (Fig. 6).

Figure 6
figure 6

Muscle activity. Muscle activity of TA (Tibialis anterior), GM (Gastrocnemius medialis), VM (Vastus medialis), RF (Rectus femoris), and BF (Biceps femoris) muscles as predicted by the linear mixed models (left column). Comparison of mean muscle activity under the different conditions (right column).


Active participation

Basic neuroscience studies have shown that motor learning is more effective when human subjects practice movements actively rather than being passively moved [18, 19, 48, 49]. Although the underlying mechanisms are not well understood yet, this principle is generally translated also to robotic neurorehabilitation [23], where researchers aim at making patients participate as actively as possible during training.

Our evaluation has shown that iSCI patients participated with higher muscle activity (Fig. 6) and higher cardiovascular effort (increased heart rate, Fig. 5) when they were training under the path control condition (COOP) than under the position control condition (POS). Theoretically, this increased activity could also be caused by the robot generating torques opposed to the movements of the patient. While there are studies investigating the effects of such robotic resistance [50], our goal was to obtain active, unobstructed participation of the patients. The fact that interaction torques did not increase under the path control conditions (Fig. 4) shows that the patients were indeed contributing actively to the movements and not working against robotic resistance.

We have included a condition with soft impedance control (SOFT) as a benchmark for the current state-of-the-art of patient-cooperative Lokomat training in clinical practice. The impedance setting (Khip = 192 Nm/rad, Kknee = 144 Nm/rad, these values correspond to a "guidance force" setting of 40% in the commercial Lokomat software) for this condition was chosen based on discussions with the physical therapy staff at University Hospital Balgrist (Zurich, Switzerland) about the lowest impedance settings they use during clinical trainings on a regular basis. Interestingly, it appears that the remaining temporal guidance (Fig. 3, right) in this compliant control mode still kept the patients in a rather passive state: Only the vastus medialis muscle was significantly more active in compliant control mode than in position control mode. All other parameters did not differ significantly (Fig. 5, Fig. 6). This observation is in line with theoretical models of human-robot interactions which predict that the human motor system will "slack" whenever possible to reduce its effort [5154]. Apparently, the free timing of movements provided by the path control strategy which requires patients to actively propel their legs through the gait pattern makes patients less likely to "slack" than the timing-based soft impedance control mode used under condition SOFT.

Thus, the iSCI patients in our experiment participated more actively during training only with the patient-cooperative path control strategy.

Modulation of activity by additional support

Unlike in our study with healthy volunteers [37], we were not able to modulate activity by adjusting the amount of additional support. Apparently, subjects reacted very inconsistently to the increased support in condition COOP+. While for some subjects the additional support was actually helpful, others felt "pushed forward" and had to put more effort in actively canceling this "perturbation". This effect may be the reason for the large variability of heart rate increase under the condition COOP+ (Fig. 5).

As already seen in the feasibility experiment with iSCI subjects in [37], iSCI patients have diverse needs for support, usually limited to specific gait phases. Therefore, the "global" support parameter ks which determines the intensity of the supportive "flow" for the whole gait cycle appears to be not sufficient to adapt the support for iSCI patients. For an impedance controller based on a reference pattern with fixed timing, gait-phase dependent adaptation of controller impedance has been demonstrated by Emken et al. [33]. For the path control strategy evaluated in this paper, which allows free timing of movements, an automatic adaptation algorithm that identifies the individual deficits of a patient as implemented for the upper extremity by Wolbrecht et al. [55] could possibly improve the training mode by providing support that is better tailored to the individual patients.

Movement variability

Variability and the possibility to make errors is considered an essential component of practice for motor learning. Bernstein's demand that training should be "repetition without repetition" [21] is still considered a crucial requirement, which is also supported by recent advances in computational models describing motor learning [23]. More specifically, a recent study by Lewek et al. [22] has shown that intralimb coordination after stroke was improved by manual training after stroke, which allowed kinematic variability, but not by position-controlled Lokomat training, which reduced kinematic variability to a minimum.

The analysis of spatiotemporal variability shows that while spatial variability is significantly increased in all three compliant modes SOFT, COOP+, and COOP compared to the stiff position control condition POS, temporal variability is only significantly increased in the path control modes COOP+ and COOP.

The virtual tunnel of the path control strategy allowed spatial variability to an extent that still ensured a functional gait pattern, therefore, it did not substantially increase the patients' risk of stumbling.

Thus, the path control strategy does not only technically provide free timing of movements, but iSCI patients also showed more temporal variability in their movements than with position control (POS) or with the compliant, but timing-controlled impedance control (SOFT).


Limitations of the path control strategy

It should be noted that a constant treadmill speed was used throughout the presented experiments. Thus, the temporal freedom of the path control mode were limited to the swing phase. Nevertheless, a substantial increase in temporal variability could be detected. To increase patient interactivity during training, we will combine the path control strategy with approaches which adapt the treadmill speed according to the intention of patients [56].

The fixed walking pattern that defines the spatial movement path may not be ideal for every patient. As in position-controlled Lokomat training, the pattern can be adapted manually by the therapist. However, it is not guaranteed that a pattern close to the "healthy" pattern of an individual patient can be achieved. For hemiparetic patients, it would be possible to derive a desired path for the affected leg from observing the unaffected leg, as proposed by Vallery et al. [32]. For iSCI patients, an adaptive re-shaping of the path, similar to the approach by Jezernik et al. [25], may improve the applicability of the path control strategy.

Limitations of the study

The present study only investigated the reactions of iSCI patients to different controllers during a single training session with short exposure to the different training modes. Clearly, the long term effects of the different training modes are much more important and should be investigated in future work. However, we believe that verifying the intended, presumably beneficial effects in a single training session was an important first step in preparation of a long term trial.

We deliberately included patients with a wide range of ambulatory skills to gain insights into the feasibility of path control training with patients at different skill levels. The distribution of walking skills comprised four fully ambulatory patients with a WISCI score of 20, indicating that they were able to independently ambulate 10 m without any walking aids. Furthermore, six patients had reduced, but good ambulatory skills (WISCI score between 12 and 19) and were able to independently ambulate 10 m using appropriate walking aids (crutches and braces). Finally, there was one patient in the transition range between non-ambulatory and ambulatory, indicated by a WISCI score of 5. As we expect the most practical benefits of patient-cooperative control strategies for patients in the transition range between non-ambulatory and ambulatory, more data regarding the feasibility with functionally more restricted patients would be desirable. Thus, future studies with the path control strategy should more explicitly focus on patients within this functional range.

As we planned to include patients with very different walking skills, we decided that it would have been very difficult to reliably standardize a control condition where patients would have walked without assistance or manual assistance of a therapist. Therefore, we performed our experiments without such a condition which would of course have allowed very interesting further analyses. Future studies which will be focusing on patients from a more narrow functional range. As these patients will have similar--and thus standardizable--needs for support during manual assisted treadmill training, it will then be feasible to include such a condition.

The limited number of patients included in the study does not provide sufficient statistical power to stratify patients according to their disability levels, which might reduce the variability in the results and provide further insights into the different effects of the evaluated control strategies on different groups of patients. The focus of the study on iSCI patients leaves it an open question whether similar results can be expected for patients with stroke or other pathologies. The feasibility of patient-cooperative training and the immediate effects for such patients needs to be investigated separately.

The choice of heart rate as a measure of effort was made because it did put a relatively low additional burden on the patients during the experiment. As discussed by Pennycott et al., heart rate may be influenced by emotional state, pain and hydration level, whereas oxygen uptake would be a more robust measure of effort during robot-aided gait training [57]. These factors may explain the large variability under the condition COOP+ where some patients may have been irritated by the increased amount of robotic support. However, as the general trend of the heart rate results is consistent with the results regarding the muscle activity of the patients, we believe that the method has captured the patients' effort in a sufficiently robust way for the sake of our research questions.


Patients with incomplete spinal cord injury participated more actively and with larger kinematic variability in patient-cooperative robot-aided gait training than in non-cooperative, position-controlled robot-aided gait training. Free timing of movements appears to be an important feature of patient-cooperativeness, as a compliant impedance control mode with fixed timing did not significantly increase active participation, but the path control strategy with free timing did.

Future development should focus on providing adaptive, patient-specific support to make training with patient-cooperative control strategies feasible for a larger population of patients. Future clinical evaluation should compare the effects of patient-cooperative robot-aided training versus non-cooperative robot-aided training and manual BWSTT in a long term randomized clinical trial.

Foot Notes

1The following notation is used throughout this paper: all vectors of joint angles and torques consist of two elements, the first one for the hip joint and the second one for the knee joint, e.g. q = (q(1), q(2))T = (qhip, qknee)T. The control algorithms discussed in this paper are always defined for a single leg. The second leg is controlled by an independent second instance of the respective control algorithm.

2The equivalent end-point stiffness of the exoskeleton depends on the joint angles and the direction of force application and, thus, can not be reflected in a single, representative number. The relationship between end-point stiffness and joint stiffness in a lower-limb exoskeleton is discussed in [41].

3The therapist was instructed to adjust ks to the minimal value that enabled the patient to walk in the path control mode. The individual support gains which were used under this condition are listed in Tab. 1.


  1. Visintin M, Barbeau H, Korner-Bitensky N, Mayo NE: A new approach to retrain gait in stroke patients through body weight support and treadmill stimulation. Stroke; a journal of cerebral circulation 1998,29(6):1122-1128.

    Article  CAS  PubMed  Google Scholar 

  2. Barbeau H, Visintin M: Optimal outcomes obtained with body-Weight support combined with treadmill training in stroke subjects. Archives of Physical Medicine and Rehabilitation 2003,84(10):1458-1465. 10.1016/S0003-9993(03)00361-7

    Article  PubMed  Google Scholar 

  3. Sullivan KJ, Brown DA, Klassen T, Mulroy S, Ge T, Azen SP, Winstein CJ: Effects of Task-Specific Locomotor and Strength Training in Adults Who Were Ambulatory After Stroke: Results of the STEPS Randomized Clinical Trial. Phys Therapy 2007,87(12):1580-1602. 10.2522/ptj.20060310

    Article  Google Scholar 

  4. Dobkin B, Barbeau H, Deforge D, Ditunno J, Elashoff R, Apple D, Basso M, Behrman A, Harkema S, Saulino M, Scott M: The evolution of walking-related outcomes over the first 12 weeks of rehabilitation for incomplete traumatic spinal cord injury: the multicenter randomized Spinal Cord Injury Locomotor Trial. Neurorehabil Neural Repair 2007, 21: 25-35. 10.1177/1545968306295556

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Winchester P, Querry R: Robotic orthoses for body weight-supported treadmill training. Physical medicine and rehabilitation clinics of North America 2006, 17: 159-172. 10.1016/j.pmr.2005.10.008

    Article  PubMed  Google Scholar 

  6. Colombo G, Wirz M, Dietz V: Driven gait orthosis for improvement of locomotor training in paraplegic patients. Spinal Cord 2001,39(5):252-255. 10.1038/

    Article  CAS  PubMed  Google Scholar 

  7. 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 

  8. 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(4):672-680. 10.1016/j.apmr.2004.08.004

    Article  PubMed  Google Scholar 

  9. Mayr A, Kofler M, Quirbach E, Matzak H, Frohlich K, Saltuari L: Prospective, Blinded, Randomized Crossover Study of Gait Rehabilitation in Stroke Patients Using the Lokomat Gait Orthosis. Neurorehabil Neural Repair 2007,21(4):307-314. 10.1177/1545968307300697

    Article  PubMed  Google Scholar 

  10. Krewer C, Müller F, Husemann B, Heller S, Quintern J, Koenig E: The influence of different Lokomat walking conditions on the energy expenditure of hemiparetic patients and healthy subjects. Gait Posture 2007, 26: 372-377. 10.1016/j.gaitpost.2006.10.003

    Article  PubMed  Google Scholar 

  11. Pohl M, Werner C, Holzgraefe M, Kroczek G, Mehrholz J, Wingendorf I, Hoölig G, Koch R, Hesse S: Repetitive locomotor training and physiotherapy improve walking and basic activities of daily living after stroke: a single-blind, randomized multicentre trial (DEutsche GAngtrainerStudie, DEGAS). Clin Rehabil 2007, 21: 17-27. 10.1177/0269215506071281

    Article  CAS  PubMed  Google Scholar 

  12. Hornby TG, Campbell DD, Kahn JH, Demott T, Moore JL, Roth HR: Enhanced Gait-Related Improvements After Therapist- Versus Robotic-Assisted Locomotor Training in Subjects With Chronic Stroke: A Randomized Controlled Study. Stroke 2008,39(6):1786-1792. 10.1161/STROKEAHA.107.504779

    Article  PubMed  Google Scholar 

  13. Hidler J, Nichols D, Pelliccio M, Brady K, Campbell DD, Kahn JH, Hornby TG: Multicenter Randomized Clinical Trial Evaluating the Effectiveness of the Lokomat in Subacute Stroke. Neurorehabil Neural Repair 2009, 23: 5-13.

    Article  PubMed  Google Scholar 

  14. Westlake K, Patten C: Pilot study of Lokomat versus manual-assisted treadmill training for locomotor recovery post-stroke. Journal of NeuroEngineering and Rehabilitation 2009, 6: 18. 10.1186/1743-0003-6-18

    Article  PubMed Central  PubMed  Google Scholar 

  15. Nooijen C, ter Hoeve N, Field-Fote E: Gait quality is improved by locomotor training in individuals with SCI regardless of training approach. Journal of NeuroEngineering and Rehabilitation 2009, 6: 36. 10.1186/1743-0003-6-36

    Article  PubMed Central  PubMed  Google Scholar 

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

    Article  Google Scholar 

  17. Hidler J, Wisman W, Neckel N: Kinematic trajectories while walking within the Lokomat robotic gait-orthosis. Clinical Biomechanics 2008,23(10):1251-1259. 10.1016/j.clinbiomech.2008.08.004

    Article  PubMed  Google Scholar 

  18. Lotze M, Braun C, Birbaumer N, Anders S, Cohen LG: Motor learning elicited by voluntary drive. Brain 2003,126(4):866-872. 10.1093/brain/awg079

    Article  PubMed  Google Scholar 

  19. Kaelin-Lang A, Sawaki L, Cohen LG: Role of voluntary drive in encoding an elementary motor memory. J Neurophysiol 2005,93(2):1099-1103. 10.1152/jn.00143.2004

    Article  PubMed  Google Scholar 

  20. Dromerick A, Lum P, Hidler J: Activity-Based Therapies. NeuroRX 2006,3(4):428-438. 10.1016/j.nurx.2006.07.004

    Article  PubMed Central  PubMed  Google Scholar 

  21. Bernstein NA: The Co-ordination and regulation of movements. first edition. Pergamon Press Ltd.; 1967.

    Google Scholar 

  22. Lewek MD, Cruz TH, Moore JL, Roth HR, Dhaher YY, Hornby TG: Allowing Intralimb Kinematic Variability During Locomotor Training Poststroke Improves Kinematic Consistency: A Subgroup Analysis From a Randomized Clinical Trial. Physical therapy 2009.

    Google Scholar 

  23. Huang VS, Krakauer JW: Robotic neurorehabilitation: a computational motor learning perspective. Journal of NeuroEngineering and Rehabilitation 2009, 6: 5. 10.1186/1743-0003-6-5

    Article  PubMed Central  PubMed  Google Scholar 

  24. Riener R, Fuhr T: Patient-driven control of FES-supported standing up: a simulation study. Rehabilitation Engineering, IEEE Transactions on 1998,6(2):113-124. 10.1109/86.681177

    Article  CAS  Google Scholar 

  25. Jezernik S, Colombo G, Morari M: Automatic gait-pattern adaptation algorithms for rehabilitation with a 4-DOF robotic orthosis. IEEE Trans Robot Autom 2004,20(3):574-582. 10.1109/TRA.2004.825515

    Article  Google Scholar 

  26. Riener R, Frey M, Bernhardt M, Nef T, Colombo G: Human-centered rehabilitation robotics. Rehabilitation Robotics, 2005. ICORR 2005. 9th International Conference on 2005, 319-322. full_text

    Chapter  Google Scholar 

  27. Emken JL, Bobrow JE, Reinkensmeyer DJ: Robotic movement training as an optimization problem: designing a controller that assists only as needed. IEEE Int. Conf. on Rehabilitation Robotics (ICORR), Chicago 2005, 307-312. full_text

    Google Scholar 

  28. Reinkensmeyer DJ, Aoyagi D, Emken JL, Galvez JA, Ichinose W, Kerdanyan G, Maneekobkunwong S, Minakata K, Nessler JA, Weber R, Roy RR, de Leon R, Bobrow JE, Harkema SJ, Edgerton VR: Tools for understanding and optimizing robotic gait training. J Rehabil Res Dev 2006,43(5):657-670. 10.1682/JRRD.2005.04.0073

    Article  PubMed  Google Scholar 

  29. Banala SK, Agrawal SK, Scholz JP: Active Leg Exoskeleton (ALEX) for Gait Rehabilitation of Motor-Impaired Patients. In Proc. IEEE 10th Int. Conf. Rehabil. Robot. Noordwijk; 2007:401-407. full_text

    Google Scholar 

  30. van Asseldonk EHF, Veneman JF, Ekkelenkamp R, Buurke JH, van der Helm FCT, van der Kooij H: The Effects on Kinematics and Muscle Activity of Walking in a Robotic Gait Trainer During Zero-Force Control. IEEE Trans. Neural Syst. Rehabil. Eng 2008,16(4):360-370. 10.1109/TNSRE.2008.925074

    Article  PubMed  Google Scholar 

  31. Duschau-Wicke A, von Zitzewitz J, Banz R, Riener R: Iterative Learning Synchronization of Robotic Rehabilitation Tasks. Proc. IEEE 10th Int. Conf. Rehabil. Robot 2007, 335-340. full_text

    Google Scholar 

  32. Vallery H, Guidali M, Duschau-Wicke A, Riener R: Patient-Cooperative Control: Providing Safe Support without Restricting Movement. In World Congress on Medical Physics and Biomedical Engineering, September 7 - 12, 2009, Munich, Germany. Volume 25/9. Edited by: Dössel O, Schlegel WC. Berlin, Heidelberg: Springer Berlin Heidelberg; 2009:166-169. full_text

    Chapter  Google Scholar 

  33. Emken JL, Harkema SJ, Beres-Jones JA, Ferreira CK, Reinkensmeyer DJ: Feasibility of Manual Teach-and-Replay and Continuous Impedance Shaping for Robotic Locomotor Training Following Spinal Cord Injury. IEEE Trans. Biomed. Eng 2008, 55: 322-334. 10.1109/TBME.2007.910683

    Article  PubMed  Google Scholar 

  34. Marchal Crespo L, Reinkensmeyer D: Review of control strategies for robotic movement training after neurologic injury. Journal of neuroengineering and rehabilitation 2009, 6: 20. 10.1186/1743-0003-6-20

    Article  PubMed Central  PubMed  Google Scholar 

  35. Krebs H, Palazzolo JJ, Dipietro L, Ferraro M, Krol J, Rannekleiv K, Hogan N: Rehabilitation robotics: Performance-based progressive robot-assisted therapy. Autonomous Robots 2003, 15: 7-20. 10.1023/A:1024494031121

    Article  Google Scholar 

  36. Reinkensmeyer DJ, Maier MA, Guigon E, Chan V, Akoner O, Wolbrecht ET, Cramer SC, Bobrow JE: Do robotic and non-robotic arm movement training drive motor recovery after stroke by a common neural mechanism? Experimental evidence and a computational model. Conf. IEEE Engineering in Medicine and Biology Society (EMBS) 2009, 2439-2441. full_text

    Google Scholar 

  37. Duschau-Wicke A, von Zitzewitz J, Caprez A, Lünenburger L, Riener R: Path Control: A Method for Patient-Cooperative Robot-Aided Gait Rehabilitation. IEEE Transactions on Neural Systems and Rehabilitation Engineering 2010, 18: 38-48. 10.1109/TNSRE.2009.2033061

    Article  PubMed  Google Scholar 

  38. Colombo G, Joerg M, Schreier R, Dietz V: Treadmill training of paraplegic patients using a robotic orthosis. J Rehabil Res Dev 2000,37(6):693-700.

    CAS  PubMed  Google Scholar 

  39. Cai LL, Fong AJ, Otoshi CK, Liang Y, Burdick JW, Roy RR, Edgerton VR: Implications of assist-as-needed robotic step training after a complete spinal cord injury on intrinsic strategies of motor learning. J Neurosci 2006,26(41):10564-10568. 10.1523/JNEUROSCI.2266-06.2006

    Article  CAS  PubMed  Google Scholar 

  40. Banala SK, Kim SH, Agrawal SK, Scholz JP: Robot Assisted Gait Training With Active Leg Exoskeleton (ALEX). IEEE Trans. Neural Syst. Rehabil. Eng 2009, 17: 2-8. 10.1109/TNSRE.2008.2008280

    Article  PubMed  Google Scholar 

  41. Vallery H, Veneman J, van Asseldonk E, Ekkelenkamp R, Buss M, van Der Kooij H: Compliant actuation of rehabilitation robots. Robotics & Automation Magazine, IEEE 2008,15(3):60-69.

    Article  Google Scholar 

  42. Hermens HJ, Freriks B, Merletti R, Stegeman D, Blok J, Rau G, Disselhorst-Klug C, Haegg G: European Recommendations for Surface Electromyography. 1999.

    Google Scholar 

  43. Martin A: Specific influences of cerebellar dysfunctions on gait. Brain 2007,130(3):786-798. 10.1093/brain/awl376

    Article  Google Scholar 

  44. Giese MA, Poggio T: Morphable Models for the Analysis and Synthesis of Complex Motion Patterns. Int J Comp Vis 2000, 38: 59-73. 10.1023/A:1008118801668

    Article  Google Scholar 

  45. Gibbons JD: Nonparametric Statistical Inference. Marcel Dekker Ltd; 1985.

    Google Scholar 

  46. Perry J: Gait Analysis: Normal and Pathological Function. Delmar Learning; 1992.

    Google Scholar 

  47. West BT, Welch KB, Galecki AT: Linear Mixed Models: A Practical Guide Using Statistical Software. Chapmann Hall/CRC Press; 2006.

    Google Scholar 

  48. van Asseldonk EH, Wessels M, Stienen AH, van der Helm FC, van der Kooij H: Influence of haptic guidance in learning a novel visuomotor task. Journal of physiology, Paris 2009,103(3-5):276-285. 10.1016/j.jphysparis.2009.08.010

    Article  PubMed  Google Scholar 

  49. Hu XL, Tong Ky, Song R, Zheng XJ, Leung WWF: A Comparison Between Electromyography-Driven Robot and Passive Motion Device on Wrist Rehabilitation for Chronic Stroke. Neurorehabilitation and Neural Repair 2009,23(8):837-846. 10.1177/1545968309338191

    Article  CAS  PubMed  Google Scholar 

  50. Lam T, Wirz M, Lünenburger L, Dietz V: Swing Phase Resistance Enhances Flexor Muscle Activity During Treadmill Locomotion in Incomplete Spinal Cord Injury. Neurorehabil Neural Repair 2008,22(5):438-446. 10.1177/1545968308315595

    Article  PubMed  Google Scholar 

  51. Emken JL, Benitez R, Sideris A, Bobrow JE, Reinkensmeyer DJ: Motor adaptation as a greedy optimization of error and effort. Journal of neurophysiology 2007,97(6):3997-4006. 10.1152/jn.01095.2006

    Article  PubMed  Google Scholar 

  52. Reinkensmeyer DJ, Wolbrecht E, Bobrow J: A computational model of human-robot load sharing during robot-assisted arm movement training after stroke. Conference proceedings: ... Annual International Conference of the IEEE Engineering in Medicine and Biology Society. IEEE Engineering in Medicine and Biology Society. Conference 2007, 2007: 4019-4023.

    Google Scholar 

  53. Franklin DW, Burdet E, Peng Tee K, Osu R, Chew CM, Milner TE, Kawato M: CNS Learns Stable, Accurate, and Efficient Movements Using a Simple Algorithm. J Neurosci 2008,28(44):11165-11173. 10.1523/JNEUROSCI.3099-08.2008

    Article  CAS  PubMed  Google Scholar 

  54. Reinkensmeyer DJ, Akoner O, Ferris DP, Gordon KE: Slacking by the human motor system: computational models and implications for robotic orthoses. 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society, IEEE 2009, 2009: 2129-2132. full_text

    Article  Google Scholar 

  55. Wolbrecht ET, Chan V, Reinkensmeyer DJ, Bobrow JE: Optimizing Compliant, Model-Based Robotic Assistance to Promote Neurorehabilitation. IEEE Trans. Neural Syst. Rehabil. Eng 2008,16(3):286-297. 10.1109/TNSRE.2008.918389

    Article  PubMed  Google Scholar 

  56. von Zitzewitz J, Bernhardt M, Riener R: A Novel Method for Automatic Treadmill Speed Adaptation. IEEE Trans. Neural Syst. Rehabil. Eng 2007,15(3):401-409. 10.1109/TNSRE.2007.903926

    Article  PubMed  Google Scholar 

  57. Pennycott A, Hunt KJ, Coupaud S, Allan DB, Kakebeeke TH: Feedback Control of Oxygen Uptake During Robot-Assisted Gait. IEEE Transactions on Control Systems Technology 2010, 18: 136-142. 10.1109/TCST.2008.2009465

    Article  Google Scholar 

  58. Maynard FM, Bracken MB, Creasey G, Ditunno JF, Donovan WH, Ducker TB, Garber SL, Marino RJ, Stover SL, Tator CH, Others: International Standards for Neurological and Functional Classification of Spinal Cord Injury. Spinal Cord 1997, 35: 266-274. 10.1038/

    Article  PubMed  Google Scholar 

  59. Itzkovich M, Gelernter I, Biering-Sorensen F, Weeks C, Laramee MT, Craven BC, Tonack M, Hitzig SL, Glaser E, Zeilig G, Aito S, Scivoletto G, Mecci M, Chadwick RJ, Masry WS, Osman A, Glass CA, Silva P, Soni BM, Gardner BP, Savic G, Bergström EM, Bluvshtein V, Ronen J, Catz A: The Spinal Cord Independence Measure (SCIM) version III: reliability and validity in a multi-center international study. Disabil Rehabil 2007,29(24):1926-1933. 10.1080/09638280601046302

    Article  CAS  PubMed  Google Scholar 

  60. Dittuno PL, Dittuno JF: Walking index for spinal cord injury (WISCI II): scale revision. Spinal Cord 2001, 39: 654-656. 10.1038/

    Article  CAS  PubMed  Google Scholar 

Download references


The authors would like to thank all patients who participated in this study. Furthermore, we are grateful to Heike Vallery, Markus Wirz, Marc Bolliger, and Huub van Hedel for their support. This work was supported by NIDRR grant H133E070013.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Alexander Duschau-Wicke.

Additional information

Authors' contributions

AD and AC contributed equally to this work. AD and AC performed the measurements of all patients, data analysis, statistical analysis, and drafted the manuscript. RR participated in the design and coordination of the study and assisted with drafting the manuscript. All authors read and approved the final manuscript.

Alexander Duschau-Wicke, Andrea Caprez contributed equally to this work.

Authors’ original submitted files for images

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 cited.

Reprints and permissions

About this article

Cite this article

Duschau-Wicke, A., Caprez, A. & Riener, R. Patient-cooperative control increases active participation of individuals with SCI during robot-aided gait training. J NeuroEngineering Rehabil 7, 43 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: