PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Comput Biol Med. Author manuscript; available in PMC 2010 December 1.
Published in final edited form as:
PMCID: PMC2784179
NIHMSID: NIHMS147097

An EMG-driven Model to Estimate Muscle Forces and Joint Moments in Stroke Patients

Abstract

Individuals following stroke exhibit altered muscle activation and movement patterns. Improving the efficiency of gait can be facilitated by knowing which muscles are affected and how they contribute to the pathological pattern. In this paper we present an electromyographically (EMG) driven musculoskeletal model to estimate muscle forces and joint moments. Subject specific EMG for the primary ankle plantar and dorsiflexor muscles, and joint kinematics during walking for four subjects following stroke were used as inputs to the model to predict ankle joint moments during stance. The model’s ability to predict the joint moment was evaluated by comparing the model output with the moment computed using inverse dynamics. The model did predict the ankle moment with acceptable accuracy, exhibiting an average R2 value ranging between 0.87 and 0.92, with RMS errors between 9.7% and 14.7%. The values are in line with previous results for healthy subjects, suggesting that EMG-driven modeling in this population of patients is feasible. It is our hope that such models can provide clinical insight into developing more effective rehabilitation therapies and to assess the effects of an intervention.

Keywords: EMG, Musculoskeletal model, Stroke, Hill-type muscle model, joint moment

INTRODUCTION

Abnormal movement patterns are associated with disability following stroke [1, 2]. Knowledge of internal forces and moments during movement is important for developing better rehabilitation programs for this population. Unfortunately, in vivo muscle force measurements are invasive and therefore not practical. For this reason, optimization techniques have been employed to predict muscle forces using a variety of cost functions [3, 4]. Other studies using forward dynamics simulations optimize muscle excitation (i.e., activation) patterns to reproduce recorded movement kinematics and to estimate muscle forces. In these studies, the cost function was the error between measured and simulated joint kinematics and ground reaction forces [5, 6]. These cost functions may have difficulty predicting co-contraction and cannot account for different muscle synergy patterns during various tasks [7]. EMG-driven biomechanical models are not subject to the same limitations since the models use recorded EMG to determine muscle force [711]. That is, these models use EMG as inputs rather than attempting to predict how muscles are activated to produce a given movement.

Although EMG-driven models have been used to investigate movement dynamics for healthy individuals during a variety of tasks, they have not been used to model muscle forcesproduced by people following stroke. Post-stroke patients often have difficulty generating an effective push-off during late stance and they exhibit abnormal activation patterns [12]. Abnormal and more variable muscle activity presents a challenging paradigm to test the robustness of our EMG-driven musculoskeletal model.

In this paper we describe our EMG-driven model and how it is implemented to estimate muscle forces and joint moments. Of specific interest was: (1) to provide a comprehensive description of the model structure, and (2) evaluate the model’s ability to predict the ankle joint moment for patients following stroke. The ability to predict joint dynamics and muscle forces in patients following stroke would benefit treatment planning and evaluation of therapies geared at improving patient function.

METHODS

Our EMG-driven model uses subject recorded joint kinematics and EMGs as inputs to estimate individual muscle forces and joint moments [10, 11]. The EMG-driven model is comprised of three separate but related sub-models: (1) anatomical, (2) muscle activation dynamics, and (3) muscle contraction dynamics. Sub-models (2) and (3) require the specification of several muscle specific parameters. These parameters characterize the force generating potential of each muscle and define how the recorded EMGs are transformed to muscle activations. Choice of parameters greatly affects the model’s ability to predict joint moments. Stochastic optimizationis used to determine ideal values for the muscle parameters during a model-tuning process which is described in the section entitled “Model Tuning”.

Anatomical model

The lower limb anatomical model was developed using SIMM (Motion Analysis Corp.) [13], and extended by Lloyd and Buchanan [14]. This model includes the four major plantar and dorsiflexor muscles: tibialis anterior (TA), lateral gastrocnemius (LG), medial gastrocnemius (MG) and soleus (Sol). The anatomical model was used to determine individual muscle tendon lengths and moment arms for the primary ankle plantar and dorsiflexors.

The muscle activation dynamics model transforms raw EMG to muscle activation. Firstly, the raw EMG were high-pass filtered to remove movement artifact, then full wave rectified, normalized by peak rectified EMG values obtained during maximum voluntary contraction (MVC), and then low-pass filtered to account for the natural filtering property of muscle [11]. A recursive filter was then used to determine neural activation, u(t), from the rectified, normalized and filtered EMG input, e(t). This process characterized the time varying features of EMG and was represented in the following discrete form [14, 15]:

