Improved Transductive Support Vector Machine for a Small Labelled Set in Motor Imagery-Based Brain-Computer Interface

Xu, Yilu; Hua, Jing; Zhang, Hua; Hu, Ronghua; Huang, Xin; Liu, Jizhong; Guo, Fumin

doi:https://doi.org/10.1155/2019/2087132

Computational Intelligence and Neuroscience

On this page

Abstract Introduction Methods Results Discussion Conclusion Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 2087132 | https://doi.org/10.1155/2019/2087132

Improved Transductive Support Vector Machine for a Small Labelled Set in Motor Imagery-Based Brain-Computer Interface

Yilu Xu,^1,2Jing Hua,²Hua Zhang,¹Ronghua Hu,¹Xin Huang,^2,3Jizhong Liu,¹and Fumin Guo¹

Academic Editor: Paolo Gastaldo

Received29 Jun 2019

Revised20 Sept 2019

Accepted10 Oct 2019

Published25 Nov 2019

Abstract

Long and tedious calibration time hinders the development of motor imagery- (MI-) based brain-computer interface (BCI). To tackle this problem, we use a limited labelled set and a relatively large unlabelled set from the same subject for training based on the transductive support vector machine (TSVM) framework. We first introduce an improved TSVM (ITSVM) method, in which a comprehensive feature of each sample consists of its common spatial patterns (CSP) feature and its geometric feature. Moreover, we use the concave-convex procedure (CCCP) to solve the optimization problem of TSVM under a new balancing constraint that can address the unknown distribution of the unlabelled set by considering various possible distributions. In addition, we propose an improved self-training TSVM (IST-TSVM) method that can iteratively perform CSP feature extraction and ITSVM classification using an expanded labelled set. Extensive experimental results on dataset IV-a from BCI competition III and dataset II-a from BCI competition IV show that our algorithms outperform the other competing algorithms, where the sizes and distributions of the labelled sets are variable. In particular, IST-TSVM provides average accuracies of 63.25% and 69.43% with the abovementioned two datasets, respectively, where only four positive labelled samples and sixteen negative labelled samples are used. Therefore, our algorithms can provide an alternative way to reduce the calibration time.

1. Introduction

A brain-computer interface (BCI) system can allow people to communicate directly with electronic equipment using their brain activity and without using their peripheral nerves and muscles [1]. In a noninvasive BCI system, electroencephalogram (EEG) signals are used to measure brain activity due to their safety and convenience [2]. In this paper, we focus on EEG signals of motor imagery (MI), which are invoked by either real or imagined movements of feet, hands, or tongue [3]. An MI-based BCI system is suitable for use in military, entertainment, and rehabilitation engineering systems.

However, due to the inherent nonstationarity of EEG signals, long and tedious calibration time is one of the key issues preventing broad use of MI-based BCI [4, 5]. Reducing the calibration time without loss of accuracy is a major challenge. To solve this problem, semisupervised learning (SSL) classifiers can use a small labelled set and a relatively large unlabelled set from the same subject for training.

In general, SSL classifiers can be categorized into generative, self-training, cotraining, graph-based, and transductive support vector machine (TSVM) models. A generative model iteratively uses the expectation maximization (EM) technique to build a probabilistic model with the aid of labelled and unlabelled data. Nevertheless, a generative model emphasizes that the labelled data must follow the Gaussian distribution [6]. A self-training model selects a supervised learning classifier as the base learner, which is retrained continually using the initial labelled data and the unlabelled data with high confidence [7–10]. Likewise, in the cotraining model, two supervised learning classifiers are iteratively trained using the other classifier’s previous classification results [11, 12]. The accuracies of the self-training and cotraining models decrease when the unlabelled data are assigned incorrect labels. A graph-based model constructs a weighted graph to explore the manifold structure behind the labelled and unlabelled data [13–16]. However, it is difficult to develop a good graph in general situations. The TSVM model learns the decision boundary going through low-density regions and maximizes the margin between different clusters using the labelled and unlabelled data [17]. Nevertheless, the TSVM model may converge to a local optimum because of the nonconvex optimization problem. Thus, each SSL model has clear disadvantages.

In a BCI system, support vector machine (SVM) has been commonly used with small, nonlinear, high-dimensional EEG-labelled sets [7, 8]. Therefore, we pay more attention to the TSVM model, which originated from SVM [18, 19]. TSVM-light was an early implementation of the TSVM model, which was used to determine the maximum margin by switching different labels for a pair of unlabelled data during each iteration [18]. However, there is a nonconvex optimization problem in TSVM-light due to the nondifferentiability of the Hinge loss function on the unlabelled samples. To tackle this drawback, concave-convex procedure (CCCP) was used to decompose the optimization problem into its concave and convex parts [20]. However, CCCP could not scale well with larger datasets. Robust TSVM (RTSVM) provided higher computational efficiency for millions of samples by using the stochastic gradient (SG) method to solve the primal optimization problem [21]. Due to insufficient domain knowledge, it remains challenging for the TSVM model to provide high accuracy when used with obscure unlabelled data [22]. Based on manifold assumption, the graph-based model can be used with a large unlabelled set to describe the global distribution of the data. Recently, many graph-based semisupervised SVM (S3VM) classifiers were studied extensively in the literature, such as spatial-spectral label propagation based on SVM (SS-LPSVM) [23] and TSVM based on active learning (AL) and graph (TSVM_AL+graph) [24]. SS-LPSVM and TSVM_AL+graph formulated information on the manifold structure using the Laplacian regularization term, which was added to the objective function in SVM. However, it was difficult to determine the optimal parameters using cross validation under the condition of small labelled sets. Consequently, such important parameters were always defined empirically. Semisupervised classification with low-density separation (LDS) was used to transform the original features of all samples into the geometric features [25]. Despite this advancement, the transformation procedure may omit important original information.

Moreover, it may be unreasonable to preset the ratio of positive to negative samples in the unlabelled set to be equal to the ratio in the labelled set in many TSVM methods [20, 21, 24, 25], especially when the small labelled set is extremely unbalanced. Incorrect estimation of this ratio may decrease the classification accuracy. To address this problem, Zhang designed a robust S3VM method via ensemble learning, where various distributions of the unlabelled set were considered [26].

