Angkoon Phinyomark

Franck Quaine

Yann Laurillau

The Relationship between Anthropometric Variables and Features of Electromyography Signal for Human-Computer Interface

The Relationship between Anthropometric Variables and Features of Electromyography Signal for Human-Computer Interface

Angkoon Phinyomark, Franck Quaine, Yann Laurillau

INTRODUCTION

Surface electromyography (EMG) signals are measured by surface electrodes that are placed on the target muscles. During muscle contractions, a compound of the whole motor unit action potentials (MUAPs) occurred in the muscles. These MUAPs are the useful information for numerous fields, e.g. rehabilitation engineering, biomechanics, ergonomics, and human-computer interfaces (HCIs) (Merletti & Parker, 2004). Surface EMG signals can also be used in a medical decision support system, e.g. the diagnosis of neuromuscular disorders [START_REF] Subasi | Medical decision support system for diagnosis for neuromuscular disorders using DWT and fuzzy support vector machines[END_REF][START_REF] Subasi | Classification of EMG signals using PSO optimized SVM for diagnosis of neuromuscular disorders[END_REF]. In this chapter, we focus on the development of the HCIs based on surface EMG signals, as called "muscle-computer interfaces or MCIs" [START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF]. These interfaces can widely use in controlling many external devices, e.g. prosthetic limbs, electric-power wheelchairs, interactive surfaces, a virtual computer mouse or keyboard, a portable music player, and in-car electronic equipment (e.g. [START_REF] Barreto | A practical EMG-based human-computer interface for users with motor disabilities[END_REF][START_REF] Benko | Enhancing input on and above the interactive surface with muscle sensing[END_REF][START_REF] Khushaba | Muscle computer interfaces for driver distraction reduction[END_REF][START_REF] Saponas | Enabling always-available input with muscle-computer interfaces[END_REF][START_REF] Shenoy | Online electromyographic control of a robotic prosthesis[END_REF][START_REF] Wei | Fusing EMG and visual data for hands-free control of an intelligent wheelchair[END_REF]Yang, Lin, Lin, & Lee, in press).

MCIs generally consist of three main modules, as shown in Fig. 1. The first module composes of two submodules: surface EMG signal acquisition (hardware) and data pre-processing (software). In the data acquisition sub-module, surface EMG signals are firstly amplified with an amplifier due to small EMG amplitude, and may be filtered by hardware filters, i.e., notch, high pass and/or low pass filters. Then, continuous surface EMG signals are sampled using an analog-to-digital converter. In case of no hardware filters, the raw EMG data can be filtered by software filters, before the EMG data go to the next module, the pattern recognition module. Different patterns of surface EMG signals are classified and matched to the control commands in this module. The second module can be divided into three sub-modules: feature extraction, dimensionality reduction, and classification algorithms. All sub-modules are in the software part. The pre-processed surface EMG data from the first module are segmented into small time slot length using an adjacent or an overlapped windowing technique, and then some features are extracted in order to emphasize the relevant structures in the EMG signals and remove noises/irrelevant parts. A feature vector is formed and can be sent directly to a classifier, or the dimension of a feature vector is reduced by the dimensionality reduction technique before sending it to the classifier. The classifier maps the extracted features (the representative of the actions associated with surface EMG signals) to the target classes (the control commands for an external device). After the control commands were generated based on the mathematical functions in the second module, the third module is a control system, which serves as an interface between the software and hardware. In other words, an output command is converted from a digital code to an analog signal for controlling an external device. In this chapter, we focus on the second module, particularly the extraction of EMG features. Nearly all previous studies on MCIs based on EMG pattern recognition concentrate on improving recognition rate together with increasing the number of discriminated motions. Currently, the recognition rate is more than 90% in discriminating 4-12 finger, wrist, hand, and forearm motions. It should be noted that the recognition rate or the classification accuracy is calculated as a ratio of the number of correct classifications to the total number of classifications. More details about many previous studies on MCIs based on EMG pattern recognition concentrate on improving recognition rate can be found in additional readings in the Additional Reading Section. However, the high recognition rates reported are usually based in single-session experiments conducted in research laboratories [START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF].

On the other hand, in clinics when the context of the real-world requirements have been given attention, there are many issues that have to find the solutions, e.g. the robustness against a variety of noises [START_REF] Boostani | Evaluation of the forearm EMG signal features for the control of a prosthetic hand[END_REF][START_REF] Zardoshti-Kermani | EMG feature evaluation for movement control of upper extremity prostheses[END_REF], the robustness over time-between sessions [START_REF] Oskouei | Intra-session and inter-day reliability of forearm surface EMG during varying hand grip forces[END_REF][START_REF] Sensinger | Adaptive pattern recognition of myoelectric signals: exploration of conceptual framework and practical algorithms[END_REF] and between days (Boschmann, Kaufmann, Platzner, & Winkler, 2010;[START_REF] Zhang | An adaptive feature extractor for gesture SEMG recognition[END_REF]-due to EMG electrode location shift [START_REF] Young | The effects of electrode size and orientation on the sensitivity of myoelectric pattern recognition systems to electrode shift[END_REF][START_REF] Subasi | Medical decision support system for diagnosis for neuromuscular disorders using DWT and fuzzy support vector machines[END_REF], muscle fatigue, and/or variation in muscle contraction effort [START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF], or the robustness in limb position variations [START_REF] Chen | Effect of upper-limb positions on motion pattern recognition using electromyography[END_REF][START_REF] Fougner | Resolving the limb position effect in myoelectric pattern recognition[END_REF][START_REF] Jiang | Effect of arm position on the prediction of kinematics from EMG in amputees[END_REF]).

Among such issues, the problem of the cross-user classification is far from being a practical one [START_REF] Cannan | Automatic circumference measurement for aiding in the estimation of maximum voluntary contraction (MVC) in EMG systems[END_REF][START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF]. The proposed systems in the literature were trained and tested independently on the EMG signals from each user. In other words, EMG-based MCIs have to design and create for each user. Based on a preliminary study of [START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF], when the classifier is trained on the data from all subjects in the database except the testing subject, the classification accuracy computed from the testing subject data is about 50%, reducing from about 95% using the independent training and testing system. It should be noted that their system classifies five single finger motions based on three time-domain and frequency-domain features extracted from eight EMG channels placed around the upper forearm [START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF]. Moreover, based on our preliminary experiments, the selection of the robust EMG features is not an efficient solution to reduce the variability of surface EMG signals between users, like other previous mentioned problems.

