# Effect of EOG Signal Filtering on the Removal of Ocular Artifacts and EEG-Based Brain-Computer Interface: A Comprehensive Study

Authors: Malik M. Naeem Mannan; M. Ahmad Kamran; Shinil Kang; Myung Yung Jeong
Journal: Complexity
 (2018)

---

## Abstract

It is a fact that contamination of EEG by ocular artifacts reduces the classification accuracy of a brain-computer interface (BCI) and diagnosis of brain diseases in clinical research. Therefore, for BCI and clinical applications, it is very important to remove/reduce these artifacts before EEG signal analysis. Although, EOG-based methods are simple and fast for removing artifacts but their performance, meanwhile, is highly affected by the bidirectional contamination process. Some studies emphasized that the solution to this problem is low-pass filtering EOG signals before using them in artifact removal algorithm but there is still no evidence on the optimal low-pass frequency limits of EOG signals. In this study, we investigated the optimal EOG signal filtering limits using state-of-the-art artifact removal techniques with fifteen artificially contaminated EEG and EOG datasets. In this comprehensive analysis, unfiltered and twelve different low-pass filtering of EOG signals were used with five different algorithms, namely, simple regression, least mean squares, recursive least squares, REGICA, and AIR. Results from statistical testing of time and frequency domain metrics suggested that a low-pass frequency between 6 and 8 Hz could be used as the most optimal filtering frequency of EOG signals, both to maximally overcome/minimize the effect of bidirectional contamination and to achieve good results from artifact removal algorithms. Furthermore, we also used BCI competition IV datasets to show the efficacy of the proposed framework on real EEG signals. The motor-imagery-based BCI achieved statistically significant high-classification accuracies when artifacts from EEG were removed by using 7 Hz low-pass filtering as compared to all other filterings of EOG signals. These results also validated our hypothesis that low-pass filtering should be applied to EOG signals for enhancing the performance of each algorithm before using them for artifact removal process. Moreover, the comparison results indicated that the hybrid algorithms outperformed the performance of single algorithms for both simulated and experimental EEG datasets.


---

## Body

## 1. Introduction

The functional dynamics of the brain have been thoroughly investigated over the course of many years using noninvasive brain imaging techniques [1–5]. Electroencephalography (EEG), for example, is a portable neuroimaging system that can be used to assess different functional brain states [6–9]. However, a recorded EEG signal is highly contaminated with nonneuronal activities from different sources including eye blinking, eye movements, muscle movements, and electrocardiography (ECG) [10–15]. Eye movements and blinking generate high-magnitude artifacts as compared with the pure neuronal activity present in EEG data [16–18]. Such interferences are commonly known as ocular artifacts [19, 20].It is widely accepted within the BCI research community that in any BCI system, neurological phenomena are the only source of control [21, 22]. Artifacts, unwanted electrical signals that arise from sources other than the brain, can interfere with neurological phenomena. Such artifacts might alter the characteristics of neurological phenomena or even be mistakenly used as the source(s) of control in BCI systems [23]. Among the different artifacts, eye movement and blinks are the most important and major sources of physiological artifacts in BCI systems [24–26]. If not removed, these artifacts could, as indicated above, be mistakenly used to control the BCI system, which is the most significant artifact-related problem [27]. As failing to deal with artifacts can result in deterioration of BCI system performance during practical applications, it is necessary to develop automatic methods to handle artifacts or to design BCI systems robust to them. Bashashati et al. showed that dealing with eye artifacts in EEG data can enhance the performance of a self-paced BCI system [24]. Erfanian and Mahmoudi used recurrent neural networks to automatically suppress ocular artifacts for improved EEG-based BCI performance [25]. Recently, Yong et al. combined stationary wavelet analysis with adaptive thresholding to automatically remove ocular artifacts from EEG data in an EEG- and eye tracker-based self-paced BCI system [26]. They showed that their system can achieve higher BCI performance than can BCIs in which artifacts are not removed. Furthermore, artifacts can also affect diagnosis and analysis in clinical research such as on sleep disorders, Alzheimer disease, and schizophrenia [28–32]. It is therefore mandatory, in either clinical or practical research, to deal with these artifacts prior to the analysis of EEG signals.Several manual and automated methods have been developed to deal with this challenging task. One straightforward approach to the reduction of ocular artifacts is to prevent eye movements as much as possible, though requiring this and achieving it are two very different things. Also, the specific request of avoiding blinks could affect the investigated states and cognitive process of the subject [33]. Another commonly employed solution is to discard those epochs of EEG data that contain ocular artifacts, though this can incur the loss of neuronal activity-related EEG data. Alternatively, several automated methods for detection and removal/reduction of ocular artifacts have been proposed such as blind source separation-based methods, wavelet transforms, regression-based analysis, and empirical mode decomposition. Among these, the most commonly employed, which are known as regression-based algorithms, are based on the removal of electrooculography (EOG) contamination from EEG data. The simplest and most common procedure for removal of ocular artifacts from EEG data entails the subtraction, from each EEG channel, of reference channel signals containing proper artifactual interference. Techniques of this kind were widely applied until the mid-1990s, due to their low computational costs and simplicity [20, 33]. Subsequently, researchers used EOG channels to record eye movement and blinking to efficiently remove ocular artifacts from EEG data. EOG-based methodologies assume that the true neuronal activity and ocular artifacts are present in a linear combination in acquired EEG signals. These methods employ regression-based analysis [10, 18, 34–46], by which, in the time domain, the contamination coefficients of EOG signals are estimated and subtracted from each EEG channel to obtain clean EEG signals. Since being very low computation demanding, they are a great tool for real-time/online BCI applications. Although these methods have been proved to be more efficient than simple reference channel regression, their performance is highly affected by many factors. For example, neuronal activity from the frontal brain area, which EOG additionally measures, might be eliminated during the subtraction process, resulting in loss of true EEG signals [47]. Furthermore, these techniques are based on the assumption that EOG signals and neuronal activity recorded in EEG signals have no correlation, which has been found to be totally invalid [48, 49]. In order to overcome these issues, recently, regression-based algorithms have been combined with blind source separation techniques in the development of automated methodologies for removal of ocular artifacts [47, 50]. These methodologies have been shown to be more effective than either only regression- or blind source separation-based techniques, but they still lack the best results. The reason might be due to the fact that the outcomes are still affected by bidirectional contamination.The main problem with regression-based techniques is that they are always affected by bidirectional contamination; for example, EEG recordings are contaminated as the results of eye movement and blinking, while EOG recordings are contaminated by neuronal activities (originating mostly from the frontal and lateral frontal areas) [18, 47, 50, 51]. Therefore, removal of ocular artifacts using EOG signals would also remove common neuronal activity present in both EEG and EOG data. In a modified version of the regression method, called filtered regression, the effects of bidirectional contamination are reduced by low-pass filtering of EOG signals prior to regression analysis [39, 43, 44, 52]. This idea is based on studies that have shown that high-frequency components in EOG channels are generated from brain activity [43, 52]. In Table 1 list of the different low-pass frequencies used by researchers, it can be seen that there is no consensus on any particular frequency of EOG signal, though there is agreement on the fact that the most of the low-frequency components in EOG signals belongs to the ocular artifacts [43, 44, 47, 51, 52]. Determining the optimal low-pass frequency for EOG signals is very important, as the outcomes for regression-based correction methods can be affected by the selected filtering frequency. To the best of our knowledge, no study has investigated the optimal low-pass EOG signal filtering limits for use with regression-based algorithms. Our hypothesis is that EOG signal filtering will enhance the process of artifact removal to reduce bidirectional contamination and if so, then what are the optimal low-pass frequency limits which will give better results from all other low-pass filtering.Table 1
Different low-pass EOG signal filtering frequencies used in literature.
Study
Low-pass filtering (Hz)
Journal
Maddirala and Shaik [67]
~5
IEEE Sensors Journal
Yang et al. [68]
~100
Neurocomputing
Kanoga et al. [69]
~60
Neurocomputing
Mannan et al. [51]
~5
Front. Hum. Neuro.
Wang et al. [70]
~8
Biomed. Signal Process. Control
Zeng et al. [71]
~15
The Sci. World J.
Sameni and Gouy-Pailler [72]
~100
J. Neurosci. Methods
Murthy and Khan [73]
~30
Research J. Biotech.
Klados et al. [47]
~5
Biomed. Signal Process. Control
Pham et al. [42]
~200
Int. J. Psychophysiology
Ghandeharion and Erfanian [61]
~45
Medical Eng. Phy.
Chan et al. [74]
~30
Annals of Biomed. Eng.
Romero et al. [43]
~7.5
Annals of Biomed. Eng.
Romero et al. [44]
~7.5
Computer in Bio. Med.
Schlogl et al. [55]
~100
Clin. Neurophy.
Puthusserypady and Ratnarajah [75]
~11.5
Signal Process.
Erfanian and Mahmoudi [25]
~40
Med. Biol. Eng. Comput.
Croft et al. [38]
~100
Phychophy.
He et al. [45]
~30
Med. Bio. Eng. Comput.
Wallstrom et al. [18]
~100
Int. J. Psych.
Moretti et al. [76]
~100
Int. J. Psych.
Croft and Barry [66]
~100
Int. J. Psych.
Croft and Barry [35]
~35
Electroen. Clin. Neuro.
Sadasivan and Dutt [65]
~10
Signal Process.In this study, we used simulated contaminated EEG and EOG datasets and motor-imagery-based experimental BCI datasets to investigate the effect of different EOG filtering on the removal of ocular artifacts from EEG data. We used 12 different low-pass EOG signal filtering and unfiltered EOG data along with five different methods from the literature, namely, simple regression [44], least mean square-based regression [53], recursive least squares based regression [45], REGICA [47], and the method developed in [51] (we reference it hereafter as automatic independent component analysis and regression (AIR)) to determine the effect of different EOG signal filtering on artifact removal from EEG data. Since the underlying artifact-free EEG (true EEG) in artificially contaminated EEG data is known; therefore, it is possible to evaluate the effect of each EOG filtering using different performance metrics. The performance evaluation indexes employed were the mean square error and the mean absolute error in the time and frequency domains, respectively. Additionally, mutual information was utilized to estimate the common information between the reconstructed EEG signal and the artifact-free EEG signal. The improvements in reconstructed EEG is also evaluated using signal-to-artifact ratio before and after the artifact removal process for all algorithms. For real EEG datasets, we evaluated the classification accuracies of each subject and each method after the artifact removal using different low-pass EOG filtering. Finally, paired t-test was employed to the results of simulated and experimental datasets to find out the optimal EOG filtering with highest statistical significance. The results of this statistical testing revealed that best results from EOG-based algorithms could be achieved if low-pass frequency is used from 6–8 Hz. Furthermore, the results of both simulated and real EEG signals indicate that hybrid algorithms performed better than simple regression and adaptive filtering. A schematic diagram of the study is shown in Figure 1.Figure 1
Schematic plan of the study.


## 2. Materials and Methods

### 2.1. Materials

This section describes the detailed procedure used to simulate contaminated EEG and EOG datasets and real EEG datasets used from BCI competition IV.


#### 2.1.1. Simulated Datasets