Feature learning is as important as classifier learning in a BCI system. The common spatial patterns (CSP) method is commonly used with EEG signals because CSP can provide efficient feature extraction and dimension reduction [27, 28]. However, CSP is a supervised feature learning method. A limited labelled set may result in an unreliable CSP transformation matrix, which can directly affect the accuracy of feature vectors in all samples and consequently decrease the classification accuracy. To solve this problem, Li introduced an S3VM method based on the self-training model, in which feature learning and classifier learning were performed jointly and iteratively. In this method, the CSP transformation matrix and SVM classifier were successively updated by exploiting the initial labelled data and all or part of the unlabelled data with new labels learned during the previous iteration [7]. Similarly, many self-training and cotraining methods classify EEG signals using different supervised algorithms as the base learners, such as linear discriminant analysis (LDA), Bayesian LDA (BLDA), biomimetic pattern recognition (BPR), or sparse representation (SR) [8–12].

Motivated by the aforementioned studies, we formulate an improved TSVM (ITSVM) method by combining the TSVM model with a graph-based model. In this method, we construct the variation of a weighted graph as proposed by Chapelle [25] in order to explore the potential distribution of all samples in a semisupervised way. Then, we introduce a comprehensive feature for each sample, which consists of its CSP feature and its geometric feature. In addition, we use CCCP to solve the nonconvex optimization problem. Inspired by Zhang [26], in order to determine the unknown distribution of the unlabelled set, we impose a new balancing constraint that considers various possible distributions of the unlabelled set. As mentioned above, feature learning is critical for the BCI system. Thus, we develop an improved self-training TSVM (IST-TSVM) method that can execute CSP and our proposed ITSVM method jointly and iteratively. The contributions of our work are summarized as follows:(1)We propose an ITSVM method that can maximize the margin between different clusters and provide different views of all samples based on their CSP and geometric features.(2)In contrast to the traditional definition, we impose a new balancing constraint on the optimization problem in TSVM to address the unknown distribution of the unlabelled set.(3)Most existing self-training methods adopt supervised methods as the base learners. Here, we present an IST-TSVM method based on our confidence criterion and semisupervised ITSVM approach to utilize the unlabelled data in feature and classifier learning.(4)We performed extensive experiments to evaluate the efficiency of our proposed algorithms using small labelled sets with balanced or unbalanced classes. In particular, IST-TSVM outperforms the competing TSVM methods when used with extremely unbalanced labelled sets.

The remainder of this paper is structured as follows. In Section 2, the TSVM model is briefly reviewed and the details of our two improved TSVM methods are described. The effectiveness of our proposed methods, using two famous MI-based BCI competition datasets, is evaluated in Section 3. A discussion of the experimental results is presented in Section 4. Finally, our conclusions are drawn in Section 5.

2. Methods

2.1. TSVM Model

Consider a dataset with L labelled samples and U unlabelled samples . The labelled samples are initially assigned binary class labels .

TSVM aims to identify the optimal hyperplane that separates the labelled and unlabelled samples with maximum margin. The linear hyperplane can be characterized by , where is the normal of the hyperplane and is a bias term. Compared with SVM, TSVM minimizes the cost function by adding an “effect term” for each unlabelled sample as follows [18]:where () is a user-specified parameter that can punish misclassified labelled or unlabelled samples. The slack variables are defined to handle inseparable data. Equation (1) can be rewritten as an unconstrained minimization problem:where is the loss function. Earlier implementations of the TSVM model adopted different loss functions. TSVM-light used the classical Hinge loss function . However, the nondifferentiability of on the unlabelled samples produces a nonconvex optimization problem. Thus, was replaced by in the LDS method [25]. However, most TSVM methods employed a symmetric ramp loss function for unlabelled samples [20, 21].

2.2. ITSVM Algorithm

Classification of EEG signals is a difficult task, especially when the labelled samples are sparse and unbalanced. In this paper, we first propose an ITSVM method that involves two stages. Specifically, in the first stage, we generate the comprehensive features for all samples based on their CSP features and geometric features to provide different views of the data. In the second stage, we use CCCP to solve the nonconvex loss function and define a new balancing constraint to adapt to the unknown distribution of the unlabelled set. Figure 1 shows a flowchart illustrating the signal processing of the EEG signals, where CSP and ITSVM are successively employed for feature learning and classifier learning.

2.2.1. Comprehensive Feature

CSP is an effective feature extraction method in a BCI system. Given an initial labelled set and the remaining unlabelled set , is the ith EEG sample that was already preprocessed. As shown in Figure 1, CSP uses the labelled samples to calculate the CSP matrix, which can be used to maximize discrimination of the two classes of EEG signals. For a given CSP matrix , the mapping of is defined as the new time series . Note that consists of m pairs of spatial filters. Then, element in the CSP feature vector for is defined as follows:where is the variance of the pth row of . Although CSP is robust against noise, a small labelled dataset may produce an unreliable CSP matrix, which directly influences the correctness of CSP features for all samples. Thus, it is valuable to explore the inherent spatial distribution with the assistance of unlabelled samples using graph-based SSL approaches.

In this paper, the original CSP features are converted into geometric features based on LDS [25]. LDS is used to build the nearest neighbour graph, and multidimensional scaling (MDS) is used to produce a new graphic representation of the data in a small number of dimensions [29]. In contrast to LDS, we replace the Euclidean distance with the cosine distance to measure the pairwise distance between two samples. The cosine distance can be used to correct inconsistencies in measurement standards that may be caused by high intersession variability among EEG signals. Moreover, it is assumed for LDS that two samples lying close to each other might belong to the same class. Then, LDS is used to calculate the shortest path between two samples in an unsupervised way. However, it is difficult to assess the classes of two samples if their shortest path has different classes of labelled samples. To overcome this problem, we build the nearest neighbour graph in a semisupervised way by maximizing the edge length between two labelled samples with different classes. The process for computing each geometric feature for the ith sample can be described as follows: Step 1: the pairwise distance between the ith and jth samples is initially weighted as follows: where and are the lengths of and , respectively. The last term in equation (4) is the cosine distance. Step 2: a fully connected graph with edge lengths is constructed as follows: Step 3: Before using to compute the shortest path length based on Dijkstra’s algorithm [30], we manually set . Therefore, it is impossible for labelled samples with different classes to exist along the shortest path. Step 4: the matrix G of minimal squared -path distances is defined as follows: Step 5: the positive eigenvalues and corresponding eigenvectors of are calculated using MDS, where . is the element of identity matrix . Both H and I are matrices [29]. Step 6: element in the new graph-based representation is defined as where and .

