SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Advances in Nonstationary Electrophysiological Signal Analysis and Processing.

Open Access Open Badges Research

Nonlinear analysis of electromyogram following gait training with myoelectrically triggered neuromuscular electrical stimulation in stroke survivors

Anirban Dutta1*, Bhawna Khattar2 and Alakananda Banerjee2

Author Affiliations

1 Department of Clinical Neurophysiology, Georg-August-University, Goettingen, Germany

2 Department of Physiotherapy and Rehabilitation, Max Super Specialty Hospital Saket, New Delhi, India

For all author emails, please log on.

EURASIP Journal on Advances in Signal Processing 2012, 2012:153  doi:10.1186/1687-6180-2012-153

The electronic version of this article is the complete one and can be found online at: http://asp.eurasipjournals.com/content/2012/1/153

Received:5 November 2011
Accepted:15 June 2012
Published:20 July 2012

© 2012 Dutta et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Neuromuscular electrical stimulation (NMES) facilitates ambulatory function after paralysis by activating the muscles of the lower extremities. The NMES-assisted stepping can either be triggered by a heel-switch (switch-trigger), or by an electromyogram (EMG)-based gait event detector (EMG-trigger). The command sources—switch-trigger or EMG-trigger—were presented to each group of six chronic (>6 months post-stroke) hemiplegic stroke survivors. The switch-trigger group underwent transcutaneous NMES-assisted gait training for 1 h, five times a week for 2 weeks, where the stimulation of the tibialis anterior muscle of the paretic limb was triggered with a heel-switch detecting heel-rise of the same limb. The EMG-trigger group underwent transcutaneous NMES-assisted gait training of the same duration and frequency where the stimulation was triggered with surface EMG from medial gastrocnemius (MG) of the paretic limb in conjunction with a heel-switch detecting heel-rise of the same limb. During the baseline and post-intervention surface EMG assessment, a total of 10 s of surface EMG was recorded from bilateral MG muscle while the subjects tried to stand steady on their toes. A nonlinear tool—recurrence quantification analysis (RQA)—was used to analyze the surface EMG. The objective of this study was to find the effect of NMES-assisted gait training with switch-trigger or EMG-trigger on two RQA parameters—the percentage of recurrence (%Rec) and determinism (%Det), which were extracted from surface EMG during fatiguing contractions of the paretic muscle. The experimental results showed that during fatiguing contractions, (1) %Rec and %Det have a higher initial value for paretic muscle than the non-paretic muscle, (2) the rate of change in %Rec and %Det was negative for the paretic muscle but positive for the non-paretic muscle, (3) the rate of change in %Rec and %Det significantly increased from baseline for the paretic muscle after EMG-triggered NMES-assisted gait training. Therefore, the study showed an improvement in paretic muscle function during a fatiguing task following gait training with EMG-triggered NMES. This study also showed that RQA parameters—%Rec and %Det—were sensitive to changes in paretic/non-paretic muscle properties due to gait training and can be used for non-invasive muscle monitoring in stroke survivors undergoing rehabilitation.

Neuromuscular electrical stimulation; Electromyogram; Foot-drop; Stroke; Recurrent quantification analysis


The World Health Organization defines stroke as “the rapid development of clinical signs and symptoms of a focal neurological disturbance lasting more than 24 h or leading to death with no apparent cause other than vascular origin” [1]. Worldwide, 15 million people suffer a stroke every year, and stroke is the second leading cause of disability after dementia. India-centered studies have presented a prevalence rate of 27–34/100,000 in the 35–44 age group to 822–1116/100,000 in the 75+ age group [2-4]. The prevalence of stroke in younger individuals is higher in India as compared with high-income nations [5]. The burden of stroke on the quality of life is much greater in this younger age group of stroke survivors [5-7]. Moreover, stroke drastically affects the economical productivity for that younger age-group adding further to the overall disease burden [5] and therefore early rehabilitation to alleviate disability has a special significance in developing countries.

The disability due to poor ability to walk in stroke patients is frequently caused by the “foot-drop” symptom that prevents the patient from being able to raise the foot during the swing phase of walking [8]. The foot slaps down on the ground after heel-strike and the toe drags during the swing phase. In addition to that the patients often lack from push-off power due to weakness in ankle plantarflexors (like, gastrocnemius muscle), which not only reduces walking speed, but also reduce ankle dorsiflexion angles during the early swing phase of gait [9]. Goldberg et al. [10] have shown that the ankle plantarflexors are important for knee flexion velocity during push-off. This leads to insufficient foot clearance that puts the patient at risk for stumbling and falling [9]. Consequently, impairment of walking has been mentioned most frequently (39–90%) as the most important disabling condition in community-dwelling stroke survivors and improving walking ability has been found to be the best way to reduce dependency [11,12]. Almost 30% of patients who survive stroke are affected by foot-drop, where the ankle dorsiflexors (e.g., tibialis anterior muscle) are functionally impaired [8,9].

