Research Article  Open Access
Md. Mostafizur Rahman, Shaikh Anowarul Fattah, "Mental Task Classification Scheme Utilizing Correlation Coefficient Extracted from Interchannel Intrinsic Mode Function", BioMed Research International, vol. 2017, Article ID 3720589, 11 pages, 2017. https://doi.org/10.1155/2017/3720589
Mental Task Classification Scheme Utilizing Correlation Coefficient Extracted from Interchannel Intrinsic Mode Function
Abstract
In view of recent increase of brain computer interface (BCI) based applications, the importance of efficient classification of various mental tasks has increased prodigiously nowadays. In order to obtain effective classification, efficient feature extraction scheme is necessary, for which, in the proposed method, the interchannel relationship among electroencephalogram (EEG) data is utilized. It is expected that the correlation obtained from different combination of channels will be different for different mental tasks, which can be exploited to extract distinctive feature. The empirical mode decomposition (EMD) technique is employed on a test EEG signal obtained from a channel, which provides a number of intrinsic mode functions (IMFs), and correlation coefficient is extracted from interchannel IMF data. Simultaneously, different statistical features are also obtained from each IMF. Finally, the feature matrix is formed utilizing interchannel correlation features and intrachannel statistical features of the selected IMFs of EEG signal. Different kernels of the support vector machine (SVM) classifier are used to carry out the classification task. An EEG dataset containing ten different combinations of five different mental tasks is utilized to demonstrate the classification performance and a very high level of accuracy is achieved by the proposed scheme compared to existing methods.
1. Introduction
Electroencephalogram (EEG) signal is used extensively nowadays by the researchers to handle different applications of braincomputer interface (BCI). EEGbased BCI systems employ electrical activity of brain to classify different EEG signals corresponding to various mental tasks precisely. One way to classify the signals effectively is to acquire discriminative features from that signal. As a matter of fact, different schemes to extract distinctive features are available in literature. For example, in [1], spectral power and asymmetry ratio based feature extraction scheme is proposed where an additional band (24–37 Hz) is used along with conventional lower spectral bands for mental task classification. This method offers comparatively satisfactory classification performance but lacks consistency for all cases. In [2], similar feature extraction scheme used in [1] is proposed; however, the difference is that it utilizes an additional high frequency band (40–100 Hz) to obtain those features. In [3], a dictionary consisting of power spectral density and common spatial pattern (CSP) algorithm is introduced to classify various mental tasks. Autoregressive (AR) model based feature extraction scheme is reported in [4] where sixthorder AR system is considered to extract feature. Moreover, in [5], multivariate AR models are taken into consideration and four different representations of AR coefficients are tested to classify mental task. In [6], feature extraction scheme based on sparse autoregressive model is investigated, which involves complex computation to exclude autoregressive coefficients that are useless in the prediction stage. In [7], a feature extraction method based on generalized Higuchi fractal dimension spectrum along with AR parameters is proposed. Wavelet transform and empirical mode decomposition (EMD) based classification methods are proposed in [8], where feature selection method is utilized for better classification performance. In [9], Stockwell transform based algorithm is proposed and mean square root of standard deviation of signal after transformation is utilized as distinctive feature. In [10], cross correlation based feature extraction scheme is introduced where cross correlation is computed between two channels keeping one of them always a fixed reference channel. Since one channel is kept fixed, the effect of considering cross correlation between all channels on overall feature quality has not been investigated. Moreover, unique choice of a fixed channel depends on various reasoning. Most of the reported algorithms except [10] are performed on the data taken from various channels while the interchannel relationship has not been utilized. It is considered that, for different types of task, different channels corresponding to different parts of the brain are stimulated. Measuring interchannel relationship in some efficient spectrotemporal domains may play a significant role to cover the spatial and temporal relationship between different channels. Thus, development of a proficient method capable of detecting and classifying different types of mental tasks utilizing the interchannel relationship is still undiscovered.
One of the main objectives of this paper is to extract robust feature by utilizing interchannel relationship. Instead of directly dealing with temporal EEG data, it is expected that features extracted from decomposed EEG data will provide more consistent characteristics. In particular, in this paper, widely used EMD is utilized to obtain intrinsic mode functions (IMFs) and first three IMFs are selected. Each IMF is utilized to compute the correlation coefficient from interchannel IMF data, which is referred to as interIMFCC method in this paper. Moreover, in the proposed method, unlike [10], only the channel information of the test frame is utilized to extract correlation coefficient and no previously defined reference signal is required for that purpose. Another objective of this paper is to observe whether the classification accuracy improves if different statistical features obtained from respective IMF are used along with interIMFCC. The support vector machine (SVM) classifier is used to carry out classification process. The effect of channel selection and that of using different kernels is investigated. Simulation details are reported considering a publicly available EEG dataset on various mental tasks.
2. Data Acquisition
A widely used EEG dataset collected by Keirn and Aunon is utilized [11] in this paper. EEG signals are acquired from the locations C3, C4, P3, P4, O1, and O2 which are denoted as the 1020 international system of electrode placement. Measurements are made considering A1 and A2 as reference. Data are bandpass filtered using an analog filter with band limit of 0.1–100 Hz and sampled at 250 Hz with 12bit quantizer. The recording is carried out for ten seconds during each session. EEG signals from seven subjects performing five different mental tasks, namely, geometrical figure rotation (R), mathematical multiplication (M), mental letter composing (L), visual counting (C), and baselineresting (B) are investigated. For notational convenience, hereafter, each task is abbreviated with an alphabet as shown in the parentheses. However, data obtained from three subjects contain fewer than ten sessions or have some recording errors. Hence, like some other existing research works [1], in this paper, data from four subjects, each having ten or more sessions, are taken into consideration.
For the purpose of analysis of each tensecond session, a number of frames with shorter time interval are investigated as EEG signal is assumed to be nonstationary. In this case, onesecond frame duration is considered with 0.5second frame shift (i.e., 50% overlap between successive frames) [2], which provides reasonable number of samples (250 samples) in each frame.
3. Proposed Method
The proposed mental task classification scheme can be divided into four major steps: empirical mode decomposition, interchannel relation, feature extraction, and classification. These steps are described in detail in the following subsections.
3.1. Empirical Mode Decomposition
Due to random nature of recordings of EEG data, it is very difficult to obtain discriminative characteristic from the time domain EEG data. Therefore, instead of directly utilizing EEG data, it may be easier to extract distinctive characteristic if decomposition is imposed on EEG data. Empirical mode decomposition is found very effective as it decomposes the signals in particular patterns preserving the originality of the signal. EMD is intuitive and adaptive, with IMFs directly derived from the signal under test without changing their domains. Moreover, each IMF contains information about how the frequency of the original signal changes in time. In Figure 1, a sample EEG signal and its four IMFs obtained from counting task are plotted. It is observed that the four IMFs, obtained after employing EMD on the test EEG signal, are lesser irregular and complex in nature than the original signal and have particular patterns. As a result, it is expected that classification performance will improve if IMFs are utilized to obtain distinctive features instead of main signal.
(a)
(b)
(c)
(d)
(e)
An IMF can be defined as a function which has equal number of maxima and minima or the difference between them is at most one. Moreover, the mean value of the envelope defined by the local maxima and the local minima is zero. In what follows, a brief description of obtaining IMFs by employing EMD on the EEG signal is described.
First, all the local maxima points of EEG data , obtained from channel , are connected to define the upper envelope and all the local minima points are connected to define the lower envelope. The new signal is reconstructed aswhere is the mean value of the envelopes and is the number of samples of EEG signal, . The whole process is iterated times until an IMF signal is generated according to the definition. The first IMF, , is defined bywhere is the mean value of the envelopes at th iteration. The residue signal is found by subtracting the constructed IMF from the main signal; that is,This residue signal is considered as the main signal to estimate the next IMF. The process continues until the residue signal is either a signal consisting of a single maxima or minima or a constant value. Finally, IMFs and a residue signal are generated after performing the whole decomposition process. For level of decomposition, can be reconstructed asHere represent the IMFs.
In the proposed method, it is observed that the number of IMFs can be extracted considering any frame is four or more. It is expected that higher order IMFs which contain low frequency information may not be necessarily required while mental task is evaluated. Alternatively, it is expected to find more distinguishable characteristics in the IMFs that contain relatively high frequency information. As shown in Figure 1(e), it is clearly observed that fourth IMF contains very low frequency information. Considering all these facts, for the sake of consistency, each channel data of a given frame is decomposed into only three IMFs. However, effect of varying the number of IMFs on the classification performance is delineated next in the result section.
3.2. Interchannel Relation
In general, it is considered that, for different types of task, different spatial locations of brain, such as central, parietal, or occipital are stimulated. It is expected that data obtained from locations of the brain that are highly stimulated due to a specific type of task will be less correlated with data obtained from other less stimulated locations. For example, tasks involving visual effects are most likely to stimulate occipital regions predominantly. Therefore, EEG data obtained from the channels located in the occipital region will be significantly different from the data obtained from other less stimulated regions. Measuring interchannel relationship may play a significant role to cover this spatial and temporal relationship between different channels for a particular type of task. In the proposed method, correlation coefficient is utilized to measure interchannel relationship.
Correlation coefficient is a kind of statistical measure to quantify relationship between two or more signals. In this paper, it is utilized as a measuring tool to obtain interchannel correlation of th and th channel. Instead of directly using EEG data, correlation coefficient is obtained from the th IMF, decomposed from EEG signal. The correlation coefficient extracted from interchannel IMF data is referred to as interIMFCC in this paper. The interIMFCC obtained from th and th channel can be estimated aswhere is the th component of the covariance matrix of the th and th channel IMFs and ; each consists of samples. It is expressed asThe covariance of and denoted by is calculated considering the following formula:Here and indicate the mean of IMF data obtained from th and th channels, respectively, and denotes the complex conjugate. In the proposed method, all possible pair of th and th channels are taken into consideration to obtain interIMFCC which is expected to provide maximum utilization of channel information. However, effect of choosing lesser pairs of channels is also investigated and presented in the result section.
One of the major advantages of utilizing interIMFCC as feature is that its values are bounded, which is . If the IMF data obtained from the channels are the same, interIMFCC is one; otherwise if there is no relationship, it is zero. To investigate the differentiating quality of interIMFCC as feature, a sample experiment considering multiplication and rotation task is performed. All fifteen different combinations of six channels denoted as C3C4, C3P3, C3P4, and so forth are universally taken into consideration to measure interIMFCC. In Figure 2, the box plot corresponding to interIMFCC obtained for fifteen different combinations of channels is presented. The boxplot indicates various statistical information, such as median, 25th and 75th percentile, and outliers of interIMFCC. There are thirty boxplots in each subfigure; each boxplot represents interIMFCC measured from a particular combination of channel for a particular type of task performed by subject 1. In comparison to the boxplots presented in Figures 2(a), 2(b), and 2(c), the presence of outliers in boxplot presented in Figure 2(d) is much higher. As discussed before, higher order IMF contains very low frequency information which is less relevant to mental tasks considered here and hence poor distinctive features are expected to be extracted if 4th IMF is used. This fact is also reflected in boxplot presented in Figure 2(d). Therefore, in what follows, our discussions are restricted only for the first three IMFs.
(a)
(b)
(c)
(d)
It is observed that the values of interIMFCCs obtained for three combinations of channels, namely, C4O1, P4O1, and O1O2, are found significantly higher in case of multiplication task than that in case of rotation task. It is to be noted that for these three combinations, O1 is considered as the reference channel obtained from left hemisphere and other nonreference channels, namely, C4, P4, and O2, are from right hemisphere. One possible reason behind this observation is that, in case of rotation task, O1 channel may get more stimulated than other channels due to the fact that, in rotation task, visually observed objects are required to rotate around their axis mentally. As a result, data obtained from O1 channel is less correlated with data obtained from other channels, specially the channels located in right hemisphere in case of rotation task. This observation corroborates the hypothesis that during performing a particular type of task if any location of the brain becomes more excited than other locations, data obtained from stimulated location will be significantly different from data obtained from comparatively less stimulated locations resulting in lower interIMFCC. Moreover, it is also found that interIMFCC values are comparatively higher in case of multiplication task than that obtained in case of rotation task. It is expected that, being an arithmetic task, multiplication involves more complexity in comparison to rotation task. As a result, all locations of brain get more excited while performing multiplication task than rotation task, which in turn leads to more correlation between channels in case of multiplication task. However, location of stimulation may vary from person to person depending on the nature of task.
3.3. Feature Extraction
In the proposed method, for the purpose of feature extraction, interIMFCCs are utilized to exploit the relationship among various channels. Moreover, statistical parameters such as root mean square (RMS), standard deviation, and entropy are also included in the feature vector to represent statistical measure of IMF data obtained from various channels. RMS depicts statistical measure of numerical values of varying quantity of the data obtained from channel of corresponding IMF. For IMF data consisting of samples, RMS can be expressed asTo measure the dispersion of the IMF data around its mean value , standard deviation is proposed as a distinctive feature. Standard deviation of IMF data obtained from channel is given byFor the purpose of measuring uncertainty of the IMF data , entropy is introduced in the feature vector. Entropy is a statistical measure of randomness that is defined aswhere indicates the probability of occurrence of a particular value of IMF data of th channel and is denoted byand indicates the number of occurrences of among the number of samples of ; that is,
In brief, for the purpose of feature extraction, at first, the raw EEG signal is preprocessed with a 60 Hz notch filter. After that, the eeg data corresponding to a channel is decomposed utilizing EMD where from each channel data three IMFs are extracted. Finally, the feature vector is formed utilizing interIMFCC and statistical parameters of IMFs, such as RMS, standard deviation, and entropy obtained from each channel. For number of selected IMFs and number of channels for each IMF, number of interIMFCCs obtained is . The number of features obtained from statistical parameters of IMFs for a test frame is . Finally the total feature dimension of the proposed method is .
3.4. Classification
Classifier selection is essential to obtain satisfactory result while performing test validation of the proposed method. In the proposed method, kernel based SVM classifier is chosen to effectively classify mental tasks due to its effectiveness and acceptability in supervised classification. To generate an dimensional decision vector , features extracted from the IMF data are provided into the classifier instead of raw EEG data. The extracted features from the training dataset consisting of frames are converted from the original space to a new representative vector space to discriminate different classes more efficiently. A class label is provided for each dimensional th frame . For two class problems with two class labels and , each frame fulfills the following inequalities considering the threshold [12]:In kernel based SVM classifier, to match with class label of the training dataset, the following discriminant function is utilized to form the decision vector, which can be expressed as [12]Here is an empirical value and kernel matrix is given by
For the purpose of classification, the performance of different kernel functions in SVM classifier is observed considering various feature extraction methods. It is found that polynomial kernel based classification outperforms other kernels in terms of classification accuracy. In all calculations of the reported classification accuracies, leaveoneout cross validation scheme is employed to generate classification result. In this scheme, each frame is tested one by one; that is, when a frame is left out for testing, remaining frames are used for training. Let us consider total number of frames with number of frames belonging to class A and number of frames belonging to class B. In the leaveoneout cross validation scheme, when one of those frames is left out for testing, remaining frames are used for training. This process is repeated times. Finally, classification accuracy is defined as the percentage of correctly identifying the class of each frame. Among total number of frames if number of frames are correctly classified, the classification accuracy can be expressed as
4. Simulation and Results
In this section, performance of various feature extraction methods is investigated considering classification accuracy obtained under different conditions, such as varying the feature dimension, utilizing different statistical parameters as feature, and use of various EEG channel locations. Moreover, effect of utilizing different kernel functions of SVM classifier on classification accuracy is also analyzed. A comparative analysis on classification performance between the proposed method and some other methods is also presented.
In the proposed method, instead of directly using channel data, corresponding IMFs are used to extract interIMFCC and statistical parameters using (4)–(10). Unless otherwise specified, polynomial kernel of SVM classifier is employed in leaveoneout cross validation manner to obtain classification accuracy. The classification task is carried out considering two types of mental tasks at a time, as conventionally done by other researchers [1, 2]. In this way, ten different combinations of the five types of tasks, as mentioned in Section 2, are possible. Here, for notational convenience, each combination of tasks is denoted with two alphabets from two different tasks. For example, MC refers to a twoclass (multiplication and counting) classification problem, BL corresponds to another twoclass (baselineresting and mental letter composing tasks) classification problem. In what follows, detailed results and analyses are presented.
4.1. Effect of Variation of Number of IMFs
The number of IMFs to be used in the feature matrix directly dictates the feature dimension. It is already mentioned that higher order IMFs which contain very low frequency information are not necessary to be considered. The distinctive quality of the proposed interIMFCC feature deteriorates for 4th IMF as shown in Figure 2(d). Hence, in the proposed method, only first three IMFs are considered. In this subsection, effect of variation of number of IMFs is demonstrated on overall classification accuracy for four subjects. Here the number of IMFs is varied from 1 to 4 and different cases like extracting only one IMF (1IMF), two IMFs (2IMFs), and so forth are considered.
In Figure 3, the box plot corresponding to classification performance obtained by varying number of IMFs is presented. The sixteen boxplots indicate various statistical information, such as median, 25th and 75th percentile, and outliers of classification accuracy. Each boxplot represents classification accuracy of ten different combinations of tasks for a subject considering particular number of IMFs to be used for feature extraction. It is found that, with the increase in number of IMFs, classification accuracy becomes more consistent for each subject until number of IMF is three. Moreover, it is observed that, for all subjects, features extracted considering the first three IMFs offer the best classification accuracy with respect to all other combinations of IMFs. That is why, although all channel data of a frame can be decomposed into four or more IMFs, only three IMFs are considered to extract feature. Meanwhile, considering three IMFs rather than four or higher number of IMFs offers a reduced feature dimension.
4.2. Effect of Different Statistical Feature
In the proposed method, as mentioned in Section 3.3, some statistical parameters are used as features, which are extracted from the channel IMF data. Effect of using conventional statistical features on classification accuracy is investigated considering ten widely used higher and lower order statistical parameters, namely, average (avg), median (med), mode (mod), maxima (max), minima (min), standard deviation (std), root mean square (RMS), entropy (ent), skewness (skew), and kurtosis (kurt). For notational convenience, hereafter, each statistical feature is abbreviated as shown in the parentheses. It is to be noted that the main objective of this paper is to demonstrate the efficacy of proposed correlation feature (interIMFCC) obtained from interchannel IMFs. It is expected that the use of proposed interIMFCC feature along with the conventional statistical features of IMFs will offer better classification performance. In this regard, two different cases are considered:(1)Use of only statistical features: each statistical feature is extracted from each of three IMFs of a channel, that is, for number of channels with number of IMFs extracted from each channel, feature dimension is .(2)Use of proposed interIMFCC feature along with statistical feature: in this case, number of interchannel correlation coefficients (interIMFCC) to be obtained from channels for each IMF is . Hence, for number of channels with number of IMFs extracted from each channel, total feature dimension is .
In Figure 4, classification accuracies considering the previously discussed two cases for the ten statistical features obtained for all subjects are shown. It is observed that classification accuracy increases if interIMFCC is combined with channel statistical information of each IMF. Statistical parameters such as std, RMS, and ent of IMFs offer better classification performance compared to some higher order statistical feature, namely, skew and kurt. Moreover, features like max and min which are likely to be more biased because of the presence of noise are avoided. First order statistical parameters such as avg, med, and mod are also excluded as EEG signals are very random in nature. Due to distinctive nature of std, RMS, and ent, these three statistical parameters are finally chosen for the feature vector along with proposed interIMFCC feature to classify mental tasks.
4.3. Effect of Utilizing Kernel of SVM Classifier
The effect of using different kernels in SVM classifier on overall classification performance of the proposed method is thoroughly investigated. In order to demonstrate the performance variation due to change in kernels, three widely used kernels are considered, namely, linear, quadratic, and polynomial kernel. To observe the variation of classification accuracies for different kernels, all 10 different combinations of tasks, namely, MC, MB, ML, MR, CB, CL, CR, BL, BR, and LR, from each subject are considered and average classification accuracy of those combination of tasks are measured from four subjects. In Figure 5, average classification accuracies for 10 different combinations of tasks by using three different kernels are plotted.
It is found that, between linear and quadratic kernel, the latter offers better classification performance. However, it is observed that the classification performances of polynomial kernel are consistently better in comparison to those obtained by linear and quadratic kernels in all cases. For that purpose, polynomial kernel of SVM classifier is chosen to classify the tasks in the proposed method.
4.4. Effect of Variation of Number of Channel Pairs
In the proposed method, all possible pairs of channels are taken into consideration to obtain interIMFCC so that maximum channel information can be utilized. However, choosing lesser pairs of channels reduce feature size effectively. Reduction in feature size definitely helps in reducing computation time. Hence, effect of variation of the number of channel pairs is presented in this subsection. It is to be noted that, in [1, 2], asymmetry ratio of a pair of channels is computed considering one channel from left hemisphere and the other channel from right hemisphere. Similarly, in this paper, the effect of measuring interIMFCC considering one channel from left hemisphere and the other from right hemisphere is investigated. This investigation is performed considering counting and baselineresting task and denoted as Experiment 1 in Figure 6. Moreover, the effect of measuring interIMFCC with respect to a specific region, denoted as Experiment 2, is also observed.
Experiment 1. For three channels located in left hemisphere C3, P3, and O1 and three channels located in right hemisphere C4, P4, and O2, possible nine combinations of computing interIMFCCs are (C3, C4), (C3, P4), (C3, O2), (P3, C4), (P3, P4), (P3, O2), (O1, C4), (O1, P4), and (O1, O2).
Experiment 2. Depending on the choice of region, such as parietal, central, or occipital, to obtain reference signals, three different investigations can be performed: (a)Considering signals of parietal region as reference, eight combinations of channels for computing interIMFCCs are possible, such as (P3, C3), (P3, C4), (P3, O1), (P3, O2), (P4, C3), (P4, C4), (P4, O1), and (P4, O2).(b)Considering signals of central region as reference, eight combinations of channels for computing interIMFCCs are possible, such as (C3, P3), (C3, P4), (C3, O1), (C3, O2), (C4, P3), (C4, P4), (C4, O1), and (C4, O2).(c)Considering signals of occipital region as reference, eight combinations of channels for computing interIMFCCs are possible, such as (O1, C3), (O1, C4), (O1, P3), (O1, P4), (O2, C3), (O2, C4), (O2, P3), and (O2, P4).
In Figure 6, a comparative analysis among these experiments is presented in terms of classification accuracy. In the above two experiments, a reduced number of channel pairs are utilized and lower classification accuracy compared to the proposed method is achieved. As a result, it is not possible to select any one particular choice of reduced number of channels to obtain acceptable classification performance in all subjects.
4.5. Performance Comparison among Various Methods
With a view to comparing the classification performance, five methods referred to as PAR4, PAR5, PAR6, EF8, and EF3 have been considered. Among these five methods, three methods are based on power asymmetry ratio (namely, PAR) computed from number of spectral bands [1, 2]. Remaining two methods are based on EMD feature (namely, EF) where corresponds to number of features to be extracted from each IMF obtained by EMD decomposition [8].
In PAR methods, depending on the number of frequency bands utilized, the methods are referred to as PAR4, PAR5, and PAR6. For example, in PAR4 method, features are extracted from the four traditionally used bands, namely, delta (<4 Hz), theta (4–7 Hz), alpha (8–13 Hz), and beta (14–20 Hz) while PAR5 utilizes an additional gamma band (23–37 Hz). In PAR6, one more additional band (40–100 Hz), along with these five bands, is proposed to compute power of spectral bands and asymmetry ratios. For one pair of channels, the asymmetry ratio for each spectral band is computed as [1]where two indices and are used to correspond to channel pairs placed in the left and right hemispheres, respectively. For example, corresponds to the spectral band power of the th channel placed in the left hemisphere and corresponds to that obtained from the th channel placed in the right hemisphere. Depending on the number of channels ( and ) in each hemisphere, total number of asymmetry ratios, denoted by , can be computed for each spectral band. As a result, the feature dimension for PAR4, PAR5, and PAR6 method is + where denotes number of spectral bands considered for these methods.
On the other hand, in EF8 method, eight features are extracted from each IMF, namely, RMS, variance, Shannon entropy, LempelZiv complexity measure, central frequency, maximum frequency, skewness, and kurtosis. However, in the proposed method, the first three of these eight statistical features are employed along with the proposed interIMFCC feature. In order to better demonstrate the effect of incorporating the interIMFCC feature, another method EF3 is considered where only the first three features are used without the proposed interIMFCC feature and classification performance of the EF3 method is also compared with that of the proposed method.
For the purpose of performance evaluation, leaveoneout cross validation technique is carried out in all methods. In Tables 1–4, the classification accuracies obtained by using four different subjects are separately reported for six methods. It is found that the classification accuracies obtained from different subjects are or more in the proposed method. In all cases, it is observed that the proposed feature extraction method outperforms other existing methods reported in this paper in terms of classification accuracy. However, in some combinations of mental tasks, existing methods offer competitive classification performance with respect to proposed method. For example, in case of BR combination of subject 1 reported in Table 1, both EF8 and proposed method achieve classification accuracy. In Table 4, it is observed that the average classification accuracies obtained by PAR6 and EF8 are very comparable with those obtained by the proposed method. For all subjects, it is found that the average classification accuracy obtained for EF8 is very close to EF3 despite having a larger feature dimension. However, after adding interIMFCC along with the three parameters used in EF3, the average classification accuracy increases drastically and for subject 2 and subject 3, it increases around 7.5% from EF3. In each reported existing method, it is observed that, for various combination of mental tasks, classification accuracy varies a lot. For example, in PAR4 method, for subject 1 and subject 4, the standard deviation of classification accuracies for various subjects is found to be 8.91% and 7.75% compared to 3.21% and 2.52% of the proposed method. It is found that the classification performance obtained by the proposed method varies from subject to subject, but not at a very large scale. For subject 2, the standard deviation obtained from different combination of mental tasks is found to be 1.16 which is the least among all four subjects. It is clearly observed that the proposed method offers consistently satisfactory classification accuracy in all cases irrespective of subjects and combination of mental tasks.