This is due to the different compositions of the body across users, which result in the different muscle types with varying architecture characteristics. This factor prevents the development of standard MCIs that are compatible with almost any user, and MCIs based on EMG pattern recognition have never really reached the general population. Being able to solve this problem, in our previous work (Phinyomark et al., 2013a) a feasibility to use anthropometric variables, a measurable characteristic of the human body, to calibrate the EMG pattern recognition systems is presented. The anthropometric variables are used to compensate for discrepancies between user body compositions [START_REF] Cannan | Automatic circumference measurement for aiding in the estimation of maximum voluntary contraction (MVC) in EMG systems[END_REF], which are based on the relationship between the EMG signal, muscle force, and muscle size consisting of cross-sectional area and length [START_REF] Hof | Comparison between EMG to force processing and kinetic analysis for the calf muscle moment in walking and stepping[END_REF][START_REF] Marras | A three-dimensional motion model of loads on the lumbar spine: I. Model structure[END_REF] [START_REF] Cannan | Automatic circumference measurement for aiding in the estimation of maximum voluntary contraction (MVC) in EMG systems[END_REF][START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF]. It is also possible to obtain estimated anthropometric variables from the anthropometric table, which can be found in many sources including academic books (e.g. De [START_REF] De Leva | Adjustments to Zatsiorsky-Seluyanov's segment inertia parameters[END_REF][START_REF] Winter | Biomechanics and Motor Control of Human Movement[END_REF]. It means that there is no need to measure all anthropometric variables and may need only a few basic variables, i.e., gender, age, body mass, or standing height. However, at the beginning all anthropometric variables proposed in this chapter are measured directly from the volunteers.

For simple MCI systems that use a thresholding technique, a maximum voluntary contraction (MVC) can be used as an additional input to calibrate and adapt the classification system from one user to another user. In order to use MVC normalization technique, it is important that the linear relationship between force generated and EMG level should be held for the muscle of interest up to the maximum force/EMG level [START_REF] Bolgla | Reliability of electromyographic normalization methods for evaluating the hip musculature[END_REF][START_REF] Vera-Garcia | MVC techniques to normalize trunk muscle EMG in healthy women[END_REF]. Although the linear relationship has been found in many previous works (e.g. [START_REF] Kamavuako | Relationship between grasping force and features of single-channel intramuscular EMG signals[END_REF][START_REF] Woods | Linear and non-linear surface EMG/force relationships in human muscles. An anatomical/functional argument for the existence of both[END_REF], several researchers have found a non-linear relationship between force generated and EMG muscle activity (e.g. [START_REF] Lawrence | Myoelectric signal versus force relationship in different human muscles[END_REF][START_REF] Woods | Linear and non-linear surface EMG/force relationships in human muscles. An anatomical/functional argument for the existence of both[END_REF]). The relationship is more complicated in case of a dynamic muscle contraction, i.e., a muscle is free to change length and a joint is free to move [START_REF] Oatis | Kinesiology: The mechanics and pathomechanics of human movement[END_REF]. Several anthropometric variables have been found to have a strong relationship with the MVC and other techniques of maximum force/EMG measurement, such as hand circumference has the strong correlation with maximal grip strength (MGS) and can be used to predict the MGS using a simple regression equation in the study of [START_REF] Li | Predicting maximal grip strength using hand circumference[END_REF], or in the study of Marras & David (2001), several linear regression equations developed from five of the fifteen anthropometric variables are used to predict the EMC (Expected Maximum Contraction) of trunk muscles.

However, for advanced MCI systems, features extracted from surface EMG signals are the important stage, and only MVC (related only with muscle force) is not enough to be used alone. Other additional useful information is needed, which should have strong relationships with many EMG features (extracted not only force level information but also other useful information, e.g. complexity and frequency). So some additional inputs should be used to calibrate the machine learning and pattern recognition instead of the simple thresholding technique. Forearm circumference is an anthropometric variable that has already been proposed to be used as a calibrating variable in an EMG recognition system in previous works [START_REF] Anakwe | Grip strength and forearm circumference in a healthy population[END_REF][START_REF] Cannan | Automatic circumference measurement for aiding in the estimation of maximum voluntary contraction (MVC) in EMG systems[END_REF]. However, the relationship between EMG and forearm circumference is not strong enough, so further anthropometric variables were evaluated their relationships with surface EMG features in this chapter. All relationships between EMG features and anthropometric variables that are strong and significant could benefit a further design of the EMG-based

EMG FEATURE EXTRACTION Definition of Feature Extraction

Feature extraction is a technique to transform raw EMG signals into a reduced representation set of features as called "a feature vector". A transformation is usually based on mathematical functions, which highlights relevant structures in surface EMG signals and rejects irrelevant parts. If the extracted features are carefully selected, the feature vector will contain effective and relevant information drawn from the whole set of raw EMG data, which can be used to represent the desired actions. Training an EMG classifier on raw EMG signal patterns requires high computational cost and their testing classification performances are very poor (Phinyomark, Phukpattaranont, & Limsakul, 2012e). Therefore, the extraction of EMG features is very important step in EMG-based gesture recognition. It is noted in the literature that the success of the EMG pattern recognition system almost entirely depends on the choice of features [START_REF] Kendall | A novel approach to surface electromyography: An exploratory study of electrode-pair selection based on signal characteristics[END_REF][START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF].

EMG feature extraction can be computed in three domains: 1) time domain, 2) frequency domain, and 3) time-frequency or time-scale domain [START_REF] Boostani | Evaluation of the forearm EMG signal features for the control of a prosthetic hand[END_REF]. Time domain features are widely used in EMG pattern recognition due to the low complexity and computational cost, which can be implemented and computed using a simple hardware processor, i.e., a mobile device or a microprocessor [START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF][START_REF] Kundu | Design of wearable, low power, single supply surface EMG extractor unit for wireless monitoring[END_REF]. Moreover, time domain features can be successfully used as a dimensionality reduction technique for time-frequency/ time-scale features, i.e., the coefficients of short-time Fourier transform, discrete wavelet transform and wavelet packet transform [START_REF] Phinyomark | Feature extraction and reduction of wavelet transform coefficients for EMG pattern classification[END_REF]. Hence, in this chapter, features in time domain are mainly focused, and time-frequency/time-scale features are not considered. A number of frequency domain features are considered to be used in developing the robust MCI systems and to provide additional information.

Computation of EMG Feature Extraction

A basic idea about the computation of EMG feature extraction is illustrated in Fig. 2. Following is a brief description of all robust EMG features, which were investigated and discussed in the present chapter. In each, x i is the i th sample of surface EMG signal amplitude and L is the length of the analysis window for computing the features. Willison Amplitude (WAMP) is a number of times resulting from the difference between two consecutive EMG amplitudes in a time segment becomes more than a predefined threshold thr. In other words, WAMP is a number of times that the length of the EMG waveform in a time segment exceeds a predefined threshold thr, which can be expressed as .

(1)

Myopulse Percentage Rate (MYOP) is an average number of times that the absolute value of EMG amplitude, i.e., myopulse output, exceeds a predefined threshold thr, and the calculation is defined as .

(2)

Slope Sign Change (SSC) or the number of turns (NT) is a number of times that the slope of the EMG waveform changes sign. This method uses three consecutive EMG signal amplitudes to detect the EMG waveform slope, which is calculated as .

(3)

Zero Crossing (ZC) is a number of times that the EMG waveform crosses the zero amplitude axes. Several works used this value without implementing threshold condition [START_REF] Boostani | Evaluation of the forearm EMG signal features for the control of a prosthetic hand[END_REF][START_REF] Zardoshti-Kermani | EMG feature evaluation for movement control of upper extremity prostheses[END_REF]. However, a predefined threshold thr should be included in the ZC computation to reduce the effect of noise [START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF]. ZC can be defined as .

(4)

A suitable predefined threshold to reduce noise is dependent on a system gain value and a background noise level. So, we can roughly estimate a predefined threshold if we know both values [START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF]. For example, assuming a system gain of 1000 and a background noise of 4 µV peak to peak, the threshold can be calculated to be ± 2 mV. It is normally chosen between 50 µV and 100 mV [START_REF] Philipson | The electromyographic signal as a measure of muscular force: A comparison of detection and quantification techniques[END_REF]. However, to obtain a high quality feature space, the suitable threshold could be defined based on taking some offline trial-and-error classifications.

Detrended Fluctuation Analysis (DFA) is a modified root mean square analysis of a random walk [START_REF] Peng | Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series[END_REF]. A slope of the line relating the logarithm of the RMS fluctuation of the profiles y k (F n ) to the logarithm of the predefined box size (n) is used as a feature value (α). This value indicates the presence of the power law (fractal) scaling, as can be expressed as

(5)

Based on our previous work [START_REF] Phinyomark | Electromyography (EMG) signal classification based on detrended fluctuation analysis[END_REF], the minimum box size n min is set at four, the maximum box size n max is set at one-tenth of the signal length, and the box size increment is based on a power of two. A least-square fit, which is applied to the profiles y k , is the quadratic polynomial fit.

Higuchi method (HG) is another fractal dimension estimator [START_REF] Higuchi | Approach to an irregular time series on the basis of the fractal theory[END_REF], which performs well in the classification of surface EMG signals and the simulated signals [START_REF] Arjunan | Decoding subtle forearm flexions using fractal features of surface electromyogram from single and multiple sensors[END_REF][START_REF] Esteller | A comparison of waveform fractal dimension algorithms[END_REF]. For m = 1, 2, …, k, a negative slope of the line relating the logarithm of the length of the curve (L mk ) to the logarithm of the discrete time interval between points (k) is used as a feature value (D). Based on the finding of [START_REF] Phothisonothai | Fractal-based EEG data analysis of body parts movement imagery tasks[END_REF], the maximum time interval k max is set at 128. This value is from the relationship L k α k D , which can be given by ( 6)

Integrated EMG (IEMG) is a summation of the absolute value of the EMG signal amplitude over the time segment, which can be expressed as .

(7)

Waveform Length (WL) is a cumulative length of the EMG waveform over the time segment. It is the IEMG of the EMG wavelength. A number of previous studies call this feature as the "wavelength", which provides simple waveform complexity information. It can be calculated by .

Maximum Fractal Length (MFL) is proposed as a feature, which can measure EMG signal patterns at low level muscle contraction [START_REF] Arjunan | Fractal features of surface electromyogram: A new measure for low level muscle activation[END_REF]. It is defined as the average fractal length of the signal measured at the smallest scale from the HG method in Eq. ( 6). If the smallest scale is set at one, MFL can be defined as the modification of the WL using a logarithm function as follows .

(9)

Mean Absolute Value (MAV) is an average of the absolute value of the EMG signal amplitude. There are many ways to call this feature, e.g. averaged absolute value, average rectified value, integral of absolute value, integrated absolute value, and the first order of the v-Order feature. It can be given by .

(10)