Neuromuscular electrical stimulation (NMES)-based orthosis has been shown to enhance walking abilities, increasing gait speed while lowering effort and has only recently developed into a therapeutic intervention for stroke rehabilitation [8,13-22]. NMES involves electrical stimulation of nerves and muscles with continuous short pulses of electrical current [23]. Hausdorff and Ring [20] showed that the number of falls reduced significantly in hemiparetic patients who used FES to correct foot-drop. Although FES provides orthotic benefit but the neurophysiological therapeutic benefit has not been investigated thoroughly. Recently, different command modalities are available where NMES-assisted stepping can either be triggered with a heel-switch (switch-trigger) or with an electromyogram (EMG)-based gait event detector (EMG-trigger). Although both the modalities provide similar orthotic benefit [24,25] where they assist in lifting the foot to clear the ground during stepping action but the relative long-term therapeutic benefit of one modality over the other in facilitating recovery of neurophysiological function has not rigorously been shown. Preliminary research demonstrated the feasibility of EMG control of NMES in partial paralysis as well as its therapeutic benefits [26-31]. Dutta et al. [31] have shown therapeutic benefits in terms of incremental improvement in the classification of the gait event using a constant EMG-based classifier while the command muscle was undergoing NMES-assisted gait training. However, capturing the neurophysiological therapeutic effects in terms of changes in the properties of the surface EMG from paretic muscles will require more advanced analysis with tools like recurrence quantification analysis (RQA), which is the focus of the current study.

RQA is a technique for detecting changes in the states of a dynamical systems that does not need any a priori constraint on data size, stationarity, and statistical distribution [32]. RQA parameters are very sensitive in detecting changes in surface EMG due to fatigue [33]. Webber et al. [34] showed that small changes in surface EMG can be detected by RQA where the parameter—%Det reflects the deterministic attributes of the signal whereas the parameter—%Rec reflects the recurrence of the signal. RQA parameters have been validated as means for fatigue assessment with potential advantages such as the higher sensitivity to neurophysiological changes in the muscle status [35]. Therefore, the objective of this study was to investigate the effects of NMES-assisted gait training with two different command sources—switch-trigger or EMG-trigger—on two RQA parameters—the percentage of recurrence (%Rec) and determinism (%Det), which were extracted from the surface EMG during fatiguing contractions of the paretic muscle in 12 chronic stroke survivors.



Ten male and two female hemiplegic chronic (>6 months post-stroke) stroke survivors (age: 56–75 years) volunteered for this study. Informed consent was obtained from all the subjects before their participation. The study was approved by the Institutional Review Board of Max Superspeciality Hospitals Saket, New Delhi, India. The clinical inclusion and exclusion criteria are tabulated in Table 1, which were evaluated by a consulting neurologist.

Table 1. Inclusion and exclusion criteria for the clinical study

Experimental setup

Experimental setup and the signal flow are shown in the top panel of Figure 1. Surface EMG signals were collected from medial gastrocnemius (MG) muscle using Ag/AgCl electrodes with 2-cm inter-electrode distance following the SENIAM guidelines [36]. The surface EMG signals were amplified and low-pass filtered (anti-aliasing, frequency cutoff = 1000 Hz) with custom-made amplifier before being sampled at 2400 Hz by 16-bit data acquisition system (NI USB-6215, National Instruments, USA). The data acquisition system also sampled the heel-switch signal to monitor heel-rise and heel-strike of the paretic foot. The digital output of the data acquisition system activated a switching circuit (clamp) by a trigger pulse that disconnected the EMG electrode inputs from the amplifier and connected them to the common ground electrode from heel-rise to heel-strike when the stimulation was delivered. The gain of the custom-made amplifier was set to prevent saturation at the maximum muscle activity. The control signal for triggering NMES—the linear envelope (LE)—was computed after digitally band-pass filtering (5th order Butterworth, 3 dB bandwidth = 10–500 Hz), then full-wave rectification and then low-pass filtering (5th order Butterworth, 3 dB frequency cutoff = 3 Hz) the sampled EMG signal. Visual biofeedback of the control signal was provided to the subject on a computer monitor (bottom panel of Figure 1) while the subjects learned to trigger the stimulation with surface EMG from their paretic MG muscle.

