Journal of Medical Signals & Sensors

: 2020  |  Volume : 10  |  Issue : 3  |  Page : 185--195

A novel method for individual age group determination based on the hand muscle synergy

Arash Maghsoudi, Fereidoon Nowshiravan Rahatabad, Parisa Rangraz 
 Department of Biomedical Engineering, Science and Research Branch, Islamic Azad University, Tehran, Iran

Correspondence Address:
Dr. Parisa Rangraz
Department of Biomedical Engineering, Islamic Azad University, Science and Research Branch, Tehran


Background: As people get older, muscles become more synchronized and cooperate to accomplish an activity, so the main purpose of this research is to determine the relationship between changes in age and the amount of muscle synergy. The presence of muscle synergies has been long considered in the movement control as a mechanism for reducing the degree of freedom of the motor system. Methods: By combining these synergies, a wide range of complex movements can be produced. Muscle synergies are often extracted from the electromyogram (EMG) signal. One of the most common methods for extracting synergies is the nonnegative matrix factorization. In this research, the EMG signal is obtained from individuals from different age groups (namely 15–20 years, 25–30 years, and 35–40 years), and after preprocessing, the muscular synergies are extracted. By processing and studying these synergies. Results: It was observed that there is a significant difference between the muscular synergy of different age groups. Furthermore, there was a significant difference in the mean value of synergy coefficients in each group, especially in motions that were accompanied by force. Conclusion: This result candidates this parameter as a biomarker to differentiate and recognize the effects of age on the individual's muscular signal. In the best case, using the synergy tool, classification of the age of persons can be done by 77%.

How to cite this article:
Maghsoudi A, Rahatabad FN, Rangraz P. A novel method for individual age group determination based on the hand muscle synergy.J Med Signals Sens 2020;10:185-195

How to cite this URL:
Maghsoudi A, Rahatabad FN, Rangraz P. A novel method for individual age group determination based on the hand muscle synergy. J Med Signals Sens [serial online] 2020 [cited 2020 Aug 3 ];10:185-195
Available from:

Full Text


In recent years, the study on muscle synergy has been considered by many researchers and physicians, and due to the high potential, it provides for treatment of complications. However, it is not completely known how the central nervous system coordinates all the muscles involved in a particular movement. It has been suggested that muscle synergy provides the nervous system with a strategy for control of the motor system. In this regard, several different models have been presented to describe various human movements, such as standing, walking, running, and pedaling.[1] For most scholars in the field of motor studies, the synergy term is associated with Bernstein and its famous problem of “motion redundancy.”[2],[3] Cheung et al. have compared the muscle synergies of a person with stroke and a healthy person while performing movements and claimed that, since synergies are related to the spinal cord, brain damage does not change the individual's inability to perform different tasks.[4],[5] Azari Pasand et al. suggested that individual's inability to perform movements is related to the inappropriate stimulation factor of synergies.[6] They postulated the electromyogram (EMG) signal as a linear combination of a number of basic vector models. The activity of the muscles is measured as the rate of ignition of the motor neurons and is limited to nonnegative values. There are, in general, different methods for testing but generally matrix-decomposition algorithms are used such as principal component analysis, independent component analysis, and nonnegative matrix factorization (NMF).[7]

Eskandari et al. studied the reliability and coordination of muscle synergy during isometric trunk insertion in sagittal and transverse synergies of healthy individuals.[8] In their study, the lumbar spine model was used in conjunction with 18 muscles. One of the limitations of this study is that only the equilibrium conditions have been taken into account.

The lifestyle of individuals and their physical and recreational activities may affect muscle strength. Low-antigravity muscle weakness plays an important role in losing the balance and hence injuries to people.[9] In a study, muscle neuromuscular adaptation has been studied with age. The results indicate that the muscular nervous system in older people is influenced by the choice of fast entropy fibers and the difference in firing motor unit.[10] Ahamed et al.[11] recorded EMG from the biceps muscle to identify the age group and they concluded that, in larger time windows, the resolution of the root-mean-square (RMS) signal is increased. Of course, no synergy has been used in their study. Khayat et al. characterized prolonged standing and its effect on postural control in old and adult individuals and showed that the signal of the center of pressure can be studied in the chaotic domain and characteristic features can be extracted from the signal in the frequency domain and chaotic domain too.[12]

