Discrete Dynamics in Nature and Society
Volume 2018 (2018), Article ID 2143873, 19 pages
https://doi.org/10.1155/2018/2143873

Research Article

Feature Extraction of EEG Signal upon BCI Systems Based on Steady-State Visual Evoked Potentials Using the Ant Colony Optimization Algorithm

S. M. Fernandez-Fraga,¹ M. A. Aceves-Fernandez,² J. C. Pedraza-Ortega,² and S. Tovar-Arriaga²

¹Department of Computer Systems Instituto Tecnológico de Querétaro, Av. Tecnológico s/n, Centro, CP 76000, Santiago de Querétaro, Mexico
²Department of Engineering, Universidad Autónoma de Querétaro, Cerro de las Campanas s/n, Centro Universitario, Las Campanas, CP 76010, Querétaro, Mexico

Correspondence should be addressed to M. A. Aceves-Fernandez; marco.aceves@gmail.com

Received 9 December 2017; Revised 16 April 2018; Accepted 5 May 2018; Published 26 June 2018

Academic Editor: Seenith Sivasundaram

Copyright © 2018 S. M. Fernandez-Fraga et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This work presents the use of swarm intelligence algorithms as a reliable method for the optimization of electroencephalogram signals for the improvement of the performance of the brain interfaces based on stable states visual events. The preprocessing of brain signals for the extraction of characteristics and the detection of events is of paramount importance for the improvement of brain interfaces. The proposed ant colony optimization algorithm presents an improvement in obtaining the key features of the signals and the detection of events based on visual stimuli. As a reference model, we used the Independent Component Analysis method, which has been used in recent research for the removal of nonrelevant and detection of relevant data from the brain’s electrical signals and also allows the collection of information in response to a stimulus and separates the signals that were generated independently in certain zones of the brain.

1. Introduction

A brain computer interface (BCI) enables direct communication between a brain and a computer translating brain activity into computer commands, thus providing no muscular interaction with the environment [1]. The development of brain computer interface (BCI) systems is mainly based on the information obtained by processing brain signals obtained by electroencephalogram (EEG). Further to obtaining signals, extraction and feature selection steps are vital for the development of applications; the BCI systems based on steady-state visual evoked potentials (SSVEP) have the advantage over other BCI systems because they have a better signal-to-noise ratio (SNR) and faster transfer rate information (ITR) [2]. Additionally it does not require intensive training and requires fewer EEG channels for application development and immunity from artifacts (flicker, movement joints, etc.) [3, 4]. BCI-SSVEP systems are characterized by using stimuli at different frequencies sensitive to the dynamic processes of the occipital part of the cerebral cortex.

In recent years, algorithms based on swarm intelligence (SI) have recently emerged as a family-based wild populations of ant colony, bees, and swarm individuals algorithms which are capable of producing low give robust computational cost solutions to various complex and optimization problems [5, 6]. The ant colony optimization (ACO) algorithm proposed by Dorigo et al. (1996) is inspired by the foraging behavior of some species of ants. Ants lay pheromones on the ground in order to make a shorter path which eventually will be followed by other members of the colony which increases the probability of other ants following the same path resulting in an optimal way from the nest to the food source. Biologists have shown that many of the behaviors observed at the level of colony of social insects can be explained by fairly simple models in which only communication (called stigmergy) is present. The idea behind the ant algorithms is to use some form of artificial stigmergy to coordinate the societies of artificial agents [7].

Deneubourg et al. (1990) proposed the original model Ant Systems (AS). Several different aspects of the behavior of ant colonies have inspired different types of ant’s algorithms, AS Dorigo, et al. (1991), ELITIST AS Colorni et al. (1992), ANT-Q Gambardella and Dorigo (1995), MAX-MIN AS Stützle and Hoos (1996), RANK-BASED AS Bullnheimer et al. (1997), HYPER-CUBE AS Blum et al. (2001) [8].

Recently, ACO has been used in biomedical data processing: Behnam and Pourghassem (2015), brain rhythms detection with periodogram patters; Erguzel et al. (2014), feature selection method for QEEG data classification; Huang et al. (2012), feature selection method to classify hand motion surface electromyography; Korürek and Nizam (2010), clustering analysis of ECG arrhythmias; and Burša and Lhotska (2007), biomedical clustering process.