In our opinion, it is important to combine CSP features with geometric features simultaneously by considering the consistency and complement of different features. First, discriminative information in the CSP features may be reduced using a small labelled set. However, the global distribution of geometric features may not be sufficiently reliable, because it is obtained from a small labelled set and a large unlabelled set. Therefore, we define a new comprehensive feature , which is a combination of the CSP feature and geometric feature with equal weight. Both and are column vectors.

2.2.2. A New Balancing Constraint

To prevent all unlabelled samples from being assigned to the same class, it is assumed in LDS that the labelled and unlabelled samples have the same ratio of positive to negative samples by adding the following balancing constraint on the minimization problem in equation (1):

Many TSVM methods follow this idea in LDS [20, 21, 24]. However, one problem is that the distribution of a limited labelled set cannot always represent that of a large unlabelled set, especially when the existing labelled set is unbalanced.

To address this problem, Zhang trained diverse base learners based on different hypotheses regarding the distribution of positive and negative unlabelled samples; an ensemble method based on clustering evaluation means was proposed [26]. Figure 2 shows a binary classification problem and two possible classification consequences.

(a)

(b)

(c)

As illustrated in Figure 2(a), the larger solid circle and square denote labelled samples in two different classes. The extra dots are the unlabelled samples. Zhang constructed a set of base learners based on various disturbance factors that were correlated with the ratio of positive to negative unlabelled samples; the ratio ranged from 1 : 9 to 9 : 1. Figures 2(b) and 2(c) illustrate different results for the two base learners. As shown in Figure 2, neither of these two classification results is satisfactory. As a result, Zhang used k-means to cluster the diverse base learners and employed the clustering evaluation index to evaluate the clustering effect [26].

In general, training multiple base learners is time-consuming. Thus, we attempt to exploit a simple method that considers all possible distributions of the unlabelled set. We assume that is the average ratio of the positive samples to all samples in the unlabelled set:

The number of positive unlabelled samples varies from 0 to , which could cover all cases. If each ratio is equally weighted, the value of is 0.5. Thus, the average ratio of positive to negative unlabelled data is 1 : 1, regardless of the distribution of the labelled samples. Therefore, we modify the balancing constraint using comprehensive features as follows:

2.2.3. Description of ITSVM

As depicted in Figure 1, in ITSVM, we use CCCP to solve the nonconvex optimization problem under a new balancing constraint after generating the comprehensive features for all samples.

In CCCP, the cost function given in equation (2) can be decomposed into convex and concave parts: . In addition, the concave part is approximated by its tangent . CCCP employs a ramp loss function for the labelled samples and a symmetric ramp loss function for the unlabelled samples, where . It is clear that is a clipped version of . is a user-defined parameter, which defines where is clipped.

Each unlabelled sample is duplicated when using . Each original unlabelled sample and the corresponding duplicated sample are assigned a positive or negative label, respectively, as follows:

By using and , the convex and concave parts of can be written as

The minimization problem can be reformulated as follows by calculating the derivative of the concave part with respect to :where

In CCCP, the balancing constraint in equation (8) can be applied to the minimization problem by introducing an extra sample and a variable [20, 21].

In ITSVM, the comprehensive features are used as the input of CCCP. Therefore, we replace the original features in equations (13) and (14) with the comprehensive features . The new balancing constraint in equation (10) can be achieved by defining a sample and a variable .

ITSVM can converge quickly after at most five iterations. In the th iteration, the hyperplane parameter group can be updated using a dual quadratic programming (QP) solver based on the generalized sequential minimal optimization (SMO) algorithm [31]. We define a linear kernel matrix such that The pseudocode of ITSVM can be seen in Algorithm 1. More details are shown in Appendix A.

	Input: the initial labelled set: .
	Input: the unlabelled set: .
	Output: the optimal hyperplane parameter group .
	Initialize with a traditional SVM on the initial labelled set;
	Compute
	Set for and ;
	for = 1 to 5 do
	Solve the following minimization problem by using SMO to find coefficients:
	,
	subject to: 0
	;
	Update and by using equations (A.6) and (A.12) as shown in Appendix A;
	Compute
	end for

2.3. IST-TSVM Algorithm

As shown in Figure 1, the CSP matrix may be the bottleneck in ITSVM when the number of labelled samples is small. Hence, we propose an IST-TSVM method that can update the CSP matrix using the expanded labelled set.

2.3.1. Confidence Criterion

Generally, the base learners in self-training methods are supervised. Here, we use a semisupervised ITSVM as the base learner, which provides higher classification accuracy than supervised methods when the labelled set is small. We use the following confidence criterion to choose some unlabelled samples with high confidence to include in our labelled set, as this allows a tradeoff between the smallest distance to the class centre and the largest distance to the hyperplane. Note that the initial labelled set and the remaining unlabelled set are defined in Section 2.2.1.

First, the comprehensive features and their decision scores for all samples are updated by CSP and ITSVM in each iteration, where is the comprehensive feature of . Then, the positive labelled set and the negative labelled set are selected from the initial labelled set based on the signs of the decision scores. and are obtained from in the same manner.

Second, the positive class-centre and the negative class-centre are calculated using and , respectively:

Third, we define the following function for all unlabelled sets while considering the distance to the class centre and the distance to the hyperplane simultaneously:

The corresponding unlabelled sample is included as a labelled sample with higher confidence when is smaller. Therefore, we, respectively, rearrange the positive unlabelled set and the negative unlabelled set according to the values of , which are sorted in the ascending order as follows:where and are the sizes of and , respectively. To maintain the distribution of the labelled set and avoid mislabelling unlabelled samples, the first N unlabelled samples are, respectively, selected from the two reordered sets and , where . These unlabelled samples with their predicted labels are used to construct the selected unlabelled set , which yields the expanded labelled set .