In this study, considering the importance of balance and power in the health of the elderly, the effect of age on muscle strength on women aged 21–80 years was determined. As it was expected, the results concurred that muscle strength decreases significantly with age. According to the results, it can be predicted that as age increases and muscle strength decreases, the standing and standing balance of individuals decreases. Furthermore, a biomarker for classifying age group based on the EMG synergy is proposed.

In the following, first, data acquisition is described, and then, the proposed algorithm including signal preprocessing, synergy calculation, time segmentation, RMS extraction from time segments, and statistical analysis of the difference between the RMSs of different age groups are described.


[Figure 1] shows the flowchart of the proposed method. According to this chart, first, data were acquired and preprocessed. Then, synergy was extracted and the similarity test was done. Furthermore, after time segmentation of synergies, RMS was calculated and compared in different age groups.{Figure 1}


Data used in this research were acquired from 18 women in three age groups (six people aged 15–20 years, six people aged 25–30 years, and six people aged 35–40 years). To minimize the effect of ECG interference in the records, all right-handed people were selected and the trials were done with the right hand.

The Biopac-MP100 and surface electrodes were used to record surface electromyogram (SEMG) data. Bipolar SEMG electrodes were placed on the surface of the skin (with 2-cm inter-electrode distance) over the brachioradialis muscle, flexor carpi, triceps, and biceps. A ninth electrode was mounted to the opposite (nonpreferred) acromion process in the function of an “electrical common” for data recording. An elastic bandage was wrapped around the EMG electrodes to secure the devices from extraneous movement; it did not impeding muscular function or movement about the shoulder and elbow joints. An amplifier with 60 dB gain was used.[13] Sampling frequency was set at 1000 Hz.[14]

For each individual, the EMG signal is recorded during six activities:

Hand in restBiceps flexion with a weight of 2 kgBiceps muscle at a maximum forceTriceps muscle at a maximum forceThe flexor carpi muscle at a maximum forceBrachioradialis muscle at a maximum force.


The initial 5 s of the signals were saved and the remaining was deleted. A fourth Butterworth bandpass filter with a cutoff frequency of 50–150 Hz was used to eliminate artifacts and noise.

Synergy extraction

Given that a four-channel EMG was recorded, a maximum of four synergies were extracted from each instance of the signal.[1],[4] A nonnegative matrix factorization (NNMF) analysis method was used to extract synergies. Given that the matrix in this method should be nonnegative, the magnitude of the signal was first calculated and then fed to the matrix analysis algorithm. [Figure 2] shows the error value in terms of the number of synergies to determine the appropriate number of synergies. As can be seen, in the case of three synergies, the error has dropped to a satisfactory level, and therefore, this number of synergy was the most appropriate choice.{Figure 2}

Then, the variance accounted for VAF is calculated using the following equation:


Where y is the main data, ˆ is extracted data of NNMF. [Figure 3] shows the amount of VAF in terms of synergy to determine the appropriate number of synergies. For each data, the results of this chart were averaged over 18 subjects and 4 reconstructed outputs. As it can be seen, in the case of three synergies, the VAF has reached about 98%, and therefore, three synergies were considered.{Figure 3}

Synergy similarity verification

To compare the synergies of various subjects, their synergies were clustered in three different categories, so that in each activity, in each cluster, only one synergy of three synergies of each subject was considered. To do this, three synergies from a random subject were selected as the center of three clusters. Then, the synergies of the subsequent random subject were assigned to each cluster based on the distance from the cluster's center, and the center of each cluster was updated at the end of each epoch. This process was continued to the last subject.

After clustering of synergies, the similarity between the synergies of each cluster was calculated using the cosine similarity criterion, which was calculated for the two vectors X and Y as following:


Comparison of root-mean-square rates at different age ranges

The 5-s signal of each subject was divided into ten windows, and energies of each of these windows were calculated. To determine the intensity of the signal in each window, the RMS of signal is used as:


Where j is the window number and N is the number of samples in each window. It should be noted that for comparing the performance of synergies in classification, the RMS value was calculated for both synergies and different channels.


Synergy clustering results

Synergies of different individuals were clustered into three clusters, as it was described in the previous section. [Figure 4] shows the average values of these similarities for different activities. As can be seen, intracluster similarity for synergy 1 and 2 was always more than 0.7, indicating that there was a high similarity between the synergies placed in these clusters.[15]{Figure 4}

Comparison of mean root-mean-squares in synergies

Using the method described in the previous section, it was observed that seventh window yielded the maximum discrepancy between the age groups. [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10] compare the dispersion values of RMSs in activities one through six in the seventh window between different age ranges. In these figures, the three upper plots represent the RMSs dispersion, and the three lower plots represent P values from t-test. As it can be seen, in the figures relating to the second, third, and sixth activities, some differences are found in the mean value of RMSs, which indicate that these activities are more appropriate than other activities (3, 1, and 5) to distinguish between the age groups.{Figure 5}{Figure 6}{Figure 7}{Figure 8}{Figure 9}{Figure 10}

According to [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10], it can be concluded that the first synergy value in the second activity, the second synergy value in the third activity, and the first synergy value in the sixth activity can be used to identify the age ranges.

Comparison of mean channel root-mean square

RMS of each channel also was calculated to be compared with the results attained from the previous section. Given that in the previous section it has been determined that activities 2, 3, and 6 can differentiate between different age categories, RMS calculation was performed only for these three activities. [Figure 11], [Figure 12], [Figure 13] show the dispersion of the RMS values of the main signals in these three activities and in the seventh window. As expected, the classification of the 3 age groups can be achieved in some of the channels. For instance, in the second activity, the third channel could partially divide the first age from the second and third ones.{Figure 11}{Figure 12}{Figure 13}

Classification using synergies

The RMS of the three most significant synergies (Section 3–2) was used for classifying using the k-nearest neighbor method. K-fold cross-validation with K = 10 has been used.

[Table 1] shows the results of this classification using the most significant synergies. K =3 was set on training data, and as it can be seen, using that K, the best accuracy was about 77%. [Figure 14] shows the related confusion matrix.{Table 1}{Figure 14}

Classification using the channels

[Table 2] shows the classification results using the channels. K =4 was set on training data and as it can be seen, using that K, the best accuracy is about 50%. [Figure 15] shows the related confusion matrix.{Table 2}{Figure 15}

Comparison of this value with the value obtained by synergies shows that synergy extraction has been able to provide better characteristics for the classification of age ranges.

Classification using a combination of synergy and channel

[Table 3] shows the results of the classification using both synergy and channel features. [Figure 16] shows the related confusion matrix. As it is seen, a slight improvement was obtained in contrast to using only synergies or channel information.{Table 3}{Figure 16}

 Discussion and Conclusion

In this study, it was demonstrated that using the EMG signal, the age group can be determined. Considering that in only three of the six activities, the classification of age groups was possible, and it can be postulated that activities which assert force and time on muscle and thus leading to fatigue can be used for the classification of age groups. Furthermore, the accuracy obtained in this work is 75% (for male samples and for the detection of a specific range of 23–29 years [11]) which is higher compared to similar works. For precision of 74%, the research was conducted for female samples and only for the detection of a specific range of 23–29 years.[9] In the present study, the testing and recording of signals from female samples have been taken, but three age groups have been investigated, and ultimately, the obtained accuracy is much better than previous works.

Financial support and sponsorship


Conflicts of interest

There are no conflicts of interest.