Resulting research ICA has been used in order to reduce the calculation time by removing and constructing irrelevant and relevant data from electrical brain signals obtained by EEG to improve the performance of BCI systems [11]. In the present work, AS model will be used as an improvement in the EEG signals obtained by visual events and ICA will be used as reference model to evaluate the performance of the proposed method.

2. Background

2.1. Independent Components Analysis (ICA)

ICA was introduced first by Comon (1994), with further developments by Bell and Sejnowski (1995). Hyvärinen (2012) and Acharya and Panda (2009) provide an overview of some recent developments in the ICA theory since 2000 on many fields as neural computing, signal processing, and statistics. In the field of biomedical sciences, Isomura and Toyoizumi (2016) used a local learning rule for ICA to understand the processing of neuronal information through natural images and their applications in neuromorphic engineering. Gramfort et al. (2014) presented a medical-software to help reduce EEG signal artifacts. Mahajan and Morshed (2015); Akhtar et al., Mammone at al., and Safieddine et al (2012) used ICA removal artifact from multichannel EEG data. Oostenveld et al. (2011) presented open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data based on ICA. Kong et al. (2008) used ICA in microarray gene expression data applications. Hung et al. (2005) used it in recognition EEG motor imagery signals. James and Hesse (2004) used it to separate measured multichannel EEG signals. Vigário et al. (2000) and Tang et al. (2002) used it to identify the artifacts and signals of interest obtained by magnetoencephalograms (MEG). Jung et al. (2000 and 2002) used it on EEG signals. Liebermeister (2002) used it to analyze gene expression patterns during the yeast cell cycle and in the identification of human lymphocytes. Martoglio et al. (2002) applied it to ovarian cancer data. Saidi et al. (2004) showed that clustering in the ICA components obtained clusters biologically more reasonably clustered than clusters obtained by PCA (Principal Component Analysis).

It was found that ICA extracts more significant and better interpretable components than PCA, and even an unexpected experimental artifact was discovered. In addition, when applied to enzymatic activities, ICA was able to provide components of greater discrimination and with greater significance than PCA components [10].

Resulting research ICA has been used in order to reduce the calculation time by removing and constructing irrelevant and relevant data from electrical brain signals obtained by EEG to improve the performance of BCI systems [11].

ICA is a method based on the statistical independence of the components that form a signal. Statistical independence is the main idea for the entire construction of the method, since we want to find within our mixed signals those that are statistically more independent than the others [9]. Statistical independence is defined as follows: let be a set of random variables with probability density function (PDF) ; these variables are mutually independent if

where is the marginal density function of (Figure 1).

Figure 1: It shows a set of independent signals with non-Gaussian PDF [9].

Let a vector of components; we assume that this vector is produced as a linear combination (mixture) of independent signals denoted by . We will assume for the implementation of the algorithm that each of the signals is statistically independent of the other signals; then

where the coefficients determine a matrix that is known as a mixture matrix, the vector as the vector of mixtures, and the vector of independent components, so we can define the system as

Figure 2 shows the relationship of (3).

Figure 2: represents the mixture of the components of mixture and the set of independent signals [10].

Consider that the mixing components take values that allow to be invertible (Figure 2). There exists a matrix , with coefficients such that it will allow us to separate the signals as

which we can define as

where

ICA method consists of applying an algorithm that allows us to find the matrix in such a way as to be a good approximation of the vector [9, 12].

2.2. Ant Colony Optimization (ACO)

The ant colony optimization algorithms are part of swarm intelligence (SI), which is the field of research that studies algorithms inspired by the social behavior of the insects and other animals. Algorithms based on swarms have recently emerged as a family of algorithms based on nature populations that are able to produce, at low cost and fast robust solutions to various complex problems and optimization [5, 13, 14], as well as some applications of grouping of data [15–18].

Ant colonies and, in general, groups of social insects, despite the simplicity of their individuals, have a highly structured social organization. Because of these organizations, ant colonies can perform complex tasks that in some cases far outweigh the individual capacities of a single ant. The fields of “ant algorithms” are models of studies derived from the observation of the behavior of real ants and these models are used as a source of inspiration for the design of new algorithms for optimization solution and distributed control problems [7].