4.6. Computation Time
Average computational time is measured to extract features from one test signal for six methods, namely, PAR4, PAR5, PAR6, EF8, EF3, and proposed method. The whole process of computation is performed using Intel(R) Core(TM) i54200M processor with 2.50 GHz clock speed and 4 GB ram. The feature dimension and the feature extraction time for six methods are listed in Table 5.
It is found that the proposed method uses a very small computation time for feature extraction compared to recently reported EF8 method. One of the reasons for such a small computation time for the proposed method is its feature dimension compared to EF8. For three selected IMFs and six channels for each IMF, the feature dimension of the proposed method is . On the contrary, for four selected IMFs and similar number of channels for each IMF, the feature dimension of the EF8 and EF3 method is and , respectively. In case of PAR4, PAR5, and PAR6, feature dimension is , , and , respectively. The PAR4, PAR5, and PAR6 method utilizes lesser time and features but the classification accuracies are lesser in these methods than proposed method.
5. Conclusion
In the proposed mental task classification scheme, interchannel correlation coefficient of each IMF is utilized to explore the relationship between channels, which is referred to as interIMFCC method. Moreover, intrachannel features, such as standard deviation, RMS, and entropy of each IMF, are also measured. Finally, both interchannel features and intrachannel features of each IMF are utilized to form feature vector and a quite satisfactory classification performance is achieved. It is observed that increase in feature dimension by considering more IMFs not necessarily provides better classification performance and thus only three IMFs from each channel are found sufficient. Effect of selecting different combinations of channels is also investigated and it is observed that considering all combinations of channels provides the best classification performance irrespective of the tasks or the subjects. Classification performance for various feature extraction methods is listed considering polynomial kernel and it is observed that the proposed method outperforms other methods in terms of classification accuracy. Results obtained from various types of investigation verify that the proposed mental task classification scheme is capable of classifying EEG signals with high classification accuracy.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 R. Palaniappan, “Utilizing gamma band to improve mental task based braincomputer interface design,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 14, no. 3, pp. 299–303, 2006. View at: Publisher Site  Google Scholar
 L. Zhang, W. He, C. He, and P. Wang, “Improving mental task classification by adding high frequency band information,” Journal of Medical Systems, vol. 34, no. 1, pp. 51–60, 2010. View at: Publisher Site  Google Scholar
 R. Ameri, A. Pouyan, and V. Abolghasemi, “EEG signal classification based on sparse representation in brain computer interface applications,” in Proceedings of the 22nd Iranian Conference on Biomedical Engineering (ICBME '15), pp. 21–24, November 2015. View at: Publisher Site  Google Scholar
 N. Y. Liang, P. Saratchandran, G. B. Huang, and N. Sundararajan, “Classification of mental tasks from EEG signals using extreme learning machine,” International Journal of Neural Systems, vol. 16, no. 1, pp. 29–38, 2006. View at: Publisher Site  Google Scholar
 C. W. Anderson, E. A. Stolz, and S. Shamsunder, “Multivariate autoregressive models for classification of spontaneous electroencephalographic signals during mental tasks,” IEEE Transactions on Biomedical Engineering, vol. 45, no. 3, pp. 277–286, 1998. View at: Publisher Site  Google Scholar
 D. Vidaurre, C. Bielza, and P. Larrañaga, “Classification of neural signals from sparse autoregressive features,” Neurocomputing, vol. 111, pp. 21–26, 2013. View at: Publisher Site  Google Scholar
 Q. Wang and O. Sourina, “Realtime mental arithmetic task recognition from EEG signals,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 21, no. 2, pp. 225–232, 2013. View at: Google Scholar
 A. Gupta, R. K. Agrawal, and B. Kaur, “Performance enhancement of mental task classification using EEG signal: a study of multivariate feature selection methods,” Soft Computing, vol. 19, no. 10, pp. 2799–2812, 2015. View at: Publisher Site  Google Scholar
 M. Hariharan, V. Vijean, R. Sindhu, P. Divakar, A. Saidatul, and S. Yaacob, “Classification of mental tasks using stockwell transform,” Computers & Electrical Engineering, vol. 40, no. 5, pp. 1741–1749, 2014. View at: Publisher Site  Google Scholar
 S. Siuly and Y. Li, “Improving the separability of motor imagery EEG signals using a cross correlationbased least square support vector machine for braincomputer interface,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 20, no. 4, pp. 526–538, 2012. View at: Publisher Site  Google Scholar
 Z. A. Keirn and J. I. Aunon, “A new mode of communication between man and his surroundings,” IEEE Transactions on Biomedical Engineering, vol. 37, no. 12, pp. 1209–1214, 1990. View at: Publisher Site  Google Scholar
 S. Y. Kung, Kernel Methods and Machine Learning, Cambridge University Press, Cambridge, 2014. View at: Publisher Site
Copyright
Copyright © 2017 Md. Mostafizur Rahman and Shaikh Anowarul Fattah. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.