Difference Absolute Mean Value (DAMV) is an average of the absolute value of the difference between the adjacent EMG amplitudes over the time segment [START_REF] Kim | Comparison of k-nearest neighbor, quadratic discriminant and linear discriminant analysis in classification of electromyogram signals based on the wrist-motion directions[END_REF][START_REF] Park | EMG pattern recognition based on artificial intelligence techniques[END_REF][START_REF] Yu | Classification of nine directions using the maximum likelihood estimation based on electromyogram of both forearms[END_REF]. It can be seen as the MAV of the EMG wavelength and sometimes is called as the "Average Amplitude Change (AAC)" [START_REF] Fougner | Proportional myoelectric control of a multifunction upper-limb prosthesis[END_REF]. It is given by .

(11)

Variance of EMG (VAR) is an average of the square value of the EMG signal amplitude. This is due to a mean value of surface EMG signals is close to zero, in which the equation is defined as .

(

) 12
Root Mean Square (RMS) is a square root of the average of the square of the EMG signal amplitude values. In other words, it has a similar value to the standard deviation of EMG signal amplitude. It can be defined as the modification of the VAR using a square root function as follows .

(

) 13
Difference Absolute Standard Deviation Value (DASDV) is a square root of the average of the square of the difference between the adjacent EMG amplitudes over the time segment. It is the RMS of the EMG wavelength, as can be formulated as .

(

) 14
Approximation Entropy (ApEn) and Sample Entropy (SampEn) are techniques used to quantify the unpredictability of fluctuations in EMG signal amplitudes over the time segment (Zhao et al., 2006a[START_REF] Zhao | EMG control for a five-fingered underactuated prosthetic hand based on wavelet transform and sample entropy[END_REF]. They provide the waveform similarity or complexity information. SampEn is developed from ApEn to avoid the bias caused by self-matching. Both methods use two input parameters: m and r, where m is the pattern length and r is the criterion of similarity [START_REF] Zhang | Sample entropy analysis of surface EMG for improved muscle activity onset detection against spurious background spikes[END_REF]. The basic idea is about the estimation of the conditional probability that the patterns of the EMG signal amplitude, which are similar to each other within a predefined tolerance r, will remain similar for the next comparison point [START_REF] Richman | Physiological time series analysis using approximate entropy and sample entropy[END_REF].

Mean Frequency (MNF) is an average frequency, which is calculated as a sum of the product of the EMG power spectrum and the frequency divided by a total sum of the EMG power spectrum. There are other ways to call this feature, such as central frequency and the spectral center of gravity. It can be defined as

, ( 15
)
where f j is the frequency of the spectrum at a frequency bin j, P j is the EMG power spectrum at a frequency bin j, and M is the length of whole frequency bin.

Median Frequency (MDF) is a frequency at which the EMG power spectrum is divided into two areas with an equal total power, as can be expressed as . ( 16)

Experimental Data for Extracting EMG Features

All features are computed from the EMG data, which are recorded from 20 subjects (10 males and 10 females). The subjects were asked to perform 8 motions consisting of forearm pronation (FP), forearm supination (FS), wrist flexion (WF), wrist extension (WE), wrist radial deviation (WR), wrist ulnar deviation (WU), hand open (HO), and hand close (HC), and maintain for 2 s. Each motion was repeated 15 trails per day for four separate days. EMG data were measured from 4 forearm muscles and an upper arm muscle: extensor carpi radialis longus (ECRL), extensor carpi ulnaris (ECU), extensor digitorum communis (EDC), flexor carpi radialis (FCR), and biceps brachii (BB). More details about data acquisition and experiments can be found in [START_REF] Phinyomark | Electromyography (EMG) signal classification based on detrended fluctuation analysis[END_REF].

To briefly explain the relationship between actions associated with surface EMG signals and EMG features, the amplitude shape of surface EMG signals acquired from all the muscles and motions from a trial in time domain is shown in Fig. 3, together with one extracted feature: MAV. It is clearly shown in the figure that the EMG amplitude shapes from all muscles are significantly different according to the direction of eight motions. For instance, we can observe that the EMG magnitudes of WE are very high in two extensor muscles, i.e., ECRL and ECU. These two muscles are the corresponding muscles to produce wrist extension, i.e., high contraction level. On the other hand, the EMG magnitudes of WF are low in both extensor muscles but are very high in the flexor muscle, FCR. As we know that the MAV feature is used to estimate the EMG amplitude/magnitude level, so MAV features of WE extracted from ECRL and ECU are higher than the other muscles. In the same way, the MAV feature of WF extracted from FCR is relatively higher than the other muscles.

ANTHROPOMETRIC VARIABLES Definition of Anthropometric Variables

Anthropometry is a measurement of the dimensions of the different parts of human body. There are two types of measurement: 1) static dimension and 2) dynamic dimension. In order to use anthropometric variables to calibrate the EMG-based MCI systems, only static dimensions are considered. There are many anthropometric variables, e.g. overhead reach height and breadth, shoulder circumference and length, and waist front and back length. However, only twelve related variables are chosen for the experiments in this chapter.

Computation of Anthropometric Variables

Following is a brief description of all anthropometric variables, which were investigated and discussed for this chapter.

Body mass is measured in kilograms (kg). It is defined as a subject's weight or mass, which can be measured by the balance type scales. A subject stands on the center of the scale platform.

Stature or standing height is measured in centimeters (cm). It is defined as the distance from the bottom of the feet to the top of the head of a subject, which can be measured by a stadiometer. A subject stands erect on the center of the base plate with heels together, and head in the Frankfort plane.

BMI or body mass index is a roughly estimation of human body fat based a subject's weight and height in kilograms/meters 2 (kg/m 2 ). It can be defined as .

(17)

The remaining anthropometric variables are measured in centimeters (cm), as shown in Fig. 4, using a tape for biceps and forearm circumferences, a sliding or small bone caliper for hand breadth and length, a wide sliding torso caliper for elbow-hand grip length, elbow-fingertip length, shoulder-elbow length and bi-deltoid breadth, and a measuring block with tape measure for forward grip reach. The specific locations are defined as the following: Biceps circumference is a linear distance around the upper arm when subject stands. The upper arm is extended forward horizontally and the elbow is flexed about 90º (the arm is abducted). It is measured at the level of the drawn biceps point landmark, which is the point of maximum quiet inspiration.

Forearm circumference is a linear distance around the lower arm with the same posture as measured the biceps circumference. It is usually measured at the level of maximum forearm circumference.

For hand breadth, subject sits with right hand flats on a table and the fingers are together and straight.

Hand breadth is measured from metacarpalphalangeal joint II to metacarpalphalangeal joint V.

Hand length is measured from the wrist landmark to dactylion with the same posture as the hand breadth.

Elbow-hand grip length is measured from the posterior tip of the olecranon process to the center of grip during holding a pencil. It can be called as the elbow to center of grip.

Elbow-fingertip length is measured from the posterior tip of the olecranon process to dactylion.

Shoulder-elbow length is measured from the right acromion landmark to the inferior tip of the olecranon process of the right elbow.

Forward grip reach or functional reach is measured from the back wall to the tip of the thumb.

Bi-deltoid breadth is measured across the body at the level of the deltoid landmarks.

Experimental Anthropometric Data

It should be noted that the anthropometric variables are measured in the same day for all twenty subjects from the right arm of the subject using the standard instruments (Centurion Kit, Rosscraft). Due to the significant difference of muscle size and force between male and female subjects, anthropometric variables have usually been investigated and discussed by gender: male and female [START_REF] Anakwe | Grip strength and forearm circumference in a healthy population[END_REF][START_REF] Holzbaur | Upper limb muscle volumes in adult subjects[END_REF]. In our experiments, the difference of anthropometric variables between male and female subjects is statistically significant at p < 0.01 for 8 parameters (body weight, stature, biceps circumference, forearm circumference, hand breadth, hand length, elbow-fingertip length, and forward grip reach) and significant at p < 0.10 for the remainders (BMI, elbow-hand grip length, shoulder-elbow length, and bi-deltoid breadth). In addition, to avoid the effect of ages we used the subjects from the same age, around 21 years old ( n o s i g n i f i c a n t d i f f e r e n c e , p = 0 . 4 5 8 ) . Twelve anthropometric variables from 10 males and 10 females are reported respectively in Table 1 andTable 2

STATISTICAL ANALYSIS

Correlation analysis is used as an evaluating function in the study. It measures a relationship between two parameters and provides a statistic known as the correlation r coefficient. This coefficient shows the degree of a linear relationship between two measured variables. Correlation coefficients can interpret as the weak or low correlation when r ≤ 0.35, the modest or moderate correlation when 0.35 < r ≤ 0.67, and the strong or high correlation when 0.67 < r ≤ 1 [START_REF] Taylor | Interpretation of the correlation coefficient: A basic review[END_REF]. Actually, the r value contains both a magnitude and a direction (positive and negative) of the relationship. However, in our experiment only the magnitude of correlation coefficient (the absolute value of average r coefficient) is used. Moreover, due to a small number of subjects n (less than 20), the significant level is set at p < 0.05 and is tested using t-test, as can be defined by .