2.3.2. Description of IST-TSVM

In our proposed IST-TSVM method, we iteratively use the CSP feature extraction method and a semisupervised ITSVM classifier. We define at most five iterations. More details are presented in Algorithm 2.

	Input: the initial labelled set: .
	Input: the remaining unlabelled set: .
	Output: the labels assigned to the unlabelled samples: .
	Initialize the expanded labelled set: ;
	for = 1 to 5 do
	Compute the CSP matrix with the expanded labelled set ;
	Calculate the CSP features of all samples using ;
	Generate the comprehensive features of all samples, as described in Section 2.2.1;
	Perform the ITSVM classifier given in Algorithm 1 to obtain the optimal parameter group using the comprehensive features of all samples;
	Calculate temporal labels of all samples based on and ;
	if then
	or
	then
	break;
	end if
	end if
	Update the labels of the unlabelled samples
	Construct the selected unlabelled set based on the confidence criterion in Section 2.3.1
	Expand the labelled set: ;
	end for

In Algorithm 2, if the predicted labels of the unlabelled samples do not change or the classification error rate of the initial labelled set increases by more than in the current iteration compared to the previous iteration, then the loop will be terminated in advance.

3. Experiments and Results

In this section, two well-known BCI competition datasets for MI are used to evaluate and compare the accuracies of our proposed approaches with SVM and classical TSVM classifiers.(i)SVM is a traditional supervised learning classifier. A Gaussian kernel is often used in BCI systems for nonlinear SVM [8, 19]. Considering the computational load of the optimal kernel parameters, a linear form is chosen for SVM and the following TSVM algorithms in our study. The version of SVM used is SVM^light [32].(ii)TSVM-light is an efficient transductive learning method. This method switches the different labels of a pair of unlabelled data and solves the optimization problem in equation (2) with a dual solver during each training iteration [18].(iii)RTSVM is used to solve the primal optimization problem with an SG method. There are three parameters to be preset. Parameter s is used in the ramp loss function. Parameters C₁ and C₂ denote the punishment factors for labelled and unlabelled samples, respectively. We use the default selection of s = −0.2 and C₁ = C₂ = 4, as suggested by the authors [21].(iv)LDS is a graph-based TSVM approach. Like RTSVM, LDS performs a gradient descent on the primal formulation. The parameters for and are empirically set to 1 and , respectively [25].(v)CCCP minimizes the cost function in the dual space by using a dual QP solver [20]. Like RTSVM, CCCP and our algorithms contain parameters s, C₁, and C₂, which are preset to −0.2, 2, and 2, respectively.(vi)To compare the balancing constraint in equation (8) used by CCCP with the one in equation (10) used by ITSVM, we propose a method named as CCCP1, which is equivalent to CCCP, except that the original features are applied to the new balancing constraint as follows: . Therefore, is set to 0 in CCCP1.

The purpose of our experiments is threefold. First, the small labelled sets are used to verify the effectiveness of SVM and all TSVM approaches. Second, the balanced and unbalanced labelled sets are used to evaluate the robustness of different classifiers. Because CCCP1 is entirely equivalent to CCCP under the condition of the balanced labelled sets, we only discuss the classification performance of CCCP1 under the condition of the unbalanced labelled sets. Third, we analyse their performance in terms of computation time.

3.1. EEG Datasets

(i)Dataset IV-a in BCI competition III: the dataset was recorded from five healthy subjects (aa, al, av, aw, and ay) with a total of 118 electrodes [33]. The dataset only contained data from four initial sessions without feedback. Each subject was shown visual cues for 3.5 s and performed three MI tasks: moving the left hand, right hand, or right foot. Only the latter two MI tasks were provided in the competition. For each subject, each MI task consisted of 140 trials. The presentation of target cues was interrupted by periods of random length ranging from 1.75 to 2.25 s, in which the subject could relax. The EEG signals were band-pass filtered between 0.05 and 200 Hz and downsampled from 1000 Hz to 100 Hz.(ii)Dataset II-a in BCI competition IV: data in this set were collected from nine subjects [34]. At the beginning of a trial, a fixation cross was shown on a black screen. Each subject then executed the desired MI tasks as directed by the visual cue in the form of an arrow pointing either to the left, right, down, or up (corresponding to moving the left hand, right hand, foot, or tongue). No feedback was provided. Twenty-two Ag/AgCl electrodes were used to record EEG signals, which were then sampled with 250 Hz and band-pass filtered between 0.5 and 100 Hz. In total, 72 trials per MI task were gathered from each subject on different days. To focus on the problem of binary classification, only MI EEG signals from the left and right hands were extracted for analysis.

3.2. Preprocessing

The two BCI competition datasets were preprocessed using the same methods. All raw EEG signals were band-pass filtered between 8 and 30 Hz using a fifth-order Butterworth filter. Then, the filtered signals were extracted from nonoverlapping time segments ranging from 0.5 to 2.5 s.

All classifiers in our experiments used CSP to generate their CSP features with three pairs of spatial filters. Our proposed algorithms added geometric features during classifier learning.

Data from every subject was randomly partitioned into two parts over ten repetitions. The first portion was used as the labelled set to train the classifier, while the second portion was used as the unlabelled set to verify the effectiveness of the classifier. To investigate the robustness of all algorithms with small labelled sets, we set M and R equal to the size of the labelled set and ratio of positive to negative labelled trials, respectively. We selected M from the set [10, 15, 20, 25, 30, 35, 40, 45, 50] and R from the set [1 : 4, 2 : 3, 1 :1, 3 : 2, 4 : 1].

3.3. Experiments with Balanced Labelled Sets

The ratio of positive to negative samples in the labelled set has a great effect on the performance of the classifier. Balanced and adequate labelled samples can provide higher classification accuracy, and vice versa, for unbalanced and sparse labelled samples.

For most semisupervised algorithms, more consideration is given to the number of labelled samples rather than the ratio of positive to negative labelled samples. In reality, both balanced and unbalanced labelled sets are common in classification problems. Therefore, we first conducted experiments with small balanced labelled sets.

