Next Article in Journal
Development of an Integrated Evaluation System for a Stretchable Strain Sensor
Next Article in Special Issue
A Review of Classification Techniques of EMG Signals during Isotonic and Isometric Contractions
Previous Article in Journal
Modeling and Compensating Temperature-Dependent Non-Uniformity Noise in IR Microbolometer Cameras
Previous Article in Special Issue
Multi-Section Sensing and Vibrotactile Perception for Walking Guide of Visually Impaired Person
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression

by
Morufu Olusola Ibitoye
1,2,*,
Nur Azah Hamzaid
1,*,
Ahmad Khairi Abdul Wahab
1,
Nazirah Hasnan
3,
Sunday Olusanya Olatunji
4 and
Glen M. Davis
1,5
1
Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur 50603, Malaysia
2
Department of Biomedical Engineering, Faculty of Engineering and Technology, University of Ilorin, P.M.B 1515, Ilorin 24003, Kwara State, Nigeria
3
Department of Rehabilitation Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur 50603, Malaysia
4
Computer Science Department, College of Computer Science & Information Technology, University of Dammam, Dammam 34212, Saudi Arabia
5
Clinical Exercise and Rehabilitation Unit, Discipline of Exercise and Sports Sciences, Faculty of Health Sciences, The University of Sydney, Sydney, 2006 NSW, Australia
*
Authors to whom correspondence should be addressed.
Sensors 2016, 16(7), 1115; https://doi.org/10.3390/s16071115
Submission received: 14 March 2016 / Revised: 6 June 2016 / Accepted: 8 June 2016 / Published: 19 July 2016
(This article belongs to the Special Issue Sensing Technology for Healthcare System)

Abstract

:
The difficulty of real-time muscle force or joint torque estimation during neuromuscular electrical stimulation (NMES) in physical therapy and exercise science has motivated recent research interest in torque estimation from other muscle characteristics. This study investigated the accuracy of a computational intelligence technique for estimating NMES-evoked knee extension torque based on the Mechanomyographic signals (MMG) of contracting muscles that were recorded from eight healthy males. Simulation of the knee torque was modelled via Support Vector Regression (SVR) due to its good generalization ability in related fields. Inputs to the proposed model were MMG amplitude characteristics, the level of electrical stimulation or contraction intensity, and knee angle. Gaussian kernel function, as well as its optimal parameters were identified with the best performance measure and were applied as the SVR kernel function to build an effective knee torque estimation model. To train and test the model, the data were partitioned into training (70%) and testing (30%) subsets, respectively. The SVR estimation accuracy, based on the coefficient of determination (R2) between the actual and the estimated torque values was up to 94% and 89% during the training and testing cases, with root mean square errors (RMSE) of 9.48 and 12.95, respectively. The knee torque estimations obtained using SVR modelling agreed well with the experimental data from an isokinetic dynamometer. These findings support the realization of a closed-loop NMES system for functional tasks using MMG as the feedback signal source and an SVR algorithm for joint torque estimation.

1. Introduction

The magnitude of the muscle force or joint torque generated during neuromuscular electrical stimulation-evoked contractions has been used as a marker of physical performance in healthy individuals [1,2], as well as a benchmark of functional recovery in individuals with neurological conditions [3,4]. To optimize neuromuscular electrical stimulation (NMES) technology in therapeutic and functional applications, real-time information about the generated muscle force or joint torque, of the controlled limb, is vital [3,5]. Such information is required; (i) to automate the neuromuscular stimulation characteristics based on the muscle state during the onset of fatigue; and (ii) to modulate muscle forces based on the requirements of the task (therapeutic or functional) to be performed, for example during sit-to-stand and sustained standing perturbations. However, joint torque is often impractical or impossible to quantify directly during real-time application of NMES [5]. Estimation of joint torque from readily available muscle characteristics (e.g., biopotentials of nerve and/or muscle activation), particularly, from physical sensors has recently become both viable and attractive [5].
One such neuromuscular biopotential is the Mechanomyographic signal (MMG), which quantifies the mechanical equivalent of an electromyographic output generated during muscle contractions [6,7]. The signal originates from the skeletal muscle contractions due principally to the shortening of the muscle fiber length and increase in its diameter [7,8]. The activation of muscle fibres and their dimensional changes during muscle contraction creates pressure waves that can be detected on the skin surface and translated into an acceleration obtained by physical sensors, such as an accelerometer [9]. The signal can represent a proxy for neuromuscular contractions [10] and has gained recent popularity due to its close relationship with muscle force [11]. Specifically, the signal is directly related to the two main force-generating mechanisms of human skeletal muscle—magnitude and pattern of motor unit recruitment and their firing rates/frequency [12,13].
Moreover, due to the convenience of MMG signal collection, its insusceptibility to skin impedance [14], flexibility of its sensing technology [15,16], and immunity from electrical stimulation artifacts associated with NMES [17], the signal has been successfully used to classify muscle activity for specific application in controlling prostheses [18], and as a control signal for muscle machine interfaces [16,19]. In addition, during NMES-evoked muscle contractions, MMG has been used to track muscle fatigue in healthy volunteers [20]. Thus, the signal may be used to estimate joint torque during voluntary and/or NMES-evoked muscle contractions [7]. However, relating MMG signals as a direct proxy for NMES-evoked muscle effort/force can be practically challenging due to the complexity and diversity of the recruitment of muscle’s motor units (MUs) [6,18]. Accordingly, the application of computational intelligence techniques for quantification of joint torque from MMG signals has been proposed through statistical predictive modelling, and then validated during voluntary contractions [13,21].
The use of machine-learning techniques has recently shown promise, subverting the dual problems of non-linearity and non-stationarity in estimation, prediction and classification tasks. For example, Youn and Kim [13,22] used an artificial neural network model to estimate elbow flexion force from MMG during voluntary isometric contractions. The investigators obtained an estimation accuracy of up to 0.892 [13] and 0.883 [22] in terms of cross-correlation coefficient, and concluded that their model is subject dependent, while suggesting the future application of other machine learning techniques including Support Vector Regression (SVR) to improve the estimation accuracy of the model [22]. However, due to the advancement in the field of signal processing, several other computational intelligence statistical regression techniques have been proposed with SVR yielding a good predictive and estimation accuracy, with often low Root Mean Square Errors (RMSE) [23] and outstanding performance [21]. Being a category of support vector machine learning technique, SVR is based on the principles of computational intelligence that is built on the kernel method (that maps data into higher dimensional space where the training sample may be linearly separable to facilitate linear regression analysis) [24]. SVR algorithms take into account the error approximation to a dataset with the ability to adapt and improve the estimation capability of a model [23], particularly when the model is used to evaluate an additional dataset for the purpose of generalization [25,26]. Moreover, SVR is robust in handling multivariate processes and offsets the limitation of traditional regression methods [27]—which cannot solve problems with high dimensional input dataset [24]. Additionally, the SVR modelling only involves a solution to a “convex optimization problems”, and unlike ANN model, it is not influenced by the “local minimal problem” [28]. Thus, the SVR algorithms could be used to build a generalized model and well suited for regression tasks [24]. Based on this strength, the technique has been successfully deployed in several fields of applications including physical therapy and exercise science during voluntary muscle activation [21], medical diagnosis [29], and a host of other related fields. However, to our knowledge, SVR modelling has not been previously used to construct a joint torque estimation model, particularly, during electrically stimulated contraction.
The purpose of this study was, therefore, to use SVR modelling to predict knee extensor joint torques from MMG signal characteristics during NMES-evoked incremental muscle contraction intensities. Since it has been suggested [13] that a combination of muscle contraction signals and related characteristics could compliment the estimation accuracy of joint torques, three input parameters (related to muscle contractions) to the SVR model were chosen (MMG signals, level of electrical stimulation or contraction intensity, and knee angle) to estimate knee torque accurately. This information is particularly applicable to research areas where a real-time proxy of muscle force is sought.