thumbnailFigure 1. Top panel: Experimental setup for EMG-triggered NMES training. The signal flow is overlaid where the raw surface EMG was collected from MG muscle, which was amplified and low-pass filtered (anti-aliasing, frequency cutoff = 1000 Hz) before being sampled at 2400 Hz by 16-bit data acquisition system (NI USB-6215, National Instruments, USA) in a personal computer (PC). The LE was computed after digitally band-pass filtering (5th order Butterworth, 3 dB bandwidth = 10–500 Hz), then full-wave rectification and the low-pass filtering (5th order Butterworth, 3 dB frequency cutoff = 3 Hz) the sampled EMG signal. Visual biofeedback of the control signal (150 ms sliding window) was provided to the subject on a computer monitor (bottom panel) while the subjects learned to trigger the stimulation by activating their paretic MG muscle by exceeding a subject-specific threshold of 75% of their maximal voluntary contraction (MVC) in conjunction with heel-rise detected by the heel-switch of the ODFS Pace (Odstock, UK) electrical stimulator. Bottom panel: EMG visual biofeedback for EMG-triggered NMES training.

After achieving proficiency in EMG-triggered FES, the subjects in the EMG-triggered group (see the next section) walked with single-channel transcutaneous NMES device (ODFS Pace, Odstock, UK) which delivered electrical stimulation using saline-soaked sponge electrodes to activate tibialis anterior of the paretic leg during swing-phase of the gait cycle. The stimulation and recording electrodes were attached to a customized sleeve as shown in Figure 2, which was worn by the subject during each visit for gait training.

thumbnailFigure 2. Photo of the subject-specific textile compression sleeve with the EMG recording electrodes and the electrical stimulation electrodes.

NMES intervention

The subjects were randomly allocated to one of the two groups: GROUP EMG-triggered or GROUP Switch-triggered, as explained below.

Group EMG-triggered

In this group, the subjects received EMG-triggered transcutaneous NMES-assisted gait training on treadmill for 1 h, five times a week for 2 weeks. The NMES was triggered at the heel-rise of the paretic side that was detected with a heel-switch only if the LE of the MG crossed 75% of the MVC during preceding double-support (push-off) phase of the gait. Before walking with EMG-triggered transcutaneous NMES, the subjects were trained to activate paretic MG muscle for push-off using biofeedback of the LE control signal in order to trigger the NMES device.

Group switch-triggered

In this group, the subjects received heel-switch-triggered transcutaneous NMES-assisted gait training on treadmill for 1 h, five times a week for 2 weeks. The NMES was triggered at the heel-rise of the paretic side that was detected with a heel-switch.

In both the groups, NMES produced muscle contractions that mimic normal voluntary ankle movement, lifting the foot off the ground and improving gait, by applying electrical pulses to the common peroneal nerve through the skin surface. Stimulation electrodes were placed over the nerve using a textile sleeve (Figure 2), which were connected by leads to the NMES device (ODFS Pace, Odstock, UK).

Baseline and post-intervention RQA

During the baseline and post-intervention surface EMG assessment, subjects were asked to stand steady on their toes for a total of 10 s with heel-rise of both their feet while the surface EMG was recorded from fatiguing bilateral MG muscles. The EMG signals were amplified and low-pass filtered (anti-aliasing, frequency cutoff = 1000 Hz) with custom-made amplifier before being sampled at 2400 Hz by 16-bit data acquisition system (NI USB-6215, National Instruments, USA). Then the sampled EMG signals were digitally band-pass filtered (5th order Butterworth, 3 dB bandwidth = 10–500 Hz) and re-sampled at 1024 Hz to be stored for offline analysis. The data acquisition system also sampled bilateral heel-switch signals to monitor heel-rise and heel-strike of both the feet.

The RQA was conducted offline from 9 s of the surface EMG signal (half seconds of leading and trailing data were discarded) from paretic and non-paretic MG to estimate the parameters—%Det and %Rec from adjacent non-overlapping epochs of 1 s (1,024 data points) for each subject. The first step in RQA was embedding procedure where vectors representing the states of the system were extracted. The phase space was assumed to be 15-dimensional [37]. Therefore, D = 15-dimensional vectors were extracted from 1 s epoch of the surface EMG signal (i.e., data = {s (1), s(2),, s (1024)}) by lagging the signal by an integer number of samples (d = 3 [35]) as shown below,

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M1">View MathML</a>