The main idea is that the principles of self-organization that allow the highly coordinated behavior of real ants can be exploited to coordinate populations of artificial agents that collaborate to solve computational problems. Several different aspects of the behavior of ant colonies have inspired different types of ant algorithms: search for food, division of labor, classification of breeding, and cooperative transport. In all these examples, the ants coordinate their activities through stigmergy, a form of indirect communication by a slight modification of the environment. In the search for food, the ant deposits a chemical (called pheromone) in the soil, which increases the likelihood that other ants will follow the same path. Biologists have shown that many of the behaviors observed at the social insect colony level can be explained through simple models in which only stigmergy communication is present. The idea behind ant algorithms is to use some form of artificial stigmergy to coordinate artificial agent societies [7].

Recently, SI algorithms have been used in biomedical data processing of biomedical data such as EEG and electrocardiogram (ECG) [19–22]. For instance, Cheng et al. (2016) used ACO for reconstructing ECG signals; Salem et al. (2014) used ant colony (AC) for the detection of epileptic seizures on EEG signals. Erguzel et al. (2014) used ACO as a method of classifying data in quantitative EEG signals. Korürek and Nizam (2010) used ACO for clustering MIT-BIH arrhythmias; Burša and Lhotska (2006) used AC and ACO as ECG interpretation.

The first proposed algorithm within the ACO metaheuristic was the Ant System (AS) to solve the Travel Agent Problem (TSP) problem in which we search for the optimal route between a set of cities that the traveler must visit only once and return to the starting point.

2.3. Ant System (AS)

The first proposed algorithm within the ACO metaheuristic was the Ant System (AS) to solve the Travel Agent Problem (TSP) problem in which we search for the optimal route between a set of cities that the traveler must visit only once and return to the starting point.

Formally the TSP can be represented by a weighted complete graph , where is a set of nodes representing cities and the set of arcs between cities. Each arc is assigned a weight (length) which is the distance between city and city , with . The objective in the TSP is to find a Hamiltonian circuit of minimum length of the graph [7]. An optimal solution for the TSP is a permutation of the indexes of nodes such that is minimal, where is defined by

In the AS algorithm, ants are placed in different random cities that simultaneously construct a path for the TSP. Each ant constructs a feasible solution to the problem by applying an iterative probabilistic transition rule called proportional random rule to decide which city to visit next. The probability with which the ant that is currently in the city chooses to go to the city is

where , and , determine the relative influence of the pheromone and the heuristic information, respectively, is the set of neighboring cities the ant k, being in city , has not visited, and the probability of choosing a city outside is 0.

Based on (7) the probability of choosing a particular arc increases with the value of the pheromone associated with the trail. Once all the ants have built a complete tour, the pheromone trail is updated; this is done first by reducing the value of the pheromone in all arcs by a constant factor, and then add pheromone in the arcs that the ants have crossed on their path. The pheromone evaporation is calculated according to

where and is the evaporation rate of the pheromone. The parameter is used to avoid the unlimited accumulation of pheromones in the trail and allows the algorithm to “forget” the bad decisions made previously. If associated value of pheromones decreases, the ants do not choose an arc exponentially based on the number of iterations. After evaporation, all the ants deposit pheromones in the bows they have crossed on their tour:

where is the amount of pheromone deposited in the arc by the ant . is calculated as

represents the length of the trail constructed and is a constant of generation of pheromone [6, 23, 24] (Figure 3).

Figure 3: Ants behavior.

3. Materials and Methods

3.1. Statement of Ethics

For the development of the evidence contained herein, they asked all participants to read and sign an informed consent form and a written letter of confidentiality for information management before participating in this study. The study protocol was based on the manual “Ethics of Scientific Research” at the Autonomous University of Queretaro [25].

3.2. Experimental Protocol

Twenty-nine healthy volunteers, 17 men and 12 women, aged 20-29 and 48-50, participated in the experiment. 17 of them with normal vision and 12 corrected vision. The EEG signals were obtained by a portable, high-resolution multibrand equipment, Emotiv Epoc+ EEG model. The equipment has 16 electrodes, 14 data channels and 2 reference channels, which were placed in the standard positions of the international 10-20 system (Figure 4). The transfer rate is 128 bps of a frequency response from 0.16 to 43 Hz with 14 to 16 bit resolution and dynamic range of +/- 4.17mV per channel.

Figure 4: International 10-20 system for Emotiv Epoc+.