(

) 18
It should be noted that degrees of freedom for entering the t-distribution is defined as n-2. Only the combinations between EMG features extracted from a muscle and anthropometric variables that have the strong and significant correlations in at least 4 motions are considered for the next stage.

RELATIONSHIP BETWEEN EMG FEATURE AND ANTHROPOMETRIC VARIABLE

To reach the high potential of EMG-based MCIs in the context of real-world requirements, the proposed anthropometric variables have been investigated the relationship with robust EMG features, which are divided into 4 groups as follows:

(1) Noise tolerance: WAMP, MYOP, SSC, and ZC;

( It should be noted that the relationship between anthropometric variables and the Hudgins' time domain feature set, i.e. SSC, ZC, WL, MAV and mean absolute value slope, is presented and discussed in detail in our previous study [START_REF] Phinyomark | EMG Feature Evaluation for Improving Myoelectric Pattern Recognition Robustness[END_REF].

Feature Set 1: Noise Tolerance

Noise is one of the major problems in the analysis of surface EMG signals. [START_REF] Zardoshti-Kermani | EMG feature evaluation for movement control of upper extremity prostheses[END_REF] proposed that the robustness of EMG features is one of the three major properties for evaluating high-quality features. In their study, the definition of robustness is limited to the tolerance of noise. Artificial white noise at different noise levels was used to evaluate the robustness of features. The cluster separability of feature spaces should be preserved as much as possible in a noisy environment. This property is important because both biological and environmental resources can generate noise and artifact. It is difficult to record only a pure signal component without noise. For this reason, noise removal or noise reduction is an important step before performing feature extraction in most of MCI systems. Many noise removal tools have been proposed, e.g. notch filter, band pass filter, adaptive filter, and wavelet filter [START_REF] Jindapetch | FPGA implementations of an ADALINE adaptive filter for power-line noise cancellation in surface electromyography signals[END_REF][START_REF] Phinyomark | Wavelet-based denoising algorithm for robust EMG pattern recognition[END_REF]. Unfortunately, there are not any filters that can remove one hundred percent of noise. Sometimes it removes some important parts of real EMG signals [START_REF] De Luca | Filtering the surface EMG signal: Movement artifact and baseline noise contamination[END_REF][START_REF] Li | Conditioning and sampling issues of EMG signals in motion recognition of multifunctional myoelectric prostheses[END_REF]. So EMG features that have high tolerance for biological and environmental noises should increase the ability of EMG-based gesture recognition.

Noises contaminated in surface EMG signal can be categorized into four major types: 1) ambient noise, 2) motion artifact, 3) the inherent instability of surface EMG signals, and 4) the inherence in electronic components in the detection and recording equipment [START_REF] De Luca | Surface electromyography: Detection and recording[END_REF][START_REF] Reaz | Techniques of EMG signal analysis: Detection, processing, classification and applications[END_REF]. The first three types have a specific band of frequencies, but the last noise type has a broaden frequency band and falls in a usable energy band of surface EMG signals. In the literature, four noisy EMG signals that consist of baseline noise, movement artifact, power-line interference, and random noise, i.e., white Gaussian noise, have been simulated and used to develop robust EMG methods [START_REF] Phinyomark | The usefulness of wavelet transform to reduce noise in the SEMG signal[END_REF]. First two simulated noises, i.e., baseline noise and movement artifact have a narrow frequency band below 20 Hz, which ranges outside the EMG energy band. We can use a high-pass filter at a cut-off frequency of 20 Hz to remove these kinds of noise [START_REF] De Luca | Filtering the surface EMG signal: Movement artifact and baseline noise contamination[END_REF]. On the other hand, power-line interference (50/60 Hz) and random noise (0-1000 Hz) have a frequency band in a range of usable EMG energy band (20-500 Hz), therefore such noises can be only reduced to a predictable level or cannot be entirely removed.

In our previous works [START_REF] Phinyomark | EMG feature extraction for tolerance of white Gaussian noise[END_REF], 2009), nine time domain and frequency domain features: RMS, WL, ZC, SSC, WAMP, HIST, AR, MDF, and MNF, were evaluated the robust performance based on two noise types: 50-Hz interference and white Gaussian noise, at different signal-to-noise ratios (0-20 dB SNRs). The experimental results showed that, on average WAMP is the most robust feature, followed by ZC and SSC features. These features are in the same feature group, together with MYOP, which compute in time domain and provide frequency information [START_REF] Phinyomark | Feature reduction and selection for EMG signal classification[END_REF]. All features in this group have a threshold parameter, which can be used to avoid some noises [START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF]. In addition, the distribution of feature spaces in this group is similar to each other, particularly for WAMP and ZC. 3, all four features extracted from the BB muscle have a strong relationship with the bideltoid breadth for female subjects. For male subjects, there are two relationships: 1) WAMP and MYOP from the ECU muscle and bi-deltoid breadth, and 2) ZC from the BB muscle and biceps circumference. All relationships share the motions of FP and WF.

Although the SSC feature has a higher average |r| value than other three features for female subjects, it does not have a strong relationship with any anthropometric variables for male subjects. SSC and ZC are the members of the Hudgins' time domain feature set [START_REF] Hudgins | A new strategy for multifunction myoelectric control[END_REF] and their classification performances are not influenced by the variation in muscle contraction effort [START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF]. Other features in Hudgins' time domain feature set are MAV, WL, and mean absolute value slope. [START_REF] Chan | Fuzzy EMG classification for prosthesis control[END_REF] re-evaluated the classification performance of the Hudgins' feature set using fuzzy approach. The SSC feature does not improve classification performance and in some subjects the classification accuracy decreases. So SSC is not recommended for developing robust MCI system.

Based on the similar distribution of feature spaces of ZC and WAMP [START_REF] Phinyomark | Feature reduction and selection for EMG signal classification[END_REF], one of them should be selected in order to avoid the information redundancy in an EMG feature vector. The selection depends on the choice of the system. If the automatically calibrated system is preferred, the ZC feature is recommended. On the other hand, if the semi-automatically calibrated system is acceptable, WAMP is a suitable feature. WAMP and MYOP have the same associations with the anthropometric variables, whereas WAMP is the most robust feature against noises. Moreover, WAMP has higher classification accuracy than other three features in the classification of upper-limb motions based on a robust linear discriminant analysis (LDA) c l a s s i f i e r with and without periodic retraining scheme (Phinyomark et al., 2011a, in press b). WAMP is also successful in the classification of lower-limb exercise activities for the elderly (Phinyomark et al., 2012a) and in the estimation of muscle force levels [START_REF] Kamavuako | Influence of the feature space on the estimation of hand grasping force from intramuscular EMG[END_REF].

Feature Set 2: Low-level and High-level Surface EMG Signals

If the low-noise EMG signals are obtained from an EMG recording system with/without noise removal algorithms, the next challenge in EMG-based gesture recognition is about weak measured EMG signals.

Simple MCI systems use the potential of a large difference between a high EMG magnitude at strong muscle contractions and a very low magnitude (close to zero) at the rest of the contractions. EMG amplitude estimators, i.e., RMS and MAV are sufficient enough for classifying the proposed activities [START_REF] Clancy | Probability density of the surface electromyogram and its relation to amplitude detectors[END_REF]. However, such systems often offer a few control schemes with a single speed of actuation based on a few EMG channels. To increase the number of control commands, a variety of motions and EMG channels have been proposed for the recognition system. In that case, it is not possible to get only strong EMG signals. Weak or low-level EMG signals (low SNR signals) must be measured.

There are two cases: 1) a weak EMG signal is measured from a major corresponding muscle when a user performs a low force motion or a little movement [START_REF] Arjunan | Fractal features of surface electromyogram: A new measure for low level muscle activation[END_REF] and 2) a weak EMG signal is measured from a minor corresponding muscle when a user performs a high force motion or a strong movement [START_REF] Phinyomark | Fractal analysis features for weak and single-channel upper-limb EMG signal[END_REF].

To deal with the analysis of low-level EMG signals, a fractal analysis is used. Fractals refer to signal patterns, which exhibit self-similarity, that have fractal dimension. Fractal dimension of surface EMG signals is found under different types and levels of muscle contraction using several fractal methods, e.g. box-counting method [START_REF] Gupta | Fractal analysis of surface EMG signals from the biceps[END_REF], correlation dimension method [START_REF] Hu | Classification of surface EMG signal with fractal dimension[END_REF], critical exponent analysis [START_REF] Phinyomark | Critical exponent analysis applied to surface electromyography (EMG) signals for gesture recognition[END_REF], and Katz method [START_REF] Gitter | Fractal analysis of the electromyographic interference pattern[END_REF]. Based on the finding in previous studies, a fractal dimension of EMG signals depends on the level of muscle contraction during strong or high level activities [START_REF] Gupta | Fractal analysis of surface EMG signals from the biceps[END_REF][START_REF] Hu | Classification of surface EMG signal with fractal dimension[END_REF]. For the low level of muscle contraction, fractal dimension does not change with change in the level of muscle contraction. It is a measure of the size and complexity of the muscles [START_REF] Arjunan | Fractal features of surface electromyogram: A new measure for low level muscle activation[END_REF].