(1) Participants. Fifteen (15) healthy subjects (all male) participated in this study. All had normal or corrected-to-normal vision. The experimental protocol was approved by the Institutional Review Board of Pusan National University. Experiment was conducted in accordance with the ethical guidelines established by the Institutional Review Board of Pusan National University and the Declaration of Helsinki. Each participant was asked to sign an informed consent form after being thoroughly informed about the nature and purpose of the study. The experiments were performed in a quiet room with dim lighting to prevent environmental disturbances. Each participant was seated in an armchair at a distance of about 1 m from a 24″ LCD monitor (ASUS; resolution: 1366 × 768).(2) Experiment 1. In this experiment, the participants were asked to sit relaxed and calm while keeping their eyes closed for 30 s. Each also was instructed to avoid moving their eyes during the experiment so as to avoid or at least minimize artifacts. After the experiment, the subjects’ data were carefully inspected for any presence of major artifacts; none was found. These datasets were then used as “clean EEG” signals for the purposes of further analysis. The data from this experiment will be referred as the “neuronal group” throughout in this paper.(3) Experiment 2. The experimental protocol was as follows. At the start of the experiment, the subject was instructed to sit relaxed and calm for 3 s. Three different word cues (blink, move horizontally, and move vertically) were used. The subjects were asked to blink their eyes or move them vertically or horizontally according to any of six visual cues (2 cues for each word) that appeared for 2 s each at the center of the screen. The interval between the cues was 2 s. At the end of the experiment, each subject was again asked to relax for 3 s. The total duration of the experiment was 30 s. The data from this experiment will be referred as “artifactual group” throughout in this paper.(4) EEG Recordings. The EEG data were acquired using an ActiCap 32-channel active electrode system with a BrainAmp DC amplifier (Brain Products GmbH, Gilching, Germany). The data sampling rate was 250 Hz. Nineteen (19) electrodes positioned according to the international 10–20 system (Fp1, Fp2, F3, F4, F7, F8, Fz, C3, C4, Cz, T7, T8, P7, P8, Pz, P3, P4, O1, and O2) were used for acquisition of EEG signals. AFz and FCz were used as the ground and reference electrodes, respectively. The impedance of all of the electrodes was reduced to below 5kΩ. The data were high-pass filtered at 0.5 Hz.(5) EOG Recordings. The EOG data were acquired using the BrainAmp ExG system (Brain Products GmbH, Gilching, Germany). Four electrodes were placed around the left and right eye to record ocular activities. All of the data were sampled at a rate of 250 Hz. Table 2 lists the class-wise low-pass frequencies used in this study.Table 2
Class-wise list of frequencies for low-pass EOG signal filtering used in this study.
Class
Low-pass frequency (Hz)
I
4
II
5
6
7
8
III
9
10
12
15
IV
20
25
35
Unfiltered(6) Simulated Datasets. In the present study, in order to investigate the optimal frequency limits, we simulated 15 artificially contaminated EEG and EOG datasets. As the underlying true EEG signal in artificially contaminated EEG data is known, such data can be used as a primary tool to determine the optimal filtering for EOG signals. We utilized EEG data recorded in an eyes-closed session to simulate contaminated EEG. Whereas such data might contain low-frequency eye movement contamination, they are notwithstanding preferred, as they tend to contain minimal overall artifacts. The alternative to this is eyes-open data acquisition. Note though that human eyes produce much-higher-amplitude signals in light than in darkness [54]. In this sense, recording of EEG signals in an eyes-closed session is preferred. However, EOG signals were acquired in an eyes-open session with different eye movements. Simple linear models were estimated to calculate the parameters of contamination for both EEG and EOG. By doing so, we can simulate signals which are bidirectionally contaminated, that is, EEG data is contaminated with EOG signals and EOG data is contaminated with EEG signals. In this sense, we can obtain simulated signals as close as possible to real signals [43].(6.1) Simulated EEG Signals. It is known that the recorded EEG signal contains pure activity from neurons, ocular artifacts, and measurement noise (artifacts from all other sources), as shown in Figure 2. In this paper, neuronal sources (EEG signals from the neuronal group) were artificially contaminated with ocular sources (EOG signals from the artifactual group) to simulate contaminated EEG signals [43]. These interferences were calculated by estimation of simple linear models between EEG and EOG recordings [44]. The detailed procedure is explained below.Figure 2
General approach for generation of EEG data.Four 2 s epochs with high EOG actives were selected from each one of the fifteen subjects from the artifactual group. A simple linear model was estimated for each 2 s epoch and for each channel. These models have two inputs corresponding toVEOGA and HEOGA signals (where XA refers to signals from the artifactual group) and one output that was each one of the 19 EEGA channels. For a better estimation of models, neuronal activity in EOG due to bidirectional contamination is reduced in order to calculate accurate parameters of contamination. For this purpose, in each model, EEGA and EOGA channels, corresponding to outputs and inputs, respectively, were low-pass filtered with the cutoff frequency corresponding to the highest value of the 99% (f99) of the total energy of these VEOGA and HEOGA signals [43]. Thus, the remainder 1% of signal energy was not considered as ocular activity (neural activity, power line interference, electrode noise, etc.) in the EOGA recordings [43]. The 99% cutoff frequencies obtained were 6.44 ± 2.43 and 7.61 ± 3.82 Hz for VEOGA and HEOGA, respectively, as mean and standard deviation for all epochs and subjects. This idea was based on and supported by the consideration that most components in the EOG signals related to the high frequencies are of neuronal source [52].A linear model used to estimate parameters of EOG interferences for each EEG channeli (i=1,2,3,…,19) and each epoch p (p=1,2,3,4) is evaluated as follows.
(1)EEGAip=αipVEOGAp+βipHEOGAp+e,where α and β are the unknown model parameters, and e is unknown error mapping. By this procedure, we got four α and four β parameters corresponding to all four epochs for each channel. These parameters were averaged (∑p=14αip/4 and ∑p=14βip/4) to obtain ocular contamination coefficients αi and βi for each channel. Finally, simulated EEG signals were generated according to Elbert’s contamination model [54] as
(2)cEEGNi=pEEGNi+αiVEOGA+βiHEOGA+v,where cEEG and pEEG are the artificially contaminated and pure EEG signals, subindex N refers to the neuronal group, and v is modeled as white Gaussian noise compensating for other noise sources. An example of simulated contaminated EEG data for electrode Fp1 is shown in Figure 3(a).Figure 3
Example of simulated signals. (a) Simulated contaminated EEG signal. (b) Simulated contaminated EOG signal.
(a)
(b)(6.2) Simulated EOG Signals. Ocular sources (EOG signals from the artifactual group) were artificially contaminated with neuronal sources (EEG signals from the neuronal group) to simulate contaminated EEG signals. Similar to simulated EEG signals, these interferences were calculated by approximation of simple linear models between EOG and EEG recordings. The detailed procedure is described below.Four 2 s epochs with no apparent EOG actives were selected from all of the fifteen subjects from the neuronal group. Neuronal contamination of EOG channels was obtained from the frontal electrodes (Fp1, Fp2, F7, and F8), which are the nearest ones to the eyes [44]. A linear model for both VHOGN and HEOGN was estimated for each 2 s epoch. These models have four inputs corresponding to EEGN signals from Fp1, Fp2, F7, and F8 and one output that was each one of the VHOGN and HEOGN.Two linear models used to estimate parameters of EEG interferences forVEOG and HEOG, and each epoch p p=1,2,3,4 is evaluated as follows:
(3)VEOGNip=a1pEEGNpFp1+a2pEEGNpFp2+a3pEEGNpF7+a4pEEGNpF8+e,(4)HEOGNip=b1pEEGNpFp1+b2pEEGNpFp2+b3pEEGNpF7+b4pEEGNpF8+e,where aj and bj j=1,2,3,4 are the unknown model parameters.By this procedure, we got fouraj and four bj parameters corresponding to all four epochs for VEOGN and HEOGN, respectively. These parameters were averaged (∑p=14ajp/4 and ∑p=14βjp/4) to obtain neuronal contamination coefficients aj and bj for each VEOGN and HEOGN, respectively. Finally, simulated EOG signals were generated as follows:
(5)cVEOGA=VEOGA+a1EEGNFp1+a2EEGNFp2+a3EEGNF7+a4EEGNF8,cHEOGA=HEOGA+b1EEGNFp1+b2EEGNFp2+b3EEGNF7+b4EEGNF8,where cVEOG and cHEOG are the artificially contaminated VEOG and HEOG, respectively. An example of simulated contaminated EOG data can be visualized in Figure 3(b).


#### 2.1.2. Experimental Datasets

In this study, the datasets from nine healthy subjects were sampled from publically available MI-based BCI signals of BCI competition IV (datasets 2a). Four different MI classes including left hand, right hand, both feet, and tongue were performed by all subjects. The experiment consisted of training sessions and evaluation sessions. For all subjects, each session consisted of six runs with short breaks. For each subject, two classes from evaluation sessions corresponding to left- and right-hand MI were selected. One run consisted of twelve trails of each class, resulting in 144 trails for each subject. Twenty-two EEG and three EOG channels [55] were used to record the data with a sampling frequency of 250 Hz and 50 Hz notch filter. All EEG channels were recorded monopolarly with the left mastoid serving as reference and the right mastoid as ground. The signals were band-pass filtered between 0.5 and 100 Hz. More details can be found in [56]. The datasets were highly contaminated with ocular artifacts which is a challenging problem in practical BCI systems [57].


### 2.2. Methods

In literature, different artifact removal algorithms have been developed to deal with the ocular contamination present in EEG signals. Broadly speaking, these algorithms can be divided into three main categories, namely, EOG-based, non-EOG-based, and hybrid algorithms. The most commonly employed among those are EOG-based regression algorithms [18, 33]. Although these algorithms are simple and perform well as compared with manual rejection, they nonetheless cause EEG data distortion due to bidirectional contamination [47]. To improve the performance of simple regression-based algorithms, researchers developed adaptive filter-based regression algorithms. Most commonly used adaptive filters are least mean squares and recursive least square-based filters, and these methods proved to be more effective as compared to simple regression. On the other hand, non-EOG-based algorithms, for example, the ICA-based algorithms, do not require any EOG signals, though their removal of artifactual ICs might cause the loss of substantial neuronal data, which is their major drawback [58, 59]. Urigüen and Garcia-Zapirain [10] suggested that multiple-combination artifact removal methods can be developed to efficiently remove artifacts from recorded EEG signals. Recently, EOG-based algorithms have been combined with non-EOG-based algorithms (ICA) to more effectively deal with ocular artifacts present in EEG data [47, 50, 51]. These algorithms were shown to outperform all of the algorithms with which they were compared in terms of artifact removal and maintenance of neuronal activity present in EEG data; however, their performance can be further improved if optimal low-pass EOG filtering is used. In this study, we used all three kinds of EOG-based algorithms (simple regression, adaptive regression, and hybrid methods) to investigate the effect of different low-pass filtering on the removal of ocular artifacts from EEG data. One may argue that there are other methods like ICA that can be used to remove ocular artifacts from EEG data without the need of EOG signals, but irrespective of their other disadvantages these methods cannot be used for real-time/online BCI applications whereas regression-based methods (simple and adaptive) are simple and fast; therefore, it can be used as an optimal option for BCI applications if their performance is enhanced [60].Next, we will briefly describe the implementation steps of the methods used in this study.


#### 2.2.1. Simple Regression Method

The simple regression method is implemented as follows [44]:
(1)
Equation (3) in [44] was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