During the experimental trials, the study subjects were seated in a comfortable chair to 70 cm apart from a standard 15-inch LCD monitor (with a screen resolution of 1024 x 768 and a refresh rate of 60 Hz).

3.3. Experiment

The experiment is based on tests performed visual discrimination by Makeig et al. (1999) for the study of BCI systems based SSVEP. The objective of the tests is to obtain data related stimuli in the brain signals during an exercise of simple attention; search to find the difference between attention-related and unrelated stimuli. The test subjects must discriminate between different types of stimuli that are presented at high speed.

Study subjects set their gaze on a cross, which is on top of one five boxes at a time that are constantly displayed (Figure 5). Each test block last for 76 s with approximately 5 s initial preparation. Stimuli presented as a series of red circles are briefly presented in any of the five boxes in a random order. The boxes are blue and one is green. The location of the latter is random during test periods. The stimuli presented within the green box are attended stimuli (events); stimuli within the blue boxes are known as unattended stimuli.

Figure 5: Five box task configuration.

Before starting the test, the subject is asked to concentrate whenever a red circle appears in the green box [26]. The outlines of the pictures are of 1.6 cm, displayed on a black background horizontal viewing angles 0° ± 2.7° to ± 5.5°, and fixed. Each stimulus is displayed for 117 ms to one of five empty spaces inside the boxes with a pseudo-random sequence. Stimuli are presented to the subject in a random position of 250 ms, 500 ms, 750 ms, and 1 s, intervals between stimuli (ISI) between 250 and 1000 ms (Figure 6). At the beginning of the test, a black screen is presented for 5 s for the stabilization of the signals; the test begins with the configuration shown in Figure 7(a). In each test, a total of 100 stimuli are presented in total: 20 attended (events) 80 unattended (Figures 7(b) and 7(c), respectively), the last 200 ms stimulus remains on the screen before disappearing [27].

Figure 6: Intervals configuring.

(a)

(b)

(c)

Figure 7: Test 5 boxes. (a) Initial setup. (b) Stimulus unattended. (c) Attended stimulus (event).

3.4. Proposed Model

Figure 8 shows the classical model ICA, where two signals, in this case O1 and O2, are mixed (Figure 8(a)), the independent components of the mixture are obtained, and then a reconstruction of the signal is made based only on the found independent components (Figure 8(b)).

(a)

(b)

Figure 8: ICA model. (a) Independent components obtained from the mixture of signals O1 and O2. (b) Signal reconstructed from the combination of independent components.

Figure 9 shows block diagram of proposed model based on the implementation of the ant colony optimization algorithm, described in Section 2.3 to the signals derived from EEG, based on an SSVEP-based BCI system.

(a)

(b)

Figure 9: Representation of the ACO-ICA implemented model. (a) Signal optimized by the ACO model from the electrodes O1 and O2 obtained by EEG. (b) Signal reconstructed from the independent components obtained of the ICA model.

After acquisition, signals from EEG, the electrodes O1 and O2, which are directly related to the visual processing of the brain signals, are selected. These signals are postprocessed using the proposed ACO algorithm to obtain an optimized signal (ACO-O1O2), Figure 9(a). Subsequently a mixture of the original signals (O1/O2) is performed separately, with the optimized signal ACO-O1O2, its mean mixture O1 with ACO-O1O2, and mixture O2 with ACO-O1O2. The ICA model, described in Section 2.1, is applied to the mixed signals to obtain its independent components.

Finally, the signal is reconstructed based on the independent components obtained, as shown on Figure 9(b). To establish the performance of the proposed model, a comparison of the original reconstructed signals (O1/O2) with the reconstructed signals using the ACO model is performed. Figure 10 showed flowchart proposed model.

Figure 10: Proposed method.

3.5. Optimization Model

The optimization of the EEG signals is based on the solution of Dorigo and Gambardella (1997), Ant System (AS), to solve the TSP. The optimization model of the present work is shown in Figure 11, as well as the algorithm implemented in Figure 12.

Figure 11: Optimization model.

Figure 12: ACO flowchart optimization process.

Illustratively, consider a pair of segments of the O1 and O2 signals (Figure 13) with 20 data points each, which represent the cities for a total of 40 cities. Consider the set of cities as and the distance (the cost of the edge) between the cities . Explicitly, cities can be seen as points of a 2-dimensional space, so we can represent them as and . The distance between cities is calculated by means of the Euclidean distance defined as