u(t)=α·e(td)β1·u(t1)β2·u(t2),
(1)

where d is the electromechanical delay, and α, β1, and β2 are coefficients that define the second-order dynamics. To realize a positive stable solution of Equation (1), the following must hold true:

β1=γ1+γ2β1=γ1×γ2αβ1β2=1
(2)

where |γ1| < 1, |γ2| < 1, as discussed by Buchanan et al. [11].

Woods and Bigland-Ritchie [16] reported a nonlinear force-EMG relationship at low levels of force for some muscles and a linear relationship at higher levels. To account for this, we used a logarithmic function for low levels of EMG and a linear function for high values, as shown below:

a(t)=dln(c·u(t)+1)0u(t)<0.3a(t)=m·u(t)+b0.3u(t)<1
(3)

where u(t) is the neural activation and a(t) is the muscle activation. The coefficients c, d, m, and b can be solved and reduced to a single parameter, A, to characterize the curvature of the relationship [17]. A varies from 0.0 to approximately 0.12.

Muscle contraction dynamics model

The muscle activations and joint kinematics were used as inputs to a modified Hill-type muscle model to calculate individual muscle forces. The muscle-tendon unit was modeled as a muscle fiber in series with a tendon. The muscle fiber had a contractile element in parallel with an elastic element and a damping element (Figure 1). Pennation angle (3) is the angle between the lines of action of the tendons and the muscle fiber. lm is muscle fiber length, lt is total length of the tendons (each one is half of the length), and lmt is musculotendon length. Muscle-tendon force (F) was calculated as below:

F=Ft=Fmax[FA(lm)·FV(vm)·a(t)+FP(lm)+bm·vm]·cos(φ)
(4)
Figure 1
Schematic of the Hill-type muscle model. The muscle-tendon unit has a muscle fiber in series with the tendons, with muscle-tendon force (F) passing through both components.

Where, Ft is tendon force, and Fm is muscle fiber force that is comprised of the three parallel elements. Fmax is maximum isometric muscle force, lm is normalized muscle fiber length, vm is normalized muscle fiber velocity, a(t) is muscle activation. FA (lm) represents the active force-length relationship [18], FV (vm) represents the force-velocity relationship [1921], FP (lm) represents the passive elastic force-length relationship [22], and bm is the damping factor [23]. These parameters are normalized to maximum isometric muscle force, optimal fiber length (lom) and maximum muscle contraction velocity (vmax).

Pennation angle ([var phi](t)) was calculated as described by Scott and Winter [24], and is shown below:

φ(t)=sin1(lmosinφolm(t))
(5)

where lm(t) is the muscle fiber length at time t, and [var phi]o is the pennation angle at muscle optimal fiber length lmo.

The coupling relationship between optimal fiber length and muscle activation (i.e., optimal fiber length increases with decreasing activation) reported by Guimaraes et al. [25] and Huijing [26] is implemented in our model using the mathematical description of Lloyd and Besier [10]:

lmo(t)=lmo(λ(1a(t))+1)
(6)

where λ is the percentage change in optimal fiber length, lmo(t) is the optimal fiber length at time t.

Tendon force varies with tendon strain ε only when tendon length lt is greater than tendon slack length lst; otherwise the tendon force is zero as described by Zajac [23]. Subsequently the normalized tendon force, Ft, is given by

Ft=0ε0Ft=1480.3ε20<ε<0.0127Ft=37.5ε0.2375ε0.0127
(7)

where tendon strain was defined as

ε=(ltlst)/lst
(8)

Muscle tendon length and activation data were used as input to the muscle model, and then muscle fiber lengths were calculated by forward integration of the fiber velocities using a Runge-Kutta–Fehlberg algorithm. After the muscle fiber lengths and velocities were determined, the muscle force was calculated using Equation (4). Once individual muscle forces were computed, they were multiplied by their respective muscle moment arms derived from the anatomical model to calculate the net joint moment.

Model Tuning