3.3.1. Classification Performance with Small Balanced Labelled Sets

For the two BCI competition datasets, the complete set for each subject consists of an equal number of positive and negative trials. Hence, the unlabelled set is also balanced after randomly selecting the same number of positive and negative labelled trials. First, we evaluated the recognition rates for the unlabelled sets using different classifiers learned from very small and balanced labelled sets (M = 10). For each subject, the classification accuracy was taken as an average from ten repetitions. Detailed results using two datasets are given in Tables 1 and 2. The highest classification performance is written in bold.

In Table 1, IST-TSVM performs better than the others. A paired t-test shows that the result of IST-TSVM (68.07 17.62) is statistically higher than that of SVM (66.36 15.53), TSVM-light (67.09 16.13), RTSVM (64.12 18.08), and LDS (63.27 18.86) (). ITSVM provides slightly higher accuracy over CCCP. In addition, all TSVM methods are superior to SVM, except for RTSVM and LDS. Previous research with the same dataset led to the following categorization: strong: al; normal: aa, aw, and ay; weak: av [35]. As shown in Table 1, the accuracies of all algorithms for the strong subject (al) are greater than 90%. For the normal subject ay, IST-TSVM provides higher accuracy than the other algorithms. However, for the normal subject aa and the weak subject av, all classifiers provide poor results.

In Table 2, IST-TSVM stands out prominently on average for these nine subjects. IST-TSVM (70.22 19.74) exhibits a significant improvement over RTSVM (68.47 19.08) and LDS (68.14 20.23) (). ITSVM provides a 0.47% improvement over CCCP. All TSVM methods outperform SVM. Furthermore, according to the accuracy data in Table 2, the subjects can be categorized as follows: strong: A03, A08, and A09; normal: A01; weak: A02, A04, A05, A06, and A07. Likewise, for the strong subjects, the recognition rates of all methods remain considerably high. TSVM-light exhibits the highest performance for the normal subject (A01). Finally, all classifiers yield accuracies at the chance level for the weak subjects.

3.3.2. Computation Time with Small Balanced Labelled Sets

In Table 3, we list the average computation time per subject for all TSVM classifiers in order to compare their operating speeds with small balanced labelled sets as mentioned above. The lowest computation time is highlighted in bold. The algorithms were implemented with a PC running Windows 7 Professional and Matlab R2015a. This PC contained an Intel (R) Core (TM) i3-6100 CPU @ 3.70 GHz and 8 GB RAM.

In Table 3, the time spent by CCCP is close to that by RTSVM. LDS is slower than RTSVM, while IST-TSVM requires much more time than ITSVM. TSVM-light is the most time-consuming algorithm. The following reasons may lead to the different running times. First, the framework of RTSVM is similar to that of CCCP. However, an SG method is used in RTSVM, while a dual solver is used in CCCP [20, 21]. Like RTSVM, a similar optimization strategy is pursued in LDS. However, LDS requires more time to compute the shortest paths for all pairs of samples. Based on CCCP, our proposed ITSVM spends more time calculating geometric features in a semisupervised way. By iteratively performing CSP and ITSVM, IST-TSVM exhibits higher accuracy at the cost of longer computation time. TSVM-light requires much more time because only one pair of unlabelled samples is switched to retrain the SVM in each iteration.

3.3.3. Classification Performance with Varying Sizes of the Balanced Labelled Sets

We also selected balanced labelled sets with different sizes to search for more convincing results. The average classification accuracies for all subjects are plotted in Figures 3(a) and 3(b), where the numbers of labelled trials on dataset IV-a from BCI competition III and dataset II-a from BCI competition IV, respectively, are variable. The horizontal axis presents different values of M in intervals of ten trials.

(a)

(b)

As shown in Figure 3(a), IST-TSVM outperforms the others when the number of labelled trials is less than 30. However, TSVM-light provides high results as the number of labelled trials increases. As shown in Figure 3(b), IST-TSVM performs better than the other algorithms in most instances. The accuracy of TSVM-light is less than that of SVM. As shown in Figures 3(a) and 3(b), SVM is superior to RTSVM and LDS in terms of accuracy. In addition, ITSVM provides slightly higher recognition rates than CCCP when the number of labelled trials is less than 20. For all classifiers, the classification accuracies improve as the number of labelled trials increases if the labelled sets are balanced.

3.4. Experiments with Unbalanced Labelled Sets

3.4.1. Classification Performance with Small Unbalanced Labelled Sets

In most datasets, the number of positive labelled samples is often equal or similar to the number of negative labelled samples. However, we do not rule out some cases. For example, the labelled set is not always balanced in the process of online training. In Tables 4 and 5, for each subject in the two datasets, we set M and R to 20 and 1 : 4, respectively. To compare different balancing constraints, the results of CCCP1 are also shown in Tables 4 and 5.

Table 4 shows that our proposed algorithms perform better than the other algorithms, even using extremely small unbalanced labelled sets. Paired t-test results show that IST-TSVM (63.25 18.02) provides higher accuracy than that of SVM (50.82 5.45), RTSVM (49.25 3.01), and LDS (48.180.67) (). Compared to SVM, TSVM-light provides higher accuracy in most instances. For the strong subject al, IST-TSVM exhibits the highest accuracy. For the normal subject aa and the weak subject av, all classifiers produce results with low accuracy due to the unbalanced labelled sets and inherent characteristics of the subjects.

Similarly, one can see the advantage of IST-TSVM in Table 5. Paired t-test results reveal a clear difference in the accuracy of IST-TSVM (69.43 20.57) and that of RTSVM (53.67 10.07) (). ITSVM performs moderately better than CCCP for seven out of nine subjects. Compared to Table 2, the average accuracies for SVM, RTSVM, and LDS decrease abruptly to approximately 55%. For most subjects, TSVM-light provides relatively higher recognition rates than SVM. Regarding the strong and normal subjects, the accuracies of CCCP and our methods are very high with the small and extremely unbalanced labelled sets. All classifiers produce an accuracy near 50% for the weak subjects.