#### 2.2.2. LMS Regression Method

The LMS regression method is implemented as follows [53]:
(1)
Least mean square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


#### 2.2.3. RLS Regression Method

RLS regression method is applied as follows [45]:
(1)
Recursive least square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


#### 2.2.4. REGICA

REGICA is implemented as follows [47]:
(1)
EEG signals are independent component analysis (ICA) decomposed.
(2)
Independent components (ICs) are filtered using recursive least square estimation with reference EOG signal.
(3)
ICs are backprojected to reconstruct EEG signal.


#### 2.2.5. AIR

AIR is implemented as follows [51]:
(1)
Contaminated EEG data are decomposed using ICA to obtain ICs.
(2)
Composite multiscale entropy and kurtosis are calculated to identify ocular-artifact-related ICs.
(3)
ICs are filtered using the linear regression model and extended recursive least mean squares.
(4)
Median absolute deviation is applied to remove any high-magnitude ocular artifacts left.
(5)
Artifact-free EEG data are obtained by backprojecting all ICs using inverse ICA.


### 2.3. Evaluation Indexes

#### 2.3.1. Mean Square Error

In this study, the performance of each algorithm for each of the EOG signal filtering ranges was evaluated using the mean square error. It was defined as [41].
(6)MSE=∑n=1NEEGout−EEGin2N,where EEGout is the reconstructed EEG and EEGin is the artifact-free EEG (EEG from the neuronal group).


#### 2.3.2. Mutual Information

The amount of mutual information between the reconstructed EEG signal and the artifact-free EEG was calculated in order to analyze the utility of each method for recovery of neuronal activity-related EEG signals. Mathematically, the calculation proceeds are as follows [61]:
(7)MI=∫−∞∞∫−∞∞fa,blogfa,bfafbdadb,where fa,b represents the joint pdf and fa and fb represent the marginal pdfs. The artifact-free EEG and reconstructed EEG were deemed to be closely related if and only if the mutual information values between them were large.


#### 2.3.3. Signal-to-Artifact Ratio

Signal-to-artifact ratio is the metric commonly used to evaluate the improvements in the corrected EEG signal as compared to the contaminated EEG signal. We calculated signal-to-artifact ratio for contaminated EEG signals as follows [62].
(8)SARB=1/N∑n=1NEEGin21/N∑n=1NEEGcn−EEGin2,where SARB is signal-to-artifact ratio before artifact removal, and EEGcn is the contaminated EEG signal. We also calculated signal-to-artifact ratio for corrected EEG as follows.
(9)SARA=1/N∑n=1NEEGin21/N∑n=1NEEGcr−EEGin2,where SARA is signal-to-artifact ratio after artifact removal and EEGcr is the corrected EEG signal. An effective artifact removal algorithm will remove all the artifacts and will have higher SARA values and consequently SARA>SARB. The gain in signal-to-artifact ratio γ can be calculated as follows.
(10)γ=10logSARASARB.Theγ value is positive if signal-to-artifact ratio is improved, negative if signal-to-artifact ratio is decreased and zero if there is no improvement.


#### 2.3.4. Mean Absolute Error

In order to measure the percentage distortion across the different frequency bands delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), and gamma (30–100 Hz), mean absolute error was defined as [51]
(11)MAE=PinEEG−PoutEEG,where P denotes the power spectrum density (PSD). PSD was estimated using the Welch method according to the following parameters: 200 sample points as the window length and 5 sample points as the overlap. The average PSD for each frequency band was calculated for all subjects.Also, mean absolute percentage error to estimate the percentage distortion in each frequency band was defined as [43].
(12)MAPE=100×PinEEG−PoutEEGPinEEG.


## 2.1. Materials

This section describes the detailed procedure used to simulate contaminated EEG and EOG datasets and real EEG datasets used from BCI competition IV.


### 2.1.1. Simulated Datasets

(1) Participants. Fifteen (15) healthy subjects (all male) participated in this study. All had normal or corrected-to-normal vision. The experimental protocol was approved by the Institutional Review Board of Pusan National University. Experiment was conducted in accordance with the ethical guidelines established by the Institutional Review Board of Pusan National University and the Declaration of Helsinki. Each participant was asked to sign an informed consent form after being thoroughly informed about the nature and purpose of the study. The experiments were performed in a quiet room with dim lighting to prevent environmental disturbances. Each participant was seated in an armchair at a distance of about 1 m from a 24″ LCD monitor (ASUS; resolution: 1366 × 768).(2) Experiment 1. In this experiment, the participants were asked to sit relaxed and calm while keeping their eyes closed for 30 s. Each also was instructed to avoid moving their eyes during the experiment so as to avoid or at least minimize artifacts. After the experiment, the subjects’ data were carefully inspected for any presence of major artifacts; none was found. These datasets were then used as “clean EEG” signals for the purposes of further analysis. The data from this experiment will be referred as the “neuronal group” throughout in this paper.(3) Experiment 2. The experimental protocol was as follows. At the start of the experiment, the subject was instructed to sit relaxed and calm for 3 s. Three different word cues (blink, move horizontally, and move vertically) were used. The subjects were asked to blink their eyes or move them vertically or horizontally according to any of six visual cues (2 cues for each word) that appeared for 2 s each at the center of the screen. The interval between the cues was 2 s. At the end of the experiment, each subject was again asked to relax for 3 s. The total duration of the experiment was 30 s. The data from this experiment will be referred as “artifactual group” throughout in this paper.(4) EEG Recordings. The EEG data were acquired using an ActiCap 32-channel active electrode system with a BrainAmp DC amplifier (Brain Products GmbH, Gilching, Germany). The data sampling rate was 250 Hz. Nineteen (19) electrodes positioned according to the international 10–20 system (Fp1, Fp2, F3, F4, F7, F8, Fz, C3, C4, Cz, T7, T8, P7, P8, Pz, P3, P4, O1, and O2) were used for acquisition of EEG signals. AFz and FCz were used as the ground and reference electrodes, respectively. The impedance of all of the electrodes was reduced to below 5kΩ. The data were high-pass filtered at 0.5 Hz.(5) EOG Recordings. The EOG data were acquired using the BrainAmp ExG system (Brain Products GmbH, Gilching, Germany). Four electrodes were placed around the left and right eye to record ocular activities. All of the data were sampled at a rate of 250 Hz. Table 2 lists the class-wise low-pass frequencies used in this study.Table 2
Class-wise list of frequencies for low-pass EOG signal filtering used in this study.
Class
Low-pass frequency (Hz)
I
4
II
5
6
7
8
III
9
10
12
15
IV
20
25
35
Unfiltered(6) Simulated Datasets. In the present study, in order to investigate the optimal frequency limits, we simulated 15 artificially contaminated EEG and EOG datasets. As the underlying true EEG signal in artificially contaminated EEG data is known, such data can be used as a primary tool to determine the optimal filtering for EOG signals. We utilized EEG data recorded in an eyes-closed session to simulate contaminated EEG. Whereas such data might contain low-frequency eye movement contamination, they are notwithstanding preferred, as they tend to contain minimal overall artifacts. The alternative to this is eyes-open data acquisition. Note though that human eyes produce much-higher-amplitude signals in light than in darkness [54]. In this sense, recording of EEG signals in an eyes-closed session is preferred. However, EOG signals were acquired in an eyes-open session with different eye movements. Simple linear models were estimated to calculate the parameters of contamination for both EEG and EOG. By doing so, we can simulate signals which are bidirectionally contaminated, that is, EEG data is contaminated with EOG signals and EOG data is contaminated with EEG signals. In this sense, we can obtain simulated signals as close as possible to real signals [43].(6.1) Simulated EEG Signals. It is known that the recorded EEG signal contains pure activity from neurons, ocular artifacts, and measurement noise (artifacts from all other sources), as shown in Figure 2. In this paper, neuronal sources (EEG signals from the neuronal group) were artificially contaminated with ocular sources (EOG signals from the artifactual group) to simulate contaminated EEG signals [43]. These interferences were calculated by estimation of simple linear models between EEG and EOG recordings [44]. The detailed procedure is explained below.Figure 2
General approach for generation of EEG data.Four 2 s epochs with high EOG actives were selected from each one of the fifteen subjects from the artifactual group. A simple linear model was estimated for each 2 s epoch and for each channel. These models have two inputs corresponding toVEOGA and HEOGA signals (where XA refers to signals from the artifactual group) and one output that was each one of the 19 EEGA channels. For a better estimation of models, neuronal activity in EOG due to bidirectional contamination is reduced in order to calculate accurate parameters of contamination. For this purpose, in each model, EEGA and EOGA channels, corresponding to outputs and inputs, respectively, were low-pass filtered with the cutoff frequency corresponding to the highest value of the 99% (f99) of the total energy of these VEOGA and HEOGA signals [43]. Thus, the remainder 1% of signal energy was not considered as ocular activity (neural activity, power line interference, electrode noise, etc.) in the EOGA recordings [43]. The 99% cutoff frequencies obtained were 6.44 ± 2.43 and 7.61 ± 3.82 Hz for VEOGA and HEOGA, respectively, as mean and standard deviation for all epochs and subjects. This idea was based on and supported by the consideration that most components in the EOG signals related to the high frequencies are of neuronal source [52].A linear model used to estimate parameters of EOG interferences for each EEG channeli (i=1,2,3,…,19) and each epoch p (p=1,2,3,4) is evaluated as follows.
(1)EEGAip=αipVEOGAp+βipHEOGAp+e,where α and β are the unknown model parameters, and e is unknown error mapping. By this procedure, we got four α and four β parameters corresponding to all four epochs for each channel. These parameters were averaged (∑p=14αip/4 and ∑p=14βip/4) to obtain ocular contamination coefficients αi and βi for each channel. Finally, simulated EEG signals were generated according to Elbert’s contamination model [54] as
(2)cEEGNi=pEEGNi+αiVEOGA+βiHEOGA+v,where cEEG and pEEG are the artificially contaminated and pure EEG signals, subindex N refers to the neuronal group, and v is modeled as white Gaussian noise compensating for other noise sources. An example of simulated contaminated EEG data for electrode Fp1 is shown in Figure 3(a).Figure 3
Example of simulated signals. (a) Simulated contaminated EEG signal. (b) Simulated contaminated EOG signal.
(a)
(b)(6.2) Simulated EOG Signals. Ocular sources (EOG signals from the artifactual group) were artificially contaminated with neuronal sources (EEG signals from the neuronal group) to simulate contaminated EEG signals. Similar to simulated EEG signals, these interferences were calculated by approximation of simple linear models between EOG and EEG recordings. The detailed procedure is described below.Four 2 s epochs with no apparent EOG actives were selected from all of the fifteen subjects from the neuronal group. Neuronal contamination of EOG channels was obtained from the frontal electrodes (Fp1, Fp2, F7, and F8), which are the nearest ones to the eyes [44]. A linear model for both VHOGN and HEOGN was estimated for each 2 s epoch. These models have four inputs corresponding to EEGN signals from Fp1, Fp2, F7, and F8 and one output that was each one of the VHOGN and HEOGN.Two linear models used to estimate parameters of EEG interferences forVEOG and HEOG, and each epoch p p=1,2,3,4 is evaluated as follows:
(3)VEOGNip=a1pEEGNpFp1+a2pEEGNpFp2+a3pEEGNpF7+a4pEEGNpF8+e,(4)HEOGNip=b1pEEGNpFp1+b2pEEGNpFp2+b3pEEGNpF7+b4pEEGNpF8+e,where aj and bj j=1,2,3,4 are the unknown model parameters.By this procedure, we got fouraj and four bj parameters corresponding to all four epochs for VEOGN and HEOGN, respectively. These parameters were averaged (∑p=14ajp/4 and ∑p=14βjp/4) to obtain neuronal contamination coefficients aj and bj for each VEOGN and HEOGN, respectively. Finally, simulated EOG signals were generated as follows:
(5)cVEOGA=VEOGA+a1EEGNFp1+a2EEGNFp2+a3EEGNF7+a4EEGNF8,cHEOGA=HEOGA+b1EEGNFp1+b2EEGNFp2+b3EEGNF7+b4EEGNF8,where cVEOG and cHEOG are the artificially contaminated VEOG and HEOG, respectively. An example of simulated contaminated EOG data can be visualized in Figure 3(b).