In our muscle activation dynamics model, γ1, γ2, d, A are four parameters that account for the relationship between EMG and muscle activation. In our muscle contraction dynamics model, optimal fiber length (lom), tendon slack length ( lst) are anatomical parameters for each muscle, and λ is the global percentage change in lom as a function of muscle activation. In addition, to account for differences in muscle strength due to variations in PCSA between people, we scale the maximum isometric force, Fmax, for the flexor and extensor muscles [27]. Because these parameters are likely subject-specific, we begin by assigning each parameter an average value obtained from the literature. Initial values are allowed to vary within physiolgoical limits based on data reported in the literature Yamaguchi et al. [28].

Our calibration process uses a hybrid forward-inverse dynamics approach, which assumesthe joint moment estimated using the EMG-driven model should match the moment computed using inverse dynamics. A simulated parallel annealing algorithm [29] was used during the calibration process to tune the model parameters to minimize the difference between the EMG-driven model joint moment and the moment computed using inverse dynamics. The initial values of lst and lom for each muscle were constrained to change within ±15% and ±5% of their initial values respectively. Ranges for γ1, γ2 and A were described earlier; the range for electromechanical delay, d, was between 10 – 100 ms [30]. Maximum isometric force was allowed to vary between 50% and 200% of the initial value for Fmax. It is important to note that Fmax for the LG, MG and Sol were constrained to scale by the same amount to prevent unrealistic scenarios in which plantarflexion is dominated by a single muscle. Initial values for γ1, γ2, d, A and Fmax were set at the midpoint of their upper and lower limits. Once the model was tuned by adjusting the parameters, it was used to predict the ankle moment using EMG and kinematics for other trials on which it was not tuned.

DATA COLLECTION AND MODEL CALCULATIONS

Four stroke patients who gave informed consent were included in this experiment (Table 1). The subjects walked at self-selected speeds. Motion analysis and force plate data were used to compute joint kinematics and moments for four walking trials. In addition, EMG was recorded from the tibialis anterior (TA), medial gastrocnemius (MG), lateral gastrocnemius (LG) and soleus (Sol) muscles. Maximum voluntary isometric contraction (MVIC) trials were collected fornormalization of the EMG. data. The raw EMG were first high-pass filtered using a forward and reverse pass 4th order Butterworth filter (cutoff frequency 50 Hz), then full wave rectified, normalized by peak rectified EMG values obtained from MVC trials, and then low-pass filtered using a forward and reverse pass 4th-order Butterworth filter (cutoff frequency 4 Hz).

Table 1
An overview of clinical information regarding the four subjects.

The EMG-driven model was tuned to each of the four walking trial for each subject. The coefficient of determination (R2) between the model (i.e., EMG-driven) and calculated (i.e., Inverse dynamics) and the normalized root mean square (RMS) error (normalized to peak-to-peak joint moment) were calculated for each of the tuning trials. Once the subject’s muscle parameters were tuned to a particular trial, the model was given EMG and joint kinematics for the subject’s other three walking trials and the tuned-model was used to predict the ankle moment. For example, the model was tuned to trial 1 and then used to predict trials 2, 3 and 4. Similarly, the model was tuned to trial 2 and used to predict trials 1, 3 and 4.

RESULTS

Gait in people who have had strokes is markedly different than in healthy subjects. For impaired subjects in this study, circumduction was present, likely to compensate for dorsiflexion weakness and prevent toe-stubbing during swing. The EMG patterns were also notably different than normal [31]. For example, the triceps surae (gastrocnemii and soleus) were active during initial contact, a time when the plantarflexors are normally silent. In addition, an increase in tibialis anterior was noted following mid-stance, with peak values occurring at roughly the same time as triceps surae during late stance (Figure 2). In addition, there was noticable variability in the activation profiles.

Figure 2
Examplar EMG for the primary ankle plantar and dorsiflexor muscles during stance for one subject. Muscle activations depicted in bold are the resulting profiles after tuning the model. The dashed line represents the muscle activity used when predicting ...

Despite these differences, the model was able to be tuned to the subjects to capture the ankle joint moments accurately (Figure 3). When average values obtained from the literature were used to characterize the muscle parameters, poor agreement was observed between the model output and the moment computed using inverse dynamics (i.e., measured). However, once the model parameters were tuned to fit the subjects, the model output was in close agreement to the measured moment. The R2 value for all the four subjects was high (Table 2), indicating the model did an excellent job of following the shape of the joint moment computed using inverse dynamics. Excellent agreement in shape and similarity in magnitude and times to peak were observed when comparing ankle moments estimated from the model with the experimentally derived values (Figure 3). The normalized RMS error for the tuning trials ranged between 1.9% – 6.9%. This improvement implies that initial values for the muscle parameters obtained from the literature were clearly not appropriate for the subject.