Overall, the results in Tables 4 and 5 show that IST-TSVM can be used to differentiate strong and weak subjects with extremely small unbalanced labelled sets. In addition, ITSVM provides greater accuracy than CCCP and CCCP1. For these two datasets, the accuracy of CCCP is close to that of CCCP1. CCCP performs moderately better than CCCP1 for two out of five subjects in dataset IV-a and for three out of nine subjects in dataset II-a. The average accuracy of CCCP1 is slightly lower than that of CCCP.

3.4.2. Classification Performance with Varying Sizes and Distributions of Labelled Sets

We randomly selected labelled sets with different sizes and distributions (M [10, 15, 20, 25, 30, 35, 40, 45, 50] and R [1 : 4, 2 : 3, 3 : 2, 4 : 1]) from the two BCI datasets. First, to further evaluate the effect of different balancing constraints, the accuracies of CCCP and CCCP1 are shown in Tables 6 and 7.

In Table 6, for each R, CCCP1 performs slightly better than CCCP no less than four times. In Table 7, CCCP1 provides slightly better accuracies than CCCP no less than five times when R is 3 : 2 or 4 : 1. In total, for each R, the average accuracy of CCCP1 is equal to or slightly higher than that of CCCP after averaging nine values in the corresponding column. To evaluate the performance of more algorithms, the classification accuracies of all classifiers except for CCCP1, with varying numbers of labelled trials (M) and ratios of positive to negative labelled trials (R), are plotted in Figures 4(a)–4(h).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

As illustrated in Figures 4(a)–4(h), IST-TSVM shows compelling validity in most cases. The differences between ITSVM and CCCP, with the extremely unbalanced labelled sets (R = 1 : 4 and 4 : 1), are more apparent than those with comparatively unbalanced labelled sets (R = 2 : 3 and 3 : 2). The recognition rates provided by TSVM-light are always higher than those of SVM. However, RTSVM and LDS have lower accuracies than SVM with the extremely unbalanced labelled sets. In contrast to Figure 3, the performances of RTSVM with R = 1 : 4 and 4 : 1 are close to 50%. Nevertheless, the accuracies of RTSVM with R = 2 : 3 and 3 : 2 are much higher. In general, most TSVM methods are more suitable for the small unbalanced labelled sets as compared to the supervised SVM, except for RTSVM and LDS. Moreover, our algorithms are comparatively insensitive to the distribution of labelled trials.

4. Discussion

In this section, we discuss how various factors affect the classification performance of our proposed algorithms.

4.1. Impact of a Dual QP Solver

All TSVM methods mentioned above can be divided into two groups. The first group solves the primal optimization problem including RTSVM and LDS. The first group is suitable for large-scale datasets that contain millions of samples. However, in MI-based BCI systems, it is difficult to collect many samples for each subject due to the large invariability between sessions. According to the experimental results, it is clear that RTSVM and LDS cannot make full use of their merits with small-scale EEG datasets. In contrast, the second group (TSVM-light, CCCP, and our proposed algorithms) minimizes the cost function using a dual QP solver. CCCP can be used to overcome the nonconvex problem in TSVM-light. Thus, in most cases, CCCP performs better than TSVM-light. Moreover, because we use CCCP to solve the optimization problem, our algorithms and CCCP exhibit similar classification accuracies, as shown in Figures 3 and 4. Consequently, a dual QP solver plays an important role in enhancing the recognition rates of small-sized EEG sets.

4.2. Impact of the Comprehensive Features

In our proposed algorithms, we generate the comprehensive features for all samples by combining the CSP features with the geometric features. Under the condition of balanced labelled sets, equation (8) in CCCP is nearly equivalent to equation (10) in ITSVM, except for the use of different features. Thus, the results in Tables 1 and 2, as well as the results in Figure 3, show that the improvement provided by ITSVM compared to CCCP is attributed to the comprehensive features. Compared to CCCP, ITSVM adds the geometric features that can provide an inherent distribution of the data based on the labelled and unlabelled data. However, because the geometric features are transformed from the CSP features, they may not be sufficiently correct when the labelled set is small. Therefore, ITSVM provides a slight improvement over CCCP.

4.3. Impact of a New Balancing Constraint

To address the unknown distribution of the unlabelled set, we consider the various distributions of the unlabelled set and create a new balancing constraint. CCCP1 is equivalent to CCCP except for different values of . is set to 0 to achieve the new constraint in CCCP1. However, is set to to achieve the traditional constraint in CCCP. Therefore, as shown in Tables 4–7, the results of CCCP1 are close to those of CCCP. For each subject in the two BCI datasets, the number of positive samples is equal to the number of negative samples. Thus, the real ratio of positive to negative unlabelled samples will be 4 : 1, if the value of R is 1 : 4. However, the assumed ratio of positive to negative unlabelled samples is 1 : 4 for CCCP and 1 : 1 for CCCP1. It is clear that these two assumptions are quite different from the real distribution of the unlabelled set. As shown in Tables 6 and 7, the average accuracy of CCCP1 is equal to or slightly higher than that of CCCP under the condition of different unbalanced labelled sets. Therefore, it is feasible that we consider all possible distributions of the unlabelled set with equal weight. Moreover, following the experimental results shown in Tables 4 and 5, as well as the results in Figures 4(a), 4(b), 4(g), and 4(h), one can see that ITSVM provides higher accuracy compared to CCCP for extremely unbalanced labelled sets. We suggest that this is due to the new constraint and comprehensive features used in ITSVM.

4.4. Impact of the Confidence Criterion and Self-Training Model

For ITSVM, the unlabelled samples are only used in the classification phase. If the labelled set is small, the CSP transformation matrix may not be very reliable. Therefore, IST-TSVM uses the unlabelled samples from feature extraction to classifier learning. Overall, IST-TSVM exhibits its superiority using small labelled sets with balanced or unbalanced classes as depicted in Figures 3 and 4. In addition, IST-TSVM can be used to distinguish strong and weak subjects, as shown in Tables 4 and 5. We postulate that the combination of the confidence criterion and self-training model effectively improves the classification accuracy of IST-TSVM. Our confidence criterion selects the most useful unlabelled samples that are close to the class centre and far from the hyperplane simultaneously. However, if these unlabelled samples lead to convergence of the classification results for unlabelled samples or sharp degeneration of recognition rates of labelled samples, our self-training model will terminate the current iteration.

