This article has Open Peer Review reports available.
Joint torque variability and repeatability during cyclic flexion-extension of the elbow
© Ballaz et al. 2016
Received: 26 April 2015
Accepted: 1 April 2016
Published: 11 April 2016
Joint torques are generally of primary importance for clinicians to analyze the effect of a surgery and to obtain an indicator of functional capability to perform a motion. Given the current need to standardize the functional evaluation of the upper limb, the aim of this paper is to assess (1) the variability of the calculated maximal elbow joint torque during cyclic elbow flexion-extension movements and (2) participant test-retest repeatability in healthy young adults. Calculations were based on an existing non-invasive method including kinematic identification and inverse dynamics processes.
Twelve healthy young adults (male n = 6) performed 10 elbow flexion-extension movement carrying five different dumbbells (0, 1, 2, 3 and 4 kg) with several flexion-extension frequencies (½, 1/3, ¼ Hz) to evaluate peak elbow joint torques.
Whatever the condition, the variability coefficient of trial peak torques remained under 4 %. Bland and Altman plot also showed good test-retest, whatever the frequency conditions for the 0, 1, 2, and 3 kg conditions.
The good repeatability of the flexion-extension peak torques represents a key step to standardize the functional evaluation of the upper limb.
In many musculoskeletal diseases muscular weakness leads to functional disability and decreased quality of life. For therapists, it is important to assess and quantify muscle strength in order to choose the most appropriate treatment or to evaluate therapy effects [1, 2]. Joint torques are generally of primary importance for clinicians to analyze the effect of a surgery on symmetry and comfort, and to obtain an indicator of functional capability to perform a motion. Joint torques are very often analyzed in patients with osteoarthritis (e.g.: [3, 4]) or scoliosis (e.g.: [5, 6]). Especially at the elbow, the change in elbow torque is an indicator of incremental release of the brachioradialis insertion footprint, for surgeons performing open reduction or internal fixation of distal radius fractures . For physio/ergo-therapists, the elbow torque is an indicator of functional capability to perform a motion, e.g. in stroke patients, and a control variable for assistive devices developed for these patients . In the rehabilitation field, strength is assessed though the measurement of the maximal joint torque [9–11], which represents the resultant action of all muscles crossing the joint, but do not provide each muscle force contribution. Studies have shown the potential of musculoskeletal simulation tools to determine the contribution of each muscle crossing a joint during movement which was otherwise impractical or impossible to obtain experimentally . According to the clinical relevance and accuracy of the used method, such quantification would help clinicians to target the best therapeutic solution. Indeed, computational model could give the opportunity to predict the effect of the muscle property modifications on joint torque production . For example, the effect of antagonist muscle release (e.g.: spasticity treatment) on joint torque production could be anticipated.
The upper limb function is of utmost importance in improving the quality of life and enhancing functional independence. Especially, elbow flexion movement has been related to motor impairment and performance . Thus, accurate modeling of elbow muscle involvement could provide an interesting tool to better understand the movement limitation. Within this process of calculating the muscle forces, joint torque is an essential intermediate variable [15–17]. Moreover, precise and repeatable quantification of the upper limb joint torque is of major importance for numerous applications (e.g. [18–20]) including exoskeletons and interactive rehabilitation devices development (e.g. [18, 21]), the understanding of the mechanisms resulting in joint rigidity (e.g. [22, 23]), or the impact of joint co-contraction on joint constraint (e.g. [17, 24]).
However, it is not always obvious to obtain accurate joint torque results that could be usefully exploited in model [25–27]. Applied to human motion analysis, several parameters can be a source of error. The major problems are linked to the inverse dynamic solution repeatability, which is affected by both the data processing and the experimental procedure. More specifically, in a top down approach, inaccuracy in movement coordinate data, joint centre of rotation location, and kinematic data processing can impact on inverse dynamics solution . Indeed, using marker-based optical motion capture systems, marker misallocation and skin movement greatly influence joint centre localisation [28, 29]. The inertia parameters of the body segments can also influence inverse dynamic solution . Lastly, the estimate of internal efforts, i.e. joint torques and muscle forces, is particularly sensitive to accelerations [31–33]. As a result, kinematic data analysis is also of greatest importance and mainly impact inverse dynamic results. Riemer et al. found that these various inaccuracies can result in uncertainties of estimated joint torques ranging from 6 % to 232 % of the peak torque during gait. As suggested in the literature however, more accurate results can be obtained with corrected kinematics based on a kinematic identification process, named solidification procedure , compared to inverse dynamics using either raw kinematic data, smoothing or low-pass filtering [35–37].
Additionally, in order to use inverse dynamics to follow patient progress, the experimental procedure should (1) allow the spontaneous adaptation of the participant to perform the task (e.g.: minimally constraint movement) and (2) result in within-subject test-retest task repeatability, according to the kinematic and dynamic movement parameters used in the model.
In light of this information, we have developed a model which quantifies the contribution of muscles crossing the elbow joint during flexion and extension movements  in order to use it as a clinical tool. The model-based process includes two consecutive steps: a kinematic identification based on procedure of solidification , combined with inverse kinematics and an inverse dynamics process that provides the elbow joint net torque (for more details, see ). As a first step to test the accuracy of the model, the aim of the present study is (1) to assess the maximal elbow joint torque variability during cyclic elbow flexion extension movements and (2) to assess participant test-retest repeatability in healthy young adults.
Twelve healthy young adults (age = 23 ± 2; male n = 6) were included in the present study. Exclusion criteria were known musculoskeletal or orthopaedic pathology, on the basis of a questionnaire in participants. The study was approved by the Research Ethics Board of Ste-Justine Hospital, Montreal, Canada (Ethics case #3362). A written informed consent was obtained from participants. The research was in compliance with the Helsinki Declaration.
During experimentation, the participant sat on the chair. The participants were asked to perform 10 cycles of flexion-extension, following the rhythm of a given metronome, with and without dumbbells. Participants had to keep the shoulder and elbow joint as motionless as possible and the dumbbell axis horizontal. Participants were involved a few minutes with the dumbbells, before beginning the experiments.
The participants had to perform ten elbow flexion-extension movements with five different masses: 0, 1, 2, 3 and 4 kg, and at three motion frequencies, 0.5 Hz (i.e. a cycle in 2 seconds), 0.33 Hz (1 cycle in 3 seconds) and 0.25 Hz (1 cycle in 4 seconds). The order of the masses and frequencies was drawn randomly by the operator. Each male participant performed the whole experimental protocol twice in order to assess test and retest reproducibility of the joint torques. The retests were performed approximately 20 min after the tests, without removing the kinematic sensor.
Joint torque quantification process
The full model joint kinematics: the system is modeled as a constrained multibody system, using kinematic loops.
The joint kinematic identification: the joint coordinates q, velocities q̇ and accelerations q̈ are numerically determined by an optimization process that estimates the joint coordinates of the multibody model that best fit the experimental joint positions.
The inverse dynamics: using recursive Newton-Euler formalism, a 3D multibody model  provides the vector Q inv of joint forces and torques during movement as follows:
These equations were symbolically generated by the ROBOTRAN software , UCL, which allows us to straightforwardly interface these equations with any numerical process, such as the optimization process presented above and the time simulation of the trials.
Data was reported as mean (standard deviation) (SD). Normality of the distributions was determined using the Kolmogorov–Smirnov test. For each frequency (0.5, 0.3, and 0.25 Hz) and mass (0, 1, 2, 3, 4 kg), the peak torque variability within each trial was assessed by computing the coefficient of variation (%CV). The aim of this intra-test variability analysis was to enable to average the peak torques of each trial for the repeatability analysis. Paired t-tests were performed to detect possible systematic bias between test and retest trial. The possibility of heteroscedasticity was examined on the basis of the Pearson product-moment correlation (r) between the mean and the absolute differences. If the correlation coefficient was significant the data were considered as heteroscedastic . Bland and Altman plots and limits of agreement analyses were also calculated to determine whether peak torque is in agreement between tests and retest trial . This method (Bland & Altman, 1986) was extensively used in different research fields in test-retest studies [43–46] and is suitable in the case of the present study . A corrected standard deviation of differences for repeated measurements, SD corrected = √(2●SD2), was used based on Bland and Altman (1986) . Statistical analysis was performed using SPSS 17.0 (IBM, Chicago, USA).
In each condition, the peak torque values were normally distributed (Kolmogorov–Smirnov test, p > 0.05).
Peak torque coefficients of variation within trial
Peak torque CV (n = 12)
Test conditions mass (Kg)/Frequency (Hz)
95 % confidence interval
This study showed that the data processing and the experimental procedure implemented in the present study resulted in a low within-trial variability, i.e. a low variability inside each trial, and a good within-participant test-retest repeatability, i.e. a good repeatability between tests of the same participant, of the elbow peak torque in typically developing young adults. As shown by the limit of agreements, expressed as a percentage of the averaged peak torque, the result repeatability was equivalent whatever the frequency, amongst 0.25, 0.33, and 0.5 Hz, or the load, amongst 0 1, 2, and 3 kg, imposed during the movement.
This study highlighted that the 4 kg resulted in a more important variability compared to the lower masses. Based on this observation, it can be assumed that increasing the mass higher than 4 kg would result in a more important variability that would not be appropriated to evaluate the joint torques. On the contrary, using lower masses, such as 0 kg, are recommended for the good repeatability, and certainly do not imply fatigue, especially in female participants.
In summary, to evaluate muscle efforts in the rehabilitation field, the repeatability of the model at low frequencies and with light loads was a key result. In patients with neurological disorder, muscular strength and movement velocity is potentially very low depending on their functional capacity. As supported by the Bland and Altman analysis (Fig. 3), at low frequency (0.25Hz) the limit of agreement represented 9.1 % of the averaged peak torque, and considering the condition without dumbbells, the limit of agreement represented 9.3 % of the averaged peak torque. Even if the literature still has no consensus on the clinically important difference in elbow torque for humans, because this torque relates to each joint and each motion, Laitenberger et al. (2015)  reported an elbow torque variability up to 24 % in healthy subjects, which confirms that the obtained repeatability of 8.5 % when all test conditions are viewed together (Fig. 2) is relevant compared to the magnitude of this measurement. As described earlier, kinematic data processing, marker misallocation and skin movement could greatly influence joint centre localisation [28, 29] and in turn greatly impact inverse dynamic solution repeatability . Riemer et al. found that these various inaccuracies can result in uncertainties of estimated joint torques ranging from 6 % to 232 % of the peak torque during gait. The methodology used in the present study in terms of kinematic data processing, based on solidification procedure , was adequate to result in a good within-participant test-retest repeatability. At the same time, these results showed that the device used (Fig. 4) was adequate to obtain a repeatable elbow flexion-extension maximal torque.
Several limits were inherent with this study. First the repeatability of the data processing and the experimental procedure was tested with a limited number of participants. Nevertheless, many conditions were tested (frequency*mass), resulting in a test-retest repeatability analysis based on 90 trials. The test-retest repeatability analysis was performed only in male because fatigue could be more present in female compared to male participants. Secondly, the present study included healthy participants, the repeatability of the data processing and the experimental procedure implemented in the present study should be tested for each targeted disease. Thirdly, the repeatability of the model has been tested without removing the maker. A Further study is required to test the reproducibility with markers replacement because markers location could impact the inverse dynamic solution. Fourthly, the gender effect on repeatability was not studied and could be a further perspective. Fifthly, the trials were randomised and a 3 min. rest period was allocated between the trials, as recommended by Kollmitzer et al. (1999)  to avoid the muscle fatigue effect in the context of a sub-maximal effort of the upper limb. However, it is never excluded that either a fatigue or a time effect may have influenced the CV results. Especially concerning the higher CV reported for the 1 kg condition, we believe that this results from the method variation that might also be seen in other conditions. Nevertheless, the CV reported for the 1 kg condition remains low, even if it represents twice the CV reported in the other conditions. Sixthly, a method with a good reproducibility does not necessarily guarantee an accurate estimation of joint torques. Reproducibility is a necessary feature that is complementary to the accuracy, guaranteeing that the results will be similar for any trial. Let us remind that it is still not possible today to check the accuracy of the joint torques in a non-invasive way, even if this information is of primary importance for clinicians to analyze the effect of a surgery and to obtain an indicator of functional capability. Being aware of these limitations, the incentive of this paper was to analyse the quality of our joint torque quantification. The present inverse dynamical model of the human body, necessarily preceded by a kinematic identification of the model configurations, is proposed as a satisfying method to estimate the joint efforts in dynamical context. This problem being deterministic, Q inv becomes a sufficiently accurate result that can be exploited as a reference for the optimization process that attempts to solve the muscle force redundancy. These results represent the first step leading to the development of an accurate assessment of elbow muscle strength in clinical environment. The ability of giving accurate elbow joint net torques during motion, without requiring an important computational cost, is the main benefit of this method. Based on these results, multibody model refinement and clinical analysis will be implemented in further studies.
The aim of this study was to assess the peak torque elbow variability and repeatability. Whatever the flexion-extension movement conditions imposed, within-trial peak torque variability was low and within-participant test-retest repeatability of the elbow joint torques resulted in good agreement. This method is promising for potential clinical applications and can be used as a basis for further comparison between efforts quantification methods or refined multibody models in the human body during motion.
Availability of data and materials
The authors’ Research Ethics Board did not allow to publicly share the data and materials of this study. However, these data and materials will be available upon request to the corresponding author and in accordance with the Research Ethics Board.
The authors wish to thank MÉDITIS program supported by FONCER-CSNG for financial support.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Damiano DL, Abel MF, Pannunzio M, Romano JP. Interrelationships of strength and gait before and after hamstrings lengthening. J Pediatr Orthop. 1999;19:352–8.PubMedGoogle Scholar
- Hinderer KA, Hinderer SR. Muscle strength development and assessment in children and adolescents, in: H.-R. K. (Ed.), Muscle strength. Edinburgh: Churchill Livingstone 1993;p.93-140Google Scholar
- Schache AG, Fregly BJ, Crossley KM, Hinman RS, Pandy MG. The effect of gait modification on the external knee adduction moment is reference frame dependent. Clin Biomech (Bristol, Avon). 2008;23(5):601–8.View ArticleGoogle Scholar
- Schwachmeyer V, Damm P, Bender A, Dymke J, Graichen F, Bergmann G. In vivo hip joint loading during post-operative physiotherapeutic exercises. PLoS One. 2013;8(10), e77807.View ArticlePubMedPubMed CentralGoogle Scholar
- D'Amico M, D'Amico G, Paniccia M, Roncoletta P, Vallasciani M. An integrated procedure for spine and full skeleton multi-sensor biomechanical analysis & averaging in posture gait and cyclic movement tasks. Stud Health Technol Inform. 2010;158:118–26.PubMedGoogle Scholar
- Raison M, Ballaz L, Detrembleur C, Mahaudens P, Lebleu J, Fisette P, Mousny M. Lombo-sacral joint efforts during gait: comparison between healthy and scoliotic subjects. Stud Health Technol Inform. 2012;176:113–6.PubMedGoogle Scholar
- Tirrell TF, Franko OI, Bhola S, Hentzen ER, Abrams RA, Lieber RL. Functional consequence of distal brachioradialis tendon release: a biomechanical study. J Hand Surg Am. 2013;38(5):920–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Cheng HS, Ju MS, Lin CC. Improving elbow torque output of stroke patients with assistive torque controlled by EMG signals. J Biomech Eng. 2003;125(6):881–6.View ArticlePubMedGoogle Scholar
- Florence JM, Pandya S, King WM, Robison JD, Baty J, Miller JP, Schierbecker J, Signore LC. Intrarater reliability of manual muscle test (Medical Research Council scale) grades in Duchenne's muscular dystrophy. Phys Ther. 1992;72:115–22.View ArticlePubMedGoogle Scholar
- Cuthbert SC, Goodheart Jr GJ. On the reliability and validity of manual muscle testing: a literature review. Chiropr Osteopat. 2007;15:4.View ArticlePubMedPubMed CentralGoogle Scholar
- Stark T, Walker B, Phillips JK, Fejer R, Beck R. Hand-held dynamometry correlation with the gold standard isokinetic dynamometry: a systematic review. PM R. 2011;3:472–9.View ArticlePubMedGoogle Scholar
- Bessonet G, Sardain P, Chéssé S. Optimal motion synthesis – dynamic modeling and numerical solving aspects. Multibody System Dynamics. 2002;8:257–78.View ArticleGoogle Scholar
- Hoy MG, Zajac FE, Gordon ME. A musculoskeletal model of the human lower extremity: the effect of muscle, tendon, and moment arm on the moment-angle relationship of musculotendon actuators at the hip, knee, and ankle. J Biomech. 1990;23(2):157–69.View ArticlePubMedGoogle Scholar
- Massie CL, Fritz S, Malcolm MP. Elbow extension predicts motor impairment and performance after stroke. Rehabil Res Pract. 2011;381978.Google Scholar
- Amarantini D, Martin L. A method to combine numerical optimization and EMG data for the estimation of joint moments under dynamic conditions. J Biomech. 2004;37:1393–404.View ArticlePubMedGoogle Scholar
- De Groote F, Pipeleers G, Jonkers I, Demeulenaere B, Patten C, Swevers J, De Schutter J. A physiology based inverse dynamic analysis of human gait: potential and perspectives. Comput Methods Biomech Biomed Engin. 2009;12:563–74.View ArticlePubMedGoogle Scholar
- Raison M, Detrembleur C, Fisette P, Samin JC. Assessment of antagonistic muscle forces during forearm flexion/extension. Comput Methods Appl Sci. 2011;23:215–38.View ArticleGoogle Scholar
- Gilliaux M, Lejeune T, Detrembleur C, Sapin J, Dehez B, Stoquart G. A robotic device as a sensitive quantitative tool to assess upper limb impairments in stroke patients: a preliminary prospective cohort study. J Rehabil Med. 2012;44(3):210–7.View ArticlePubMedGoogle Scholar
- Zampagni ML, Casino D, Zaffagnini S, Visani A, Marcacci M. Trend of the carrying angle during flexion-extension of the elbow joint: a pilot study. Orthopedics. 2008;31(1):76.View ArticlePubMedGoogle Scholar
- Lin JH, McGorry RW, Banks JJ. Exposures and physiological responses in power tool operations: fastening vs. unfastening threaded hardware. J Occup Environ Hyg. 2010;7(5):290–7.View ArticlePubMedGoogle Scholar
- Malosio M, Pedrocchi N, Vincenti F. Tosatti LM. Rehabilitation Robotics. 2011;5975393.Google Scholar
- Endo T, Okuno R, Yokoe M, Akazawa K, Sakoda S. A novel method for systematic analysis of rigidity in Parkinson's disease. Mov Disord. 2009;24:2218–24.View ArticlePubMedGoogle Scholar
- Park BK, Kwon Y, Kim JW, Lee JH, Eom GM, Koh SB, Jun JH, Hong J. Analysis of viscoelastic properties of wrist joint for quantification of parkinsonian rigidity. IEEE Trans Neural Syst Rehabil Eng. 2011;19(2):167–76.View ArticlePubMedGoogle Scholar
- Gross R, Leboeuf F, Rémy-Néris O, Perrouin-Verbe B. Unstable gait due to spasticity of the rectus femoris: gait analysis and motor nerve block. Ann Phys Rehabil Med. 2012;55(9-10):609–22.View ArticlePubMedGoogle Scholar
- Riemer R, Hsiao-Wecksler ET, Zhang X. Uncertainties in inverse dynamics solutions: a comprehensive analysis and an application to gait. Gait Posture. 2008;27:578–88.View ArticlePubMedGoogle Scholar
- Vaughan CL, Davis BL, O'Connor JC. Dynamics of Human Gait. 2nd ed. Cape Town, South Africa: Kiboho; 1992. p. 137.Google Scholar
- Winter DA. Biomechanics and motor control of human movement. Thirdth ed. Hoboken, New-Jersey: John Wiley & Sons, Inc; 2005. p. 720.Google Scholar
- Kuo YL, Tully EA, Galea MP. Skin movement errors in measurement of sagittal lumbar and hip angles in young and elderly subjects. Gait Posture. 2008;27:264–70.View ArticlePubMedGoogle Scholar
- Stagni R, Leardini A, Cappozzo A, Grazia Benedetti M, Cappello A. Effects of hip joint centre mislocation on gait analysis results. J Biomech. 2000;33:1479–87.View ArticlePubMedGoogle Scholar
- Rao G, Amarantini D, Berton E, Favier D. Influence of body segments' parameters estimation models on inverse dynamics solutions during gait. J Biomech. 2006;39:1531–6.View ArticlePubMedGoogle Scholar
- Cahouet V, Luc M, David A. Static optimal estimation of joint accelerations for inverse dynamics problem solution. J Biomech. 2002;35:1507–13.View ArticlePubMedGoogle Scholar
- Cappozzo A, Leo T, Pedotti A. A general computing method for the analysis of human locomotion. J Biomech. 1975;8:307–20.View ArticlePubMedGoogle Scholar
- Challis JH, Kerwin DG. Quantification of the uncertainties in resultant joint moments computed in a dynamic activity. J Sports Sci. 1996;14:219–31.View ArticlePubMedGoogle Scholar
- Cheze L, Fregly BJ, Dimnet J. A solidification procedure to facilitate kinematic analyses based on video system data. J Biomech. 1995;28:879–84.View ArticlePubMedGoogle Scholar
- Cappozzo A, Catani F, Leardini A, Benedetti MG, Croce UD. Position and orientation in space of bones during movement: experimental artefacts. Clin Biomech (Bristol, Avon). 1996;11:90–100.View ArticleGoogle Scholar
- Chiari L, Della Croce U, Leardini A, Cappozzo A. Human movement analysis using stereophotogrammetry. Part 2: instrumental errors. Gait Posture. 2005;21:197–211.View ArticlePubMedGoogle Scholar
- Lu TW, O'Connor JJ. Bone position estimation from skin marker co-ordinates using global optimisation with joint constraints. J Biomech. 1999;32:129–34.View ArticlePubMedGoogle Scholar
- Kim K, Song WK, Lee J, Lee HY, Park DS, Ko BW, Kim J. Kinematic analysis of upper extremity movement during drinking in hemiplegic subjects. Clin Biomech (Bristol, Avon). 2014;29(3):248–56.View ArticleGoogle Scholar
- de Leva P. Adjustments to Zatsiorsky-Seluyanov's segment inertia parameters. J Biomech. 1996;29(9):1223–30.View ArticlePubMedGoogle Scholar
- Samin JC, Fisette P. Symbolic modeling of multibody systems. Netherlands: Kluwer Academic Publishers; 2003. p. 479.View ArticleGoogle Scholar
- Atkinson G, Nevill AM. Statistical methods for assessing measurement error (reliability) in variables relevant to sports medicine. Sports Med. 1998;26(4):217–38.View ArticlePubMedGoogle Scholar
- Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.View ArticlePubMedGoogle Scholar
- Høyer E, Opheim A, Strand LI, Moe-Nilssen R. Temporal and spatial gait parameters in patients dependent on walking assistance after stroke: reliability and agreement between simple and advanced methods of assessment. Gait Posture. 2014;40(1):101–6.View ArticlePubMedGoogle Scholar
- Thompson P, Beath T, Bell J, Jacobson G, Phair T, Salbach NM, Wright FV. Test-retest reliability of the 10-metre fast walk test and 6-minute walk test in ambulatory school-aged children withcerebral palsy. Dev Med Child Neurol. 2008;50(5):370–6.View ArticlePubMedGoogle Scholar
- Birmingham TB, Hunt MA, Jones IC, Jenkyn TR, Giffin JR. Test-retest reliability of the peak knee adduction moment during walking in patients with medial compartment knee osteoarthritis. Arthritis Rheum. 2007;57(6):1012–7.View ArticlePubMedGoogle Scholar
- de Zwart BC, Frings-Dresen MH, van Duivenbooden JC. Test-retest reliability of the Work Ability Index questionnaire. Occup Med (Lond). 2002;52(4):177–81.View ArticleGoogle Scholar
- Laitenberger M, Raison M, Périé D, Begon M. Refinement of the upper limb joint kinematics and dynamics using a subject-specific closed-loop forearm model. Multibody System Dynamics. 2015;33(4):413–38.View ArticleGoogle Scholar
- Kollmitzer J, Ebenbichler GR, Kopf A. Reliability of surface electromyographic measurements. Clin Neurophysiol. 1999;110(4):725–34.View ArticlePubMedGoogle Scholar