### 2.1.2. Experimental Datasets

In this study, the datasets from nine healthy subjects were sampled from publically available MI-based BCI signals of BCI competition IV (datasets 2a). Four different MI classes including left hand, right hand, both feet, and tongue were performed by all subjects. The experiment consisted of training sessions and evaluation sessions. For all subjects, each session consisted of six runs with short breaks. For each subject, two classes from evaluation sessions corresponding to left- and right-hand MI were selected. One run consisted of twelve trails of each class, resulting in 144 trails for each subject. Twenty-two EEG and three EOG channels [55] were used to record the data with a sampling frequency of 250 Hz and 50 Hz notch filter. All EEG channels were recorded monopolarly with the left mastoid serving as reference and the right mastoid as ground. The signals were band-pass filtered between 0.5 and 100 Hz. More details can be found in [56]. The datasets were highly contaminated with ocular artifacts which is a challenging problem in practical BCI systems [57].


## 2.1.1. Simulated Datasets

(1) Participants. Fifteen (15) healthy subjects (all male) participated in this study. All had normal or corrected-to-normal vision. The experimental protocol was approved by the Institutional Review Board of Pusan National University. Experiment was conducted in accordance with the ethical guidelines established by the Institutional Review Board of Pusan National University and the Declaration of Helsinki. Each participant was asked to sign an informed consent form after being thoroughly informed about the nature and purpose of the study. The experiments were performed in a quiet room with dim lighting to prevent environmental disturbances. Each participant was seated in an armchair at a distance of about 1 m from a 24″ LCD monitor (ASUS; resolution: 1366 × 768).(2) Experiment 1. In this experiment, the participants were asked to sit relaxed and calm while keeping their eyes closed for 30 s. Each also was instructed to avoid moving their eyes during the experiment so as to avoid or at least minimize artifacts. After the experiment, the subjects’ data were carefully inspected for any presence of major artifacts; none was found. These datasets were then used as “clean EEG” signals for the purposes of further analysis. The data from this experiment will be referred as the “neuronal group” throughout in this paper.(3) Experiment 2. The experimental protocol was as follows. At the start of the experiment, the subject was instructed to sit relaxed and calm for 3 s. Three different word cues (blink, move horizontally, and move vertically) were used. The subjects were asked to blink their eyes or move them vertically or horizontally according to any of six visual cues (2 cues for each word) that appeared for 2 s each at the center of the screen. The interval between the cues was 2 s. At the end of the experiment, each subject was again asked to relax for 3 s. The total duration of the experiment was 30 s. The data from this experiment will be referred as “artifactual group” throughout in this paper.(4) EEG Recordings. The EEG data were acquired using an ActiCap 32-channel active electrode system with a BrainAmp DC amplifier (Brain Products GmbH, Gilching, Germany). The data sampling rate was 250 Hz. Nineteen (19) electrodes positioned according to the international 10–20 system (Fp1, Fp2, F3, F4, F7, F8, Fz, C3, C4, Cz, T7, T8, P7, P8, Pz, P3, P4, O1, and O2) were used for acquisition of EEG signals. AFz and FCz were used as the ground and reference electrodes, respectively. The impedance of all of the electrodes was reduced to below 5kΩ. The data were high-pass filtered at 0.5 Hz.(5) EOG Recordings. The EOG data were acquired using the BrainAmp ExG system (Brain Products GmbH, Gilching, Germany). Four electrodes were placed around the left and right eye to record ocular activities. All of the data were sampled at a rate of 250 Hz. Table 2 lists the class-wise low-pass frequencies used in this study.Table 2
Class-wise list of frequencies for low-pass EOG signal filtering used in this study.
Class
Low-pass frequency (Hz)
I
4
II
5
6
7
8
III
9
10
12
15
IV
20
25
35
Unfiltered(6) Simulated Datasets. In the present study, in order to investigate the optimal frequency limits, we simulated 15 artificially contaminated EEG and EOG datasets. As the underlying true EEG signal in artificially contaminated EEG data is known, such data can be used as a primary tool to determine the optimal filtering for EOG signals. We utilized EEG data recorded in an eyes-closed session to simulate contaminated EEG. Whereas such data might contain low-frequency eye movement contamination, they are notwithstanding preferred, as they tend to contain minimal overall artifacts. The alternative to this is eyes-open data acquisition. Note though that human eyes produce much-higher-amplitude signals in light than in darkness [54]. In this sense, recording of EEG signals in an eyes-closed session is preferred. However, EOG signals were acquired in an eyes-open session with different eye movements. Simple linear models were estimated to calculate the parameters of contamination for both EEG and EOG. By doing so, we can simulate signals which are bidirectionally contaminated, that is, EEG data is contaminated with EOG signals and EOG data is contaminated with EEG signals. In this sense, we can obtain simulated signals as close as possible to real signals [43].(6.1) Simulated EEG Signals. It is known that the recorded EEG signal contains pure activity from neurons, ocular artifacts, and measurement noise (artifacts from all other sources), as shown in Figure 2. In this paper, neuronal sources (EEG signals from the neuronal group) were artificially contaminated with ocular sources (EOG signals from the artifactual group) to simulate contaminated EEG signals [43]. These interferences were calculated by estimation of simple linear models between EEG and EOG recordings [44]. The detailed procedure is explained below.Figure 2
General approach for generation of EEG data.Four 2 s epochs with high EOG actives were selected from each one of the fifteen subjects from the artifactual group. A simple linear model was estimated for each 2 s epoch and for each channel. These models have two inputs corresponding toVEOGA and HEOGA signals (where XA refers to signals from the artifactual group) and one output that was each one of the 19 EEGA channels. For a better estimation of models, neuronal activity in EOG due to bidirectional contamination is reduced in order to calculate accurate parameters of contamination. For this purpose, in each model, EEGA and EOGA channels, corresponding to outputs and inputs, respectively, were low-pass filtered with the cutoff frequency corresponding to the highest value of the 99% (f99) of the total energy of these VEOGA and HEOGA signals [43]. Thus, the remainder 1% of signal energy was not considered as ocular activity (neural activity, power line interference, electrode noise, etc.) in the EOGA recordings [43]. The 99% cutoff frequencies obtained were 6.44 ± 2.43 and 7.61 ± 3.82 Hz for VEOGA and HEOGA, respectively, as mean and standard deviation for all epochs and subjects. This idea was based on and supported by the consideration that most components in the EOG signals related to the high frequencies are of neuronal source [52].A linear model used to estimate parameters of EOG interferences for each EEG channeli (i=1,2,3,…,19) and each epoch p (p=1,2,3,4) is evaluated as follows.
(1)EEGAip=αipVEOGAp+βipHEOGAp+e,where α and β are the unknown model parameters, and e is unknown error mapping. By this procedure, we got four α and four β parameters corresponding to all four epochs for each channel. These parameters were averaged (∑p=14αip/4 and ∑p=14βip/4) to obtain ocular contamination coefficients αi and βi for each channel. Finally, simulated EEG signals were generated according to Elbert’s contamination model [54] as
(2)cEEGNi=pEEGNi+αiVEOGA+βiHEOGA+v,where cEEG and pEEG are the artificially contaminated and pure EEG signals, subindex N refers to the neuronal group, and v is modeled as white Gaussian noise compensating for other noise sources. An example of simulated contaminated EEG data for electrode Fp1 is shown in Figure 3(a).Figure 3
Example of simulated signals. (a) Simulated contaminated EEG signal. (b) Simulated contaminated EOG signal.
(a)
(b)(6.2) Simulated EOG Signals. Ocular sources (EOG signals from the artifactual group) were artificially contaminated with neuronal sources (EEG signals from the neuronal group) to simulate contaminated EEG signals. Similar to simulated EEG signals, these interferences were calculated by approximation of simple linear models between EOG and EEG recordings. The detailed procedure is described below.Four 2 s epochs with no apparent EOG actives were selected from all of the fifteen subjects from the neuronal group. Neuronal contamination of EOG channels was obtained from the frontal electrodes (Fp1, Fp2, F7, and F8), which are the nearest ones to the eyes [44]. A linear model for both VHOGN and HEOGN was estimated for each 2 s epoch. These models have four inputs corresponding to EEGN signals from Fp1, Fp2, F7, and F8 and one output that was each one of the VHOGN and HEOGN.Two linear models used to estimate parameters of EEG interferences forVEOG and HEOG, and each epoch p p=1,2,3,4 is evaluated as follows:
(3)VEOGNip=a1pEEGNpFp1+a2pEEGNpFp2+a3pEEGNpF7+a4pEEGNpF8+e,(4)HEOGNip=b1pEEGNpFp1+b2pEEGNpFp2+b3pEEGNpF7+b4pEEGNpF8+e,where aj and bj j=1,2,3,4 are the unknown model parameters.By this procedure, we got fouraj and four bj parameters corresponding to all four epochs for VEOGN and HEOGN, respectively. These parameters were averaged (∑p=14ajp/4 and ∑p=14βjp/4) to obtain neuronal contamination coefficients aj and bj for each VEOGN and HEOGN, respectively. Finally, simulated EOG signals were generated as follows:
(5)cVEOGA=VEOGA+a1EEGNFp1+a2EEGNFp2+a3EEGNF7+a4EEGNF8,cHEOGA=HEOGA+b1EEGNFp1+b2EEGNFp2+b3EEGNF7+b4EEGNF8,where cVEOG and cHEOG are the artificially contaminated VEOG and HEOG, respectively. An example of simulated contaminated EOG data can be visualized in Figure 3(b).


## 2.1.2. Experimental Datasets

In this study, the datasets from nine healthy subjects were sampled from publically available MI-based BCI signals of BCI competition IV (datasets 2a). Four different MI classes including left hand, right hand, both feet, and tongue were performed by all subjects. The experiment consisted of training sessions and evaluation sessions. For all subjects, each session consisted of six runs with short breaks. For each subject, two classes from evaluation sessions corresponding to left- and right-hand MI were selected. One run consisted of twelve trails of each class, resulting in 144 trails for each subject. Twenty-two EEG and three EOG channels [55] were used to record the data with a sampling frequency of 250 Hz and 50 Hz notch filter. All EEG channels were recorded monopolarly with the left mastoid serving as reference and the right mastoid as ground. The signals were band-pass filtered between 0.5 and 100 Hz. More details can be found in [56]. The datasets were highly contaminated with ocular artifacts which is a challenging problem in practical BCI systems [57].