5. Conclusion

In summary, we introduce two improved TSVM algorithms with the goal of reducing the calibration time for BCI subjects on the premise of accurate classification in MI-based BCI systems. Our algorithms effectively incorporate a graph-based model and a self-training model into the TSVM model. To capture the inherent distribution of all samples, we use a cosine distance to measure the pairwise distance between two samples and build the nearest neighbour graph by considering the influence of labelled samples with different classes. Then, to provide different views of each sample, we combine the discriminative CSP feature with a global geometric feature embedded in the nearest neighbour graph. In addition, we replace the traditional balancing constraint with a new balancing constraint in the optimization problem to address the unknown distribution of the unlabelled set. Moreover, to make full use of unlabelled samples, we develop a confidence criterion and self-training process to iteratively retrain the CSP matrix and ITSVM classifier using the initial labelled samples and the unlabelled samples with high confidence in the IST-TSVM method. Extensive experiments show that IST-TSVM is particularly powerful and outperforms all other TSVM algorithms using small labelled sets with balanced or unbalanced classes. However, there remain opportunities for improvement. For example, there is no clear difference between ITSVM and CCCP in some cases. Thus, we will further explore the geometric characteristics of all samples in future investigations. Furthermore, in order to adapt to online MI training, we plan to develop an iterative feedback strategy with fewer unlabelled samples.

Appendix

A. Derivation of the ITSVM algorithm

In ITSVM, by using the comprehensive features, the new balancing constraint, and the following definition:the minimization problem given in equations (13) and (14) can be rewritten as follows:where

We introduce the Lagrangian variables and as follows:where , for We compute the derivatives as follows:

For simplification, an extra sample is defined. Setting and the derivatives to zero yields

Then, the minimization problem in equation (A.2) can be rewritten as

If we define and for , and note the linear kernel matrix such thatand perform

The minimization problem in equation (A.7) can be rewritten as follows:

We can extend the method to the nonlinear case by defining the kernel matrix as follows:

For simplification, we only consider the linear case. In order to obtain the optimal hyperplane parameter group , five iterations are executed in Algorithm 1. In each iteration, the bounds in equation (A.10) on the are adjusted after each update of and the coefficients are found by SMO. Then, the hyperplane normal can be updated by using equation (A.6). The hyperplane bias can be obtained by using the following constraints:

Data Availability

Two datasets were employed in this study for binary classification, which are publicly available: (1) dataset IVa, BCI competition III [33]: this dataset contains EEG signals from 5 subjects, who performed 2-class MI tasks: right hand and foot. (2) dataset IIa, BCI competition IV [34]: this dataset contains EEG signals from 9 subjects, who performed 4-class MI tasks: left hand, right hand, foot, and tongue MI. In this dataset, only EEG signals from left and right hands were used. Our code and results are availableat https://github.com/xuyilu1980/tsvm.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (Grant nos. 61861021 and 61863027), Science and Technology Planning Project of Education Department of Jiangxi Province (Grant no. GJJ170272), Advantage Technological Innovation Team Construction Project of Jiangxi Province (Grant no. 20171BCB24001), and Natural Science Foundation of Jiangxi Province (Grant no. 20171BAB201013). The authors thank LetPub (http://www.letpub.com) for its linguistic assistance during the preparation of this manuscript.

References