In our previous work [START_REF] Phinyomark | Fractal analysis features for weak and single-channel upper-limb EMG signal[END_REF], a fractal dimension estimator DFA was examined the performance in the classification of low-level EMG signal patterns compared with another fractal estimator HG [START_REF] Arjunan | Decoding subtle forearm flexions using fractal features of surface electromyogram from single and multiple sensors[END_REF]. The experimental results showed that, DFA has better cluster separability than HG [START_REF] Phinyomark | Electromyography (EMG) signal classification based on detrended fluctuation analysis[END_REF], as well as classification performance [START_REF] Phinyomark | Fractal analysis features for weak and single-channel upper-limb EMG signal[END_REF]. On the other hand, features based on a magnitude detector provide better performance in the classification of high-level EMG signal patterns than fractal features. So a combination of fractal and magnitude features is recommended [START_REF] Arjunan | Fractal properties of surface electromyogram for classification of low-level hand movements from single-channel forearm muscle activity[END_REF].

The past several years, magnitude features extracted from the first-order difference of EMG time series, i.e., MFL, DAMV, DASDV, and WL, present higher accuracy than magnitude features extracted from the original time series, i.e., VAR, IEMG, MAV, and RMS [START_REF] Brzostowski | Analysis of human arm motions recognition algorithms for system to visualize virtual arm[END_REF][START_REF] Kim | Comparison of k-nearest neighbor, quadratic discriminant and linear discriminant analysis in classification of electromyogram signals based on the wrist-motion directions[END_REF][START_REF] Kim | Improvements of multi-features extraction for EMG for estimating wrist movements[END_REF][START_REF] Yu | Classification of nine directions using the maximum likelihood estimation based on electromyogram of both forearms[END_REF]. In addition, the distribution of feature spaces in both groups is similar to each other for some motions and muscles. [START_REF] Phinyomark | Fractal analysis features for weak and single-channel upper-limb EMG signal[END_REF].

For magnitude features, strong associations are found only with the bi-deltoid breadth. For female subjects, only features extracted from the BB muscle have a strong relationship with the variable. All relationships share the motions of FP, FS, and WR. For male subjects, two muscles have the strong relationships: ECU and BB muscles. Although feature space distribution of features extracted from the original and the first-order difference of EMG signals is similar to each other for some motions and muscles, the best feature from each group is usually selected in many previous studies to be a member of an optimal multiple feature set based on the feature selection scheme [START_REF] Phinyomark | EMG Feature Evaluation for Improving Myoelectric Pattern Recognition Robustness[END_REF].

For the first-order differencing magnitude features, MFL outperforms others, followed respectively by DAMV, DASDV, and WL, based on many state-of-the-art classifiers, i.e., the LDA, the quadratic discriminant analysis (QDA), the artificial neural network (ANN), the k-nearest neighbor (k-NN), and the maximum likelihood estimation (MLE), with and without the periodic retraining scheme [START_REF] Arjunan | Decoding subtle forearm flexions using fractal features of surface electromyogram from single and multiple sensors[END_REF][START_REF] Kim | Comparison of k-nearest neighbor, quadratic discriminant and linear discriminant analysis in classification of electromyogram signals based on the wrist-motion directions[END_REF][START_REF] Phinyomark | EMG Feature Evaluation for Improving Myoelectric Pattern Recognition Robustness[END_REF][START_REF] Yu | Classification of nine directions using the maximum likelihood estimation based on electromyogram of both forearms[END_REF]. Although WL performs the best in accuracy, stability, and computation load among the single features in several previous studies [START_REF] Oskoei | Support vector machine-based classification scheme for myoelectric control applied to upper limb[END_REF][START_REF] Phinyomark | Evaluation of EMG Feature Extraction for Hand Movement Recognition Based on Euclidean Distance and Standard Deviation[END_REF], it is because of its popularity and other first-order differencing features were not included in the evaluating studies. In total, MFL is recommended to be used as a representative feature in this group instead of WL.

For the original magnitude features, all features have the same relationship with the anthropometric variable. Due to its similarity of IEMG to MAV and its similarity of VAR to RMS, the IEMG and VAR features should be excluded, because in comparison IEMG and VAR result in weaker performance in classification than MAV and RMS [START_REF] Oskoei | Support vector machine-based classification scheme for myoelectric control applied to upper limb[END_REF][START_REF] Phinyomark | Feature reduction and selection for EMG signal classification[END_REF]. For MAV and RMS, both features are close together in the success rate of classification, as well as the computational cost. The selection of MAV and RMS should be based on the probability density function (PDF) of the measured EMG signals. If the EMG PDF is close to the Gaussian density, an optimal EMG amplitude estimator is RMS [START_REF] Hogan | Myoelectric signal processing: Optimal estimation applied to electromyography-Part I: Deviation of the optimal myoprocessor[END_REF] based on both theoretically and experimentally. On the other hand, if the EMG PDF is close to the Laplacian density, an optimal EMG amplitude estimator is MAV [START_REF] Clancy | Probability density of the surface electromyogram and its relation to amplitude detectors[END_REF]. Moreover, MAV and RMS performance in the estimation of muscle force is similar to the performance of WAMP and WL features [START_REF] Kamavuako | Influence of the feature space on the estimation of hand grasping force from intramuscular EMG[END_REF].

Feature Set 3: Fluctuating EMG Signals over Time

Multifunction EMG pattern recognition systems that allow high recognition rate in research laboratories are usually based in single-session experiments. Most of related works collected EMG data only from a single day. A next major problem is about the fluctuation of surface EMG signals over time. It is widely known that EMG signals measured in one day are relatively different from the EMG signals measured in another day, even on the same muscle and also same subject [START_REF] Jain | Improving long term myoelectric decoding, using an adaptive classifier with label correction[END_REF][START_REF] Kaufmann | Fluctuating EMG signals: Investigating long-term effects of pattern matching algorithms[END_REF]. This can be due to many factors, such as the EMG electrode location shift, the variation in muscle contraction effort, or change in the architecture characteristics of the subject [START_REF] Saponas | Making muscle-computer interfaces more practical[END_REF][START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF]. So it is possible that features extracted from the initial training data (one day) and the present testing data (another day) are significant difference, which will degrade the performance of the pattern matching algorithms. In order to implement the MCIs based on EMG pattern recognition in clinics, the robustness over time of the system (i.e., feature extraction and classification algorithm) is very important.

In the literature, however, this problem has rarely been evaluated. [START_REF] Sensinger | Adaptive pattern recognition of myoelectric signals: exploration of conceptual framework and practical algorithms[END_REF] and [START_REF] Jain | Improving long term myoelectric decoding, using an adaptive classifier with label correction[END_REF] found that the performance of the EMG pattern recognition systems can degrade within hours after initial classifier training. Zhang et al. (2007) studied the effect of the number of days used as the training data sets. They found that if the data used for training the classifier decrease, the classification accuracy also decreases based on EMG signals recorded from five different days. [START_REF] Kaufmann | Fluctuating EMG signals: Investigating long-term effects of pattern matching algorithms[END_REF] evaluated the performance of five state-of-the-art classifiers: support vector machine (SVM), decision tree (DT), the ANN, k-NN, and LDA, in the classification of surface EMG signals recorded from 21 days. They also found that the performance of the EMG pattern recognition systems degrades with increasing time difference between initial classifier training and testing data for all studied classifiers (about 8%-15%), except the LDA (about 3.6%) using the Hudgins' time domain feature set.

Using the same EMG data recorded during 21 days, fifty state-of-the-art EMG features in time domain and frequency domain were evaluated in one of our previous works [START_REF] Phinyomark | EMG Feature Evaluation for Improving Myoelectric Pattern Recognition Robustness[END_REF]. Among the single features, SampEn performed the highest classification accuracy (93.37%) by training the LDA classifier using EMG from only the first day without any retraining classifier again, followed by ApEn (84.68%) and MFL (82.07%). The EMG data were recorded from four surface EMG channels on the forearm for 21 days with 5-6 trials per day, and in each trial 11 wrist and finger motions were performed [START_REF] Kaufmann | Fluctuating EMG signals: Investigating long-term effects of pattern matching algorithms[END_REF].

SampEn and ApEn features extracted from the FCR muscle have a strong relationship with the forward grip reach for female subjects. Both relationships share the four motions of FP, FS, HO, and HC. In addition to the robustness over time, SampEn can use to detect the onset of EMG activity while suppressing spurious background spikes [START_REF] Zhang | Sample entropy analysis of surface EMG for improved muscle activity onset detection against spurious background spikes[END_REF].