We define the matrix of nodes formed by the order of the points in the sequence of the signal and the values of the points of each signal (Figure 14). To the graph , , is a matrix of dimension defined by (Figure 15), which represent the set of distances between the respective points and (Figure 16). Figure 17 shows the graph obtained.

Figure 13: Signals segments O1 and O2.

Figure 14: Node matrix.

Figure 15: Matrix of distances between cities . Each column stores the distances between the points corresponding, where .

Figure 16: Representation of the distances between points.

Figure 17: Resulting graph from the process of calculating the distances of the segments of signals O1 and O2.

We define the probability matrix (probs) as the matrix that represents the probabilities of traversing an arch between nodes (Figure 18(a)). The pheromone matrix is the matrix that represents the amount of pheromone deposited by each ant as it passes through a city (Figure 18(b)). Figure 18(c) shows the representation of the route matrix (tours) that indicates the nodes visited by each ant based on the probability matrix.

(a)

(b)

(c)

Figure 18: (a) Representation of probability matrix. Each box stores the probability of traversing the arch . (b) Representation of the pheromone matrix, each box represents the amount of pheromone deposited by an ant when traversing the arch . (c) Representation of the route matrix, the column on the left represents the departure node at position and the right column represents the destination node at position .

Figure 19 represents the simulation of the path of two ants during 5 iterations and Figure 20 represents the flowchart algorithm.

(a)

(b)

(c)

(d)

(e)

Figure 19: Optimization process of 2 ants during 5 iterations.

Figure 20: AS flowchart.

4. Results and Discussion

The computational simulation of the proposed model was developed in MATLAB V9.1.0441655 (R2016b), as well as EEGLAB V13.6.b. Figure 21(a) shows a segment of 5 s of the original signals O1 and O2 obtained from EEG. The ACO parameters used to experiment were pheromone generation constant , influence parameter , attractive parameter , pheromone evaporation parameter , amount of ants 20, and number of iterations for the tours 5. The optimization signal obtained with the ACO algorithm is showed in Figure 21(b). Test durations were approximately 65-72 s.

(a)

(b)

Figure 21: (a) EEG signals of the electrodes O1 and O2 and (b) representation of the signal optimized by means of the ACO algorithm based on the signals of the electrodes O1 and O2.

The representation of the mixture of the original signals as well as the histogram of each source signal, which approximates the PDF representation of the signals, is shown in Figures 22 and 23, respectively; Figure 24 showed the independent components (IC1, IC2) obtained from ICA process.

Figure 22: Representation of the mixtures of the original signals.

Figure 23: The probability density function of the original signals and the optimized signal.

Figure 24: ICA independents components.

Based on ICA model, Figure 25 shows the original signals O1, O2 and the result of mixing them to obtain the reconstructed signal. We can observe that, for the present example, the reconstruction of the electrode O1 is very similar to the original; in the case of O2, we can appreciate the loss of some elements of the original signal.

Figure 25: Comparison of the original signal and the reconstructed signal (extracted) implementing in the ICA model the mixture of signals O1 and O2.

The comparison of the original signals with the reconstructed signals generated by mixing the signals O1 and O2 with the optimized signal ACO is shown in Figures 26 and 27, respectively. We can observe in both cases less loss of elements in the reconstructed signals.

Figure 26: Comparison of the original signal and the reconstructed signal (extracted) implementing in the ACO model the mixture of signals O1 and ACO.

Figure 27: Comparison of the original signal and the reconstructed signal (extracted) implementing in the ACO model the mixture of signals O2 and ACO.

After ICA process, we can observe in the histograms of Figure 28 that in IC2 the dispersion is much smaller in the case of the mixtures with ACO (Figures 28(b), 28(d), and 28(f)), and in the case of IC1, the dispersion is similar (Figures 28(a), 28(c), and 28(e)), where we can interpret that the mixture of ACO concentrates more elements of the signal.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 28: Histogram and fitted normal PDF and IC1/IC2 component activity. (a) IC1 from mixed O1-O2. (b) IC2 from mixed O1-O2. (c) IC1 from mixed O1-ACO. (d) IC2 from mixed O1-ACO. (e) IC1 from mixed O2-ACO. (f) IC2 from mixed O2-ACO.