J. R. Wolpaw, N. Birbaumer, D. J. McFarland, G. Pfurtscheller, and T. M. Vaughan, “Brain-computer interfaces for communication and control,” Clinical Neurophysiology, vol. 113, no. 6, pp. 767–791, 2002.
View at: Publisher Site | Google Scholar
B. Blankertz, G. Dornhege, M. Krauledat, K.-R. Müller, and G. Curio, “The non-invasive Berlin Brain-Computer Interface: fast acquisition of effective performance in untrained subjects,” Neuroimage, vol. 37, no. 2, pp. 539–550, 2007.
View at: Publisher Site | Google Scholar
G. Pfurtscheller, C. Brunner, A. Schlögl, and F. H. Lopes da Silva, “Mu rhythm (de)synchronization and EEG single-trial classification of different motor imagery tasks,” Neuroimage, vol. 31, no. 1, pp. 153–159, 2006.
View at: Publisher Site | Google Scholar
I. Hossain, A. Khosravi, I. Hettiarachchi et al., “Multiclass informative instance transfer learning framework for motor imagery-based brain-computer interface,” Computational Intelligence and Neuroscience, vol. 2018, Article ID 6323414, 12 pages, 2018.
View at: Publisher Site | Google Scholar
Y. Xu, Q. Wei, H. Zhang et al., “Transfer learning based on regularized common spatial patterns using cosine similarities of spatial filters for motor-imagery BCI,” Journal of Circuits, Systems and Computers, vol. 28, no. 7, Article ID 1950123, 2019.
View at: Publisher Site | Google Scholar
M. Tkachenko and H. W. Lauw, “Comparative relation generative model,” IEEE Transactions on Knowledge and Data Engineering, vol. 29, no. 4, pp. 771–783, 2017.
View at: Publisher Site | Google Scholar
Y. Li and C. Guan, “A semi-supervised SVM learning algorithm for joint feature extraction and classification in brain computer interfaces,” in Proceedings of the International Conference of the IEEE Engineering in Medicine and Biology Society, pp. 2570–2573, New York, NY, USA, September 2006.
View at: Google Scholar
M. Chen, X. Tan, and L. Zhang, “An iterative self-training support vector machine algorithm in brain-computer interfaces,” Intelligent Data Analysis, vol. 20, no. 1, pp. 67–82, 2016.
View at: Publisher Site | Google Scholar
L. Zhang, Y. Chen, X. Tan, C. He, and L. Zhang, “An improved self-training algorithm for classifying motor imagery electroencephalography in brain-computer interface,” Journal of Medical Imaging and Health Informatics, vol. 7, no. 2, pp. 330–337, 2017.
View at: Publisher Site | Google Scholar
J. Wang, Z. Gu, Z. Yu, and Y. Li, “An online semi-supervised P300 speller based on extreme learning machine,” Neurocomputing, vol. 269, pp. 148–151, 2017.
View at: Publisher Site | Google Scholar
J. Meng, X. Sheng, D. Zhang et al., “Improved semisupervised adaptation for a small training dataset in the brain-computer interface,” IEEE Journal of Biomedical and Health Informatics, vol. 18, no. 4, pp. 1461–1472, 2014.
View at: Publisher Site | Google Scholar
Y. Ren, Y. Wu, and Y. Ge, “A co-training algorithm for EEG classification with biomimetic pattern recognition and sparse representation,” Neurocomputing, vol. 137, pp. 212–222, 2014.
View at: Publisher Site | Google Scholar
A. Blum and S. Chawla, “Learning from labeled and unlabeled data using graph mincuts,” in Proceedings of the Eighteenth International Conference on Machine Learning, pp. 19–26, San Francisco, CA, USA, June-July 2001.
View at: Google Scholar
M. Zhao, T. W. S. Chow, Z. Zhang, and B. Li, “Automatic image annotation via compact graph based semi-supervised learning,” Knowledge-Based Systems, vol. 76, pp. 148–165, 2015.
View at: Publisher Site | Google Scholar
H. Gan, Z. Li, W. Wu, Z. Luo, and R. Huang, “Safety-aware graph-based semi-supervised learning,” Expert Systems with Applications, vol. 107, pp. 243–254, 2018.
View at: Publisher Site | Google Scholar
P. Qian, C. Xi, M. Xu et al., “SSC-EKE: semi-supervised classification with extensive knowledge exploitation,” Information Sciences, vol. 422, pp. 51–76, 2018.
View at: Publisher Site | Google Scholar
V. N. Vapnik, “An overview of statistical learning theory,” IEEE Transactions on Neural Networks, vol. 10, no. 5, pp. 988–999, 1999.
View at: Publisher Site | Google Scholar
T. Joachims, “Transductive inference for text classification using support vector machines,” in Proceedings of International Conference on Machine Learning, pp. 200–209, Bled, Slovenia, June 1999.
View at: Google Scholar
X. Liao, D. Yao, and C. Li, “Transductive SVM for reducing the training effort in BCI,” Journal of Neural Engineering, vol. 4, no. 3, pp. 246–254, 2007.
View at: Publisher Site | Google Scholar
R. Collobert, F. Sinz, J. Weston et al., “Large scale transductive SVMs,” Journal of Machine Learning Research, vol. 7, no. 1, pp. 1687–1712, 2006.
View at: Google Scholar
H. Cevikalp and V. Franc, “Large-scale robust transductive support vector machines,” Neurocomputing, vol. 235, pp. 199–209, 2017.
View at: Publisher Site | Google Scholar
Y. Li and Z. Zhou, “Towards making unlabeled data never hurt,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 37, pp. 175–188, 2015.
View at: Publisher Site | Google Scholar
L. Wang, S. Hao, Q. Wang, and Y. Wang, “Semi-supervised classification for hyperspectral imagery based on spatial-spectral Label Propagation,” ISPRS Journal of Photogrammetry and Remote Sensing, vol. 97, pp. 123–137, 2014.
View at: Publisher Site | Google Scholar
X. Wang, J. Wen, S. Alam, Z. Jiang, and Y. Wu, “Semi-supervised learning combining transductive support vector machine with active learning,” Neurocomputing, vol. 173, pp. 1288–1298, 2016.
View at: Publisher Site | Google Scholar
O. Chapelle and A. Zien, “Semi-supervised classification by low density separation,” in Proceedings of the 10th International Workshop on Artifical Intelligence and Statistics, pp. 57–64, Bridgetown, Barbados, January 2005.
View at: Google Scholar
D. Zhang, L. Jiao, X. Bai, S. Wang, and B. Hou, “A robust semi-supervised SVM via ensemble learning,” Applied Soft Computing, vol. 65, pp. 632–643, 2018.
View at: Publisher Site | Google Scholar
H. Ramoser, J. Muller-Gerking, and G. Pfurtscheller, “Optimal spatial filtering of single trial EEG during imagined hand movement,” IEEE Transactions on Rehabilitation Engineering, vol. 8, no. 4, pp. 441–446, 2000.
View at: Publisher Site | Google Scholar
M. Grosse, C. Liefhold, K. Gramann et al., “Beamforming in noninvasive brain computer interfaces,” IEEE Transactions on Biomedical Engineering, vol. 56, no. 4, pp. 1209–1219, 2009.
View at: Publisher Site | Google Scholar
T. F. Cox and M. A. Cox, “Multidimensional scaling,” Journal of the Royal Statistical Society, vol. 46, no. 2, pp. 1050–1057, 2001.
View at: Google Scholar
E. W. Dijkstra, “A note on two problems in connexion with graphs,” Numerische Mathematik, vol. 1, no. 1, pp. 269–271, 1959.
View at: Publisher Site | Google Scholar
J. Platt, “Fast training of support vector machines using sequential minimal optimization,” in Advances in Kernel Methods-Support Vector Learning, B. Schlköpf and C. Burges, Eds., pp. 185–208, MIT Press, Cambridge, UK, 1999.
View at: Google Scholar
T. Joachims, “Making large-scale support vector machine learning practical,” in Advances in Kernel Methods-Support Vector Learning, B. Schlköpf and C. Burges, Eds., pp. 169–184, MIT Press, Cambridge, UK, 1999.
View at: Google Scholar
Graz University, BCI Competition III Datasets Iva, Graz University, Graz, Austria, 2004, http://www.bbci.de/competition/iii/#datasetIva.
Graz University, “BCI competition IV datasets 2a,” 2008, http://www.bbci.de/competition/iv/#dataset2a.
View at: Google Scholar
A. Atyabi, M. H. Luerssen, and D. M. W. Powers, “PSO-based dimension reduction of EEG recordings: implications for subject transfer in BCI,” Neurocomputing, vol. 119, pp. 319–331, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Yilu Xu et al. 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.

PDF Download Citation

Download other formats

Order printed copies

Views

1250

Downloads

1476

Citations