## 2.2. Methods

In literature, different artifact removal algorithms have been developed to deal with the ocular contamination present in EEG signals. Broadly speaking, these algorithms can be divided into three main categories, namely, EOG-based, non-EOG-based, and hybrid algorithms. The most commonly employed among those are EOG-based regression algorithms [18, 33]. Although these algorithms are simple and perform well as compared with manual rejection, they nonetheless cause EEG data distortion due to bidirectional contamination [47]. To improve the performance of simple regression-based algorithms, researchers developed adaptive filter-based regression algorithms. Most commonly used adaptive filters are least mean squares and recursive least square-based filters, and these methods proved to be more effective as compared to simple regression. On the other hand, non-EOG-based algorithms, for example, the ICA-based algorithms, do not require any EOG signals, though their removal of artifactual ICs might cause the loss of substantial neuronal data, which is their major drawback [58, 59]. Urigüen and Garcia-Zapirain [10] suggested that multiple-combination artifact removal methods can be developed to efficiently remove artifacts from recorded EEG signals. Recently, EOG-based algorithms have been combined with non-EOG-based algorithms (ICA) to more effectively deal with ocular artifacts present in EEG data [47, 50, 51]. These algorithms were shown to outperform all of the algorithms with which they were compared in terms of artifact removal and maintenance of neuronal activity present in EEG data; however, their performance can be further improved if optimal low-pass EOG filtering is used. In this study, we used all three kinds of EOG-based algorithms (simple regression, adaptive regression, and hybrid methods) to investigate the effect of different low-pass filtering on the removal of ocular artifacts from EEG data. One may argue that there are other methods like ICA that can be used to remove ocular artifacts from EEG data without the need of EOG signals, but irrespective of their other disadvantages these methods cannot be used for real-time/online BCI applications whereas regression-based methods (simple and adaptive) are simple and fast; therefore, it can be used as an optimal option for BCI applications if their performance is enhanced [60].Next, we will briefly describe the implementation steps of the methods used in this study.


### 2.2.1. Simple Regression Method

The simple regression method is implemented as follows [44]:
(1)
Equation (3) in [44] was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


### 2.2.2. LMS Regression Method

The LMS regression method is implemented as follows [53]:
(1)
Least mean square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


### 2.2.3. RLS Regression Method

RLS regression method is applied as follows [45]:
(1)
Recursive least square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


### 2.2.4. REGICA

REGICA is implemented as follows [47]:
(1)
EEG signals are independent component analysis (ICA) decomposed.
(2)
Independent components (ICs) are filtered using recursive least square estimation with reference EOG signal.
(3)
ICs are backprojected to reconstruct EEG signal.


### 2.2.5. AIR

AIR is implemented as follows [51]:
(1)
Contaminated EEG data are decomposed using ICA to obtain ICs.
(2)
Composite multiscale entropy and kurtosis are calculated to identify ocular-artifact-related ICs.
(3)
ICs are filtered using the linear regression model and extended recursive least mean squares.
(4)
Median absolute deviation is applied to remove any high-magnitude ocular artifacts left.
(5)
Artifact-free EEG data are obtained by backprojecting all ICs using inverse ICA.


## 2.2.1. Simple Regression Method

The simple regression method is implemented as follows [44]:
(1)
Equation (3) in [44] was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


## 2.2.2. LMS Regression Method

The LMS regression method is implemented as follows [53]:
(1)
Least mean square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


## 2.2.3. RLS Regression Method

RLS regression method is applied as follows [45]:
(1)
Recursive least square estimation was used to estimate the parameters of EOG signals.
(2)
Artifact-free EEG was reconstructed by subtracting estimated VEOG and HEOG from contaminated EEG.


## 2.2.4. REGICA

REGICA is implemented as follows [47]:
(1)
EEG signals are independent component analysis (ICA) decomposed.
(2)
Independent components (ICs) are filtered using recursive least square estimation with reference EOG signal.
(3)
ICs are backprojected to reconstruct EEG signal.


## 2.2.5. AIR

AIR is implemented as follows [51]:
(1)
Contaminated EEG data are decomposed using ICA to obtain ICs.
(2)
Composite multiscale entropy and kurtosis are calculated to identify ocular-artifact-related ICs.
(3)
ICs are filtered using the linear regression model and extended recursive least mean squares.
(4)
Median absolute deviation is applied to remove any high-magnitude ocular artifacts left.
(5)
Artifact-free EEG data are obtained by backprojecting all ICs using inverse ICA.


## 2.3. Evaluation Indexes

### 2.3.1. Mean Square Error

In this study, the performance of each algorithm for each of the EOG signal filtering ranges was evaluated using the mean square error. It was defined as [41].
(6)MSE=∑n=1NEEGout−EEGin2N,where EEGout is the reconstructed EEG and EEGin is the artifact-free EEG (EEG from the neuronal group).


### 2.3.2. Mutual Information

The amount of mutual information between the reconstructed EEG signal and the artifact-free EEG was calculated in order to analyze the utility of each method for recovery of neuronal activity-related EEG signals. Mathematically, the calculation proceeds are as follows [61]:
(7)MI=∫−∞∞∫−∞∞fa,blogfa,bfafbdadb,where fa,b represents the joint pdf and fa and fb represent the marginal pdfs. The artifact-free EEG and reconstructed EEG were deemed to be closely related if and only if the mutual information values between them were large.


### 2.3.3. Signal-to-Artifact Ratio

Signal-to-artifact ratio is the metric commonly used to evaluate the improvements in the corrected EEG signal as compared to the contaminated EEG signal. We calculated signal-to-artifact ratio for contaminated EEG signals as follows [62].
(8)SARB=1/N∑n=1NEEGin21/N∑n=1NEEGcn−EEGin2,where SARB is signal-to-artifact ratio before artifact removal, and EEGcn is the contaminated EEG signal. We also calculated signal-to-artifact ratio for corrected EEG as follows.
(9)SARA=1/N∑n=1NEEGin21/N∑n=1NEEGcr−EEGin2,where SARA is signal-to-artifact ratio after artifact removal and EEGcr is the corrected EEG signal. An effective artifact removal algorithm will remove all the artifacts and will have higher SARA values and consequently SARA>SARB. The gain in signal-to-artifact ratio γ can be calculated as follows.
(10)γ=10logSARASARB.Theγ value is positive if signal-to-artifact ratio is improved, negative if signal-to-artifact ratio is decreased and zero if there is no improvement.


### 2.3.4. Mean Absolute Error

In order to measure the percentage distortion across the different frequency bands delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), and gamma (30–100 Hz), mean absolute error was defined as [51]
(11)MAE=PinEEG−PoutEEG,where P denotes the power spectrum density (PSD). PSD was estimated using the Welch method according to the following parameters: 200 sample points as the window length and 5 sample points as the overlap. The average PSD for each frequency band was calculated for all subjects.Also, mean absolute percentage error to estimate the percentage distortion in each frequency band was defined as [43].
(12)MAPE=100×PinEEG−PoutEEGPinEEG.


## 2.3.1. Mean Square Error

In this study, the performance of each algorithm for each of the EOG signal filtering ranges was evaluated using the mean square error. It was defined as [41].
(6)MSE=∑n=1NEEGout−EEGin2N,where EEGout is the reconstructed EEG and EEGin is the artifact-free EEG (EEG from the neuronal group).


## 2.3.2. Mutual Information

The amount of mutual information between the reconstructed EEG signal and the artifact-free EEG was calculated in order to analyze the utility of each method for recovery of neuronal activity-related EEG signals. Mathematically, the calculation proceeds are as follows [61]:
(7)MI=∫−∞∞∫−∞∞fa,blogfa,bfafbdadb,where fa,b represents the joint pdf and fa and fb represent the marginal pdfs. The artifact-free EEG and reconstructed EEG were deemed to be closely related if and only if the mutual information values between them were large.


## 2.3.3. Signal-to-Artifact Ratio

Signal-to-artifact ratio is the metric commonly used to evaluate the improvements in the corrected EEG signal as compared to the contaminated EEG signal. We calculated signal-to-artifact ratio for contaminated EEG signals as follows [62].
(8)SARB=1/N∑n=1NEEGin21/N∑n=1NEEGcn−EEGin2,where SARB is signal-to-artifact ratio before artifact removal, and EEGcn is the contaminated EEG signal. We also calculated signal-to-artifact ratio for corrected EEG as follows.
(9)SARA=1/N∑n=1NEEGin21/N∑n=1NEEGcr−EEGin2,where SARA is signal-to-artifact ratio after artifact removal and EEGcr is the corrected EEG signal. An effective artifact removal algorithm will remove all the artifacts and will have higher SARA values and consequently SARA>SARB. The gain in signal-to-artifact ratio γ can be calculated as follows.
(10)γ=10logSARASARB.Theγ value is positive if signal-to-artifact ratio is improved, negative if signal-to-artifact ratio is decreased and zero if there is no improvement.


## 2.3.4. Mean Absolute Error