Feature Set 4: Muscle Fatigue

Muscle fatigue is generally resulted when a user performs the repetitive proposed motions for a long time and cannot produce the certain level of force with the muscle [START_REF] De Luca | Myoelectrical manifestations of localized muscular fatigue in humans[END_REF]. Generally, muscle fatigue results in an increase in EMG signal amplitude and a downward shift of EMG frequency spectrum [START_REF] Cifrek | Surface EMG based muscle fatigue evaluation in biomechanics[END_REF]. However, [START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF] reported that the effect of muscle fatigue on the classification performance of eleven time domain features (i.e., WAMP, ZC, SSC, WL, VAR, MAV, v-Order, log-Detector, EMG histogram, autoregressive coefficients, and cepstrum coefficients) is very weak. On the other hand, for frequency domain features, MNF and MDF have been hailed so far as the gold standard for the muscle fatigue detection with surface EMG signals [START_REF] Al-Mulla | A review of non-invasive techniques to detect and predict localised muscle fatigue[END_REF].

The experimental results showed that MDF does not have a strong relationship with any anthropometric variables, whereas MNF extracted from the FCR muscle of female subjects has a strong relationship with the hand length for four motions: FS, WU, HO, and HC at the average correlation coefficient |r| of 0.76 (in range of 0.71-0.87). Further, MNF performs better performance in the classification of upper-limb motions than MDF and other frequency domain features, i.e., peak frequency, frequency ratio, power spectrum ratio, and the variance of central frequency [START_REF] Phinyomark | Feature reduction and selection for EMG signal classification[END_REF][START_REF] Phinyomark | EMG Feature Evaluation for Improving Myoelectric Pattern Recognition Robustness[END_REF]. The modified MNF can also be used to estimate muscle force like time domain features (Phinyomark et al., 2012g;[START_REF] Thongpanja | Mean and median frequency of EMG signal to determine muscle force based on time-dependent power spectrum[END_REF]. If EMG frequency information is needed in developing EMG-based MCIs, the MNF feature is recommended to be used as an optimal frequency domain feature for both the classification of actions-based EMG signal and the assessment of muscle fatigue during actions.

SEMI-AUTOMATIC AND AUTOMATIC CALIBRATION SYSTEMS

As mentioned in the Introduction section, we can use the anthropometric variables to calibrate the systems in two different ways: a weighting factor for a classifier and a normalizing value for EMG features. Both ways can be used for adapting the system automatically. However, the difference between semi-automatic and automatic calibration systems is about the way to measure anthropometric variables. From twelve proposed anthropometric variables, only two variables: forearm and biceps circumferences can measure directly from an EMG measuring armband device, which is composed of surface EMG electrodes and the circumference measurement [START_REF] Cannan | Automatic circumference measurement for aiding in the estimation of maximum voluntary contraction (MVC) in EMG systems[END_REF]. They can be used as calibrated inputs for an automatic calibration system. On the other hand, the remaining variables cannot measure directly from the EMG measuring armband device. These variables need to be measured manually, so they can be used as calibrated inputs for a semi-automatic calibration system. From the experiments, 7 out of the 12 proposed variables have strong relationships with EMG features, which can be divided into two groups: 1) the bideltoid breadth, stature, elbow-fingertip length, BMI, forward grip reach, and hand length for a semiautomatic calibration system, and 2) the biceps circumference for an automatic calibration system.

Among the seven variables, the bi-deltoid breadth has more associations with EMG features than other variables. It can be used to calibrate many robust EMG features consisting of WAMP, MFL, MAV, and RMS. The use of a feature vector of WAMP, MFL, and MAV/RMS would result in the high classification accuracy under robust conditions such as noisy environment, variation in muscle contraction effort, and low-and-high-level muscle contractions. For instance, using this EMG multiple feature set produced 89.97±7.34% accuracy when classified the EMG data recorded from 21 days without the retraining classifier scheme [START_REF] Kaufmann | Fluctuating EMG signals: Investigating long-term effects of pattern matching algorithms[END_REF] and only the bi-deltoid breadth variable is used to calibrate the recognition system. However, if the system is not limited by a low processor, the multiple robust feature sets with respect to all studied robust conditions, which consisted of some robust features: WAMP, MFL, MAV/RMS, DFA, SampEn, and MNF, would result in high classification accuracy and could be used across users by calibrating with the correspond anthropometric variables. Hence, in future works, the stability of the possible multiple feature sets with respect to several studied robust conditions should be evaluated their performance together with the development of calibration techniques using the correspond anthropometric variables.

On the other hand, only a combination of the ZC and the biceps circumference has a strong relationship, which can be used for an automatic calibration system. However, due to a rapid increased number of EMG wearable devices and a success of classifying motions associated with surface EMG recorded from the armband, new EMG features and other useful motions/muscles should be evaluated their associations with forearm and biceps circumferences in future studies. In addition to both the circumferences, hand circumference may be a useful anthropometric variable, because it is possible to record surface EMG signals using electrode sensor ring around the hand. This variable was also found that it has a strong relationship with the maximal grip strength and can use to create a simple model to predict the maximal grip strength [START_REF] Li | Predicting maximal grip strength using hand circumference[END_REF]. In order to explore the potential use of automatic calibration systems, a review of the EMG pattern recognition using electrode arm bands is presented in the following subsection.

EMG sensor arm bands

To measure wrist and finger actions associated with surface EMG signals, traditionally, multiple surface electrodes need to be placed right above all the corresponding muscles. However, it is difficult to identify the exact position of the muscles [START_REF] Kendall | A novel approach to surface electromyography: An exploratory study of electrode-pair selection based on signal characteristics[END_REF][START_REF] Vigreux | Factors influencing quantified surface EMGs[END_REF]. Because users typically have no detailed knowledge of human anatomy hence the EMG system is limited to only medical applications, and several muscles are closely together so it is not possible to acquire EMG from some specific single muscles. Moreover, the placing electrodes in this way require more time consuming. Errors in electrode placement degrade the classification accuracy of EMG-based MCIs [START_REF] Tkach | Study of stability of time-domain features for electromyographic pattern recognition[END_REF][START_REF] Young | The effects of electrode size and orientation on the sensitivity of myoelectric pattern recognition systems to electrode shift[END_REF][START_REF] Young | Improving myoelectric pattern recognition robustness to electrode shift by changing interelectrode distance and electrode configuration[END_REF].