Figure 3
The bold solid line is the ankle moment during stance for 1 subject computed using inverse dynamics (i.e., measured). The thinsolid line depicts the modeled moment using average values obtained from the literature for the muscle parameters. Poor agreement ...
Table 2
Average R2 and normalized RMS errors expressed as a percentage of the peak plantar flexion moment. The model was tuned to trial 1 and then used to predict trials 2 and 3. Similarly, the model was tuned to trial 2 and used to predict trials 1 and 3. Thus, ...

The model also did well in predicting joint moments (without re-tuning) given novel data. Once the muscle parameters were tuned to a particular subject for one trial, the model was then used to predict the ankle moment using EMG and kinematics for a different trial for that subject (Figure 4). Note that different activation patterns were used for the trial to which the model wastuned and the muscle activity for the trial that was predicted (Figure 2). In addition, notice the difference in magnitude between the trial used to tune the model (Figure 3) and the ankle moment that was predicted (Figure 4). Normalized RMS error and R2 for the predicted trial was 5.8% and 0.95 respectively. The averageR2 and normalized RMS errors ranged between 0.87 – 0.92, and 9.7% – 14.7% respectively.

Figure 4
The bold solid line is the ankle moment during stance for 1 subject computed using inverse dynamics (i.e., measured). The thin solid line depicts the model predicted moment using the subject’s EMG and joint kinematics. Despite the variability ...

DISCUSSION

The objectives of this paper were to describe in detail our EMG-driven model and discuss how it is implemented to predict joint moments and muscle forces. The second goal of this paper was to evaluate the model’s ability to predict joint moments in patients following stroke. Joint moments and muscle activation patterns are abnormal in this population and therefore present a challenging paradigm with which to test our EMG-driven model. Model results were verified by comparing the predicted and inverse dynamic joint moments during novel walking trials.

Overall, the model did predict the ankle joint moment for the four subjects in this study, despite the variability in muscle activation patterns and joint moments between trials. This may be associated with dysfunction caused by spasticity, reduced coordinative control, muscle weakness, and impaired sensory-motor control [12, 32, 33]. Results for the four patients were comparable to data reported for healthy individuals. Bogey et al. [34] reported calibration results on the ankle of healthy subjects during walking with a joint moment R2 of 0.97 and a normalized RMS error of 4.8%, but they did not use the calibrated model to predict novel trials. Buchanan et al. [35] had success predicting ankle joint moments with an R2 of 0.94 and a normalized RMS error of approximately 9% for a healthy subject during gait. In another study, Manal and colleagues reported R2 values ranging between 0.87 and 0.95 for three healthy subjects during stance [36]. Taken together, these studies and our results suggest that EMG-driven modeling can be used to estimate muscle forces and joint kinetics in patients following stroke.

Lloyd and Besier used an EMG-driven model to estimate knee moments and muscle forces in healthy subjects during dynamic tasks [10]. The model was tuned using a variety of dynamic activities, including jogging and cutting trials. This was done to expose the model to a wide range of muscle activiation patterns and differing joint moment profiles. The model was able to predict the knee joint moment for a range of tasks with good accuracy. This is a more robust test of the model than simply tuning it to a walking trial and then predicting another walking trial for the same subject. For example, Kadaba et al. have shown that joint kinematics, kinetics and muscle activation patterns are similar within subject from trial to trial [37]. Thus, once a model is tuned to a walking trial for a healthy subject, the model will likely do well predicting another walking trial since the input data (i.e., EMG and kinematics) will be similar to the data used to tune the model.

Neurologically impaired subjects become fatigued easily and therefore testing them under a variety of experimental conditions is not feasible. In our study we tuned the model to a walking trial, and once tuned, predicted other walking trials for the subjects. This paradigm is a good test of the model because the EMG and moment data were variable from trial to trial, unlike data for healthy subjects.

Since our EMG-driven model can be used to predict novel trials, it may be possible to determine how muscle activation patterns for a stroke patient would have to be augmented to achieve a desired moment profile. These calculated corrective changes in muscle activation patterns may be used as reference to derive appropriate stimulation patterns, which can be applied in stroke patients’ gait training with the use of functional electrical stimulation.