The Euclidean distance matrix (<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M2">View MathML</a>) for the 15-dimensional vector was created where the elements were given by

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M3">View MathML</a>


The Euclidean distance matrix (<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M4">View MathML</a>) was then normalized (<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M5">View MathML</a>) by the average distance between the vectors

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M6">View MathML</a>


The recurrence map (RM) was then constructed from <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M7">View MathML</a> such that the pixel <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M8">View MathML</a> was ON if <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M9">View MathML</a> and OFF if <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M10">View MathML</a>, where the threshold was <a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M11">View MathML</a>[35]. In the RM, the pixels that were ON were defined as recurrent (Re) if they share local neighborhoods in higher dimensional space. Therefore, the percentage of the plot occupied by Re points was defined as %Rec and the percentage of Re forming lines parallel to diagonal was defined as %Det. Lines were constituted by two or more points (up to 20) that were diagonally adjacent. Therefore,

<a onClick="popup('http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.asp.eurasipjournals.com/content/2012/1/153/mathml/M12">View MathML</a>


where L was the number of points forming lines.

For each subject, the first-order regression line was computed from nine data points for %Rec and %Det during fatiguing contractions of the bilateral MG. The rate of change in %Rec and %Det was defined as the slope of the regression line and the initial value as the intercept of the regression line.


The regression lines capturing the variations in %Rec and %Det of the surface EMG during 9 s of fatiguing contractions of the paretic (affected limb) and non-paretic (unaffected limb) MG muscle before (baseline) FES-assisted gait training are shown in Figure 3 for all the subjects (N = 12). The rate of change (slope in%/second) of %Rec and %Det for the paretic MG with their 95% confidence interval (CI) were −0.411 (CI = −0.614 to −0.208) and −1.449 (CI = −1.978 to −0.92), respectively. In comparison, the rate of change (%/second) of %Rec and %Det of surface EMG for the non-paretic MG were 0.193 (CI = 0.072–0.313) and 1.149 (CI = 0.724–1.573), respectively. Therefore, %Rec and %Det decreased over 9 s of fatiguing contractions for the paretic MG while it increased for the non-paretic MG, for all the subjects pooled together (N = 12). The initial values (intercept in %) of %Rec and %Det were 5.978 (CI = 4.836–7.12) and 60.829 (CI = 57.852–63.806), respectively, for the paretic MG while 1.064 (CI = 0.387–1.741) and 41.442 (CI = 39.052–43.829), respectively, for the non-paretic MG. Therefore, there was a significant difference in the initial values of %Rec and %Det for all the subjects pooled together (N = 12), which were lower for the non-paretic MG than the paretic MG.

thumbnailFigure 3. The pre-intervention (baseline) variations in %Rec and %Det of the surface EMG during 9 s of fatiguing contractions of the paretic and non-paretic MG muscle for all the subjects pooled together (N= 12). The linear regression line with 95% CI is also shown. (a) Variation in %Rec for paretic MG, (b) variation in %Rec for non-paretic MG, (c) variation in %Det for paretic MG, (d) variation in %Det for non-paretic MG.

The subjects were randomly divided into two groups of six subjects each—EMG-triggered (N = 6) and switch-triggered (N = 6). The changes in the slope and intercept of the regression line capturing the variations in %Rec and %Det of the surface EMG during 9 s of fatiguing contractions of the paretic (affected limb) and non-paretic (unaffected limb) before (baseline) and after (post-intervention) FES-assisted gait training are tabulated in Additional file 1: Table S2 for each group. The data showed significant changes in the slope and intercept of the regression line for both the RQA parameters—%Rec and %Det of the paretic MG for the EMG-trigger group while the switch-trigger group showed no significant changes post-intervention. Specifically, the rate of change of %Rec during fatiguing contractions increased from −0.059 (CI = −0.088 to −0.031) to −0.008(CI = −0.017 to −0.0006) while that for %Det increased from −1.625 (CI = −2.371 to −0.88) to −0.4426 (CI = −0.913 to −0.028) for the EMG-trigger group. Also, the initial value of %Rec decreased from 7.727 (CI = 6.139–9.315) to 2.302 (CI = 1.807 to −2.797) while that for %Det decreased from 64.164 (CI = 59.968–68.359) to 49.901(CI = 47.252–52.55) for the EMG-trigger group. The corresponding values for the switch-trigger group remained statistically unchanged post-intervention.

Additional file 1. Table S2. The changes in the slope and intercept of the regression line for the parameters—%Det and %Rec during fatiguing contractions of paretic (affected) and non-paretic (unaffected) MG muscle before (baseline) and after (post-intervention) FES-assisted gait training, for the two groups—EMG-triggered and switch-triggered, of six subjects each. The table shows the Mean and 95% CI of the data. Significant differences between the baseline and post-intervention values are shaded.

Format: PDF Size: 6KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data


RQA has been shown to be promising for detecting neurophysiological changes in the muscle due to fatigue or other factors [35]. Farina et al. [35] showed that both the parameters—%Rec and %Det—were influenced by the conduction velocity and synchronization of the motor units of the muscle and therefore sensitive to the muscle properties. They also showed that %Det was highly correlated with the spectral variables and to some extent %Rec as well, and moreover both the %Det and %Rec were more sensitive to muscle fatigue than the spectral variables. The mean spectral frequency (MNF) and %Det were highly correlated with a negative correlation coefficient close to −1.

In this study, the parameters from RQA captured the differences in the muscle properties between paretic and non-paretic MG muscle (Figure 3) as well as the changes affected by NMES-assisted gait training interventions—EMG-triggered versus switch-triggered (Additional file 1: Table S2). From the baseline data, it was found that the parameters—%Rec and %Det—had a higher initial value for the paretic MG than the non-paretic MG. A higher value for %Det indicated a lower MNF for the paretic MG when compared to the non-paretic MG. Moreover, the rate of change of %Rec and %Det during fatiguing contraction was negative for the paretic MG and positive for the non-paretic MG. A positive rate of change in %Det indicated a decrease in MNF during fatiguing contractions which is normal [38]; however, a negative rate of change indicated an increase in MNF which may be pathological for the paretic MG. Moreover, the rate of change in %Rec and %Det during the fatiguing contractions of the paretic MG increased after NMES-assisted gait training for the EMG-trigger group. Felici et al. [39] have shown a higher rate of change in %Det in weight lifters when compared to the control group during sustained isometric contractions. A significant increase in the rate of change in %Det during fatiguing contractions may indicate exercise-effect on the paretic MG muscle brought by EMG-triggered NMES-assisted gait training. Therefore, we found that our short-duration EMG-triggered FES-assisted gait therapy improved neurophysiological properties of the muscle and a more intensive treatment with duration beyond 2 weeks might have shown additional benefits [40-43].

However, lack of accessible rehabilitation facilities and high cost of rehabilitation at the well-equipped clinics are current challenges in India which usually lead to high-dropout of individuals from rehabilitation therapy following stroke. Such de-conditioned chronic stroke survivor will need to recondition themselves with a gradual increase in the intensity (number of hours per day) and frequency (number of days per week) over the duration of the FES-assisted gait therapy. In fact, about 20% of stroke survivors suffer muscle atrophy which underlie worsening metabolic fitness in the chronic phase of stroke including gross muscular atrophy, altered muscle molecular phenotype, increased intramuscular area fat, elevated tissue inflammatory markers, and diminished peripheral blood flow dynamics [44]. EMG-triggered FES-assisted gait training may alleviate these debilitating conditions where increased intensity and frequency of rehabilitation may help [40]. In fact recent studies in India on therapeutic benefits of FES-assisted gait training in conjunction with conventional physiotherapy have shown the additive effects of FES on reducing spasticity, improving dorsiflexor strength, improving walking ability, and improving metabolic fitness [41,42]. However, stroke presents with heterogeneous patient-specific impairments in motor, sensory, tone, visual, perceptual, cognition, aphasia, apraxia, coordination, and equilibrium. Therefore, the functional limitations following stroke are varied, including gait dysfunction, fall risk, limited activities of daily living, difficulties in swallowing, reduced upper extremity function, altered communication, besides others. Based on the residual function of a stroke survivor, the number of muscles that can serve as EMG command source [31], the number of muscles that need FES assistance, intensity, frequency, and duration of the EMG-triggered FES-assisted gait training need to be decided. The ability of a stroke survivor to undergo FES-assisted gait therapy also depends on their cardiovascular and neuromuscular capacity besides psychological factors such as motivation. Also, the FES-assisted gait training should be started as soon as the patient becomes clinically stable since early intervention has shown better functional outcomes, survival rates, and reduced length of required therapy [43]. Therefore, we are currently investigating a home-based rehabilitation model where stroke survivors can use a low-cost EMG-triggered FES device daily. However, the efficacy of such home-based models in affecting therapeutic benefits needs to be validated with neurophysiological testing with advanced computational tools like RQA where parameters like %Rec and %Det can non-invasively monitor changes in muscle properties due to NMES-assisted gait training during rehabilitation.


The overarching goal of NMES-assisted gait training is to facilitate long-term therapeutic changes in muscle properties. However, a sensitive biomarker is needed to capture these neurophysiological changes during motor rehabilitation. In this study, the parameters derived from RQA captured the significant changes in muscle properties brought by EMG-triggered NMES-assisted gait training when compared to switch-triggered NMES-assisted gait training. The subjects were forced to volitionally activate the paretic MG muscle to generate the EMG-command signal during EMG-triggered NMES-assisted gait training which exercised the paretic MG muscle during a functional task. Therefore, EMG-triggered NMES-assisted gait training could be more effective in providing long-term therapeutic benefit in muscle function than switch-triggered NMES-assisted gait training. This study also provided an understanding of RQA parameters—%Rec and %Det, computed from surface EMG of the paretic and non-paretic muscles during fatiguing contractions.

Competing interests

The authors declare no existence of competing interests during the conduction of this study where the professional judgment concerning patients’ welfare and the validity of the research may be influenced by financial gain or personal opinion.


Dr. Dutta (Ph.D.) was supported by the Alexander von Humboldt Foundation, Germany. Drs. Banerjee (P.T.) and Khattar (P.T.) were supported by the Max Superspecialty Hospitals Saket, New Delhi, India. The voluntary help received from the subjects is gratefully acknowledged. The clinical study is registered with the Clinical Trials Registry of India (http://www.ctri.nic.in/Clinicaltrials/pmaindet2.php?trialid=3821).


  1. World Health Organisation, in Preventing Chronic Diseases, ed. by . A Vital Investment (World Health Organisation, Geneva, Switzerland, 2005)

  2. P Dalal, M Bhattacharjee, J Vairale, P Bhat, UN millennium development goals: can we halt the stroke epidemic in India? Ann. Indian Acad. Neurol 10, 130–136 (2007). Publisher Full Text OpenURL

  3. P Dalal, S Malik, M Bhattacharjee, ND Trivedi, J Vairale, P Bhat, Population-based stroke survey in Mumbai, India: incidence and 28-day case fatality. Neuroepidemiology 31, 254–261 (2008). PubMed Abstract | Publisher Full Text OpenURL

  4. SE Sridharan, JP Unnikrishnan, S Sukumaran, PN Sylaja, SD Nayak, PS Sarma, K Radhakrishnan, Incidence, types, risk factors and outcome of stroke in a developing country: the Trivandrum Stroke Registry. Stroke 40, 1212–18 (2009). PubMed Abstract | Publisher Full Text OpenURL

  5. M Tripathi, D Vibha, Stroke in young in India. Stroke Res. Treatment 2011, 1–6 (2011)

  6. B Shah, Development of Sentinel Health Monitoring Centers for Surveillance of Risk Factors for Non-Communicable Disease in India (Report from the Division of Non-communicable Diseases (Indian Council of Medical Research, New Delhi, 2005)

  7. B Shah, P Mathur, Workshop Report on Stroke Surveillance in India (Division of noncommunicable Diseases, Indian Council of Medical Research, New Delhi, 2006)

  8. MH Granat, DJ Maxwell, ACB Ferguson, KR Lees, JC Barbanel, Peroneal stimulator: evaluation for the correction of spastic foot-drop in hemiplegia. Arch. Phys. Med. Rehabil. 77(1), 19–24 (1996). PubMed Abstract | Publisher Full Text OpenURL

  9. V Weerdesteyn, M de Niet, HJ van Duijnhoven, AC Geurts, Falls in individuals with stroke. J Rehabil. Res. Dev. 45(8), 1195–1213 (2008). PubMed Abstract | Publisher Full Text OpenURL

  10. SR Goldberg, FC Anderson, MG Pandy, SL Delp, Muscles that influence knee flexion velocity in double support: implications for stiff-knee gait. J Biomech. 37(8), 1189–1196 (2004). PubMed Abstract | Publisher Full Text OpenURL

  11. NA Solomon, HA Glick, CJ Russo, J Lee, KA Schulman, Patient preferences for stroke outcomes. Stroke 25, 1721–1725 (1994). PubMed Abstract | Publisher Full Text OpenURL

  12. H Barbeau, J Fung, The role of rehabilitation in the recovery of walking in the neurological population. Curr. Opin. Neurol. 14(6), 735–740 (2001). PubMed Abstract | Publisher Full Text OpenURL

  13. R Merletti, A Andina, M Galante, I Furlan, Clinical experience of electronic peroneal stimulators in 50 hemiparetic patients. Scand. J Rehabil. Med. 11(3), 111–121 (1979). PubMed Abstract OpenURL

  14. U Bogataj, N Gros, M Malezic, B Kelih, M Kilajic, R Acimovic, Restoration of gait during two to three weeks of therapy with multichannel electrical stimulation. Phys. Ther. 69(5), 319–327 (1989). PubMed Abstract | Publisher Full Text OpenURL

  15. PN Taylor, JH Burridge, AL Dunkerley, DE Wood, JA Norton, C Singleton, ID Swain, Clinical use of the Odstock dropped foot stimulator: Its effect on speed and effort of walking. Arch. Phys. Med. Rehabil. 80(12), 1577–1583 (1999). PubMed Abstract | Publisher Full Text OpenURL

  16. LR Sheffler, MT Hennessey, GG Naples, J Chae, Peroneal nerve stimulation versus an ankle foot orthosis for correction of footdrop in stroke: impact on functional ambulation. Neurorehabil. Neural Repair 20(3), 355–360 (2006). PubMed Abstract | Publisher Full Text OpenURL

  17. AI Kottink, HJ Hermens, AV Nene, MJ Tenniglo, CG Groothuis-Oudshoorn, MJ IJzerman, Therapeutic effect of an implantable peroneal nerve stimulator in subjects with chronic stroke and footdrop: a randomized controlled trial. Phys. Ther 88(4), 437–448 (2008). PubMed Abstract | Publisher Full Text OpenURL

  18. JH Burridge, PN Taylor, ID Swain, The effects of common peroneal stimulation on the effort and speed of walking: a randomized controlled trial with chronic hemiplegic patients. Clin. Rehabil. 11(3), 201–210 (1997). PubMed Abstract | Publisher Full Text OpenURL

  19. JH Burridge, M Ladouceur, Clinical and therapeutic applications of neuromuscular stimulation: a review of current use and speculation into future developments. Neuromodulation 4, 147–155 (2001). PubMed Abstract | Publisher Full Text OpenURL

  20. JM Hausdorff, H Ring, Effects of a new radio frequency controlled neuroprosthesis on gait symmetry and rhythmicity in patients with chronic hemiparesis. Am. J Phys. Med. Rehabil. 87(1), 4–13 (2007)

  21. JJ Daly, J Zimbelman, KL Roenigk, JP McCabe, JM Rogers, K Butler, R Burdsall, JP Holcomb, EB Marsolais, RL Ruff, Recovery of coordinated gait: randomized controlled stroke trial of functional electrical stimulation (FES) versus no FES, with weight-supported treadmill and over-ground training. Neurorehabil. Neural Repair 25(7), 588–96 (2011). PubMed Abstract | Publisher Full Text OpenURL

  22. JJ Daly, RL Ruff, Feasibility of combining multi-channel functional neuromuscular stimulation with weight-supported treadmill training. J Neurol. Sci. 225, 105–115 (2004). PubMed Abstract | Publisher Full Text OpenURL

  23. R Kobetic, RJ Triolo, EB Marsolais, Muscle selection and walking performance of multichannel FES systems for ambulation in paraplegia. IEEE Trans. Rehabil. Eng. 5(1), 23–29 (1997). PubMed Abstract | Publisher Full Text OpenURL

  24. A Dutta, R Kobetic, RJ Triolo, Ambulation after incomplete spinal cord injury with EMG-triggered functional electrical stimulation. IEEE Trans. Biomed. Eng. 55(2 Pt 1), 791–794 (2008). PubMed Abstract | Publisher Full Text OpenURL

  25. A Dutta, R Kobetic, RJ Triolo, Gait initiation with electromyographically triggered electrical stimulation in people with partial paralysis. J Biomech. Eng. 131(8), 0810021–9 (2009)

  26. L Vodovnik, C Long, JB Reswick, A Lippay, D Starbuck, Myo-electric control of paralyzed muscles. IEEE Trans. Biomed. Eng. 12(3), 169–172 (1965). PubMed Abstract OpenURL

  27. DA Bolton, JH Cauraugh, HA Hausenblas, Electromyogram-triggered neuromuscular stimulation and stroke motor recovery of arm/hand functions: a meta-analysis. J Neurol. Sci. 223(2), 121–127 (2004). PubMed Abstract | Publisher Full Text OpenURL

  28. G Francisco, J Chae, H Chawla, S Kirshblum, R Zorowitz, G Lewis, S Pang, Electromyogram-triggered neuromuscular stimulation for improving the arm function of acute stroke survivors: a randomized pilot study. Arch. Phys. Med. Rehabil. 79(5), 570–575 (1998). PubMed Abstract | Publisher Full Text OpenURL

  29. J Cauraugh, K Light, S Kim, M Thigpen, M Beh Thigpen, A Behrman, Chronic motor dysfunction after stroke: recovering wrist and finger extension by electromyography triggered neuromuscular stimulation. Stroke 31(6), 1360–1364 (2000). PubMed Abstract | Publisher Full Text OpenURL

  30. RW Fields, Electromyographically triggered electric muscle stimulation for chronic hemiplegia. Arch. Phys. Med. Rehabil. 68(7), 407–414 (1987). PubMed Abstract OpenURL

  31. A Dutta, R Kobetic, RJ Triolo, An objective method for selecting command sources for myoelectrically controlled lower extremity neuroprostheses. J Rehabil. Res. Dev. 48(8), 935–948 (2011). PubMed Abstract | Publisher Full Text OpenURL

  32. JP Eckmann, SO Kamphorst, D Ruelle, Recurrence plots of dynamical systems. Europhys. Lett 4, 973–977 (1987). Publisher Full Text OpenURL

  33. S Ikegawa, M Shinohara, T Fukunaga, JP Zbilut, CL Webber, Nonlinear time-course of lumbar muscle fatigue using recurrence quantifications. Biol. Cybern. 82, 373–382 (2000). PubMed Abstract | Publisher Full Text OpenURL

  34. CL Webber, MA Schmidt, JM Walsh, Influence of isometric loading on biceps EMG dynamics as assessed by linear and nonlinear tools. J Appl. Physiol 78, 814–822 (1995). PubMed Abstract | Publisher Full Text OpenURL

  35. D Farina, L Fattorini, F Felici, G Filligoi, Nonlinear surface EMG analysis to detect changes of motor unit conduction velocity and synchronization. J Appl. Physiol 93, 1753–1763 (2002). PubMed Abstract | Publisher Full Text OpenURL

  36. R Merletti, H Hermens, Introduction to the special issue on the SENIAM European Concerted Action. J Electromyogr. Kinesiol. 10(5), 283–6 (2000). PubMed Abstract | Publisher Full Text OpenURL

  37. G Filligoi, F Felici, Detection of hidden rhythms in surface EMG signals with a nonlinear time-series tool. Med. Eng. Phys 21, 439–448 (1999). PubMed Abstract | Publisher Full Text OpenURL

  38. R Merletti, M Knaflitz, CJ De Luca, Myoelectric manifestations of fatigue in voluntary and electrically elicited contractions. J Appl. Physiol 69, 1810–1820 (1990). PubMed Abstract | Publisher Full Text OpenURL

  39. F Felici, A Rosponi, P Sbriccoli, GC Filligoi, L Fattorini, M Marchetti, Linear and nonlinear analysis of surface electromyograms in weightlifters. Eur. J Appl. Physiol 84, 337–342 (2001). PubMed Abstract | Publisher Full Text OpenURL

  40. G Kwakkel, RC Wagenaar, TW Koelman, GJ Lankhorst, JC Koetsier, Effects of intensity of rehabilitation after stroke. A research synthesis. Stroke 28(8), 1550–1556 (1997). PubMed Abstract | Publisher Full Text OpenURL

  41. SK Sabut, PK Lenka, R Kumar, M Mahadevappa, Effect of functional electrical stimulation on the effort and walking speed, surface electromyography activity, and metabolic responses in stroke subjects. J Electromyogr. Kinesiol. 20(6), 1170–1177 (2010). PubMed Abstract | Publisher Full Text OpenURL

  42. SK Sabut, C Sikdar, R Kumar, M Mahadevappa, Functional electrical stimulation of dorsiflexor muscle: effects on dorsiflexor strength, plantarflexor spasticity, and motor recovery in stroke patients. NeuroRehabilitation 29(4), 393–400 (2011). PubMed Abstract | Publisher Full Text OpenURL

  43. M Musicco, L Emberti, G Nappi, C Caltagirone, Early and long-term outcome of rehabilitation in stroke patients: the role of patient characteristics, time of initiation, and duration of interventions. Arch. Phys. Med. Rehabil. 84(4), 551–558 (2003). PubMed Abstract | Publisher Full Text OpenURL

  44. FM Ivey, RF Macko, AS Ryan, CE Hafer-Macko, Cardiovascular health and fitness after stroke. Topics Stroke Rehabil. Winter 12(1), 1–16 (2005). Publisher Full Text OpenURL