Recent studies have designed and developed multichannel sensor rings/armbands to solve this problem. In this way, end-users will not have expertise and time to work on the electrode placement. It can be clearly seen from a recent commercial EMG sensor armband (i.e., MYO armband, see www.getmyo.com). The MYO armband is composed of a large amount of small active electrodes in a form of sensor ring and can transfer the EMG data over a wireless network, which is feasible for use outside the laboratory. This device is easy to set up and configure by a user, and provides an always-available and highly personalized input interface with is a low cost (inexpensive There are two useful positions on the forearm that have been proposed in related works on EMG pattern recognition:

(1) Lower forearm positions, i.e., around the wrist (Rhee, You, & Shin, 2011;[START_REF] Tang | Hand motion classification using a multi-channel surface electromyography sensor[END_REF][START_REF] You | Finger motion decoding using EMG signals corresponding various arm postures[END_REF][START_REF] Li | Electromyography pattern-recognition-based control of powered multifunctional upperlimb prostheses[END_REF]; (2) Upper forearm positions, i.e., below the elbow [START_REF] Andews | Optimal electrode configurations for finger movement classification using EMG[END_REF][START_REF] Benko | Enhancing input on and above the interactive surface with muscle sensing[END_REF][START_REF] Du | Portable hand motion classifier for multi-channel surface electromyography recognition using gray relational analysis[END_REF][START_REF] Du | The effect of combining stationary wavelet transform and independent component analysis in the multichannel SEMGs hand motion identification system[END_REF][START_REF] Khushaba | Electromyogram (EMG) feature reduction using mutual components analysis for multifunction prosthetic fingers control[END_REF][START_REF] Khushaba | Muscle computer interfaces for driver distraction reduction[END_REF][START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF][START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF][START_REF] Shyu | A new electrode system for hand action discrimination[END_REF][START_REF] Smith | Real-time myoelectric decoding of individual finger movements for a virtual target task[END_REF][START_REF] Smith | Continuous decoding of finger position from surface EMG signals for the control of powered prostheses[END_REF][START_REF] Tenore | Towards the control of individual fingers of a prosthetic hand using surface EMG signals[END_REF][START_REF] Tenore | Decoding of individual finger movements using surface electromyography[END_REF].

For the placement of the EMG electrodes around the forearm, there are two possible ways:

(1) An adapted scheme, i.e., the distances between every two channels of electrodes are changed and dependent on the forearm sizes of the users [START_REF] Andews | Optimal electrode configurations for finger movement classification using EMG[END_REF][START_REF] Benko | Enhancing input on and above the interactive surface with muscle sensing[END_REF][START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF][START_REF] Tang | Hand motion classification using a multi-channel surface electromyography sensor[END_REF];

(2) A fixed scheme, i.e., the distances between every two channels of electrodes are fixed, except one pair of electrodes in which the distance is varied depending on forearm circumference, e.g. one pair across the ulnar border [START_REF] Du | The effect of combining stationary wavelet transform and independent component analysis in the multichannel SEMGs hand motion identification system[END_REF][START_REF] Du | Portable hand motion classifier for multi-channel surface electromyography recognition using gray relational analysis[END_REF][START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF][START_REF] Shyu | A new electrode system for hand action discrimination[END_REF][START_REF] Smith | Continuous decoding of finger position from surface EMG signals for the control of powered prostheses[END_REF][START_REF] Smith | Real-time myoelectric decoding of individual finger movements for a virtual target task[END_REF]. In this scheme, however, there is no guarantee that electrodes are placed over the same muscles in all users.

For lower forearm positions, the EMG electrodes are designed to place on one side of the forearm: anterior forearm or posterior forearm, as a half wristband. In [START_REF] You | Finger motion decoding using EMG signals corresponding various arm postures[END_REF][START_REF] You | Finger flexion motion inference from sEMG signals[END_REF] and [START_REF] Rhee | Recognition of finger motion with sEMG and gyrosensor signals[END_REF], four EMG channels are placed on the anterior side of the forearm (flexor muscles). The recognition system is used to classify all five single finger motions and three multi-finger motions (indexmiddle fingers, middle-ring fingers, and hand close) and provides the average accuracy of about 95%-97%. On the other hand, [START_REF] Tang | Hand motion classification using a multi-channel surface electromyography sensor[END_REF] used six EMG channels placed on the posterior side of the forearm (extensor muscles) instead of the anterior forearm to recognize three single finger motions (thumb, index finger, middle finger), and eight multi-finger motions. The success recognition rate is higher than 89%. [START_REF] Tang | Hand motion classification using a multi-channel surface electromyography sensor[END_REF] mentioned that six channels are enough to cover the circumference of the posterior side and cover all useful extensor muscles: extensor digitorum, extensor pollicis longus and brevis, extensor indicis, and extensor digiti minimi.

For upper forearm positions, the EMG electrodes are designed to place on both sides of the forearm. The number of EMG channels placed in a narrow uniform band/ring, which is proposed in the literature, is seven sensors [START_REF] Du | The effect of combining stationary wavelet transform and independent component analysis in the multichannel SEMGs hand motion identification system[END_REF][START_REF] Du | Portable hand motion classifier for multi-channel surface electromyography recognition using gray relational analysis[END_REF][START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF][START_REF] Shyu | A new electrode system for hand action discrimination[END_REF] and eight sensors (Andrews et al., 2009;[START_REF] Khushaba | Electromyogram (EMG) feature reduction using mutual components analysis for multifunction prosthetic fingers control[END_REF][START_REF] Khushaba | Muscle computer interfaces for driver distraction reduction[END_REF][START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF][START_REF] Smith | Real-time myoelectric decoding of individual finger movements for a virtual target task[END_REF]. When multiple sets of sensors were placed around the forearm, reasonable EMG signals with useful information could be acquired from the armband even the sensors are only approximately placed [START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF]. The landmark sensor armband is placed below the elbow in several conditions, e.g. approximately one third of the distance from the proximal end of a line from medial epicondyle to the head of the radius (Andrews et al., 2009;[START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF], or approximately two inches or five centimeters below the elbow [START_REF] Du | The effect of combining stationary wavelet transform and independent component analysis in the multichannel SEMGs hand motion identification system[END_REF][START_REF] Du | Portable hand motion classifier for multi-channel surface electromyography recognition using gray relational analysis[END_REF][START_REF] Shyu | A new electrode system for hand action discrimination[END_REF][START_REF] Smith | Real-time myoelectric decoding of individual finger movements for a virtual target task[END_REF]. For the landmark of the starting electrode pair, there are several ways, such as placing just superior to the ulna (Andrews et al., 2009) or over the flexor carpi radialis [START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF]. On average, the proposed systems achieve more than 90%-95% accuracy. [START_REF] Tenore | Decoding of individual finger movements using surface electromyography[END_REF][START_REF] Tenore | Towards the control of individual fingers of a prosthetic hand using surface EMG signals[END_REF] reported that there is no statistically significant difference (p < 0.05) in the classification accuracy using 32 electrodes placed on upper and lower forearm (6 sensor rings) and 19 electrodes placed on only upper forearm (3 sensor rings). The accuracy in classifying ten single flexed and extended finger motions is greater than 90% even in a transradial amputee. It means that the multichannel sensor ring may be a kind of redundant sensor. [START_REF] Mogk | Crosstalk in surface electromyography of the proximal forearm during gripping tasks[END_REF] found that the EMG signal amplitudes between adjacent electrode pairs (3 cm apart) share 40% common signal (or crosstalk), and the common signal is reduced to about 10% at 6 cm spacing and 2.5% at 9 cm, while only 2% common signal or crosstalk is found between flexor and extensor electrode pairs. On the other hand, we could place a similar sensor band on the other forearm to provide even more input possibilities [START_REF] Saponas | Demonstrating the feasibility of using forearm electromyography for muscle-computer interfaces[END_REF].

CONCLUSIONS AND FUTURE WORKS

This chapter presents the relationships between robust EMG features and twelve related anthropometric variables for MCIs. In feature extraction view point, the robust EMG feature is recommended for each feature set 1-4, which is based on robust conditions in the context of real-world requirements. It can be summarized, as the following:

(1) WAMP is the most robust feature against a variety of noises, i.e., random noise and power-line interference, followed by ZC. WAMP is also an optimal feature to estimate muscle force. It has a strong relationship with bi-deltoid breadth. (2) DFA, a fractal (complexity) feature, is the suitable feature for the classification of low-level EMG signals, and MFL and MAV/RMS are the optimal magnitude features for the classification of high-level EMG signals. DFA has a strong relationship with the stature and the elbow-fingertip length variables, whereas MFL and MAV/RMS have strong relationships with breadth.

(3) SampEn is the most robust feature for the variability of the muscle contraction over time. It is not only reliable for long-term usage but also for noise tolerance, i.e., spurious background spikes. SampEn has a strong relationship with forward grip reach. (4) MNF is used as a standard muscle fatigue detector. It is reliable for a noisy environment and can provide additional information in frequency domain. It has a strong relationship with hand length.

The optimal robust f e a t u r e s m e n t i o n e d a b o v e h a v e s t r o n g r e l a t i o n s h i p s w i t h s e v e r a l experimental anthropometric variables, i.e., biceps circumference, bi-deltoid breadth, standing height, elbow-fingertip length, forward grip reach, and hand length. These variables can be used to calibrate the MCI systems in two ways. The first way is a manual calibrated input system or a semi-automatic calibration, which can be implemented using all anthropometric variables. In this way, the system can develop using many robust features, for instance, a robust multiple EMG feature set may consist of WAMP, DFA, MFL, MAV/RMS, SampEn, and MNF, and each feature can be calibrated by its related anthropometric variables. The second way is as an auto-calibration system. Based on this condition, forearm and bicep circumferences are two out of twelve variables that can be measured directly from wearable EMG device and used automatically to adapt the recognition system. However, in our experiments (eight upper-limb motions and five muscle positions) only the ZC feature has a strong and significant relationship with the biceps circumference. In future works, new features and/or other motions and muscles should be evaluated their associations with anthropometric variables.

In addition to the measured anthropometric data, we mentioned in the Introduction Section that that it is possible to get the estimated anthropometric data from the published anthropometric tables. So in future works the estimated anthropometric data should be evaluated their relationships with EMG features too.

If there are the strong and significant relationships between them, it means that we can add the anthropometric tables into the classification system and use there information to calibrate the system without measuring the anthropometric variables directly from the user. Further, the calibration techniques using the anthropometric variables: a weighting factor for a classifier and/or a normalizing value for EMG features should be developed and evaluated in future works. A simple idea for using anthropometric variables as a normalizing value for EMG features can be found in our preliminary study (Phinyomark et al., 2013a).

ACKNOWLEDGMENT

This work is funded by the FI MSTIC University Joseph Fourier Grenoble 1 (TIGRE project).

Figure 1 .

 1 Figure1. Three main modules of MCIs: the surface EMG acquisition and pre-processing, the pattern recognition, and the control system.

Figure 2 .

 2 Figure 2. An example of EMG feature extraction procedure in time domain.

Figure 3 .

 3 Figure 3. An example of surface EMG signals in time domain from the first trail in the first day of subject 1 (5 muscles and 8 motions with rest state) and the extracted MAV features.

Figure 4 .

 4 Figure 4. Anthropometric measurements: stature, forward grip reach, biceps circumference, forearm circumference, hand breadth, hand length, shoulder-elbow length, elbow-hand grip length, elbowfingertip length, and bi-deltoid breadth.

  2) Low-level and high-level surface EMG signals: DFA, HG, MFL, DAMV, DASDV, WL, VAR, IEMG, MAV, and RMS; (3) Fluctuating EMG signals over time: SampEn and ApEn; (4) Muscle fatigue: MNF and MDF.

  see our preliminary work-Phinyomark et al., 2013a-for more details about this relationship). They can automatically or semiautomatically calibrate a system. As a result of the proposed technique, EMG-based MCIs that require little or no user-specific training EMG data are possible. It is important to note that it is easy to measure anthropometric variables, and some variables can be measured directly together with surface EMG signals via an armband device

  MCIs. It could semi-automatically/automatically adapt the setting of the EMG recognition systems to a wider population. It can be used in two different ways: 1) a normalized value for EMG features, and 2) a weighting factor for a classifier. Moreover, wearable and wireless EMG devices are already commercially available, e.g. MYO (see www.thalmic.com/myo/), a n d h a s b e e n d e v e l o p e d r a p i d l y , t h u s using anthropometric variables to calibrate the EMG recognition system would become more practical and important in the near future.