Figure 29 shows the power spectrum density (PSD) for the ICA components at each electrode for the study mixtures, more specifically for the corresponding frequency to the stimuli performed (0.895 Hz). The PSD in the IC1 component behaves similarly for O1 in all mixtures (Figures 29(a), 29(b), and 29(c)); however, for O2 the power is higher for mixtures with ACO (Figures 29(b) and 29(c)). The higher the power, the greater the probability of finding the stimulus frequency. In Figure 30 we observed the stimulus frequency (0.895 Hz), and its harmonics (1.79, 2.685, 3.58,...) are present in all the signals generated by the study, which implies that the components retain one of the main features of the signal. The spectrograms in Figure 31 show that the stimulus frequency manifests in all signals, both original and reconstructed.

(a)

(b)

(c)

Figure 29: Power spectrum density of channels and ICA components. (a) Mixed O1-O2. (b) Mixed O1-ACO. (c) Mixed O2-ACO. Notes: shown O1 left side, O2 right side.

Figure 30: Concentrate representation of PSD. Original and components signals. Green line: stimulus frequency; orange lines: harmonics of stimulus frequency.

Figure 31: Location of the stimulus frequency.

Table 1 shows the results obtained from the implementation of the ICA model. The ORIGINAL column represents the mixture of the original signals and the EXTRACT column the signals reconstructed from independent components. The values shown represent the correlation between the signals (1 would represent that a signal is identical to another).

Table 1: ICA results.

USR
O1O2
O1O2
O1ACO
O1ACO
O2 ACO
O2 ACO

ORIGINAL
EXTRACT
ORIGINAL
EXTRACT
ORIGINAL
EXTRACT

A001
0.65275
0.64599
0.75899
0.75391
0.88779
0.88507

A002
0.25834
0.25506
0.60584
0.60343
0.53443
0.53015

A003
0.41599
0.41223
0.60876
0.6058
0.75087
0.74864

A004
0.39095
0.38796
0.47605
0.47328
0.8953
0.89389

A005
0.35604
0.35225
0.66601
0.66341
0.57413
0.56765

A006
0.40944
0.40562
0.73551
0.73328
0.627
0.62342

A007
0.48767
0.48454
0.74114
0.7392
0.69802
0.69534

A008
0.97218
0.97113
0.97218
0.97113
0.56166
0.53865

A009
0.97932
0.97543
0.9797
0.97589
0.99571
0.99511

A010
0.98365
0.98232
0.97366
0.97152
0.97501
0.973

A011
0.38596
0.3816
0.70025
0.69753
0.62466
0.62172

A012
0.9068
0.90512
0.17832
0.18623
0.23504
0.24062

A013
0.62089
0.61857
0.80176
0.80034
0.76475
0.76228

A014
0.55065
0.54694
0.80358
0.80159
0.71116
0.70751

A015
0.68614
0.68214
0.82285
0.82033
0.83757
0.83452

A016
0.91164
0.90966
0.9582
0.95722
0.94837
0.94699

A017
0.81977
0.81738
0.9228
0.92169
0.88674
0.88483

A018
0.94852
0.947
0.98148
0.98091
0.96487
0.96364

A019
0.79789
0.7935
0.89979
0.89741
0.86434
0.85981

A021
0.50396
0.50081
0.73624
0.7342
0.70659
0.70329

A022
0.90516
0.90358
0.90516
0.64055
0.64671
0.64055

A023
0.73834
0.73378
0.95915
0.95833
0.76874
0.76502

A024
0.53992
0.53509
0.78506
0.78224
0.66021
0.6516

A025
0.77073
0.75807
0.89839
0.89247
0.81629
0.80984

A026
0.04135
0.035098
0.42813
0.42234
0.31358
0.30386

A027
0.93049
0.92726
0.97826
0.97722
0.95103
0.94875

A028
0.31401
0.30947
0.66291
0.65989
0.59055
0.5851

A029
0.41971
0.41529
0.73451
0.73181
0.57841
0.56956

A030
0.009597
0.014805
0.66377
0.66023
0.15246
0.61731

Comparing the original O1-O2 signals mixed with the ACO signal, we have 82.76% better correlation O1-ACO versus O1-O2, 10.34% with equal correlation, and only 6.9% better correlation between O1-O2 (Figure 32(a)). We have 86.2% better correlation between O2-ACO versus O1-O2 and a 13.8% better correlation between O1-O2 (Figure 32(b)).