2. Materials and Methods

2.1. Experimental Protocol

To validate the performance of the proposed SVR model, a calibrated commercial dynamometer (System 4; Biodex Medical System, Shirley, NY, USA) was used to record isometric knee torques produced by NMES-evoked contractions of the knee extensors (Figure 1). Eight healthy male volunteers aged 23.4 (1.3) year (mean (SD)), body mass 70.4 (5.9) kg and height (1.72 (0.05)) m participated in this experiment. All were in good physical condition and were duly informed about the study protocol prior to giving their written informed consent. The study was approved by the University of Malaya Medical Ethics Committee (Approval No: 1003.14 (1)). As portrayed in Figure 1, the participants were set-up, as has been previously described by Brown and Weir [30] for voluntary isometric knee torque measurements. The dynamometer seat was adjusted so that each participant's lateral femoral condyle was aligned with the axle of the dynamometer [31]. To ensure consistency of body position and dynamometer lever arm, for subsequent trials, notes were taken of each participant’s relevant anatomical positions.

2.1.1. NMES-Evoked Muscle Contractions and Knee Torque Measurements

A familiarization session (mimicking the actual test) preceded data collection to familiarize the participants to the study protocol and to habituate them to NMES-evoked knee extensors muscle contractions of maximally tolerable intensity. Thereafter, neuromuscular stimulation of square-wave pulses at 30 Hz frequency and 400 μs pulse duration, and incremental current amplitude from 20 mA to 80 mA (in 10 mA increments; i.e., seven different intensities of NMES or trial levels) was administered to elicit isometric torque of the knee extensors lasting 4 s [32]. Stimulation pulses were delivered through a commercially available computer-controlled neurostimulator (RehaStimTM, Hasomed GmbH, Magdeburg, Germany) using 9 × 15 cm2 self-adhesive electrodes (Hasomed GmbH, D 39114, Magdeburg, Germany) on the dominant leg [33]. To preclude voluntary effort, the participants were carefully instructed not to assist or resist NMES-evoked muscle contractions. A similar stimulation protocol has been used for strength training with tolerable discomfort [34] and without eliciting rapid muscle fatigue [35]. During each trial, the NMES-evoked torque at maximum stimulation intensity (80 mA) was taken as the NMES-evoked peak torque (PT). The PT value was used to normalize the submaximal contraction levels across participants’ data. The adopted stimulation electrode position has been recommended by Levin and colleagues [36]—the anode electrode placed at “~5 cm proximal position to the patella and the cathode electrode at ~8 cm distal to the inguinal area over the rectus femoris (RF) muscle belly near the expected location of the motor points” (Figure 1). In order to accommodate the effect of joint angle on joint torque [34,37], the experiment was conducted at three different randomized knee angles: 30°, 60°, and 90° (where 0° represented full knee extension). A duration of 48 h was allowed between each angle position, and there was a 10 min recovery between each trial to minimize potential muscle fatigue.

2.1.2. MMG Acquisition and Processing

Simultaneous with the NMES-evoked muscle contraction, MMG signals were collected using an accelerometer-based sensor (Sonostics BPS-II VMG transducer, sensitivity 30 V/g). The sensor was attached directly on the muscle belly (i.e., at the midpoint between the inguinal crease and the superior border of the patella [38] (Figure 1 and Figure 2) by means of double-sided adhesive tapes (3M Center St. Paul, MN, USA). The MMG signals were collected from the RF muscle as a simple representation of the knee extensors and a major contributor to the NMES-evoked knee torque production [39]. The signals were collected at 2 kHz sampling frequency and were digitally band-pass filtered at 20–200 Hz, amplified and stored by AcqKnowledge data acquisition and analysis software (MP150, BIOPAC Systems Inc., Santa Barbara, CA, USA) for offline analysis in the LabVIEW software environment (version 12.0, National Instruments, Austin, TX, USA) using custom written programs.
The peak torque values, MMG-root mean square (RMS) and peak to peak (PTP) amplitudes were obtained during NMES-evoked isometric contractions from 2 s epoch of the 4 s MMG and torque recordings [40] at each contraction level across the three joint angles. The selected 2 s epoch of the signals coincided with the middle position at which there was probable maximum muscle recruitment, without on-transients or off-transients of force rise at the beginning and the end of muscle contractions, respectively [40].
Thereafter, the MMG signals at each contraction level were normalized (by the equivalent value of the MMG signal at the highest stimulation intensity/contraction level (80 mA)) and fed into the proposed SVR model for training. Previous investigations [6,12] have validated the legitimacy of these MMG features for muscle force assessment, and, therefore, they were equally used as joint torque predictors in this study. The equations used to compute the MMG-RMS and MMG-PTP features [15] are as follows:
MMG RMS = 1 N i = 1 N x i 2   where   i = 1 N x i 2   = x 1 2 + x 2 2 + x 3 2 +   + x N 2 )
MMG PTP =   2 2   RMS
here, x i is the i th sample of the MMG signals and N represents the number of samples in the epoch considered.