Arash Maghsoudi received BSc from Amirkabir University of Technology, Tehran, Iran in 2004, Master's and Ph.D from Sharif university of Technology, Terhan, Iran in 2006 and 2012 respectively. He was a Researcher at Nanyang Technological University, Singapore in 2010-2011. He is currently assistant professor at department of biomedical engineering, Science and research branch, Azad University, Tehran, Iran. His research interest includes Biological signal processing, Brain Connectomics and image processing.



Fereidoon Nowshiravan Rahatabad is the Clinical Engineering Faculty Member in Azad University, Science and Research Branch Tehran, Iran. He finished his MSc at 2005 in bioelectrical engineering from Science Branch of Islamic Azad University, and he received his BSc at 2003 in bio- electrical engineering from Science Branch of Islamic Azad University. His research is mostly focus on Biological System Control and modeling.



Parisa Rangraz received the B.Sc. degree in 2004, the M.Sc. degree in 2007 and the Ph.D. degree in 2012, all in biomedical engineering from the Science and Research Branch, Islamic Azad University. Currently she is an Assistant Professor at Biomedical Engineering Department of Science and Research Branch of Islamic Azad University, Tehran, Iran. Her research interests are Ultrasound in Medicine, Medical Image processing, Medical signal processing and Ultrasound guided HIFU.



1Kristiansen M, Samani A, Madeleine P, Hansen EA. Muscle synergies during bench press are reliable across days. J Electromyogr Kinesiol 2016;30:81-8.
2Latash ML, Scholz JP, Schöner G. Toward a new theory of motor synergies. Motor Control 2007;11:276-308.
3Aruin AS, Almeida GL, Latash ML. Organization of a simple two-joint synergy in individuals with down syndrome. Am J Ment Retard 1996;101:256-68.
4Cheung VC, Piron L, Agostini M, Silvoni S, Turolla A, Bizzi E. Stability of muscle synergies for voluntary actions after cortical stroke in humans. Proc Natl Acad Sci U S A 2009;106:19563-8.
5Kaboudvand N, Tohidkhah F, Gharibzadeh SH. Extracting and Study of Synchronous Muscle Synergies During fast Arm Reaching Movements; 20th Iranian Conference on Biomedical Engineering; 2013.
6Azari Pasand P, Fallah A, Maleki A. Algorithm of myoelectric signals processing for the control of prosthetic robotic hands. J Comp Sci Technol 2018;18:28-34.
7Tresch MC, Cheung VC, d'Avella A. Matrix factorization algorithms for the identification of muscle synergies: Evaluation on simulated and experimental data sets. J Neurophysiol 2006;95:2199-212.
8Eskandari AH, Sedaghat-Nejad E, Rashedi E, Sedighi A, Arjmand N, Parnianpour M. The effect of parameters of equilibrium-based 3-D biomechanical models on extracted muscle synergies during isometric lumbar exertion. J Biomech 2016;49:967-73.
9Akbari M, Mousavi R. The relationship of lower limbs muscle strength with aging healthy women. J Mazandaran Univ Med Sci 2010;20:1-2.
10Yamada H, Masuda T, Okada M. Age-related EMG variables during maximum voluntary contraction. Percept Mot Skills 2002;95:10-4.
11Ahamed N, Alqahtani M, Altwijri O, Rahman M, Sundaraj K. Age-related EMG responses of the biceps brachii muscle of young adults. Biomed Res 2016;27:787-93.
12Khayat O, Zargarchi M, Razjouyan J, Siahi M, Nowshiravan Rahatabd F. Age effect on static balance: Center-of-pressure analysis in chaotic and frequancy domains. J Paramed Sci Rehabil 2014;3:53-60.
13Merletti R. Standards for reporting EMG data. J Electromyogr Kinesiol 1999;9:3-4.
14Zachry T, Wulf G, Mercer J, Bezodis N. Increased movement accuracy and reduced EMG activity as the result of adopting an external focus of attention. Brain Res Bull 2005;67:304-9.
15Doan A, Halevy A, Ives Z. Principles of Data Integration. Ch. 7. Elsevier; 2012.