(a)

(b)

(c)

(d)

Figure 32: Signals correlation. (a) Original signals O1-O2 versus O1-ACO. (b) Original signals O1-O2 versus O2-ACO. (c) Reconstructed signals O1-O2 versus O1-ACO. (d) Reconstructed O1-O2 versus O2-ACO.

In the case of reconstructed signals, we have an 82.76% better correlation in O1-ACO versus O1-O2, 6.9% with equal scores, and a 10.34% better correlation in O1-O2 (Figure 32(c)). We obtain an 86.2% better correlation of O2-ACO versus O1O2 and only a 13.8% better correlation between O1-O2 versus O2-ACO (Figure 32(d)).

The results of the correlation of the signals show that the signals optimized with ACO, either in the original signals or in the reconstructed signals, present better characteristics in the selection of components than in the original signals without optimizing.

In Figure 33, we can see that 24 of 29 users in O1 and 24 of 29 users in O2 obtained a better signal reconstructed using the optimization signal ACO (1 mean identical signals). Three users with the same result are shown in O1 and 2 users in O2, and only 2 in O1 and 3 in O2 where the use of the O1 and O2 signals show an improvement in the reconstruction.

Figure 33: Comparative result of the original signals versus signals with ant colony optimization.

5. Conclusions

The methods inspired by swarm intelligence can be applied in many stages of the process of biomedical data processing; in this work, we have presented techniques inspired by the ants with respect to the processing of EEG signal data. We have presented and evaluated the AS method to obtain optimized signals that can be used to extract its characteristics and improve the characteristics of the original signals.

The results obtained show a good optimization capacity when treating the set of EEG signals as a TSP problem solved with ACO; applying our proposal to real biomedical datasets allows us to search methodologies that improve the performance of the applications. Our experimental results demonstrate an improvement in the obtaining of optimized signals. The use of IS-based methodologies in biomedical signals is a new branch that must be exploited in the future.

Conflicts of Interest

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

Acknowledgments

The authors acknowledge the support received from the National Institute of Technology of México, Campus Institute of Technology of Queretaro, and the Autonomous University of Queretaro for the accomplishment of the present work.

References

C. K. Loo, A. Samraj, and G. C. Lee, “Evaluation of methods for estimating fractal dimension in motor imagery-based brain computer interface,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 724697, 8 pages, 2011.

Y.-H. Lu, C.-Y. Chen, P. G. Menon, K.-T. Liu, and H.-H. Lin, “Hemodynamic effects of endoleak formation in abdominal aortic aneurysm patients with stent-graft implants,” Journal of Medical and Biological Engineering, vol. 34, no. 6, pp. 554–558, 2014.

N. Kamel, A. Samraj, and A. Mousavi, “Whitening of background brain activity via parametric modeling,” Discrete Dynamics in Nature and Society, vol. 2007, no. 1, 2007.

J. T. Murphy, “Spehlmann’s evoked potential primer visual, auditory and somatosensory evoked potentials in clinical diagnosis. 1994. Edited by Karl E. Misulis. Published by Butterworth-Heinemann. 243 pages,” Canadian Journal of Neurological Sciences / Journal Canadien des Sciences Neurologiques, vol. 23, no. 01, p. 84, 1996.

C. Blum and D. Merkle, Swarm Intelligence: Introduction and Applications, Natural Computing Series, Springer, 2008.

M. Dorigo and L. M. Gambardella, “Ant colony system: a cooperative learning approach to the traveling salesman problem,” IEEE Transactions on Evolutionary Computation, vol. 1, no. 1, pp. 53–66, 1997.

M. Dorigo and T. Stutzle, Ant Colony Optimization, MIT Press, Combridge, MA, USA, 2004.

M. Dorigo, M. Birattari, and T. Stützle, “Ant colony optimization,” IEEE Computational Intelligence Magazine, vol. 1, no. 4, pp. 28–39, 2006.

D. J. Sandoval Salazar, Anßlisis de componentes independientes aplicado al estudio de la actividad cerebral (Doctoral dissertation), Universidad Nacional de Colombia, 2014.

M. Scholz, Approaches to analyse and interpret biological profile data, Universitat Potsdam, 2006.