2.2. Support Vector Regression Modelling Approach

SVR algorithm was proposed in this study because of its optimal predictive performance even with small dataset [41] and the ability to learn both linear and non-linear relationships between predictors and outcome (actual) variables [42]. Such relationships have been used in establishing a pattern whereby unknown outcomes could be predicted accurately [21,24]. Theoretically, SVR is derived from the statistical learning theory [43,44] and employs ε -insensitive loss function [24] which measures the flatness of the generated pattern as well as maximum allowable deviations of the targets from the predicted values for all given training datasets ( x 1 , y 1 ) , .......... , ( x k , y k ) with k number of samples [45]. However, a function used for the SVR analysis should not only approximate the training data adequately but also predicts accurately the value of y for the future data x [25]. Such a function, with w , x dot product in the space of R , is represented in linear form by Equation (3) for a set of training samples.
f ( x , α ) = w , x + b
where w R and b R
For the purpose of establishing the goal of SVR in ensuring the flatness of the Equation (3), small value of w is desired through minimization of the Euclidean norm w 2 [42] which makes the optimization problem of the regression to be represented by Equation (4):
minimize   1 2 w 2 s u b j e c t   t o { y i w , x i b ε w , x i + b y i ε }
Equation (4) holds on the assumption [46] that there exists a function that is capable of providing error which is less than ε for all training pairs of the dataset. The slack variables ( ξ i   and   ξ * i ) , which represent the upper and lower constraints on the system output, are often introduced in order to permit some errors that are associated with real life problems [44,46]. Therefore, Equation (4) is modified and presented as Equation (5).
minimize   1 2 w 2 + C i = 1 k ( ξ i + ξ * i ) s u b j e c t   t o { y i w , x i b ε + ξ i w , x i + b y i ε + ξ * i ξ i , ξ * i 0     for   all   i = 1 , 2 , ...... , k }
The optimization problem in Equation (5) is better solved, through the ε -insensitive loss function, by using Lagrangian multipliers ( η i , η i * , λ i   and   λ i * ) to transform the problem into dual space representation. Therefore, the Lagrangian for the Equation (5) is presented in Equation (6).
L = 1 2 w 2 + C i = 1 k ( ξ i + ξ i * ) i = 1 k λ i ( ε + ξ i y i + w , x i + b ) i = 1 k λ i * ( ε + ξ i * + y i w , x i b ) i = 1 k ( η i ξ i + η i * ξ i * )
It is easier to locate the saddle point of the Lagrangian function defined in Equation (6) by equating the partial derivatives of the Lagrangian ( with   respect   to   w , b , ξ i   and   ξ i * ) to zero in order to obtain the expressions presented in Equations (7)–(9):
w = i = 1 k ( λ i * λ i ) . x i
η i = C λ i
η * i = C λ i *
The optimization equation is maximized by substituting Equations (7)–(9) in Equation (6) to arrive at Equation (10):
maximize   1 2 i = 1 k j = 1 k ( λ i * λ i ) ( λ j * λ j ) ( x j . x i ) ε i = 1 k ( λ i * + λ i ) + i = 1 k y i ( λ i * λ i ) s u b j e c t   to   i = 1 k ( λ i * λ i ) = 0 , 0 λ * i   and   λ i C
The solutions ( λ i *   and   λ i ) obtained from Equation (10) are substituted in Equation (3) and presented in Equation (11):
f ( x , α ) = i = 1 k ( λ i * λ i ) x i , x + b
However, since the concept of kernel function through ‘’kernel tricks’’ allows SVR to solve non-linear problems by mapping the original non-linear data into higher dimensional feature space where a linear model could be constructed [47], a proper selection of kernel function allows optimization of SVR performance [47]. The regression function in feature space, after inserting the kernel function K x i , x , could be written as presented in Equation (12).
f ( x , α ) = i = 1 k ( λ i * λ i ) K x i , x + b
Kernel functions help in transforming datasets into hyperplane [47]. The variables of the kernel function determine the structure of high dimensional feature space which controls the complexity of the final solution. As applied in this study, Equations (13)–(16) describe several kernel functions that are obtainable in the literature [48] which include Polynomial, Linear, Gaussian (radial basis function (RBF)) and Sigmoid.
K ( x i , x j ) = ( x i x j + 1 ) d
K ( x i , x j ) = x i T . x j
K ( x i , x j ) = exp ( γ x i x j d )
K ( x i , x j ) = tanh ( γ x i T x j + r )
where γ , r , a n d   d are kernel parameters and, x i and x j represent vectors in the input space—vectors of features computed from training or test subset [23].

2.2.1. Model Development

MATLAB software environment (Version 12, The MathWorks, Inc., Natick, MA, USA) using SVR coding was used for the computational aspect of this research work. Prior to the use of the dataset, the dataset was partitioned into two components to adhere to the SVR modelling approach [23,49]—a machine-learning “training” subset and a “testing” subset in a ratio of 7:3, via stratified sampling to ensure effective random partitioning [50]. Specifically, 70% of the dataset was used for training and the remaining 30% was used for testing the SVR model via test-set cross-validation method. This allowed a regression analysis to be performed on the training dataset while estimating the future generalization accuracy, of the model, on the remaining testing subset. For further detail on the working principle of the proposed SVR model, readers are referred to Vapnik et al. (1997) [24], Lin S-W et al. (2008) [47], Shamshirband et al. (2014) [23], and Akande et al. (2015) [49].