In order to measure the percentage distortion across the different frequency bands delta (0.5–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta (12–30 Hz), and gamma (30–100 Hz), mean absolute error was defined as [51]
(11)MAE=PinEEG−PoutEEG,where P denotes the power spectrum density (PSD). PSD was estimated using the Welch method according to the following parameters: 200 sample points as the window length and 5 sample points as the overlap. The average PSD for each frequency band was calculated for all subjects.Also, mean absolute percentage error to estimate the percentage distortion in each frequency band was defined as [43].
(12)MAPE=100×PinEEG−PoutEEGPinEEG.


## 3. Results

This study investigated the effect of low-pass filtering on the removal of ocular artifacts from EEG data. According to the literature, different studies used different low-pass filtering for EOG signals ranging from 4 to 100 Hz. In this study, we used unfiltered and twelve different low-pass EOG filtering with simulated EEG datasets in an effort to find the optimal one from which best results could be achieved. For this purpose, we used five different methods from the literature with five performance evaluation metrics. Table3 lists the average mean square errors of all of the simulated datasets and for all of the electrodes with each method. It can be seen that the mean square error was lowest when the low-pass filter from 6 to 8 Hz was used. Furthermore, the mutual information scores also were calculated in a time domain analysis; the results of which are shown in Table 4. Similar to the mean square error results, the average mutual information for all of the datasets and all of the electrodes was the maximum when one of the 6 to 8 Hz low-pass filter was used with a different method. Table 5 shows the average improvement gained in reconstructing artifact-free EEG by calculating signal-to-artifact ratio before and after the artifact removal. It can be seen that corrected EEG with every EOG faltering showed improved signal-to-artifact ratio but results of this analysis also indicate similar outcomes as in the case of mean square error and mutual information. Moreover, the investigation of the optimal filtering range was carried out also in the frequency domain, by calculating the mean absolute error and mean absolute percentage values for the different frequency bands. The effect of bidirectional contamination can be best observed and analyzed in frequency domain by evaluating the distortion produce in different frequency bands using different filtering. Results of frequency domain analysis are shown in Tables 6–10. Except delta band, mean absolute error for all bands was lowest when 6 or 7 Hz low-pass EOG filtering was used. In case of delta band, least mean squares and recursive least squares showed lowest errors with 6 Hz of low-pass filter; but in case of REGICA and AIR unfiltered EOG, data showed lowest errors. Moreover, Figure 4 depicts and compare the pure EEG and output EEG with different EOG filtering to analyze the effect of bidirectional contamination for time and frequency domain. It can be seen in Figure 4(a) through highlighted regions that different levels of distortions were introduced to EEG signal with different low-pass filtering of EOG signals. Specifically, when using high value low-pass filter (e.g., 35 Hz or unfiltered, magenta, and green line), the distortion in the neuronal signal is much more as compared to other low-pass filtering’s (e.g., 7 Hz, blue line). It can be seen that the reconstructed EEG from 7 Hz low-pass filtering (blue line) follows the true EEG (black line) very closely as compared to all other outputs. Furthermore, Figure 4(b) illustrates the effect of bidirectional contamination in frequency domain. It can be seen through the highlighted box that artifact-free EEG with 7 Hz successfully recovered the frequencies similar to true EEG signal whereas all other produced distortion in the frequencies. This advocate our hypothesis that bidirectional contamination could be reduced by using optimal low-pass filtered EOG signals.Table 3
Mean square error values (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression(mean ± SD)
LMS(mean ± SD)
RLS(mean ± SD)
REGICA(mean ± SD)
AIR(mean ± SD)
4
16.42
26.35
10.90
13.03
6.73
16.49
5.45
12.90
4.64
6.55
5
14.80
23.74
9.42
11.43
5.22
12.39
3.88
8.19
3.07
4.21
6
14.19
22.74
8.94
11.09
4.68
11.22
3.27
6.59
2.46
3.37
7
14.02
22.45
8.97
11.21
4.59
11.23
3.09
6.22
2.28
3.15
8
14.06
22.46
9.24
11.48
4.72
11.67
3.11
6.22
2.30
3.16
9
14.18
22.61
9.59
11.75
4.93
12.19
3.20
6.31
2.39
3.27
10
14.31
22.79
9.89
11.97
5.13
12.62
3.30
6.43
2.49
3.41
12
14.51
23.06
10.25
12.23
5.39
13.15
3.47
6.67
2.66
3.63
15
14.66
23.28
10.49
12.40
5.58
13.53
3.61
6.93
2.80
3.81
20
14.78
23.45
10.64
12.52
5.71
13.85
3.71
7.17
2.90
3.95
25
14.83
23.52
10.71
12.57
5.77
14.00
3.76
7.29
2.95
4.01
35
14.88
23.59
10.77
12.60
5.83
14.09
3.80
7.37
3.00
4.07
Unfiltered
15.00
23.78
10.91
12.71
5.95
14.18
3.92
7.48
3.12
4.24Table 4
Mutual information scores (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression(mean ± SD)
LMS(mean ± SD)
RLS(mean ± SD)
REGICA(mean ± SD)
AIR(mean ± SD)
4
1.49
0.24
1.49
0.17
1.60
0.20
1.68
0.17
1.68
0.17
5
1.52
0.26
1.53
0.18
1.65
0.22
1.74
0.19
1.75
0.18
6
1.54
0.28
1.55
0.19
1.68
0.23
1.78
0.20
1.78
0.19
7
1.54
0.28
1.56
0.19
1.69
0.24
1.79
0.21
1.80
0.20
8
1.54
0.28
1.55
0.20
1.69
0.24
1.79
0.21
1.80
0.20
9
1.53
0.27
1.54
0.20
1.68
0.23
1.79
0.21
1.79
0.20
10
1.52
0.27
1.53
0.19
1.67
0.23
1.78
0.20
1.78
0.19
12
1.52
0.26
1.52
0.19
1.66
0.22
1.76
0.20
1.77
0.19
15
1.51
0.26
1.51
0.19
1.64
0.22
1.75
0.19
1.76
0.18
20
1.51
0.26
1.50
0.19
1.64
0.22
1.74
0.19
1.75
0.18
25
1.50
0.26
1.50
0.18
1.63
0.22
1.74
0.19
1.75
0.18
35
1.50
0.26
1.50
0.18
1.63
0.22
1.73
0.19
1.74
0.18
Unfiltered
1.49
0.25
1.49
0.18
1.62
0.21
1.73
0.18
1.73
0.18Table 5
Signal-to-artifact ratio for contaminated EEG, corrected EEG, and gain (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
SAR before artifact removal
Simple
Adaptive filtering
Hybrid methods
Regression(mean ± SD)
LMS(mean ± SD)
RLS(mean ± SD)
REGICA(mean ± SD)
AIR(mean ± SD)
After
Gain
After
Gain
After
Gain
After
Gain
After
Gain
4
78.54
112.7
1.44
116.3
1.49
120.7
1.54
123.2
1.57
124.0
1.58
5
114.5
1.46
118.3
1.51
123.9
1.58
127.3
1.63
128.3
1.64
6
115.4
1.47
119.2
1.52
125.7
1.60
129.8
1.66
130.8
1.67
7
115.7
1.47
119.2
1.52
126.2
1.61
130.7
1.67
131.8
1.68
8
115.6
1.47
118.9
1.52
125.9
1.61
130.7
1.67
131.7
1.68
9
115.4
1.47
118.5
1.51
125.4
1.60
130.2
1.66
131.2
1.68
10
115.2
1.47
118.2
1.51
124.9
1.60
129.7
1.66
130.7
1.67
12
114.9
1.46
117.7
1.50
124.2
1.59
128.9
1.65
129.9
1.66
15
114.7
1.46
117.5
1.50
123.7
1.58
128.4
1.64
129.3
1.65
20
114.5
1.46
117.3
1.50
123.4
1.58
128.0
1.63
128.9
1.65
25
114.5
1.46
117.2
1.50
123.3
1.57
127.8
1.63
128.7
1.64
35
114.4
1.46
117.1
1.50
123.1
1.57
127.7
1.63
128.6
1.64
Unfiltered
114.2
1.46
116.9
1.49
122.8
1.57
127.2
1.62
128.1
1.64Table 6
Mean absolute error (MAE) and mean absolute percentage error (MAPE) values in delta frequency band (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
4
0.613
9.249
0.399
7.721
0.319
6.774
0.332
6.750
0.329
6.734
5
0.544
5.858
0.331
4.517
0.250
3.400
0.263
3.575
0.260
3.547
6
0.539
5.968
0.327
4.013
0.245
2.816
0.258
3.053
0.255
3.027
7
0.539
6.061
0.329
4.068
0.245
2.787
0.258
3.014
0.256
2.977
8
0.539
6.080
0.335
4.150
0.247
2.828
0.258
3.045
0.256
3.010
9
0.538
6.082
0.340
4.206
0.248
2.847
0.257
3.055
0.256
3.021
10
0.538
6.080
0.344
4.248
0.249
2.857
0.257
3.058
0.256
3.024
12
0.538
6.077
0.349
4.290
0.251
2.865
0.257
3.060
0.255
3.027
15
0.538
6.076
0.351
4.312
0.251
2.870
0.257
3.060
0.255
3.028
20
0.538
6.077
0.352
4.326
0.251
2.872
0.257
3.060
0.255
3.028
25
0.538
6.078
0.352
4.333
0.251
2.874
0.257
3.060
0.255
3.028
35
0.538
6.078
0.353
4.340
0.251
2.875
0.257
3.061
0.255
3.029
Unfiltered
0.538
6.083
0.353
4.392
0.251
2.889
0.257
3.059
0.255
3.027Table 7
Mean absolute error (MAE) and mean absolute percentage error (MAPE) values in theta frequency band (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
4
0.281
10.048
0.280
10.386
0.278
10.401
0.279
10.449
0.278
10.446
5
0.145
5.536
0.141
6.120
0.138
6.030
0.137
6.042
0.137
6.040
6
0.093
3.590
0.088
3.605
0.085
3.409
0.081
3.359
0.081
3.354
7
0.094
2.020
0.092
1.879
0.087
1.589
0.082
1.484
0.083
1.468
8
0.113
3.292
0.111
2.590
0.106
2.531
0.099
2.435
0.100
2.448
9
0.130
4.134
0.128
3.395
0.122
3.317
0.114
3.205
0.116
3.231
10
0.140
4.511
0.138
3.764
0.132
3.677
0.124
3.565
0.125
3.593
12
0.149
4.753
0.147
4.009
0.140
3.910
0.131
3.794
0.133
3.824
15
0.151
4.816
0.150
4.084
0.142
3.979
0.133
3.853
0.135
3.884
20
0.152
4.827
0.150
4.108
0.142
3.999
0.134
3.863
0.135
3.894
25
0.152
4.829
0.150
4.113
0.143
4.003
0.134
3.864
0.135
3.895
35
0.152
4.829
0.150
4.116
0.143
4.006
0.134
3.864
0.135
3.895
Unfiltered
0.152
4.835
0.150
4.146
0.143
4.029
0.134
3.862
0.135
3.893Table 8
Mean absolute error (MAE) and mean absolute percentage error (MAPE) values in alpha frequency band (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
4
0.021
0.773
0.020
0.699
0.020
0.688
0.019
0.677
0.019
0.675
5
0.020
0.704
0.020
0.638
0.019
0.623
0.019
0.611
0.018
0.608
6
0.018
0.609
0.018
0.544
0.017
0.528
0.017
0.518
0.017
0.514
7
0.022
0.544
0.021
0.445
0.020
0.425
0.019
0.421
0.019
0.416
8
0.040
0.849
0.038
0.764
0.037
0.746
0.035
0.712
0.035
0.714
9
0.069
1.640
0.066
1.552
0.064
1.542
0.062
1.479
0.062
1.487
10
0.097
2.621
0.094
2.508
0.091
2.495
0.088
2.409
0.088
2.421
12
0.140
4.405
0.136
4.202
0.132
4.180
0.127
4.074
0.128
4.093
15
0.171
6.216
0.166
5.811
0.161
5.789
0.155
5.651
0.156
5.678
20
0.180
6.793
0.175
6.318
0.170
6.294
0.164
6.145
0.164
6.175
25
0.181
6.852
0.175
6.372
0.171
6.347
0.164
6.195
0.165
6.225
35
0.181
6.862
0.176
6.382
0.171
6.358
0.164
6.204
0.165
6.234
Unfiltered
0.181
6.865
0.176
6.387
0.171
6.360
0.164
6.204
0.165
6.235Table 9
Mean absolute error (MAE) and mean absolute percentage error (MAPE) values in beta frequency band (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
4
0.007
1.018
0.007
1.122
0.005
0.794
0.005
0.679
0.005
0.676
5
0.007
0.959
0.006
1.039
0.005
0.703
0.004
0.573
0.004
0.573
6
0.007
0.938
0.006
0.982
0.004
0.646
0.003
0.510
0.003
0.506
7
0.007
0.944
0.006
0.961
0.004
0.637
0.003
0.491
0.003
0.486
8
0.007
1.003
0.006
0.987
0.004
0.678
0.003
0.519
0.003
0.512
9
0.007
1.128
0.007
1.074
0.005
0.768
0.004
0.595
0.004
0.588
10
0.009
1.367
0.008
1.228
0.007
0.930
0.006
0.790
0.006
0.784
12
0.016
2.410
0.016
2.336
0.014
2.057
0.013
1.917
0.013
1.916
15
0.034
5.027
0.034
5.208
0.033
4.925
0.031
4.719
0.031
4.725
20
0.066
10.752
0.067
11.173
0.065
10.705
0.062
10.381
0.062
10.405
25
0.094
15.951
0.095
16.437
0.092
15.767
0.089
15.289
0.089
15.324
35
0.122
18.770
0.124
19.315
0.121
18.521
0.117
17.945
0.117
17.996
Unfiltered
0.127
19.029
0.129
19.581
0.125
18.771
0.121
18.184
0.122
18.237Table 10
Mean absolute error (MAE) and mean absolute percentage error (MAPE) values in gamma frequency band (mean of all EEG channels and mean of all simulated subjects).
Low-pass EOG filter (Hz)
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
MAE
MAPE
4
0.037
0.159
0.034
0.146
0.026
0.110
0.021
0.091
0.020
0.088
5
0.034
0.146
0.031
0.133
0.022
0.095
0.018
0.075
0.017
0.071
6
0.033
0.141
0.030
0.129
0.021
0.089
0.016
0.068
0.015
0.064
7
0.033
0.140
0.031
0.130
0.021
0.089
0.016
0.066
0.015
0.062
8
0.033
0.141
0.031
0.134
0.022
0.092
0.016
0.067
0.015
0.063
9
0.034
0.142
0.032
0.138
0.022
0.095
0.016
0.069
0.015
0.065
10
0.034
0.143
0.033
0.141
0.023
0.098
0.017
0.071
0.016
0.067
12
0.034
0.144
0.034
0.144
0.024
0.102
0.017
0.073
0.016
0.069
15
0.034
0.145
0.034
0.146
0.024
0.104
0.018
0.075
0.017
0.071
20
0.034
0.146
0.035
0.148
0.025
0.106
0.018
0.077
0.017
0.073
25
0.035
0.156
0.036
0.161
0.026
0.119
0.019
0.090
0.018
0.087
35
0.046
0.270
0.047
0.281
0.037
0.242
0.031
0.214
0.030
0.211
Unfiltered
0.300
3.581
0.308
3.593
0.296
3.502
0.287
3.375
0.287
3.395Figure 4
Analysis of bidirectional contamination by comparing reconstructed EEG obtained using different EOG low-pass filtering. (a) Time domain. (b) Frequency domain.
(a)
(b)We used pairedt-test to statistically compare results from all the metrics to find out if there are any differences in the outputs with different low-pass filtering. Before applying paired t-test, as listed in Table 2, all EOG filtering frequencies used in this study were split into four classes which are 4 Hz (belongs to delta band), 5–8 Hz (belongs to theta band), 9–15 Hz (belongs to alpha and low-beta band), and 20 Hz–unfiltered (belongs to high-beta and gamma band). We divided this statistical testing into two steps. In the first step, we analyzed results of paired t-test to select an optimal low-pass filtering class which showed minimum errors with significantly increased results from all four classes as listed in Figure 5. Mean square errors obtained in the range of 5–8 Hz were lowest with 6.81 ± 11.21 (averaged for all methods) when compared to 4 Hz (8.83 ± 15.06, p<0.027), 9–15 Hz (7.14 ± 11.60, nonsignificant difference), 20 Hz–unfiltered (7.65 ± 12.32, p<0.1 except for the LMS method). In case of frequency domain analysis in delta band, 5–8 Hz (MAE: 0.326 ± 0.14; MAPE: 3.89 ± 2.30%; p<0.026), 9–15 Hz (MAE: 0.329 ± 0.14; MAPE: 3.86 ± 2.16%; p<0.10), and 20 Hz–unfiltered (MAE: 0.330 ± 0.14; MAPE: 3.88 ± 2.12%; p<0.159) showed low errors when compared to 4 Hz (MAE: 0.398 ± 0.17; MAPE: 7.45 ± 8.73%). There was a nonsignificant statistical difference observed when all other ranges were compared (5–8 Hz versus 9–15 Hz, p<0.589; 5–8 Hz versus 20 Hz–unfiltered, p<0.518; and 9–15 Hz versus 20 Hz-unfiltered, p<0.965) for delta band. Five to eight Hz class (theta: MAE: 0.105 ± 0.03; MAPE: 3.44 ± 5.27%; alpha: MAE: 0.023+ 0.007; MAPE: 0.60+ 0.39%; beta: MAE: 0.004+ 0.001; MAPE: 0.73+ 0.57%) showed highly significantly increased results for theta, alpha, and beta frequency bands when compared with 4 Hz (theta: MAE: 0.279 ± 0.14; MAPE: 10.35 ± 13.66%; p<0.013; alpha: MAE: 0.020 ± 0.007; MAPE: 0.70 ± 0.60%; p<0.160; beta: MAE: 0.006 ± 0.001; MAPE: 0.86 ± 0.57%; p<0.053), 9–15 Hz (theta: MAE: 0.134 ± 0.04; MAPE: 3.86 ± 6.74%; p<0.031; alpha: MAE: 0.11 ± 0.03; MAPE: 3.51 ± 2.23%; p<0.001; beta: MAE: 0.015 ± 0.004; MAPE: 2.23 ± 1.84%; p<0.001), and 20 Hz–unfiltered (theta: MAE: 0.143 ± 0.04; MAPE: 4.14 ± 7.31%; p<0.016; alpha: MAE: 0.171 ± 0.05; MAPE: 6.39 ± 4.91%; p<0.001; beta: MAE: 0.100 ± 0.03; MAPE: 15.93 ± 11.40%; p<0.001). In case of gamma band, 4 Hz (MAE: 0.028 ± 0.008; MAPE: 0.19 ± 0.03%) showed significantly improved performance when compared with 20 Hz–unfiltered (MAE: 0.097 ± 0.03; MAPE: 0.99 ± 1.06%; p<0.001) and nonsignificant difference when compared with 9–15 Hz (MAE: 0.025 ± 0.01; MAPE: 0.11 ± 0.04%; p<0.885), but the errors were low for 5–8 Hz (MAE: 0.023 ± 0.009; MAPE: 0.10 ± 0.04%; p<0.079). Both 5–8 Hz and 9–15 Hz have nonsignificant differences for gamma band but showed statistically increased results when compared with 20 Hz–unfiltered (p<0.001). Overall, in summary, 5–8 Hz outperformed 4 Hz for all metrics except for alpha band (p<0.079) and when compared with 9–15 Hz and 20 Hz–unfiltered showed highly significantly increased results for theta, alpha, and beta bands (p<0.05) and better results with nonsignificant difference for all other metrics. From this analysis, it can be concluded that 5–8 Hz perform better as compared to all other low-pass filtering classes. Next, we analyzed results of each low-pass frequency in 5–8 Hz range to check whether there is any single optimal frequency which showed lowest errors with high-level of statistical significance (Figure 6). In this analysis, no significant differences were observed between most of the metrics when they compared statistically. However, outputs from 6 Hz (MSE: 6.70 ± 11.00; theta: MAE: 0.086 ± 0.02, MAPE: 3.46 ± 6.03%; beta: MAE: 0.004 ± 0.001, MAPE: 0.72 ± 0.56%; gamma: MAE: 0.023 ± 0.009, MAPE: 0.098 ± 0.03%) low-pass filter show significantly increased results when compared with 5 Hz (MSE: 7.28 ± 12.01; theta: MAE: 0.140 ± 0.06, MAPE: 5.95 ± 10.15%; beta: MAE: 0.005 ± 0.001, MAPE: 0.77 ± 0.56%; gamma: MAE: 0.024 ± 0.008, MAPE: 0.10 ± 0.03%) low-pass filter results with p<0.074 in mean square error, theta, beta, and gamma bands (except for LMS in gamma band with nonsignificant difference), while in delta and alpha bands low significant differences were observed (p<0.153). There was no significant difference observed between 6, 7, and 8 Hz except for theta and alpha band in 6 Hz versus 8 Hz (p<0.072) and theta, alpha, and beta bands in 7 Hz versus 8 Hz (p<0.096). It can be concluded that results from 6–8 Hz showed statistically increased results when compared with outputs from all other low-pass filtering and unfiltered EOG signals; but overall, there is no significant difference observed in 6–8 Hz except for few cases. Furthermore, these results advocate our hypothesis that unfiltered EOG signals cause high bidirectional contamination and low-pass filtering should be applied to EOG signals before using them in artifact removal algorithms.Figure 5
Statistical increases between different metrics for artifact corrections with different EOG filtering (pairedt-tests were used between class A versus class B). Bold and underline: p<0.01; Bold: p<0.05; Italic: p<0.1; Color indicate better class.Figure 6
Statistical increases between different metrics for artifact corrections with different EOG filtering of best class from Figure5 (paired t-tests were used between subclass A versus subclass B). Bold and underline: p<0.01; Bold: p<0.05; Italic: p<0.1; Color indicate better class.In this study, we also used real EEG datasets to verify the efficacy of the proposed framework. The five algorithms described above were applied separately to each subject’s data to remove ocular artifacts with four low-pass EOG filtering. These EOG filtering (4 Hz, 7 Hz, 12 Hz, and UF) were chosen such that there is one low-pass filtering from each class (Table2). We compared the classification accuracies obtained after the application of each method with each EOG low-pass filtering. In BCI studies, common spatial pattern (CSP) is the most commonly used filtering technique to extract features from EEG signals [63]. Generally, the goal of the CSP is to find spatial filters by maximizing the variance of one class while minimizing the variance of the other to discriminate the two populations of EEG signals [64]. Finally, we used linear discernment analysis (LDA) for classification of the extracted features of the two classes due to its simplicity and low computational cost. For each subject and each artifact-free EEG obtained after using different EOG filtering, 6 runs of 6-fold cross-validation were used to calculate the classification accuracies. The 6-fold cross-validation randomly divides the data into six equal partitions and use five set of partitions for training and 1 set of partition for testing. This process was repeated for six times, and the average accuracy for each subject was calculated. The average classification accuracies of each subjects for all six sessions with artifact-free EEG data using different low-pass EOG filtering, and each method are listed in Table 11. Similar to the results of simulated signals, it can be visualized from Table 11 that results from the candidate of class II low-pass filtering showed highest classification accuracies (REG: 69.52 ± 4.52; LMS: 73.45 ± 2.04; RLS: 73.99 ± 2.75; REGICA: 77.46 ± 3.59; AIR: 76.85 ± 3.27) for all subjects and all methods when compared with the results from the candidates of all other classes (class I: REG: 66.43 ± 4.11; LMS: 70.13 ± 2.45; RLS: 70.83 ± 2.59; REGICA: 73.37 ± 3.64; AIR: 73.53 ± 4.71; class III: REG: 67.20 ± 6.09; LMS: 70.37 ± 2.84; RLS: 72.14 ± 4.65; REGICA: 74.53 ± 3.00; AIR: 73.30 ± 3.65; class IV: REG: 67.12 ± 3.75; LMS: 69.90 ± 2.40; RLS: 69.67 ± 4.48; REGICA: 73.07 ± 2.21; AIR: 72.91 ± 3.06). We further validated these results statistically using paired t-test. This analysis revealed that classification accuracies obtained with 7 Hz showed highly significant results when compared with outputs from 4 Hz (p<0.028 for all methods), 12 Hz (p<0.01 for LMS, REGICA, and AIR and p<0.163 for REG and RLS), and unfiltered EOG (p<0.032 for all methods). Furthermore, it could also be noted that hybrid methods demonstrated high-classification accuracies as compared to simple regression and adaptive filtering methods. These results from experimental EEG datasets verified the results from simulated EEG signals that a low-pass EOG filtering from 6–8 Hz could be used to remove artifacts efficiently. As it can be seen from Table 11 that classification accuracies obtained from all low-pass filtered EOG are higher than those obtained when unfiltered EOG was used. From Figure 7, the t-testing also revealed that results from each low-pass filtering showed statistically increased classification accuracies as compared to the results from artifact-free EEG obtained after using unfiltered EOG signals. Therefore, these results also validate our hypothesis that low-pass filtering of EOG can be used to minimize the effect of bidirectional contamination problem.Table 11
Classification accuracies (mean of all six runs for each subject).
Low-pass EOG filter (Hz)
Sub
Simple
Adaptive filtering
Hybrid methods
Regression
LMS
RLS
REGICA
AIR
4
1
65.28
70.83
71.53
77.08
77.08
2
61.11
73.61
66.67
78.47
80.56
3
67.36
65.97
71.53
70.14
72.92
4
64.58
66.67
69.44
69.44
71.53
5
63.89
70.14
72.22
70.83
65.97
6
63.19
70.14
70.83
70.83
69.44
7
72.92
70.14
73.61
72.92
76.39
8
66.67
72.22
67.36
72.22
70.14
9
72.92
71.53
74.31
78.47
77.78
Mean ± SD
66.43 ± 4.11
70.13 ± 2.45
70.83 ± 2.59
73.37 ± 3.64
73.53 ± 4.71
7
1
69.44
72.92
79.17
81.25
81.25
2
63.19
75.69
72.22
80.56
79.17
3
73.61
72.92
72.92
75.69
75.00
4
69.44
70.14
71.53
73.61
74.31
5
64.58
72.22
75.69
73.61
71.53
6
65.28
71.53
74.31
74.31
75.00
7
73.61
76.39
76.39
76.39
79.86
8
70.14
74.31
70.14
78.47
75.69
9
76.39
75.00
73.61
83.33
79.86
Mean ± SD
69.52 ± 4.52
73.45 ± 2.04
73.99 ± 2.75
77.46 ± 3.59
76.85 ± 3.27
12
1
76.39
71.53
76.39
77.78
77.08
2
57.64
71.53
70.14
72.22
75.69
3
68.75
70.83
66.67
74.31
70.83
4
66.67
68.75
65.28
72.92
68.06
5
61.81
65.97
74.31
72.92
67.36
6
61.81
69.44
75.00
72.92
73.61
7
71.53
76.39
77.78
77.08
75.00
8
66.67
68.75
68.06
70.83
75.69
9
73.61
70.14
75.69
79.86
76.39
Mean ± SD
67.20 ± 6.09
70.37 ± 2.84
72.14 ± 4.65
74.53 ± 3.00
73.30 ± 3.65
UF
1
67.36
72.22
77.78
73.61
75.00
2
65.28
70.83
68.06
70.83
75.69
3
70.14
68.06
64.58
70.14
67.36
4
67.36
70.14
65.97
71.53
73.61
5
60.42
65.97
75.00
72.92
71.53
6
63.19
66.67
70.14
73.61
71.53
7
70.83
71.53
66.67
73.61
70.83
8
67.36
71.53
66.67
73.61
72.92
9
72.22
72.22
72.22
77.78
77.78
Mean ± SD
67.12 ± 3.75
69.90 ± 2.40
69.67 ± 4.48
73.07 ± 2.21
72.91 ± 3.06Figure 7
Statistical increases between classification accuracies after artifact corrections with different EOG filtering (pairedt-tests were used between group 1 versus group 2). Bold and underline: p<0.01; Bold: p<0.05; Italic: p<0.1; Color indicate better class.


## 4. Discussions

Many studies have shown that the performance of BCI applications can be reduced due to the presence of ocular artifacts in EEG data [23–27]. Among different methods, EOG-based algorithms are simple and fast due to which could be used as a good tool for real-time/online BCI applications if their performance is enhanced, since it is highly affected by bidirectional contamination [10, 44, 47, 51, 60]. The simplest solution to this problem is low-pass filtering EOG signals before using them in artifact removal algorithm [10, 39]. In efforts to overcome the effect of bidirectional contamination, a number of studies have applied different low-pass filtering on EOG signals ranging from 5 to 100 Hz, but there is no consensus on which low-pass frequency should be used for optimal results. The idea of low-pass filtering EOG is based on studies that have demonstrated that high-frequency components in EOG signals are generated from brain activities [52], and this is supported by some studies [10, 39, 43, 44, 47, 51]. It has been previously shown that performance of simple regression-based algorithms can be improved by using low-pass filtered EOG signals (7.5 Hz) as compared to unfiltered EOG signals [43]. Thus, the performance of EOG-based algorithms could be highly affected and dependent on the low-pass EOG signal filtering. Various studies have utilized different low-pass frequencies for removal of ocular artifacts from EEG data [35, 43–45, 47, 51, 65, 66]. Table 1 lists the different low-pass filtering used in those studies, note though that in literature, there is still no evidence on the optimal low-pass frequency of EOG signals. In this light, it is very important to investigate the optimal low-pass filtering for EOG signals before using them in artifact removal process, not only for efficient removal/reduction of artifacts but also for enhancement of the classification accuracies and communication rates of the current BCI systems.In this study, we used unfiltered and twelve different low-pass frequencies to filter simulated EOG signals before using them in artifact removal algorithm. The frequencies with their categorization used in this study are listed in Table2. EOG-based five algorithms from simple [44], adaptive [45, 53], and hybrid [47, 51] categories have been chosen to investigate the task. The performance of each algorithm was evaluated in both the time and frequency domains in order to reach a conclusion for optimal low-pass filtering of EOG signals. In the time domain, the mean square error, mutual information scores, and gain in signal-to-artifact ratio were used as evaluation metrics [41, 61, 62], whereas in the frequency domain, the mean absolute error and mean absolute percentage error were employed [43, 51]. The results for each algorithm with each evaluation metric are shown in Tables 3–10. Time and frequency results indicate that there is a reduction of ocular artifacts by using low-pass filtering on EOG signals as compared to unfiltered EOG outputs. However, for some low-pass filtering, corrected EEG signals specifically from the frontal area (e.g., Fp1) showed distortion in the neuronal component of the EEG signals. Since we argued throughout the paper that optimal filtering will reduce the bidirectional contamination and hence will result in efficient removal of artifacts from EEG data, therefore, we tried to analyze the effect of bidirectional contamination in Figure 4. It can be seen that by using filtered EOG, specifically 7 Hz (blue line), the distortion in EEG for both time and frequency domain is very less as compared to that of when using unfiltered or high values of low-pass filtering (e.g., 35 Hz, green line) of EOG signals. Moreover, the effect of bidirectional contamination can also be analyzed by observing the errors obtained using different performance metrics, that is, lower the errors means less effect of bidirectional contamination. Furthermore, statistical testing is utilized by means of paired t-test to check any differences and improvements between all low-pass filtering results. From time domain metrics, results of paired t-test were displayed only for mean square errors (mutual information and signal-to-artifact ratio also have similar results) and in the frequency domain for mean absolute errors of all frequency bands. In each low-pass filtering class, all the results were averaged before applying statistical testing. From Figure 5, the results of this statistical testing revealed that the frequency class 5–8 Hz have low errors and statistically significant results as compared to all other frequency ranges with p<0.05 in most of the cases (bold in purple color). Finally, we applied paired t-test on results of each low-pass frequency belonging to 5–8 Hz class to see whether there are any differences in the results with these low-pass filtering of EOG signals. The results of this analysis are listed in Figure 6. This testing revealed that, although in most of the cases, there are no significant differences in the results (p>0.1), but 6–8 Hz EOG filtering shows statistically increased results as compared to 5 Hz low-pass filter (p<0.1).Although, simulated signals are the primary tool to analyze the performance of algorithms, but validation with real EEG signals is the ultimate goal specifically for applications like BCIs. We also used MI-based BCI signals to further validate the results obtained through the comprehensive analysis using simulated datasets. Four low-pass EOG filtering, one from each group defined in Table2 was selected to analyze the effect of artifact removal on the classification accuracies of the MI-based BCI. Classification accuracies were used as an evaluation metric to verify the results, that is, higher the classification accuracies means that the corresponding low-pass EOG filtering could be used as an optimal frequency class to remove artifacts from EEG signals and to overcome the bidirectional contamination problem. It can be seen that all methods show best results when 7 Hz low-pass EOG filtering was used to remove ocular artifacts. Furthermore, statistical testing is utilized by means of paired t-test to check any differences and improvements between the classification accuracies obtained with each method and all low-pass filtering. The results of this analysis are shown in Figure 7. The results of this statistical testing revealed that the classification accuracies obtained after artifact removal from each method by using 7 Hz low-pass filtered EOG signals show a significant statistical increase as compared to 4 Hz, 12 Hz, and unfiltered EOG with p<0.05 in most of the cases (bold in purple color). Furthermore, results from all low-pass EOG filtering showed statistically significant results when compared to the classification accuracies obtained after the application of unfiltered EOG signals. These results from experimental EEG data not only validate the outcomes from simulated datasets but also support our hypothesis that low-pass filtering should be applied to EOG signals before using them for artifact removal to reduce the effect of bidirectional contamination.The main focus of this study was to analyze the effect of different low-pass filtering of EOG signals on the removal of ocular artifacts from EEG data, but the classification accuracies for real EEG signals could be improved by incorporating more features and by using enhanced classifier. For instance, EEG signals can be divided into different subfrequency bands to calculate more CSP features for each subband. Also, it has been shown previously that the performance of many other classifiers like support vector machine (SVM) is better than LDA but at the cost of more computations. Although the present investigation can be considered helpful in optimizing EOG signal filtering, the comprehensive comparison of the performances of other algorithms with no need of EOG signals for artifact removal still remains to be determined. Furthermore, the performances of BCI systems with EOG-based and non-EOG-based methods also should be investigated to conclude the optimal method for BCI applications. Therefore, in our future studies, we will use simulated and experimental EEG to evaluate artifact removal and BCI performance by including more methods like independent component analysis, canonical correlation analysis, empirical mode decomposition, and wavelet transform.


## 5. Conclusions

The optimal performance of a BCI depends on the effective removal/reduction of ocular artifacts from EEG recordings. Since the efficiency of an algorithm’s removal of ocular activities is highly affected by bidirectional contamination, it is very important to use the optimal low-pass filtering for EOG signals in order to overcome/minimize the effect of bidirectional contamination. In the literature, there is still no evidence on the optimal low-pass frequency of EOG signals. In this study, we investigated the optimal EOG signal filtering for efficient removal of ocular artifacts from EEG data using fifteen artificially contaminated EEG and EOG datasets. Results from statistical testing of this investigation suggest that low-pass frequency from 6–8 Hz could be used as the optimal EOG signal filtering frequency for good results in terms of artifact removal and retrieval of true EEG signals. Furthermore, MI-based BCI datasets were utilized to validate the results of simulated signals. Classification accuracies obtained by class II showed statistically increased results as compared to results from all other classes. Moreover, the performance of each algorithm was enhanced by applying low-pass filtering to EOG signals before using them for artifact removal process. Overall, hybrid algorithms (REGICA and AIR) showed better performances as compared to regression and adaptive filtering methods for both simulated and experimental signals.



---

*Source: 4853741-2018-07-04.xml*