There are several assumptions and limitations of this study. First, although model parameters are tuned for each subject, the model does not account for differences in musculoskeletal size. Second, hemiparetic muscles of post-stroke patients have been shown to have altered fiber type proportions [38] and fiber atrophy [39]. Such changes of muscle architecture may influence the muscle force-activation relationship of these patients. Since such changes are difficult to model, we included adjustable parameters, such as maximal muscle strength, optimal fiber length and tendon slack length. These parameters were adjusted to account for the impaired musculature of each subject, and it was verified that our calibrated models could be used to predict novel trials.

SUMMARY

In conclusion, our EMG-driven model did predict the ankle joint moment during stance for four patients following stroke. The ability to accurately predict the joint moment is a first step in estimating muscle forces and joint loads associated with pathological gait. We are continuing to refine this model for the study of stroke and other pathologies (e.g., osteoarthritis). The long term goal of this work is to apply such models and provide clinicians with insight regarding patient specific muscle function and joint loading.

Acknowledgments

This work supported, in part, by NIH R01-HD38582 and NIH R01-AR46386. We wish to thank Daniel Benoit, Wei-Li Hsu and Trisha Kesar for assistance with the data collection.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

References

1. Mulroy S, Gronley J, Weiss W, Newsam C, Perry J. Use of cluster analysis for gait pattern classification of patients in the early and late recovery phases following stroke. Gait Posture. 2003;18:114–25. [PubMed]
2. Neckel N, Pelliccio M, Nichols D, Hidler J. Quantification of functional weakness and abnormal synergy patterns in the lower limb of individuals with chronic stroke. J Neuroengineering Rehabil. 2006;3:17. [PMC free article] [PubMed]
3. Crowninshield RD, Brand RA. A physiologically based criterion of muscle force prediction in locomotion. J Biomech. 1981;14:793–801. [PubMed]
4. Collins JJ. The redundant nature of locomotor optimization laws. J Biomech. 1995;28:251–67. [PubMed]
5. Neptune RR, Kautz SA, Zajac FE. Contributions of the individual ankle plantar flexors to support, forward progression and swing initiation during walking. J Biomech. 2001;34:1387–98. [PubMed]
6. Thelen DG, Anderson FC. Using computed muscle control to generate forward dynamic simulations of human walking from experimental data. J Biomech. 2006;39:1107–15. [PubMed]
7. Buchanan TS, Shreeve DA. An evaluation of optimization techniques for the prediction of muscle activation patterns during isometric tasks. J Biomech Eng. 1996;118:565–74. [PubMed]
8. Buchanan TS, Moniz MJ, Dewald JP, Zev Rymer W. Estimation of muscle forces about the wrist joint during isometric tasks using an EMG coefficient method. J Biomech. 1993;26:547–60. [PubMed]
9. Manal K, Gonzalez RV, Lloyd DG, Buchanan TS. A real-time EMG-driven virtual arm. Comput Biol Med. 2002;32:25–36. [PubMed]
10. Lloyd DG, Besier TF. An EMG-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo. J Biomech. 2003;36:765–76. [PubMed]
11. Buchanan TS, Lloyd DG, Manal K, Besier TF. Neuromusculoskeletal Modeling: Estimation of Muscle Forces and Joint Moments and Movements From Measurements of Neural Command. J Appl Biomech. 2004;20:367–395. [PMC free article] [PubMed]
12. Den Otter AR, Geurts AC, Mulder T, Duysens J. Abnormalities in the temporal patterning of lower extremity muscle activity in hemiparetic gait. Gait Posture. 2007;25:342–52. [PubMed]
13. Delp SL, Loan JP, Hoy MG, Zajac FE, Topp EL, Rosen JM. An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans Biomed Eng. 1990;37:757–67. [PubMed]
14. Lloyd DG, Buchanan TS. A model of load sharing between muscles and soft tissues at the human knee during static tasks. J Biomech Eng. 1996;118:367–76. [PubMed]
15. Thelen DG, Schultz AB, Fassois SD, Ashton-Miller JA. Identification of dynamic myoelectric signal-to-force models during isometric lumbar muscle contractions. J Biomech. 1994;27:907–19. [PubMed]
16. Woods JJ, Bigland-Ritchie B. Linear and non-linear surface EMG/force relationships in human muscles. An anatomical/functional argument for the existence of both. Am J Phys Med. 1983;62:287–99. [PubMed]
17. Manal K, Buchanan TS. A one-parameter neural activation to muscle activation model: estimating isometric joint moments from electromyograms. J Biomech. 2003;36:1197–202. [PubMed]
18. Gordon AM, Huxley AF, Julian FJ. The variation in isometric tension with sarcomere length in vertebrate muscle fibres. J Physiol. 1966;184:170–92. [PubMed]
19. Hill AV. The heat of shortening and the dynamic constants of muscle. Proceedings of the Royal Society of London Series B-Biological Sciences. 1938;126:136–195.
20. Zajac FE. Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. Crit Rev Biomed Eng. 1989;17:359–411. [PubMed]
21. Epstein M, Herzog W. Theoretical models of skeletal muscle. Wiley; New York: 1998.
22. Schutte LM. PhD Thesis. Stanford University; 1992. Using Musculoskeletal Models to Explore Strategies for Improving Performance in Electrical Stimulation-Induced Leg Cycle Ergometry.
23. Schutte LM, Rodgers MM, Zajac FE, Glaser RM. Improving the efficacy of electrical stimulation induced leg cycle ergometry: an analysis based on a dynamic musculoskeletal model. IEEE Trans Rehabil Eng. 1993;1:109–125.
24. Scott SH, Winter DA. A comparison of three muscle pennation assumptions and their effect on isometric and isotonic force. J Biomech. 1991;24:163–7. [PubMed]
25. Guimaraes AC, Herzog W, Hulliger M, Zhang YT, Day S. Effects of muscle length on the EMG-force relationship of the cat soleus muscle studied using non-periodic stimulationof ventral root filaments. J Exp Biol. 1994;193:49–64. [PubMed]
26. Huijing PA. Important experimental factors for skeletal muscle modelling: non-linear changes of muscle length force characteristics as a function of degree of activity. Eur J Morphol. 1996;34:47–54. [PubMed]
27. Buchanan TS. Evidence that maximum muscle stress is not a constant: differences in specific tension in elbow flexors and extensors. Med Eng Phys. 1995;17:529–36. [PubMed]
28. Yamaguchi GT, Sawa AG-U, Moran DW, Fessler MJ, Winters JM. A survey of human musculotendonactuator parameters. In: Winters JM, Woo SL-Y, editors. Multiple Muscle Systems: Biomechanics and Movement Organization. Springer; New York: 1990. pp. 717–773.
29. Higginson JS, Neptune RR, Anderson FC. Simulated parallel annealing within a neighborhood for optimization of biomechanical systems. J Biomech. 2005;38:1938–42. [PubMed]
30. Hull ML, Hawkins D. Analysis of muscular work in multisegmental movements: application to cycling. In: Winters JM, Woo SL-Y, editors. Multiple muscle systems: biomechanics and movement organization. Springer; New York: 1990. pp. 621–638.
31. Perry J. Gait analysis: normal and pathological function. SLACK inc; Thorofare, N.J.: 1992.
32. Knutsson E. Gait control in hemiparesis. Scand J Rehabil Med. 1981;13:101–8. [PubMed]
33. Burridge JH, Wood DE, Taylor PN, McLellan DL. Indices to describe different muscle activation patterns, identified during treadmill walking, in people with spastic drop-foot. Med Eng Phys. 2001;23:427–34. [PubMed]
34. Bogey RA, Perry J, Gitter AJ. An EMG-to-force processing approach for determining ankle muscle forces during normal human gait. IEEE Trans Neural Syst Rehabil Eng. 2005;13:302–10. [PubMed]
35. Buchanan TS, Lloyd DG, Manal K, Besier TF. Estimation of muscle forces and joint moments using a forward-inverse dynamics model. Med Sci Sports Exerc. 2005;37:1911–6. [PubMed]
36. Manal K, Cohen S, Buchanan TS. Testing an EMG-Driven Model to Estimate Ankle Moments and Muscle Forces During Walking. 2nd International Ankle Symposium.2004.
37. Kadaba MP, Ramakrishnan HK, Wootten ME. Measurement of lower extremity kinematics during level walking. J Orthop Res. 1990;8:383–92. [PubMed]
38. Jakobsson F, Edstrom L, Grimby L, Thornell LE. Disuse of anterior tibial muscle during locomotion and increased proportion of type II fibres in hemiplegia. J Neurol Sci. 1991;105:49–56. [PubMed]
39. Hachisuka K, Umezu Y, Ogata H. Disuse muscle atrophy of lower limbs in hemiplegic patients. Arch Phys Med Rehabil. 1997;78:13–8. [PubMed]