2.2.2. Optimal Parameters Search Approach

The accuracy of a SVR model is dependent on the model parameters’ selection [23]. However, due to the possibility of many different combinations of SVR parameters, it is often difficult to obtain optimal SVR parameters [51]. To solve this problem systematically, and in order to obtain possible optimized parameters of SVR for an accurate estimation, a hybrid optimization search technique, which has been recommended [52], was adopted and a test-set cross-validation technique was deployed [53]. The approach is as follows: for every partitioned training and testing subsets, the performance measures were noted for the SVR parameters values including the regularization factor C (bound on the Lagrangian multiplier), λ (conditioning parameter for quadratic programming (QP) methods), ε (epsilon) and η (kernel option) as well as the related kernel functions [49]. Thereafter, this computational step was repeated for every available SVR kernel function with an incremental step of the parameters’ values. The parameters’ optimal values and the kernel function associated with the best performance measure were identified. The search procedures are presented summarily in Figure 3.
A mathematical implementation [49] of the test-set cross-validation technique is as described in Algorithm 1 as follows: K i ( j ) was defined where K contains all the available kernel functions (and i , j and k are the indexes for the kernel functions) while i y , j y and k y represent the indexes for optimal kernel function. The total number of the available kernel function is represented by n i . The maximum values of C and η were assumed to be n j and n k , respectively. The recorded performance measures were stored in p f .
Algorithm 1. Optimal parameter search algorithm
I n i t i a l i z a t i o n ; i y = 0 , j y = 0 , k y = 0 , q x = 0
  f o r   i = 1 : n i
  f o r   j = 1 : n j p f = f ( K i ( j ) )
  f o r   k = 1 : n k
  p f = f ( K i ( j ) ) {Performance measure for the present parameters combination}
  i f   p f   i s   b e t t e r   t h a n   q x   t h e n   q x = p f
  i y = i ,   j y = j ,   k y = k   { s t o r i n g   t h e   i n d e x   o f   t h e   b e s t   p a r a m e t e r }
 end
 end
 end

2.2.3. Model Statistical Performance Criteria

To evaluate the performance of the proposed model, common measures of association, between the actual and the estimated values, were employed, including correlation coefficient (r) and coefficient of determination (R2) to quantify the “goodness of fit’’, and Root Mean Square Error (RMSE) to quantify the error of estimate. For further details on their mathematical formulae, readers are referred to Youn and Kim (2011) [22] and Olatunji et al. (2014) [54].

3. Results and Discussion

Table 2 describes the actual experimental dataset used in this study. The results of the statistical analysis of the dataset are presented in Table 3. The suitability and applicability of the chosen dataset are revealed from the mean, maximum value, median, standard deviation, and minimum value. The MMG-RMS, MMG-PTP, level of electrical stimulation or contraction intensity, and knee angle obtained experimentally were the input to the SVR model to estimate the knee torque. Results of performance measures obtained from the training subset and testing subset are as shown in Table 4.
To our knowledge, this is the first attempt to use a SVR modelling technique for NMES-evoked knee torque estimation from MMG signal. The outcomes of the developed SVR model (Table 4) indicated high correlation as well as low RMSE, and the model could, therefore, be adjudged as accurate. Moreover, high accuracy of the trained system as evident by the coefficient of determination (R2 = 94%), in predicting knee torque confirmed a reliable pattern between the predictors and the outcome which might be otherwise difficult to learn using linear regression.
During the training period of the model, the estimated torques were positively correlated with the actual values drawn from the experimental data (actual vs. predicted values) for both the training (Figure 4A) and testing (Figure 4B) subsets.
In addition, the cross-plots of the “training’’ subsets (actual vs. predicted values) as shown in Figure 5 also confirmed the high accuracy of the “training” subsets. However, since the actual performance of any model is better accessed by the testing outcome [55], the accuracy of the developed SVR model was tested using 30% of the available data samples (i.e., the reserved 30% that was not used in model development). It was interesting to note that, the model also performed satisfactorily during testing phase (R2 = 89%).
This high correlation indicated that the estimated knee torque by the SVR model was very close to the actual experimentally recorded joint torque (from an isokinetic dynamometer) for each data sample. For better visualization and understanding of the outcome of this study, the cross-plot of testing sets (actual vs. predicted values) has been portrayed in Figure 6. The level of accuracy in the testing phase (R2 = 94%) of the model development indicates that the model is stable, efficient and not over-fitted. This was based on the suggestion of Tay and Cao [56] that an overfitted model could perform excellently on the training set (r > 0.90) but will perform poorly on testing set [56]. Therefore, the developed SVR model in this study achieved a good performance for both training and testing sets. These results are comparable to that of Youn and Kim [22], where an artificial neural network model has been successfully used to estimate elbow force during voluntary contractions. Meanwhile, the potential of the SVR model for NMES-evoked joint torque estimation which has not been previously documented has also been demonstrated in our study.
Moreover, Figure 5 and Figure 6 portrayed the closeness of the predicted torque by the proposed SVR model to the actual experimental values. It could be noted that almost all the predicted points fit exactly on the experimental point or at least fits very closely to the target experimental point. Taken together, it could be inferred that the real-time knee torque information which is vital for the closed-loop implementation of NMES [3,5] in physical therapy and exercise science might be reliably estimated by our proposed method. Nevertheless, we acknowledge the following limitation in our study design: The performance of the developed model is limited to torque estimation during NMES-evoked isometric knee extension in healthy volunteers. In the future studies, we will verify the performance of the model using MMG signal and torque data from participants with neurological conditions. This will allow us to examine and improve the performance of the model, and to derive clinically relevant characteristics about the muscle force recruitment in clinical populations.

4. Conclusions

