ORIGINAL RESEARCH article
Enhancing Classification Performance of Functional Near-Infrared Spectroscopy- Brain–Computer Interface Using Adaptive Estimation of General Linear Model Coefficients
- 1Department of Mechatronics Engineering, Air University, Islamabad, Pakistan
- 2Department of Electrical and Computer Engineering, Institute of Systems and Robotics, University of Coimbra, Coimbra, Portugal
- 3Faculty of Engineering and Computer Sciences, National University of Modern Languages, Islamabad, Pakistan
In this paper, a novel methodology for enhanced classification of functional near-infrared spectroscopy (fNIRS) signals utilizable in a two-class [motor imagery (MI) and rest; mental rotation (MR) and rest] brain–computer interface (BCI) is presented. First, fNIRS signals corresponding to MI and MR are acquired from the motor and prefrontal cortex, respectively, afterward, filtered to remove physiological noises. Then, the signals are modeled using the general linear model, the coefficients of which are adaptively estimated using the least squares technique. Subsequently, multiple feature combinations of estimated coefficients were used for classification. The best classification accuracies achieved for five subjects, for MI versus rest are 79.5, 83.7, 82.6, 81.4, and 84.1% whereas those for MR versus rest are 85.5, 85.2, 87.8, 83.7, and 84.8%, respectively, using support vector machine. These results are compared with the best classification accuracies obtained using the conventional hemodynamic response. By means of the proposed methodology, the average classification accuracy obtained was significantly higher (p < 0.05). These results serve to demonstrate the feasibility of developing a high-classification-performance fNIRS-BCI.
A brain–computer interface (BCI) system bypasses the peripheral nervous system and provides means of communication for patients suffering from motor disabilities or in a persistent vegetative state using devices, such as robotic arms or other prostheses (Wolpaw et al., 2002). The brain signals are acquired either invasively or non-invasively. Although the quality of brain signals acquired using invasive methods is better than those using non-invasive methods, their acquisition entails extensive surgical risk (Wester et al., 2009). With non-invasive methods, on the other hand, there is no such risk. Non-invasive techniques include electroencephalography (EEG) (Wolpaw et al., 2002; Pfurtscheller et al., 2003; Salvaris and Sepulveda, 2010; Cong et al., 2011, 2015; Jin et al., 2011, 2014, 2015; Choi, 2013; Chen et al., 2015), functional magnetic resonance imaging (fMRI) (Enzinger et al., 2008; Sorger et al., 2009), and functional near-infrared spectroscopy (fNIRS) (Ferrari et al., 1985; Kato et al., 1993; Coyle et al., 2004, 2007; Naito et al., 2007; Naseer and Hong, 2013; Naseer et al., 2014; Noori et al., 2017). Over the course of the past decade, fNIRS-based BCI systems have been the focus of considerable research interest and discussion due to their portability, affordable cost and better temporal resolution relative to fMRI. Moreover, compared with the EEG system, they offer better spatial resolution and a superior signal-to-noise ratio (Hu et al., 2012; Hong et al., 2014). In general, fNIRS has evolved into a neuroimaging technique that has contributed to ground-breaking advances in the understanding of human brain functionality (Irani et al., 2007; Aqil et al., 2012b; Ferrari and Quaresima, 2012; Hong and Nguyen, 2014; Hong and Naseer, 2016; Hong and Santosa, 2016). fNIRS utilizes near-infrared (NI) light within the 650–1000 nm wavelength range to measure changes in the concentrations of oxygenated and deoxygenated hemoglobin [ΔcHbO(t) and ΔcHbR(t)] according to the modified Beer–Lamberts Law (Delpy et al., 1988; Villringer et al., 1993; Hoshi et al., 1994; Hoshi and Tamura, 1997). Since the introduction of the principle of NI spectroscopy by Jobsis (1977), fNIRS has been used effectively for functional and structural brain imaging as well as for BCI purposes (Naseer and Hong, 2015a,b; Nguyen et al., 2016; Zafar and Hong, 2017). The first step in fNIRS-BCI is to acquire signals from a suitable mental task. Over the past decade, the mental tasks used by fNIRS-BCI researchers have been motor imagery (MI), mental rotation (MR), mental arithmetic, music imagery, and letter padding (Zhang et al., 2011a; Ayaz et al., 2013, 2014; Khan et al., 2014; Khan and Hong, 2015). In this study, we used right-hand MI and MR as the brain activity. Generation of control commands for fNIRS-based BCI systems proceeds according to the following conventional steps: first, acquisition of the desired signals; second, removal of motion artifacts and physiological noises; third, extraction of significant information (features), usually from the hemodynamic signals’ physical properties; fourth and finally, classification of the extracted features preparatory to generation of the desired control commands. Researchers have devoted considerable efforts to the improvement and enhancement of classification accuracies for fNIRS-BCI, specifically by use of different features and classifiers (Ayaz et al., 2013, 2014; Naseer and Hong, 2013; Naseer et al., 2014; Noori et al., 2016; Qureshi et al., 2016; Khan and Hong, 2017). In this paper, we propose that features be extracted from the estimated coefficients of the general linear model (GLM).
The GLM methodology was first employed by Abdelnour and Huppert (2009) in a fNIRS-based BCI study. Since that time, multiple GLM-based fNIRS studies have been performed for noise removal and brain mapping (Hu et al., 2010; Zhang et al., 2011b, 2012; Aqil et al., 2012a; Kamran and Hong, 2013). Abdelnour and Huppert (2009) have proposed the use of filter coefficients obtained by Kalman filtering as the features for classification. They assumed that different brain activities will produce different filter coefficients, using which different signals can be classified. Similarly, recursive least square estimation (Aqil et al., 2012a), and wavelet transform (Khoa and Nakagawa, 2008; Abibullaev et al., 2011; Abibullaev and An, 2012) have also been used for brain mapping using GLM. In this study, GLM is used with least square to estimate filter coefficient values. Afterward these values are used to extract features. To the best of our knowledge, this is the first work that uses filter coefficient values to extract statistical features that can be used for classification. In the proposed methodology, signals are acquired from the left motor cortex of the brain for right-hand MI (clenching of the right hand) and rest tasks, whereas MR (rotation of rectangular box) and rest signals are acquired from the prefrontal cortex; these signals are filtered to remove physiological noises and the GLM coefficients are extracted using the least squares estimation (LSE) technique; the feature values of these coefficients are then fed to support vector machine (SVM) for classification. The motivation of using GLM-based features for fNIRS data came from Abdelnour and Huppert (2009). They showed promising results using beta (β) values extracted from GLM as features. In this study, authors have used GLM with least square to estimate β values. Afterward β values are used to extract features in order to calculate classification accuracies.
Materials and Methods
A total of 10 subjects participated in the experiments. Five subjects performed MI (right-hand clenching) versus rest, whereas the other five performed MR (rotation of rectangular box) versus rest. The reason for introducing two different experiments was to establish generalization of the proposed methodology. The subjects were each seated on a comfortable chair in front of a display screen and asked to restrict their body movements as much as possible during the experiment. Verbal consent was obtained from all of the subjects after explaining the experimental paradigm in detail. The subjects had little or no previous experience of fNIRS recording. This work was approved by the Institutional Review Board of Pusan National University. All experiments were conducted in accordance with the ethical standards encoded in the latest Declaration of Helsinki. The complete details of the baseline system (conventional methodology) can be found in Naseer et al. (2016a,b).
The first 20 s was the rest period, required in order to set up the baseline condition; it was followed by 20 s of a right-hand MI task (clenching of the right hand), followed by another 20 s rest period that allowed the signals to return to their baseline values before the start of the next trial. This pattern was repeated 11 times; the total duration of experiment for each subject, therefore, was 440 s. During the MI task, the subjects were asked to imagine clenching of their right hand with a self-paced frequency of around 1 Hz; during the rest period, they were asked to relax.
Similar to MI task in MR task the first 20 s was the rest period, required in order to set up the baseline condition; it was followed by 10 s of object rotation (rotation of rectangular box), followed by another 20 s rest period that allowed the signals to return to their baseline values before the start of the next trial. This pattern was repeated 10 times; the total duration of experiment for each subject, therefore, was 330 s. During the MR task, the subjects were asked to imagine rotation of a rectangular box; during the rest period, they were asked to relax.
In order to acquire fNIRS signals from the right motor cortex of the brain, a multi-channel continuous-wave imaging system (dynamic near-infrared optical tomography; two wavelengths: 760 and 830 nm; NIRx Medical Technologies, NY, USA) with a sampling rate of 1.81 Hz was employed. The acquired signals’ light intensities were first converted to ΔcHbO(t) and ΔcHbR(t) using the modified Beer–Lamberts Law
where ΔA(t; λj) (j = 1, 2) is the unit-less absorbance (optical density) variation of a light emitter of wavelength λj, aHbX(λj) is the extinction coefficient of HbX (HbO and HbR) in μM−1mm−1, d is the unit-less differential path length factor, and l is the distance (in millimeters) between the emitter and detector. As shown in Figure 1A, four emitters and five detectors were positioned over the left motor cortex of the brain for right-hand MI task. Figure 1B shows eight emitters and three detectors placed on prefrontal cortex of the brain region in order to acquire signals for MR, the distance between each emitter–detector pair was of 3 cm. This emitter–detector distance is in accordance with the literature (McCormick et al., 1992; Gratton et al., 2006). In order to remove physiological noises (heartbeat, respiration) from the obtained signals, the Butterworth filter of order four was used with a cut-off frequency of 0.6 Hz; for removal of low-oscillation Mayer waves, a high-pass filter with a cut-off frequency of 0.01 Hz was used (Naseer and Hong, 2015a,b).
Figure 1. Optode placement and channel location. (A) 12-channel with 4 detectors and 5 emitters on the left motor cortex and (B) 12-channel with 8 detectors and 3 emitters on the prefrontal cortex of brain region.
General Linear Model
The GLM has been very widely utilized by researchers of fNIRS-BCI systems in order to identify brain-activation patterns for multiple cognitive tasks (Abdelnour and Huppert, 2009; Hu et al., 2010; Zhang et al., 2011b, 2012; Aqil et al., 2012a; Kamran and Hong, 2013). The GLM-based methods were developed initially for fMRI-based functional brain mapping. To analyze fMRI data, GLM methodology has been developed to explain the timeline blood oxygenation level dependent signal. Currently, they are frequently used in fNIRS studies. The GLM defines measured data in the form of a linear combination of several variables and an error term. The observation of hemodynamic changes can be expressed as
where the y vector represents the measured data (in fNIRS, the vector is the observed time-series of the hemodynamic response), G is the design matrix obtained by convolving the canonical hemodynamic response with the experimental box-car function (Ye et al., 2009), β is the set of coefficients for the functional response that we want to estimate, and e is the error term. The vital part of the model function of a GLM is the box-car function, which reflects the temporal structures of the experimental paradigm and is convolved with the canonical hemodynamic response function (Ye et al., 2009). As physiological noises had already been removed using the Butterworth filter, only one explanatory variable (the design matrix) was used to extract the β values.
Least Squares Estimation
Least squares estimation is used to estimate the β values from the GLM. The time-course values predicted by the model are obtained by linear combination of the predictors
In order to achieve a good fit, the β values should be close to the predicted values that are as close as possible to the measured values y. Thus, the system of equations should be rearranged as
Although the GLM methodology does not estimate β values, it can be applied to minimize the sum of squared error values by using
where e′e shows the vector notation for the sum of squares. Utilizing LSE, the β weights minimizing the square error values are obtained by
The resulting matrix (G−1G)−1 plays an important role in the calculation of the β values. The remaining term on the right side, G′y, evaluates a vector containing as many elements as predictors. Figure 2 plots the ΔcHbO(t) signals for the MI task and rest period with their corresponding adaptively estimated β values for subject 5.
Feature Extraction and Classification
In this study, the statistical properties of the β values were used as the features. Signal peak (SP), signal skewness (SSk), signal mean (SM), signal variance (SV), signal kurtosis (SK), and signal slope (SS) were extracted from the β values obtained by LSE. The SSk values were determined by measuring the asymmetry of the signal values around the mean relative to a normal distribution:
where E is the expected value of Y and σ is the SD of Y
where N is the number of observations and Yi represents the β values. The variance is calculated as follow:
where μ is the mean value of Y. The kurtosis is computed as follows:
The SS is calculated using the polyfit function in MATLAB®. The SP values, which measure the peaks of signals, were determined using MATLAB® max function. These features were calculated across all 12 channels for the MI and MR. All of the feature values were scaled between 0 and 1 using the equation
where x ∈ Rn represents the original feature values, x′ denotes the rescaled feature values between 0 and 1, max(x) is the largest value, and min(x) is the smallest value. After extracting the features from the β values, SVM was used to classify the MI and MR tasks (Naseer et al., 2016b). SVM maximizes the margins between classes by creating hyperplanes that minimize the cost function
where wT, xi ∈ R2 and b ∈ R1, ||w||2 = wTw, C is the trade-off parameter between the error and the margin, ξi is the measure of the training data, and zi is the class label for the i-th sample. The most significant advantage of SVM is that it can be used as a linear as well as a non-linear classifier; in fact, in this study, a third-degree polynomial kernel function was used with C = 0.5. Ten-fold cross-validation was utilized to extract the classification accuracies for the MI and MR tasks versus rest periods. Moreover, in order to measure classification performance, recall and precision were calculated for both paradigms as follows:
where TP, FP, and FN denote true positive, false positive, and false negative, respectively. These values were calculated from confusion matrix (Fawcett, 2006).
Multiple feature combinations were used in order to extract significant classification accuracies for proposed and conventional methodologies. The classification accuracies obtained for the five subjects using the proposed method for MI versus rest were 79.5, 83.7, 82.6, 81.4, and 84.1% using SM and SSk, whereas those for MR versus rest were and 85.5, 85.2, 87.8, 83.7, and 84.8% using SP and SSk. To establish the superiority of the proposed method over the previous methods, the classification accuracies using the conventional hemodynamic response feature also were calculated. Figures 3A,B provides a schematic of the conventional and proposed methodology for fNIRS-based BCI study. Furthermore, the classification accuracies obtained for the five subjects using the conventional method for MI versus rest were 60.4, 78.9, 70.4, 68.9, and 54.4% using SM and SP, whereas those for MR versus rest were and 66.7, 73.0, 72.2, 68.5, and 63.3% using SM and SP. Tables 1 and 2 list the classification accuracies, precisions, and recalls of all subjects using the proposed methodology and the conventional method, for all possible two-feature combinations for MI versus rest task, respectively. Tables 3 and 4 list the classification accuracies, precisions, and recalls of all subjects using the proposed methodology and the conventional method, for all possible two-feature combinations for MR versus rest task, respectively. The results show that in MI task the optimal feature combinations that yielded best classification accuracies were “SM and SSk” and “SM and SP” for beta values and conventional hemodynamic response, respectively. In MR task, the optimal feature-combination that yielded best classification accuracies were “SP and SSk” and “SM and SP” for beta values and conventional hemodynamic response, respectively. In order to ensure that the data are normally distributed Kolmogorov–Smirnov method was applied, the significant value was found to be greater than 0.05 which shows normal distribution of the data. These high classification accuracies of the proposed method relative to the conventional method were statistically verified by a statistical significance test (the Student’s t-test): the p-values obtained by performing t-test on the subject-wise accuracy scores was less than 0.05, which confirmed the statistical significance of the proposed methodology’s superior performance for both tasks.
Table 1. Classification performances of proposed methodology for motor imagery versus rest task across all feature combinations.
Table 2. Classification performances of conventional methodology for motor imagery versus rest task across all feature combinations.
Table 3. Classification performances of proposed methodology for mental rotation versus rest task across all feature combinations.
Table 4. Classification performances with proposed and conventional methodologies for mental rotation versus rest task.
In previous studies, researchers have focused their efforts on enhancing the classification performance of multiple mental tasks in order to generate commands effective for control of external devices or for communication with patients suffering from amyotrophic lateral sclerosis, locked in syndrome, or other physical disabilities. However, distinct BCI signals for a specific mental task were unsuitable for classification, even when using current advanced methods. Previously, Tai and Chau (2009); Khan and Hong (2015); Naseer et al. (2016a); Naseer et al. (2016b) have used features extracted directly from hemodynamic response in order to acquire classification accuracies. In this study, a novel methodology that proceeds by adaptive estimation of GLM coefficients and extraction of the classification performances of MI versus rest and MR versus rest task were developed and evaluated. The results indicated enhanced classification performance as compared with a conventional hemodynamic-response-based fNIRS-BCI. Moreover, the proposed methodology can enhance classification performance if a user is not able to generate distinct brain signals for a specific mental task. The GLM methodology has been frequently employed to analyze time-series fMRI data: Abdelnour and Huppert (2009) first used the GLM in an fNIRS study in order to minimize physiological noises; soon thereafter, Hu et al. (2010) developed a novel online data analysis scheme using the GLM and Kalman estimator to reduce physiological noises for finger-tapping experiments; Aqil et al. (2012a) presented an online brain-imaging framework for finger-tapping tasks using GLM and a recursive least squares estimation method; Zhang et al. (2011b, 2012) tested multiple recursive algorithms for removal of physiological noises and, thereby, extraction of better neuron-related concentration changes in observed fNIRS data. All of these studies used the GLM for the removal of physiological noises and demonstrated brain-activation mapping for multiple cognitive tasks. However, the GLM coefficients, as estimated using LSE, have not been used as features for classification. The difference in classification accuracies is possible since in the conventional method we use statistical features obtained directly from HbO signals; whereas in the proposed method, we use statistical features obtained from β values. Multiple feature combinations have been used in order to determine optimal feature combination, which yields best classification accuracies, using proposed and conventional methodologies for both mental tasks. It was found that in MI task the optimal feature-combinations that yielded best classification accuracies were “SM and SSk” and “SM and SP” for beta values and conventional hemodynamic response, respectively. In MR task, the optimal feature-combination that yielded best classification accuracy were “SP and SSk” and “SM and SP” for beta values and conventional hemodynamic response, respectively. The proposed method has shown improved overall classification accuracies as compared to conventional methodology. This study showed that there is a significant difference between the classification accuracies of the proposed and conventional methodologies: the result is improved by an average of 18.4% for MI versus rest and 16.7% for MR versus rest using the proposed method. Moreover, it was found that features extracted from proposed methodology are statistically significant from conventional methodology for both paradigms.
It should be noted that in Naseer et al. (2016a), the best two- and three-feature combinations yielded accuracies of more than 90% for a seemingly very similar classification task. In this work, the accuracies obtained are in the range of 70%. These differences in the accuracies might be attributed to different recording conditions and different mental tasks. It is observed that the signal quality in mental arithmetic tasks is better as compared to MI and MR tasks. This could be attributed to user training as well. The subjects used in Naseer et al. (2016a) were regular fNIRS-based BCI users. All subjects in this paper had little or no experience of fNIRS recording/BCI training. The effect of using current method on the data from Naseer et al. (2016a) can be evaluated in future works.
This study has some limitations. The first is that only six features were used for classification. The combination of several other statistical features acquired from β values also should be utilized as features, as, thereby, classification performance could be further enhanced. The second limitation is that only SVM was used as the classifier. The positive effects of several other classifiers, however, have been seen. As shown in Naseer et al. (2016b), classification accuracies acquired using artificial neural networks (ANNs) are better than those acquired using SVM and, therefore, ANN’s can be considered for classification in future studies. The third limitation is that the proposed methodology is complex as compared to conventional method since an extra step of calculating general linear model coefficients is involved. This will increase the computational cost as well. The fourth limitation of this study is the fact that only two mental tasks for each paradigm were considered, which restricts this study to a two-class BCI problem. Certainly in any case, it can be upgraded to a multi-class BCI problem in a further study.
In conclusion, we present a novel methodology for enhanced classification accuracy of two-class fNIRS-based BCI. The hemodynamic signals of five subjects were modeled using the GLM and the beta values estimated by LSE were used to extract the features for classification. The classification accuracies obtained using the proposed methodology were significantly higher than those obtained using conventional hemodynamic-response-based features. These results indeed show enhanced classification performance relative to the conventional methodology and represent a step forward in the important task of making fNIRS-based BCIs more accurate and reliable.
This work was approved by the Institutional Review Board of Pusan National University. All experiments were conducted in accordance with the ethical standards encoded in the latest Declaration of Helsinki.
NQ conceived this study and was involved in the experiments, data processing, and writing of the manuscript. FN was involved in the experiments and data analysis. HN and RK were involved in data analysis, rechecking of results, and revision. SS helped in revision of the manuscript. NN was involved in the writing of the manuscript and supervised the research.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer, RC, and handling editor declared their shared affiliation, and the handling editor states that the process nevertheless met the standards of a fair and objective review.
We would like to thank Prof. Keum-Shik Hong from Pusan National University for providing us an opportunity to visit his lab and use the equipment in his lab for data acquisition.
Abdelnour, A. F., and Huppert, T. (2009). Real-time imaging of human brain function by near-infrared spectroscopy using an adaptive general linear model. Neuroimage 46, 133–143. doi:10.1016/j.neuroimage.2009.01.033
Abibullaev, B., and An, J. (2012). Classification of frontal cortex hemodynamic responses during cognitive tasks using wavelet transforms and machine learning algorithms. Med. Eng. Phys. 34, 1394–1410. doi:10.1016/j.medengphy.2012.01.002
Abibullaev, B., An, J., and Moon, J. I. (2011). Neural network classification of brain hemodynamic responses from four mental tasks. Int. J. Optomechatronics 5, 340–359. doi:10.1080/15599612.2011.633209
Aqil, M., Hong, K.-S., Jeong, M.-Y., and Ge, S. S. (2012b). Detection of event-related hemodynamic response to neuroactivation by dynamic modeling of brain activity. Neuroimage 63, 553–568. doi:10.1016/j.neuroimage.2012.07.006
Ayaz, H., Onaral, B., Izzetoglu, K., Shewokis, P. A., McKendrick, R., and Parasuraman, R. (2013). Continuous monitoring of brain dynamics with functional near infrared spectroscopy as a tool for neuroergonomic research: empirical examples and a technological development. Front. Hum. Neurosci. 7:871. doi:10.3389/fnhum.2013.00871
Ayaz, H., Shewokis, P. A., Scull, L., Libon, D. J., Feldman, S., Eppig, J., et al. (2014). Assessment of prefrontal cortex activity in amyotrophic lateral sclerosis patients with functional near infrared spectroscopy. J. Neurosci. Neuroeng. 3, 41–51. doi:10.1166/jnsne.2014.1095
Chen, X., Wang, Y., Nakanishi, M., Gao, X., Jung, T. P., and Gao, S. (2015). High-speed spelling with a noninvasive brain-computer interface. Proc. Natl. Acad. Sci. U.S.A. 112, E6058–E6067. doi:10.1073/pnas.1508080112
Cong, F., Kalyakin, I., Chang, Z., and Ristaniemi, T. (2011). Analysis on subtracting projection of extracted independent components from EEG recordings. Biomed. Tech. (Berl) 56, 223–234. doi:10.1515/BMT.2011.102
Cong, F., Lin, Q. H., Kuang, L. D., Gong, X. F., Astikainen, P., and Ristaniemi, T. (2015). Tensor decomposition of EEG signals: a brief review. J. Neurosci. Methods 248, 59–69. doi:10.1016/j.jneumeth.2015.03.018
Coyle, S. M., Ward, T. E., and Markham, C. M. (2007). Brain-computer interface using a simplified functional near-infrared spectroscopy system. J. Neural Eng. 4, 219–226. doi:10.1088/1741-2560/4/3/007
Coyle, S. M., Ward, T. E., Markham, C. M., and McDarby, G. (2004). On the suitability of near-infrared (NIR) systems for next generation brain-computer interfaces. Physiol. Meas. 25, 815–822. doi:10.1088/0967-3334/25/4/003
Delpy, D. T., Cope, M., van der Zee, P., Arridge, S., Wray, S., and Wyat, J. (1988). Estimation of optical path length through tissue from direct time of flight measurement. Phys. Med. Biol. 33, 1433–1442. doi:10.1088/0031-9155/33/12/008
Enzinger, C., Ropele, S., Fazekas, F., Loitfelder, M., Gorani, F., Seifert, T., et al. (2008). Brain motor system function in a patient with complete spinal cord injury following extensive brain-computer interface training. Exp. Brain Res. 190, 215–223. doi:10.1007/s00221-008-1465-y
Ferrari, M., Giannini, I., Sideri, G., and Zanette, E. (1985). Continuous non-invasive monitoring of human brain by near infrared spectroscopy. Adv. Exp. Med. Biol. 191, 873–882. doi:10.1007/978-1-4684-3291-688
Ferrari, M., and Quaresima, V. (2012). A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application. Neuroimage 63, 921–935. doi:10.1016/j.neuroimage.2012.03.049
Gratton, G., Brumback, C. R., Gordon, B. A., Pearson, M. A., Low, K. A., and Fabiani, M. (2006). Effects of measurement method, wavelength, and source-detector distance on the fast optical signal. Neuroimage 32, 1576–1590. doi:10.1016/j.neuroimage.2006.05.030
Hong, K.-S., and Naseer, N. (2016). Reduction of delay in detecting initial dips from functional near-infrared spectroscopy signals using vector-based phase analysis. Int. J. Neural. Syst. 26, 1650012. doi:10.1142/S012906571650012X
Hong, K.-S., and Santosa, H. (2016). Decoding four different sound-categories in the auditory cortex using functional near-infrared spectroscopy. Hear. Res. 333, 157–166. doi:10.1016/j.heares.2016.01.009
Hoshi, Y., Onoe, H., Watanabe, Y., Andersson, J., Bergstrom, M., Lilja, A., et al. (1994). Non-synchronous behavior of neuronal-activity, oxidative-metabolism and blood-supply during mental tasks in man. Neurosci. Lett. 172, 129–133. doi:10.1016/0304-3940(94)90679-3
Hu, X.-S., Hong, K.-S., Ge, S. S., and Jeong, M.-Y. (2010). Kalman estimator and general linear model-based on-line brain activation mapping by near-infrared spectroscopy. Biomed. Eng. Online 9, 82. doi:10.1186/1475-925X-9-82
Irani, F., Platek, S. M., Bunce, S., Ruocco, A. C., and Chute, D. (2007). Functional near-infrared spectroscopy (fNIRS): an emerging neuroimaging technology with important applications for the study of brain disorders. Clin. Neuropsychol. 21, 9–37. doi:10.1080/138540406009100
Jin, J., Allison, B. Z., Zhang, Y., Wang, X., and Cichocki, A. (2014). An ERP-based BCI using an oddball paradigm with different faces and reduced errors in critical functions. Int. J. Neural. Syst. 24, 1450027. doi:10.1142/S0129065714500270
Jin, J., Sellers, E. W., Zhou, S., Zhang, Y., Wang, X., and Cichocki, A. (2015). A P300 brain-computer interface based on a modification of the mismatch negativity paradigm. Int. J. Neural. Syst. 25, 1550011. doi:10.1142/S0129065715500112
Kamran, M. A., and Hong, K.-S. (2013). Linear parameter-varying model and adaptive filtering technique for detecting neuronal activities: an fNIRS study. J. Neural Eng. 10, 056002. doi:10.1088/1741-2560/10/5/056002
Kato, T., Kamei, A., Takashima, S., and Ozaki, T. (1993). Human visual cortical function during photic stimulation monitoring by means of near-infrared spectroscopy. J. Cereb. Blood Flow Metab. 13, 516–520. doi:10.1038/jcbfm.1993.66
Khoa, T. Q. D., and Nakagawa, M. (2008). Functional near-infrared spectroscope for cognition brain tasks by wavelets analysis and neural networks. Int. J. Biol. Life Sci. 4, 28–33. doi:10.1186/1753-4631-2-3
Naito, M., Michioka, Y., Ozawa, K., Ito, Y., Kiguchi, M., and Kanazawa, T. (2007). A communication means for totally locked-in ALS patients based on changes in cerebral blood volume measured with near-infrared light. IEICE Trans. Inf. Syst. E90D, 1028–1037. doi:10.1093/ietisy/e90-d.7.1028
Naseer, N., and Hong, K.-S. (2013). Classification of functional near-infrared spectroscopy signals corresponding to the right- and left-wrist motor imagery for development of a brain-computer interface. Neurosci. Lett. 553, 84–89. doi:10.1016/j.neulet.2013.08.021
Naseer, N., Hong, M. J., and Hong, K.-S. (2014). Online binary decision decoding using functional near-infrared spectroscopy for the development of brain-computer interface. Exp. Brain Res. 232, 555–564. doi:10.1007/s00221-013-3764-1
Naseer, N., Noori, F. M., Qureshi, N. K., and Hong, K.-S. (2016a). Determining optimal feature-combination for LDA classification of functional near-infrared spectroscopy signals in brain-computer interface application. Front. Hum. Neurosci. 10:237. doi:10.3389/fnhum.2016.00237
Naseer, N., Qureshi, N. K., Noori, F. M., and Hong, K.-S. (2016b). Analysis of different classification techniques for two-class functional near-infrared spectroscopy based brain-computer interface. Comput. Intell. Neurosci. 2016, 11. doi:10.1155/2016/5480760
Noori, F. M., Naseer, N., Qureshi, N. K., Nazeer, H., and Khan, R. A. (2017). Optimal feature selection from fNIRS signals using genetic algorithms for BCI. Neurosci. Lett. 647, 61–66. doi:10.1016/j.neulet.2017.03.013
Noori, F. M., Qureshi, N. K., Khan, R. A., and Naseer, N. (2016). “Feature selection based on modified genetic algorithm for optimization of functional near-infrared spectroscopy (fNIRS) signals for BCI,” in Robotics and Artificial Intelligence (ICRAI), 2nd Int. Conference (Rawalpindi: IEEE), 50–53.
Pfurtscheller, G., Müller, G. R., Pfurtscheller, J., Gerner, H. J., and Rupp, R. (2003). ’Thought’ – control of functional electrical stimulation to restore hand grasp in a patient with tetraplegia. Neurosci. Lett. 351, 33–36. doi:10.1016/S0304-3940(03)00947-9
Qureshi, N. K., Noori, F. M., Abdullah, A., and Naseer, N. (2016). “Comparison of classification performance for fNIRS-BCI system,” in Robotics and Artificial Intelligence (ICRAI), 2nd Int. Conference (Rawalpindi: IEEE), 54–57.
Salvaris, M., and Sepulveda, F. (2010). Classification effects of real and imaginary movement selective attention tasks on a P300-based brain-computer interface. J. Neural Eng. 7, 056004. doi:10.1088/1741-2560/7/5/056004
Sorger, B., Dahmen, B., Reithler, J., Gosseries, O., Maudoux, A., Laureys, S., et al. (2009). Another kind of ’BOLD response’: answering multiple-choice questions via online decoded single-trial brain signals. Prog. Brain Res. 177, 275–292. doi:10.1016/S0079-6123(09)17719-1
Tai, K., and Chau, T. (2009). Single-trial classification of NIRS signals during emotional induction tasks: towards a corporeal machine interface. J. Neuroeng. Rehabil. 6, 39. doi:10.1186/1743-0003-8-34
Villringer, A., Planck, J., Hock, C., Schleinkofer, L., and Dirnagl, U. (1993). Near infrared spectroscopy (NIRS): a new tool to study hemodynamic changes during activation of brain function in human adults. Neurosci. Lett. 154, 101–104. doi:10.1016/0304-3940(93)90181-J
Wester, B. A., Lee, R. H., and LaPlaca, M. C. (2009). Development and characterization of in vivo flexible electrodes compatible with large tissue displacements. J. Neural Eng. 6, 024002. doi:10.1088/1741-2560/6/2/024002
Wolpaw, J. R., Birbaumer, N., McFarland, D. J., Pfurtscheller, G., and Vaughan, T. M. (2002). Brain-computer interfaces for communication and control. Clin. Neurophysiol. 113, 767–791. doi:10.1016/S1388-2457(02)00057-3
Zhang, H., Zhang, Y. J., Duan, L., Ma, S. Y., Lu, C. M., and Zhu, C. Z. (2011a). Is resting-state functional connectivity revealed by functional near-infrared spectroscopy test-retest reliable? J. Biomed. Opt. 16, 067008. doi:10.1117/1.3591020
Zhang, Y., Sun, J., and Rolfe, P. (2011b). Reduction of global interference in functional multi-distance near-infrared spectroscopy using empirical mode decomposition and recursive least squares: a Monte Carlo study. J. Eur. Opt. Soc. 16, 067008. doi:10.2971/jeos.2011.11033
Keywords: functional near-infrared spectroscopy, brain–computer interface, general linear model, least squares estimation, adaptive estimation, support vector machine
Citation: Qureshi NK, Naseer N, Noori FM, Nazeer H, Khan RA and Saleem S (2017) Enhancing Classification Performance of Functional Near-Infrared Spectroscopy- Brain–Computer Interface Using Adaptive Estimation of General Linear Model Coefficients. Front. Neurorobot. 11:33. doi: 10.3389/fnbot.2017.00033
Received: 10 January 2017; Accepted: 22 June 2017;
Published: 17 July 2017
Edited by:Ganesh R. Naik, University of Technology Sydney, Australia
Reviewed by:Rifai Chai, University of Technology Sydney, Australia
Jing Jin, East China University of Science and Technology, China
Felix Putze, University of Bremen, Germany
Copyright: © 2017 Qureshi, Naseer, Noori, Nazeer, Khan and Saleem. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Noman Naseer, email@example.com