M. Mastriani and A. E. Giraldez, “Enhanced directional smoothing algorithm for edge-preserving smoothing of synthetic-aperture radar images,” Journal of Measurement Science Review, vol. 4, no. 3, pp. 1–11, 2004.

A. Hyvärinen and E. Oja, “Independent component analysis: algorithms and applications,” Neural Networks, vol. 13, no. 4-5, pp. 411–430, 2000.

S. Deng, Y. Li, H. Guo, and B. Liu, “Solving a Closed-Loop Location-Inventory-Routing Problem with Mixed Quality Defects Returns in E-Commerce by Hybrid Ant Colony Optimization Algorithm,” Discrete Dynamics in Nature and Society, vol. 2016, 2016.

B. K. Panigrahi, Y. Shi, and M.-H. Lim, Eds., Handbook of swarm intelligence, vol. 8 of Adaptation, Learning, and Optimization, Springer-Verlag, Berlin, 2011.

M. Dorigo, M. Birattari, C. Blum, M. Clerc, and T. Stützle, “2008., Ant Colony Optimization and Swarm Intelligence: 6th International Conference, ANTS, 2008,” in Proceedings of the 2008., Ant Colony Optimization and Swarm Intelligence: 6th International Conference, ANTS, 2008, A. Winfield, Ed., vol. 5217, Brussels , Belgium, 2008.

L. Chen, X. H. Xu, and Y. X. Chen, “An adaptive ant colony clustering algorithm,” in Proceedings of the In Machine Learning and Cybernetics, 2004, vol. 3, pp. 1387–1392, 2004.

P. S. Shelokar, V. K. Jayaraman, and B. D. Kulkarni, “An ant colony approach for clustering,” Analytica Chimica Acta, vol. 509, no. 2, pp. 187–195, 2004.

A. Abraham and V. Ramos, “Web usage mining using artificial ant colony clustering and linear genetic programming,” in Proceedings of the 2003 Congress on Evolutionary Computation, CEC 2003, pp. 1384–1391, aus, December 2003.

M. K. Ahirwal, A. Kumar, and G. K. Singh, “Adaptive filtering of EEG/ERP through Bounded Range Artificial Bee Colony (BR-ABC) algorithm,” Digital Signal Processing, vol. 25, no. 1, pp. 164–172, 2014.

O. Salem, A. Naseem, and A. Mehaoua, “Epileptic seizure detection from EEG signal using Discrete Wavelet Transform and Ant Colony classifier,” in Proceedings of the 2014 1st IEEE International Conference on Communications, ICC 2014, pp. 3529–3534, aus, June 2014.

R. N. Khushaba, A. Alsukker, A. Al-Ani, and A. Al-Jumaily, “Intelligent artificial ants based feature extraction from wavelet packet coefficients for biomedical signal classification,” in Proceedings of the 2008 3rd International Symposium on Communications, Control, and Signal Processing, ISCCSP2008, pp. 1366–1371, mlt, March 2008.

M. Bura, L. Lhotska, and M. Burša, “Ant Colony Inspired Clustering,” in Proceedings of the in Biomedical Data Processing. In 3rd European Symposium on Nature-inspired Smart Information Systems [CD-ROM], Aachen: NiSIS, 2007.

E. A. Mancera-Galván, B. A. G. Licon, and K. Rodríguez-Vázquez, “Optimización mediante algoritmo de hormigas aplicado a la recolección de residuos sólidos en UNAM-CU,” Research in Computing Science, vol. 94, pp. 163–177, 2015.

S. A. Adubi and S. Misra, “A comparative study on the ant colony optimization algorithms,” in Proceedings of the 11th International Conference on Electronics, Computer and Computation, ICECCO 2014, nga, October 2014.

Arellano J. S, Hall R. T., and Hernández Arriaga J. (2014) Ética de la investigación científica. Universidad Autónoma de Querétaro, ISBN: 978-607-513-127-6, 2014.

E. Courchesne, J. Townsend, N. A. Akshoomoff et al., “Impairment in shifting attention in autistic and cerebellar patients,” Behavioral Neuroscience, vol. 108, no. 5, pp. 848–865, 1994.

S. Makeig, M. Westerfield, T.-P. Jung et al., “Functionally independent components of the late positive event-related potential during visual spatial attention,” The Journal of Neuroscience, vol. 19, no. 7, pp. 2665–2680, 1999.