Based on its previous estimation accuracy in relevant fields, SVR modelling was used in this study through the integration of relevant variables to predict NMES-evoked knee torque. The model was developed through training and testing via test-set cross-validation technique with available dataset partitioned into training and testing subsets. Using the SVR methodology, the predicted knee torque was positively correlated with the actual values drawn from the experimental data for the training subset. Thereafter, to check the predictive ability of the model, the trained model was tested using the reserved testing subset that was not used in model development. The model performance was measured based on the correlation coefficient and RMSE. The outcomes from the developed SVR model showed an accurate prediction of the knee torque, characterized by high correlation coefficient—up to 0.97 and 0.94; and coefficient of determination—up to 94% and 89%, and low RMSE of 9.48 and 12.95, for the training and testing cases, respectively. These results, which have not been previously reported, indicated a close similarity between the estimated joint torque by the SVR model and the actual experimental data obtained from the laboratory experiment. Additionally, the present study has uniquely shown that a SVR model could estimate NMES-evoked knee torque, generated by a synchronous modulation of muscle fibres’ motor units [57], from MMG signal. Therefore, the good performance achieved in this study will motivate further studies in a similar direction to facilitate accurate estimations of joint torque using datasets from clinical populations—in which the NMES technology is more relevant, particularly among those with spinal cord injury. Moreover, since SVR models can be adapted for classification tasks [43], in the future, the developed model will be used to classify fresh and fatiguing muscle contractions of knee extensors, from MMG signals, during standing and ambulation tasks. Such models might offset the need to contend with the stimulation artifact [3,5] often encountered with the application of surface electromyographic signal as NMES feedback source.

Acknowledgments

The research was performed with a grant from the Ministry of Higher Education, Malaysia and the University of Malaya through HIR Grant No. UM.C/625/1/HIR/MOHE/ENG/39. The authors wish to gratefully acknowledge Dr. Khameel. B. Mustapha of The University of Nottingham, Malaysia Campus for valuable suggestions.

Author Contributions

M.O.I. conceived the experiment; M.O.I., N.A.H., and A.K.A. designed the experiment; M.O.I., N.A.H., and N.N. performed the experiment; M.O.I., S.O.O., A.K.A., and G.M.D. performed data analysis; N.N and A.K.A. contributed material and analysis tools; N.N. and G.M.D. provided the clinical guidance; A.K.A. and G.M.D. checked the technical soundness, clarity and the flow of the study. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NMESNeuromuscular Electrical Stimulation
MMGMechanomyography
SVRSupport Vector Regression
SVMSupport Vector Machine
RMSERoot Mean Square Error
PTPeak Torque
RFRectus Femoris
RMSRoot Mean Square
PTPPeak to Peak