Table 1 .

 1 . Twelve anthropometric variables of 10 male subjects (M1-M10) with the mean and the standard deviation (SD) of each variable.

	Variables	Subjects	M1 M2 M3 M4 M5 M6 M7 M8 M9 M10 Mean SD
	Body weight	49.0 62.0 65.0 74.0 73.0 58.0 57.0 55.0 63.0 54.0 61.0 8.1
	Stature		166.5 170.0 173.0 177.0 172.0 170.0 170.0 167.5 167.0 164.0 169.7 3.7
	BMI		17.7 21.5 21.7 23.6 24.7 20.1 19.7 19.6 22.6 20.1 21.1 2.1
	Biceps circumference	23.1 26.2 27.5 32.6 32.4 30.7 26.1 23.6 26.7 25.6 27.5 3.4
	Forearm circumference 22.5 22.4 25.2 29.1 28.2 26.8 23.4 23.1 24.3 22.7 24.8 2.5
	Hand breadth	7.4 7.2 8.4 9.1 8.9 12.6 11.7 7.5 7.8 7.6 8.8 1.9
	Hand length	17.2 17.8 15.5 18.2 18.9 19.1 19.3 17.4 18.3 17.1 17.9 1.1
	Elbow-hand grip length 34.4 35.1 37.6 39.3 39.7 33.9 33.8 36.0 38.5 38.4 36.7 2.3
	Elbow-fingertip length	46.8 49.7 52.2 52.2 52.6 46.2 46.2 47.2 50.0 49.7 49.3 2.5
	Shoulder-elbow length	37.7 33.6 36.9 37.8 38.2 35.7 38.4 35.4 38.3 36.6 36.9 1.6
	Forward grip reach	75.6 74.5 80.7 89.0 80.4 74.3 75.7 72.6 83.5 78.1 78.4 5.0
	Bi-deltoid breadth	39.3 46.4 44.6 54.7 47.3 37.4 30.4 42.1 45.9 47.4 43.6 6.7

Table 2 .

 2 Twelve anthropometric variables of 10 female subjects (F1-F10) with the mean and the SD of each variable.

	Variables	Subjects	F1	F2	F3	F4	F5	F6	F7	F8	F9 F10 Mean SD
	Body weight	47.0 45.0 53.0 46.0 54.0 45.0 43.0 56.0 50.0 49.0 48.8 4.4
	Stature		150.0 160.0 156.0 155.0 160.0 146.0 159.0 163.0 167.0 162.0 157.8 6.3
	BMI		20.9 17.6 21.8 19.1 21.1 21.1 17.0 21.1 17.9 18.7 19.6 1.8
	Biceps circumference	24.4 19.8 27.6 24.2 24.8 25.5 22.6 24.9 22.3 21.5 23.8 2.2
	Forearm circumference 21.6 19.1 22.9 21.5 22.2 21.9 19.7 22.7 20.6 20.9 21.3 1.2
	Hand breadth	7.4 6.4 7.3 7.5 6.9 6.6 6.9 8.1 7.4 7.2 7.2 0.5
	Hand length	15.6 16.7 16.8 16.6 16.4 14.7 16.1 15.2 17.1 17.2 16.2 0.8
	Elbow-hand grip length 32.8 35.0 34.8 34.6 34.4 32.6 34.3 35.2 38.4 36.2 34.8 1.6
	Elbow-fingertip length	43.3 46.8 46.4 46.4 46.2 39.1 46.2 47.1 49.0 47.7 45.8 2.8
	Shoulder-elbow length	34.9 35.3 36.1 35.8 33.1 34.2 33.0 37.1 37.8 36.0 35.3 1.6
	Forward grip reach	65.1 72.0 60.8 68.1 75.7 65.2 70.2 63.3 74.4 70.6 68.5 4.9
	Bi-deltoid breadth	41.9 36.8 43.4 41.2 41.0 37.8 36.9 37.0 37.2 38.7 39.2 2.5

Table 3 .

 3 Correlation coefficients |r| between anthropometric variables and EMG features (WAMP, MYOP, SSC, and ZC) in cases of strong and significant relationships (at least 4 movements for a muscle) based on 10 males and/or 10 females.

	Feature Gender	Anthropometric variable	Muscle position	Movements	Average |r| (min-max)
	WAMP Male	Bi-deltoid breadth	ECU	FP, WE, WF, WU, HC	0.75 (0.72-0.78)
		Female Bi-deltoid breadth	BB	FP, FS, WF, WR	0.72 (0.69-0.77)
	MYOP Male	Bi-deltoid breadth	ECU	FP, WE, WF, WU, HC	0.77 (0.76-0.77)
		Female Bi-deltoid breadth	BB	FP, WF, WR, WU	0.72 (0.68-0.81)
	SSC	Female Bi-deltoid breadth	BB	FP, FS, WF, WU, HC	0.80 (0.70-0.86)
	ZC	Male	Biceps circumference BB	FP, FS, WF, HO	0.76 (0.69-0.85)
		Female Bi-deltoid breadth	BB	FP, FS, WF, WU, HC	0.77 (0.70-0.84)
	From Table

Table 4 .

 4 Correlation coefficients |r| between anthropometric variables and EMG features (DFA, HG, MFL, DAMV, DASDV, WL, VAR, IEMG, MAV, and RMS) in cases of strong and significant relationships (at least 4 movements for a muscle) based on 10 males and/or 10 females.Features in Table4can be divided into three feature groups based on the explanation above. Due to a large difference of feature spaces of the fractal features c o m p a r e d t o the magnitude features, anthropometric variables that have strong associations w i t h f r a c t a l f e a t u r e s a r e d i f f e r e n t t o o t h e r magnitude features. DFA computed from the FCR muscle has a strong relationship with the stature and elbow-fingertip length for female subjects. On the other hand, HG computed from the EDC muscle has a strong relationship with the BMI for male subjects. Both DFA and HG features estimate the fractal dimension of surface EMG signals, so the selection of fractal features is dependent on the gender of target users. On average DFA is more accurate than HG about 2%-5%

	Feature Gender	Anthropometric variable	Muscle position	Movements	Average |r| (min-max)
	DFA	Female Stature	FCR	FP, WF, WR, WU, HO	0.74 (0.67-0.83)
		Female Elbow-fingertip length FCR	WF, WU, HO, HC	0.73 (0.68-0.79)
	HG	Male	BMI	EDC	FP, WU, HO, HC	0.72 (0.67-0.81)
	MFL	Male	Bi-deltoid breadth	ECU	FP, WE, WF, WU, HC	0.75 (0.70-0.78)
		Both	Bi-deltoid breadth	BB	WE, WR, WU, HC	0.69 (0.68-0.69)
		Female Bi-deltoid breadth	BB	FP, FS, WF, WR, WU, HC	0.76 (0.70-0.78)
	DAMV Male	Bi-deltoid breadth	ECU	FP, WE, WF, WU, HC	0.72 (0.67-0.75)

Table 5 .

 5 Correlation coefficients |r| between anthropometric variables and EMG features (SampEn, ApEn) in cases of strong and significant relationships (at least 4 movements for a muscle) based on 10 males and/or 10 females.

	Feature Gender	Anthropometric variable	Muscle position	Movements	Average |r| (min-max)
	SampEn Female Forward grip reach	FCR	FP, FS, WU, HO, HC	0.75 (0.69-0.86)
	ApEn	Female Forward grip reach	FCR	FP, FS, HO, HC	0.76 (0.67-0.89)

KEY TERMS AND DEFINITIONS