References

  1. Brocherie, F.; Babault, N.; Cometti, G.; Maffiuletti, N.; Chatard, J.-C. Electrostimulation training effects on the physical performance of ice hockey players. Med. Sci. Sports Exerc. 2005, 37, 455–460. [Google Scholar] [CrossRef] [PubMed]
  2. Parker, M.G.; Bennett, M.J.; Hieb, M.A.; Hollar, A.C.; Roe, A.A. Strength response in human quadriceps femoris muscle during 2 neuromuscular electrical stimulation programs. J. Orthop. Sports Phy. Ther. 2003, 33, 719–726. [Google Scholar] [CrossRef] [PubMed]
  3. Braz, G.P.; Russold, M.; Davis, G.M. Functional electrical stimulation control of standing and stepping after spinal cord injury: A review of technical characteristics. Neuromodul. Technol. Neural Interface 2009, 12, 180–190. [Google Scholar] [CrossRef] [PubMed]
  4. Deley, G.; Denuziller, J.; Babault, N.; Taylor, J.A. Effects of electrical stimulation pattern on quadriceps isometric force and fatigue in individuals with spinal cord injury. Muscle Nerve 2015, 52, 260–264. [Google Scholar] [CrossRef] [PubMed]
  5. Popović, D.B. Advances in functional electrical stimulation (fes). J. Electromyogr. Kinesiol. 2014, 24, 795–802. [Google Scholar] [CrossRef] [PubMed]
  6. Orizio, C. Muscle sound: Bases for the introduction of a mechanomyographic signal in muscle studies. Crit. Rev. Biomed. Eng. 1993, 21, 201–243. [Google Scholar] [PubMed]
  7. Ibitoye, M.O.; Hamzaid, N.A.; Zuniga, J.M.; Abdul Wahab, A.K. Mechanomyography and muscle function assessment: A review of current state and prospects. Clin. Biomech. 2014, 29, 691–704. [Google Scholar] [CrossRef] [PubMed]
  8. Farina, D.; Li, X.; Madeleine, P. Motor unit acceleration maps and interference mechanomyographic distribution. J. Biomech. 2008, 41, 2843–2849. [Google Scholar] [CrossRef] [PubMed]
  9. Watakabe, M.; Itoh, Y.; Mita, K.; Akataki, K. Technical aspects of mechnomyography recording with piezoelectric contact sensor. Med. Biol. Eng. Comput. 1998, 36, 557–561. [Google Scholar] [CrossRef] [PubMed]
  10. Yuan-Ting, Z.; Frank, C.B.; Rangayyan, R.M.; Bell, G.D. A comparative study of simultaneous vibromyography and electromyography with active human quadriceps. IEEE Trans. Biomed. Eng. 1992, 39, 1045–1052. [Google Scholar]
  11. Matheson, G.O.; Maffey-Ward, L.; Mooney, M.; Ladly, K.; Fung, T.; Zhang, Y.T. Vibromyography as a quantitative measure of muscle force production. Scand. J. Rehabil. Med. 1997, 29, 29–35. [Google Scholar] [PubMed]
  12. Beck, T.W.; Housh, T.J.; Johnson, G.O.; Weir, J.P.; Cramer, J.T.; Coburn, J.W.; Malek, M.H. Mechanomyographic amplitude and mean power frequency versus torque relationships during isokinetic and isometric muscle actions of the biceps brachii. J. Electromyogr. Kinesiol. 2004, 14, 555–564. [Google Scholar] [CrossRef] [PubMed]
  13. Youn, W.; Kim, J. Estimation of elbow flexion force during isometric muscle contraction from mechanomyography and electromyography. Med. Biol. Eng. Comput. 2010, 48, 1149–1157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Alves, N.; Chau, T. The design and testing of a novel mechanomyogram-driven switch controlled by small eyebrow movements. J. NeuroEng. Rehabil. 2010, 7, 1–10. [Google Scholar] [CrossRef] [PubMed]
  15. Ibitoye, M.O.; Hamzaid, N.A.; Zuniga, J.M.; Hasnan, N.; Wahab, A.K.A. Mechanomyographic parameter extraction methods: An appraisal for clinical applications. Sensors 2014, 14, 22940–22970. [Google Scholar] [CrossRef] [PubMed]
  16. Silva, J.; Heim, W.; Chau, T. A self-contained, mechanomyography-driven externally powered prosthesis. Arch. Phys. Med. Rehabil. 2005, 86, 2066–2070. [Google Scholar] [CrossRef] [PubMed]
  17. Orizio, C.; Diemont, B.; Esposito, F.; Alfonsi, E.; Parrinello, G.; Moglia, A.; Veicsteinas, A. Surface mechanomyogram reflects the changes in the mechanical properties of muscle at fatigue. Eur. J. Appl. Physiol. Occup. Physiol. 1999, 80, 276–284. [Google Scholar] [CrossRef] [PubMed]
  18. Hong-Bo, X.; Yong-Ping, Z.; Jing-Yi, G. Classification of the mechanomyogram signal using a wavelet packet transform and singular value decomposition for multifunction prosthesis control. Physiol. Meas. 2009, 30, 441–457. [Google Scholar]
  19. Barry, D.T.; Leonard, J.A.; Gitter, A.J.; Ball, R.D. Acoustic myography as a control signal for an externally powered prosthesis. Arch. Phys. Med. Rehabil. 1986, 67, 267–269. [Google Scholar] [PubMed]
  20. Gobbo, M.; Cè, E.; Diemont, B.; Esposito, F.; Orizio, C. Torque and surface mechanomyogram parallel reduction during fatiguing stimulation in human muscles. Eur. J. Appl. Physiol. 2006, 97, 9–15. [Google Scholar] [CrossRef] [PubMed]
  21. Xie, H.-B.; Zheng, Y.-P.; Guo, J.-Y.; Chen, X.; Shi, J. Estimation of wrist angle from sonomyography using support vector machine and artificial neural network models. Med. Eng. Phys. 2009, 31, 384–391. [Google Scholar] [CrossRef] [PubMed]
  22. Youn, W.; Kim, J. Feasibility of using an artificial neural network model to estimate the elbow flexion force from mechanomyography. J. Neurosci. Methods 2011, 194, 386–393. [Google Scholar] [CrossRef] [PubMed]
  23. Shamshirband, S.; Petković, D.; Amini, A.; Anuar, N.B.; Nikolić, V.; Ćojbašić, Ž.; Mat Kiah, M.L.; Gani, A. Support vector regression methodology for wind turbine reaction torque prediction with power-split hydrostatic continuous variable transmission. Energy 2014, 67, 623–630. [Google Scholar] [CrossRef]
  24. Vapnik, V.; Golowich, S.; Smola, A. Support vector method for function approximation, regression estimation, and signal processing. In Advances in Neural Information Processing Systems; Mozer, M.C., Jordan, M.I., Petsche, T., Eds.; The MIT Press: Cambridge, MA, USA, 1997; Volume 9, pp. 281–287. [Google Scholar]
  25. Yang, H.; Huang, K.; King, I.; Lyu, M.R. Localized support vector regression for time series prediction. Neurocomputing 2009, 72, 2659–2669. [Google Scholar] [CrossRef]
  26. Jiang, H.; He, W. Grey relational grade in local support vector regression for financial time series prediction. Expert Syst. Appl. 2012, 39, 2256–2262. [Google Scholar] [CrossRef]
  27. Yu, W.; Liu, T.; Valdez, R.; Gwinn, M.; Khoury, M. Application of support vector machine modeling for prediction of common diseases: The case of diabetes and pre-diabetes. BMC Med. Inform. Decis. Mak. 2010, 10, 16. [Google Scholar] [CrossRef] [PubMed]
  28. Ziai, A.; Menon, C. Comparison of regression models for estimation of isometric wrist joint torques using surface electromyography. J. NeuroEng. Rehabil. 2011, 8, 1–12. [Google Scholar] [CrossRef] [PubMed]
  29. Güler, N.; Koçer, S. Use of support vector machines and neural network in diagnosis of neuromuscular disorders. J. Med. Syst. 2005, 29, 271–284. [Google Scholar] [CrossRef] [PubMed]
  30. Brown, L.E.; Weir, J.P. Asep procedures recommendation i: Accurate assessment of muscular strength and power. J. Exerc. Physiol. 2001, 4, 1–21. [Google Scholar]
  31. Bickel, C.S.; Slade, J.M.; VanHiel, L.R.; Warren, G.L.; Dudley, G.A. Variable-frequency-train stimulation of skeletal muscle after spinal cord injury. J. Rehabil. Res. Dev. 2004, 41, 33–40. [Google Scholar] [CrossRef] [PubMed]
  32. Orizio, C.; Solomonow, M.; Baratta, R.; Veicsteinas, A. Influence of motor units recruitment and firing rate on the soundmyogram and emg characteristics in cat gastrocnemius. J. Electromyogr. Kinesiol. 1992, 2, 232–241. [Google Scholar] [CrossRef]
  33. Adams, G.R.; Harris, R.T.; Woodard, D.; Dudley, G.A. Mapping of electrical muscle stimulation using mri. J. Appl. Physiol. 1993, 74, 532–537. [Google Scholar] [PubMed]
  34. Selkowitz, D.M. Improvement in isometric strength of the quadriceps femoris muscle after training with electrical stimulation. Phys. Ther. 1985, 65, 186–196. [Google Scholar] [PubMed]
  35. Babault, N.; Pousson, M.; Ballay, Y.; Van Hoecke, J. Activation of human quadriceps femoris during isometric, concentric, and eccentric contractions. J. Appl. Physiol. 2001, 91, 2628–2634. [Google Scholar] [PubMed]
  36. Levin, O.; Mizrahi, J.; Isakov, E. Transcutaneous fes of the paralyzed quadriceps: Is knee torque affected by unintended activation of the hamstrings? J. Electromyogr. Kinesiol. 2000, 10, 47–58. [Google Scholar] [CrossRef]
  37. Ebersole, K.T.; Housh, T.J.; Johnson, G.O.; Evetovich, T.K.; Smith, D.B.; Perry, S.R. The effect of leg flexion angle on the mechanomyographic responses to isometric muscle actions. Eur. J. Appl. Physiol. Occup. Physiol. 1998, 78, 264–269. [Google Scholar] [CrossRef] [PubMed]
  38. Ryan, E.D.; Cramer, J.T.; Egan, A.D.; Hartman, M.J.; Herda, T.J. Time and frequency domain responses of the mechanomyogram and electromyogram during isometric ramp contractions: A comparison of the short-time fourier and continuous wavelet transforms. J. Electromyogr. Kinesiol. 2008, 18, 54–67. [Google Scholar] [CrossRef] [PubMed]
  39. Shinohara, M.; Kouzaki, M.; Yoshihisa, T.; Fukunaga, T. Mechanomyogram from the different heads of the quadriceps muscle during incremental knee extension. Eur. J. Appl. Physiol. Occup. Physiol. 1998, 78, 289–295. [Google Scholar] [CrossRef] [PubMed]
  40. Katsavelis, D.; Joseph Threlkeld, A. Quantifying thigh muscle co-activation during isometric knee extension contractions: Within-and between-session reliability. J. Electromyogr. Kinesiol. 2014, 24, 502–507. [Google Scholar] [CrossRef] [PubMed]
  41. Shin, K.-S.; Lee, T.S.; Kim, H.-J. An application of support vector machines in bankruptcy prediction model. Expert Syst. Appl. 2005, 28, 127–135. [Google Scholar] [CrossRef]
  42. Owolabi, T.O.; Akande, K.O.; Olatunji, S.O. Support vector machines approach for estimating work function of semiconductors: Addressing the limitation of metallic plasma model. Appl. Phys. Res. 2014, 6, 122–132. [Google Scholar] [CrossRef]
  43. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  44. Burges, C.C. A tutorial on support vector machines for pattern recognition. Data Min. Knowl. Discov. 1998, 2, 121–167. [Google Scholar] [CrossRef]
  45. Gupta, S.M. Support vector machines based modelling of concrete strength. World Acad. Sci. Eng. Technol. 2007, 36, 305–311. [Google Scholar]
  46. Pal, M.; Goel, A. Prediction of the end-depth ratio and discharge in semi-circular and circular shaped channels using support vector machines. Flow Meas. Instrum. 2006, 17, 49–57. [Google Scholar] [CrossRef]
  47. Lin, S.-W.; Ying, K.-C.; Chen, S.-C.; Lee, Z.-J. Particle swarm optimization for parameter determination and feature selection of support vector machines. Expert Syst.Appl. 2008, 35, 1817–1824. [Google Scholar] [CrossRef]
  48. Mahmoud, S.A.; Olatunji, S.O. Automatic recognition of off-line handwritten arabic (Indian) numerals using support vector and extreme learning machines. Int. J. Imaging 2009, 2, 34–53. [Google Scholar]
  49. Akande, K.O.; Owolabi, T.O.; Olatunji, S.O. Investigating the effect of correlation-based feature selection on the performance of support vector machines in reservoir characterization. J. Nat. Gas Sci. Eng. 2015, 22, 515–522. [Google Scholar] [CrossRef]
  50. Akbani, R.; Kwek, S.; Japkowicz, N. Applying support vector machines to imbalanced datasets. In Machine Learning: Ecml 2004; Boulicaut, J.-F., Esposito, F., Giannotti, F., Pedreschi, D., Eds.; Springer: Berlin, Germany, 2004. [Google Scholar]
  51. Cherkassky, V.; Ma, Y. Practical selection of svm parameters and noise estimation for svm regression. Neural Netw. 2004, 17, 113–126. [Google Scholar] [CrossRef]
  52. Rıza Yıldız, A. A novel hybrid immune algorithm for global optimization in design and manufacturing. Robot. Comput.-Integr. Manuf. 2009, 25, 261–270. [Google Scholar] [CrossRef]
  53. Owolabi, T.O.; Akande, K.O.; Olatunji, S.O. Estimation of surface energies of hexagonal close packed metals using computational intelligence technique. Appl. Soft Comput. 2015, 31, 360–368. [Google Scholar] [CrossRef]
  54. Olatunji, S.O.; Selamat, A.; Abdul Raheem, A.A. Improved sensitivity based linear learning method for permeability prediction of carbonate reservoir using interval type-2 fuzzy logic system. Appl. Soft Comput. 2014, 14, 144–155. [Google Scholar] [CrossRef]
  55. Gencoglu, M.T.; Uyar, M. Prediction of flashover voltage of insulators using least squares support vector machines. Expert Syst. Appl. 2009, 36, 10789–10798. [Google Scholar] [CrossRef]
  56. Tay, F.E.H.; Cao, L. Application of support vector machines in financial time series forecasting. Omega 2001, 29, 309–317. [Google Scholar] [CrossRef]
  57. Gregory, C.M.; Bickel, C.S. Recruitment patterns in human skeletal muscle during electrical stimulation. Phys. Ther. 2005, 85, 358–364. [Google Scholar] [PubMed]
Figure 1. Experimental set-up at 90° knee angle showing the arrangement of stimulation electrodes (A) cathode, (B) anode Neuromuscular Electrical Stimulation (NMES) electrodes, and (C) Mechanomyographic signal (MMG) sensor in a representative participant.
Figure 1. Experimental set-up at 90° knee angle showing the arrangement of stimulation electrodes (A) cathode, (B) anode Neuromuscular Electrical Stimulation (NMES) electrodes, and (C) Mechanomyographic signal (MMG) sensor in a representative participant.
Sensors 16 01115 g001
Figure 2. Schematic representation of the experimental setup. Stimulation electrodes (A) cathode, (B) anode NMES electrodes, and (C) MMG sensor.
Figure 2. Schematic representation of the experimental setup. Stimulation electrodes (A) cathode, (B) anode NMES electrodes, and (C) MMG sensor.
Sensors 16 01115 g002
Figure 3. Flow chart of the procedure for obtaining optimal parameters (Table 1) for the proposed SVR model.
Figure 3. Flow chart of the procedure for obtaining optimal parameters (Table 1) for the proposed SVR model.
Sensors 16 01115 g003
Figure 4. Plots of the correlation coefficients for the training (A) and testing; (B) subsets.
Figure 4. Plots of the correlation coefficients for the training (A) and testing; (B) subsets.
Sensors 16 01115 g004
Figure 5. Cross plots of training sets—actual vs. predicted values: The plots show the performance of SVR with Gaussian kernel for torque prediction on the training set.
Figure 5. Cross plots of training sets—actual vs. predicted values: The plots show the performance of SVR with Gaussian kernel for torque prediction on the training set.
Sensors 16 01115 g005
Figure 6. Cross plots of testing set sets—actual vs. predicted values: The plots show the performance of SVR with Gaussian kernel for torque prediction on the testing set.
Figure 6. Cross plots of testing set sets—actual vs. predicted values: The plots show the performance of SVR with Gaussian kernel for torque prediction on the testing set.
Sensors 16 01115 g006
Table 1. Optimal parameters for the proposed Support Vector Regression model.
Table 1. Optimal parameters for the proposed Support Vector Regression model.
C879
Hyper-parameter (Lambda)2−15
Epsilon ( ε )0.1205
Kernel option54
KernelGaussian (RBF)
Table 2. Summary of the datasets: Mechanomyographic signal (MMG) characteristics at seven Neuromuscular Electrical Stimulation (NMES) intensities, at three knee angles and their respective peak torque values.
Table 2. Summary of the datasets: Mechanomyographic signal (MMG) characteristics at seven Neuromuscular Electrical Stimulation (NMES) intensities, at three knee angles and their respective peak torque values.
Stimulation Intensity (mA)Knee Angle
30°60°90°
PTRMSPTPPTRMSPTPPTRMSPTP
2013.9 (3.7)14.7 (9.9)23.6 (16.4)4.1 (0.7)17.4 (18.4)21.7 (23.0)4.3 (5.7)20.4 (20.0)22.0 (26.9)
3023.3 (19.7)51.9 (22.4)55.8 (24.7)9.7 (8.5)37.4 (21.2)38.7 (19.5)11.0 (10.2)51.3 (30.0)50.8 (30.6)
4058.2 (23.6)75.3 (29.5)73.54 (19.1)27.6 (24.2)77.7 (36.6)65.88 (19.2)21.4 (15.0)93.4 (45.4)84.0 (34.1)
5076.6 (19.3)84.2 (15.2)85.04 (14.9)51.5 (26.2)82.6 (27.3)72.7 (14.9)40.7 (18.5)115.7 (39.6)101.0 (33.8)
6086.1 (20.2)104.9 (22.5)94.86 (18.2)74.7 (19.2)94.9 (30.4)85.27 (14.7)62.1 (12.3)104.3 (29.0)101.1 (28.5)
7091.1 (21.5)100.2 (6.2)98.34 (5.7)91.0 (8.2)88.1 (9.7)90.57 (14.4)84.2 (12.3)118.5 (22.0)113.2 (10.7)
80100.0 (0)100.0 (0)100.0 (0)100.0 (0)100.0 (0)100.0 (0)100.0 (0)100.0 (0)100.0 (0)
Abbreviations: Stimulation Intensity—level of electrical stimulation or contraction intensity, PT—Peak torque, RMS—Normalized MMG-RMS%, PTP—Normalized MMG-PTP%. Values are reported in mean (standard deviation) for N = 8.
Table 3. Statistical parameters of the datasets.
Table 3. Statistical parameters of the datasets.
Input ParametersMeanMaxMedianStdevMin
Participants
Weight (kg)70.180695.963
Age (years)23.42523.51.321
Stimulation intensity (mA)5080502020
Knee angle (°)60906024.530
Normalized MMG-RMS%77.8188.186.940.04
Normalized MMG-PTP%75.2163.581.634.84.6
Peak torque53.9108.457.2380
Table 4. Performance measures that determined the accuracy of the developed model.
Table 4. Performance measures that determined the accuracy of the developed model.
Performance MeasuresTrainingTesting
r0.970.94
R294%89%
RMSE9.4812.95

Share and Cite

MDPI and ACS Style

Ibitoye, M.O.; Hamzaid, N.A.; Abdul Wahab, A.K.; Hasnan, N.; Olatunji, S.O.; Davis, G.M. Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression. Sensors 2016, 16, 1115. https://doi.org/10.3390/s16071115

AMA Style

Ibitoye MO, Hamzaid NA, Abdul Wahab AK, Hasnan N, Olatunji SO, Davis GM. Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression. Sensors. 2016; 16(7):1115. https://doi.org/10.3390/s16071115

Chicago/Turabian Style

Ibitoye, Morufu Olusola, Nur Azah Hamzaid, Ahmad Khairi Abdul Wahab, Nazirah Hasnan, Sunday Olusanya Olatunji, and Glen M. Davis. 2016. "Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression" Sensors 16, no. 7: 1115. https://doi.org/10.3390/s16071115

APA Style

Ibitoye, M. O., Hamzaid, N. A., Abdul Wahab, A. K., Hasnan, N., Olatunji, S. O., & Davis, G. M. (2016). Estimation of Electrically-Evoked Knee Torque from Mechanomyography Using Support Vector Regression. Sensors, 16(7), 1115. https://doi.org/10.3390/s16071115

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop