Decoding the Attentional Demands of Gait through EEG Gamma Band Features

Rehabilitation techniques are evolving focused on improving their performance in terms of duration and level of recovery. Current studies encourage the patient’s involvement in their rehabilitation. Brain-Computer Interfaces are capable of decoding the cognitive state of users to provide feedback to an external device. On this paper, cortical information obtained from the scalp is acquired with the goal of studying the cognitive mechanisms related to the users’ attention to the gait. Data from 10 healthy users and 3 incomplete Spinal Cord Injury patients are acquired during treadmill walking. During gait, users are asked to perform 4 attentional tasks. Data obtained are treated to reduce movement artifacts. Features from δ(1 − 4Hz), θ(4 − 8Hz), α(8 − 12Hz), β(12 − 30Hz), γlow(30 − 50Hz), γhigh(50 − 90Hz) frequency bands are extracted and analyzed to find which ones provide more information related to attention. The selected bands are tested with 5 classifiers to distinguish between tasks. Classification results are also compared with chance levels to evaluate performance. Results show success rates of ∼67% for healthy users and ∼59% for patients. These values are obtained using features from γ band suggesting that the attention mechanisms are related to selective attention mechanisms, meaning that, while the attention on gait decreases the level of attention on the environment and external visual information increases. Linear Discriminant Analysis, K-Nearest Neighbors and Support Vector Machine classifiers provide the best results for all users. Results from patients are slightly lower, but significantly different, than those obtained from healthy users supporting the idea that the patients pay more attention to gait during non-attentional tasks due to the inherent difficulties they have during normal gait. This study provides evidence of the existence of classifiable cortical information related to the attention level on the gait. This fact could allow the development of a real-time system that obtains the attention level during lower limb rehabilitation. This information could be used as feedback to adapt the rehabilitation strategy.


Introduction
According to the World Report on Disability provided by the World Health Organization (WHO), more than 1000 million people from the entire world suffer some sort of disability, which represents the 15% of the world population. Between 110 and 190 millions of adults have significant difficulties to perform daily activities. The number of disabled people is rising due to the population ageing and the increase of chronic diseases [1]. Spinal Cord Injury (SCI) is one of the most concerning diseases that lead to motor disability. According to the National Institutes of Health (NIH), among neurological disorders, the cost to society of SCI is exceeded only by the cost of mental retardation [2]. These considerable costs come from the rehabilitation devices, clinical staff and the long time periods required to perform the rehabilitation therapies needed by SCI patients. Emerging from these problems, governments are investing on technologies oriented to improve rehabilitation therapies and many research groups are focusing their studies on this topic [3][4][5].
In classical physical rehabilitation, patients are rehabilitated by therapists or devices that induce on them the movements they cannot do on their own. Over the last years, the field of neurorehabilitation has proved in multiple occasions that this process could be widely improved by involving the patients in a neurological way [6]. The brain is a learning organ, capable of restoring lost neural paths during the rehabilitation process. Using this capability, known as neuroplasticity [7], it is possible to reduce rehabilitation periods and improve the recovery results. To that end, it is necessary to provide some sort of neurological feedback to the patient during the rehabilitation. In [8][9][10], a virtual reality has been used to provide the patients with visual feedback of an avatar performing, simultaneously, the rehabilitation movements. Results show high improvements in the rehabilitation process and also in the patient's motivation and involvement. Another way of providing neurological feedback is through the use of Brain-Computer Interfaces (BCI) [11,12]. BCIs obtain neural information from the brain by acquiring the electrical signals on the scalp. This information can be used to obtain patient's intentions and mental state. Using patient intentions to trigger the rehabilitation strategy enhances the effectiveness of the therapy in terms of time and performance [13]. In [14], an upper limb rehabilitation system is tested on stroke patients showing evidence that, among this type of patients, a robotic rehabilitation based on a motor imagery BCI result in greater motor improvement than standard robotic rehabilitation. On the other hand, there are BCIs designed to obtain parameters related to cognitive mechanisms like concentration, workload and attention. These parameters play an important role during rehabilitation as the mental state of the patient has a huge influence over the therapy performance [15].
In the case of lower limb rehabilitation, attention has been proved to be an important parameter that affects the final performance of the therapy [16,17]. Slow gait and poor stability have been also associated to a low attention and cognition capabilities [18,19]. In [20,21], dual-task strategies are also used to test the gait stability and variability in elder people according to the attention paid on the gait process showing similar results. In BCI studies, attention has been a well studied parameter, but it is susceptible to be subjectively understood depending on the focus of each specific work. Attention and cognitive mechanisms have been related to different electroencephalographic (EEG) phenomena produced mainly on the alpha, beta and gamma bands. Power spectral variations on alpha and beta bands have been related to changes on brain's attentional demands [22][23][24] and the increase of theta/beta bands power ratio and alpha band peaks have been used in [25] to diagnose attention-deficit and hyperactivity disorders. Phase synchronization on gamma band has been related to visual-spatial selective attention [26][27][28][29], which is the act of focusing on a particular object, action or stimulus for a period of time, while simultaneously ignoring irrelevant information that is also occurring. Also, amplitude changes of evoked potentials like P300 have been related to the level of attention paid to an external stimulus [30,31].
In the current work, a evaluation study is performed to find a relationship between the attention level that a user pays to the gait and the EEG brain signals. The main goal of this study is to set up the basis for a future system capable of classifying in real time the attention level of SCI patients during lower limb rehabilitation. The development of such system could have important implications on the rehabilitation strategies, being a first step to a rehabilitation system capable of changing the therapy parameters in order to fit the mental state of the patient. To perform this study, cortical information from experiments where healthy users and incomplete SCI patients walk on a treadmill are extracted during several tasks related to the attentional demands on gait. This study faces two critical points. The first one is to evaluate the validity of the EEG signals measured during walking, which is a subject of controversy in this field. In [32] it is stated that it is possible to measure good quality EEG signals during movements like walking, cycling and sitting to measure attention by evaluating P300 amplitudes (1-4 Hz), on the other hand, a recent study [33] claims that, during treadmill walking, the EEG signals are polluted with movement artifacts that change depending on the subject, conditions, electrodes and strides, affecting mainly low frequencies (*1 − 8 Hz) and which cannot be removed just using the current artifact removal techniques. The second point is to correlate the results obtained with the state-of-the-art studies to find a relationship between the attention paid on the gait an the current investigations on attention and cognitive mechanisms.
To deal with these problems, all frequency bands associated to cortical signals are evaluated to find those optimal to classify the attention parameters of interest and, after that, several classifiers are tested to get the one that provides better results. The experimental conditions have been defined to avoid the maximum amount of noises and the signals are processed with artifact removal algorithms in order to reduce, as much as possible, their contributions and evaluate only the cortical contribution.

Materials and Methods
Acquisition system EEG data are acquired through 32 channels using pseudo-active electrodes located on the scalp through an elastic cap (g.GAMMAcap, g.Tec, GmbH, Austria) with the following spatial distribution: FZ, FC5, FC3, FC1, FCZ, FC2, FC4, FC6, C5, C3, C1, CZ, C2, C4, C6, CP5, CP3, CP1,  CPZ, CP2, CP4, CP6, P5, P3, PZ, P4, P6, PO7, PO3, POZ, PO4 and PO8 according to the international system 10/10 using AFz position as ground and a monoauricular reference in the right earlobe. Electrical signals are preamplified (g.GAMMAbox, g.Tec, GmbH, Austria) before their 1200 Hz digitalization using two commercial amplifiers (g.USBamp, g.Tec, GmbH, Austria). A 50 Hz Notch filter is also applied to remove the power line interference. The user is wearing the EEG cap which is connected to the amplifiers using a couple of extensor wires. A treadmill is used to create a steady gait pattern. An antistatic wrist strap connects the user's wrist to the amplifier ground to avoid noises related to the treadmill vibration. In front of the user, a screen placed on the treadmill is connected to the recording computer to provide visual guidance to the user during the experiments. The user walks on a treadmill while a screen located at the eye-level provides guidance to perform different attention-related tasks. During the experiment, EEG signals are recorded from 32 channels located over the cortex through the g.GAMMAcap. Electrical signals are preamplified through 2 g.GAMMAboxes located in the user hip and digitalized in the g.USBamplifiers. An antistatic wrist strap connects the user's wrist with the amplifiers ground to remove treadmill's electrical noise. The digitalized data are recorded in a computer system.

Experimental paradigm
The experiment performed is based on the dual task paradigm, commonly used in the literature to evaluate attention measurements [34,35]. In this work, users walk on a treadmill at 2 km/h and 0 degrees of tilt while they are asked to perform 4 different 1-minute tasks that induce changes in the attention paid to gait. In Fig 2 a graphical representation of a run is shown. During the first task the subject walks normally, looking straight ahead, without any distraction. This task represents a standard attention level on the gait as the user is not fully focused on the gait but not distracted by any other task. In the second and the third task, the subject is asked to perform several mathematical operations and to watch a video on the screen, respectively. Mathematical operations are composed by simple additions and subtractions with numbers between 1 and 9 presented to the user in a friendly interface. The video is presented soundless with subtitles to avoid the appearance of auditory potentials and to keep the interest of the user. Both tasks represent a low attention level since the user is focusing on a non-related gait task. Finally, during the fourth task, the user is asked to walk following some marks located on the treadmill trail. These marks have been consciously located following an unsteady distribution to force the user to keep a high attention level on the variable gait pattern. The order of task performance was the same shown on

Participants
Twelve healthy users have performed 2 sessions of the experiment, 4 women and 8 men, all of them right handed, with ages between 22 and 32 (26.3 ± 3.8). Also 3 incomplete SCI patients have performed 1 session of the experiment, all of them right handed men, with ages between 26 and 58 (44 ± 16.3). Healthy users are degree and Ph.D. students from the Miguel Hernández University of Elche with no known diseases, and patients are from the National Hospital for Spinal Cord Injury (Spain) and they have incomplete SCI with motor lesions between C5 and C6 level. All patients selected were able to walk by themselves or using simple assistive devices like crutches or walkers. They do not suffer from cerebral injury so their brain processes should not present huge differences from those manifested on healthy users [36]. All users have been previously informed about the experimental procedure and they have signed an informed consent according to the Helsinki declaration. The experimental procedure were approve by the ethics committee of the Miguel Hernández University of Elche (Spain).
Time processing. Although cortical processes can be described and generalized, electrocortical signal's amplitudes experience huge changes in the time domain depending on the user and the recording day [38]. Also, recording EEG signals during human walking induces several sources of noise that contaminate the cortical signals of interest [33]. For these reasons it is important to apply processing methods to remove the artifacts affecting the brain signals and also a standardization of the signals to make possible the comparison between users and sessions.
First, all data are bandpass filtered between 0.5 and 100 Hz to remove blink artifacts associated to low frequencies [39] and electromyographic (EMG) artifacts associated to high frequencies [40]. After that, each channel from every session is visually inspected to find outliers. On Fig 3A, the first session of user 1 is shown to graphically see the artifacts found on some channels during a recording. The simplest way to deal with these channels is their rejection with the consequent loss of information. In this work, to keep some information from these channels and preserve the data dimensionality and the final number of features extracted, the noisy channels are going to be removed and reconstructed using data from a spatially distributed set of recordings sites by means of spatial position as in [41]. Using this method no extra information is added to the signal (only redundant information from surrounding electrodes), allowing the recovery of some lost information from a specific scalp area using the information left on the neighbor areas due to volume conduction [42,43]. Fig 3B shows the signals after replacing  the contaminated areas. The amount of reconstructed information from noisy electrodes is the 2.70% of the total available data from both healthy users and patients. The next step is the standardization of the signals in the time domain. To do that it is necessary to find a parameter that represents the brain signals and that does not experience huge variations during a session. Also, an important aspect of this process is to reduce to its minimum the loss of information after applying the standardization processes. The selected parameter to perform this standardization has been developed during this research and for simplifications it is going to be referenced as Maximum Visual Threshold (MV Threshold). This parameter is computed by windowing an EEG channel in L-samples epochs and averaging the maximum value of each epoch. Eq (1) shows the formula to obtain the MV Threshold for the electrode e where X e ðiÀ1ÞÁLþ1:iÁÃLÞ is the L-samples epoch number i for i = 1, 2, 3. . .N being N the total number of epochs.
In Fig 4A a highpass filtered EEG signal is shown in blue. The red line represents the MV Threshold calculated for L = 1200 samples (1 second). The green dotted line represents the step-by-step computation of the MV Threshold. During the initial epochs, the average variation (green line) presents an unsteady behavior but after averaging several tens of epochs, the variation is very small and it is robust against isolated high amplitude peaks. Since the value of the MV Threshold depends on the epoch width L, it is important to select this L value taking into account the width of the epochs that are going to be used to obtain the features of our signals. Fig 4B shows the MV Threshold for different values of L. The black point shows the one selected for the current work. For L values close to 1 (0.83 milliseconds), the value of MV Threshold is really close to the signal average which is almost 0. Prior to L = 110 (92 milliseconds) the value of MV Threshold presents an accelerated increase and it becomes more stable after this point. This phenomenon is closely related to the amount of spectral information of the signal. The minimum epoch length needed to get all the spectral information desired from a signal is computed as t min = 1/f min , being f min the minimum frequency of interest in the signals analyzed. In this case, the L = 110 is associated to low frequencies (*10 Hz) corresponding to the α band, the activation of which is expected from EEG data obtained during motor tasks, implying that most of signal information appears in this band. However, current EEG data are bandpass filtered between 1-100 Hz which means that the minimum epoch length to have access to all the bandwidth is 1 seconds (L = 1200) which also in the range where MV Threshold presents a stable behaviour.
The MV Threshold is calculated for each electrode and it can be considered constant during a whole session. For each session, 32 MV Thresholds are computed, one per channel with L = 1200 samples. The signal standardization is computed according to Eq (2) where V(t) e is the time domain EEG signal of the electrode e and MVThreshold j is the MV Threshold of the channel j, with j = 1, 2, 3, . . .Ch being Ch the total number of channels.
The values of the MV Thresholds are different on each session and user, but they remain constant during a single session. The use of this standardization makes possible the comparison of data between users without the loss of information. The bandwidth of interest can be recovered from the signals and the variation of spectral power between electrodes is not lost by averaging the MV Thresholds obtained for each of the channels.
Features extraction. Once the signals are filtered and standardized, breaks from each session are removed and the remaining data is divided into one minute tasks (4 per run). This division is performed prior to epoch segmentation to avoid the existence of epoch containing cortical information of two consecutive tasks. Each task is divided into 1 second epochs (  proposed in this work related to frequency bands comparison and task classification. After that, the spectrum of every epoch is computed using the Maximum Entropy Method presented by Burg in [44]. This method is an autoregressive calculation technique introduced to compute the spectrum finite-length epochs of real time sampled data. The AR-parameters are computed by minimizing the sum of the square forward and backward prediction errors [45]. This method reduces the minimum entropy components of the signal which are associated to those components than provide less amount of information reducing this way the random components (usually associated with electronic noise) and improving the components with more information. For each epoch the spectrum is computed between 1 and 100 Hz with a spectral resolution of 1 Hz.
The computed spectrum of the epochs is used to obtain 6 different features sets according to the frequency bands analyzed. As it was mentioned at the beginning of the processing section, the frequency bands evaluated are delta δ (1-4 Hz), theta θ (4-8 Hz), alpha α (8-12 Hz), beta β (12-30 Hz), low gamma γ low (30-50 Hz) and high gamma γ high (50-90 Hz). For each frequency band and epoch, 32 features are obtained as the sum of all frequencies between the lowest and the highest frequency of each band (for each channel).
This feature extraction process is performed for each individual session, i.e twice per each healthy subject and once per each patient.
Frequency band selection. To take an objective decision regarding the optimal frequency band to classify the attention on the gait, there are several possibilities based on the use of mathematical parameters and indices that provide information related to class differentiation. In [46], Millán et al. introduce a coefficient that expresses the separability between 2 classes according to their mean and standard deviation across all the epochs. To use this coefficient it is necessary to obtain a single value of mean and standard deviation per class which is difficult when each task is represented by more than one feature per epoch. A more suitable parameter for the present study is the Battacharyya distance (Bdist) which was firstly introduced by Battacharyya in [47]. This parameter provides a measurement of the similarity of 2 sets of features and it is closely related to the Battacharyya coefficient which is a measure of the amount of overlap between 2 statistical samples or populations. According to literature, this measurement has widely proved its reliability for signal selection purposes [48][49][50]. As mentioned, Bdist is used to compare 2 classes, however in the current work, to obtain information related with the separability of 4 classes, the Bdist has been calculated for all the possible combinations of the 4 classes.
Since the Bdist can provide any value in the range 0 < Bdist < 1, it is necessary to use a selection criteria to decide which frequency bands present the higher class separability. In [51], Choi and Lee have performed a study where they represent the Bayes error [52] of 2 task classification against the Battacharyya distance between 2 classes. As a result, they obtain a logarithmic behavior being the Bayes error between *10% and 30% for Bdist = 0.5, *1% and 7% for Bdist = 2 and *0% and 2% for Bdist = 3.5. For this work, 2 classes are going to be selected as highly separable if their Bdist = > 3.5. Only features from healthy subjects are used in the computation of Bdist values. They represent most of the cortical data recorded during this work and the attention level paid from patients is expected to be higher during non-attentional tasks. Best features selected by this method are applied to patients data during the classification stage.

Classification
Five different classifiers have been considered to test the capability of distinguishing between the 4 gait attention tasks. The classifiers used are properly enumerated below: 1. Support Vector Machine (SVM) [53] with Radial Base Function kernel with C = 512 and y = 0.002 (parameters obtained in a previous work [54]).

Decision Tree Learning (DTL) [58].
To validate the success rate results obtained by each classifier for each session (2 session x 12 healthy users + 1 session x 3 patients = 27 sessions) an 8-fold cross validation has been used where each run has been used as fold.

Chance level computation
To validate the results and to select the best classifier, it is necessary to confirm the significance level between the classification results and the chance level. Applying the simplest mathematical statement, the chance level for a 4-task classification system assuming class equality (which is the case of the current study) is 25%. But for real finite data analysis the chance level presents several variations according to the specific conditions of the data studied. In [59], Müller-Putz et al. introduce a mathematical method to calculate the range of values corresponding to the chance level according to the number of tasks classified and the number of epochs used in the classification stage for a 2-tasks classification system. According to Müller-Putz's work, the confidence interval around the expected chance level can be calculated using Eq (3) where n is the number of epochs, z 1À a 2 is the 1 À a 2 quantile of the standard normal distribution N(0, 1), α is the level of confidence required (typically 0.1 and 0.05) andp is the unbiased estimator computed according to Eq (4). In this case, n is still the number of epochs and X is the averaged probability of all the individual probabilities of a correct classified epoch X i with a classifier that performs a random classification. The Eq (5) shows the mathematical formula to compute X , which according to the definition provided corresponds to the mathematical value of the chance level for an ideal random classifier (50% for 2 tasks). For more than 2 classes the model can be extended as shown in [60].

Features spatial representation
After the signal epoch segmentation and features extraction, a single 32 features vector is computed for each frequency band and task as the average of the epochs of all sessions (from both healthy and patients data). This makes 6 frequency bands per 4 tasks, a total of 24 feature vectors. The spatial distribution of these features has been represented in Fig 5 in order to visually appreciate the difference between each task depending on the frequency range evaluated. On this representation it is possible to appreciate some differentiation between tasks and frequencies. Focusing on the minimum and maximum spectral power values used to represent each frequency, it is shown that the spectral power of EEG signals gets lower when the frequency evaluated increases. These results, based on visual inspection and spectral power level, are quite subjective in order to decide which frequency bands provides more information to differentiate between the attention tasks presented.

Frequency band separability
To obtain an objective measurement related to class differentiation, the Bdist is used. In Table 1 this parameter is shown for each 2-tasks combination and frequency band. The letter labelling for the class combination is the following one: A for "Normal Walking" class, B for "Mathematical Operations" class, C for "Video" class and D for "Following Marks" class. All the tasks combinations that achieve the criteria selected (Bdist = > 3.5) have been highlighted in bold in Table 1. Also the first and second maximum values for each combination are marked with one or two asterisks respectively. Table 1 shows high Bdist values for all task combinations that include task A (Normal Walking) which means that it presents a huge separability from the rest of tasks on every frequency band. On the other hand, the remaining combinations present a low Bdist for low frequencies achieving the selection criteria only on γ low band (except for BC combination) and γ high band. Another interesting fact obtained from this table is that, for all frequencies, the combination BC ("Mathematical Operation" class and "Watching Video" class) has the lowest separability features. Taking into account that these tasks or classes have been defined as the same level of attention (low gait attention level), these results suggest that, indeed, both tasks present lots of similarities. For that reason, both γ low and γ high bands are going to be selected for the classification stage where several classifiers are going to be evaluated in the classification of the 4 tasks mentioned.

Classification results
The 5 classifiers previously described are applied to on the data. The average and the standard deviation for each user and classifier after an 8-fold cross validation are computed for γ low and γ high bands and represented on Table 2. For patients, 1 session is used to compute the success rate while results from 2 sessions are averaged in the case of healthy users. The total averaged values for each classifier and frequency band are also computed. These values are computed using all the single fold values obtained from every subject. From these results it is possible to approximate averaged classification results, e.g from high gamma band frequency features and KNN classifier, a 66.79% and a 59.08% of success rate are obtained from healthy subjects and patients respectively.
Significance against chance level. In the current work, for each session there are 3840 epochs but only 960 are used to test the classification results on each fold iteration. The chance level range of a single fold (from the 8-cross fold validation) is computed for n = 960, X ¼ 0:25 (4 tasks) and α = 0.05 (according to Eq 3), obtaining the range 23.9 < successrate < 26.3 ( Fig  6). Any iteration of the 8-cross fold validation providing a classification result within this range should be considered random.   Decoding the Attentional Demands of Gait through EEG Gamma Band Features value of the chance level for each group (a 192/24 vector, respectively, with random values from the chance level range computed). The significance of these success rates is analyzed from these results. To that end, the chance level is compared to the success rate obtained for each classifier running a Wilcoxon Sum-Rank Test with a confidence interval of 95% and then applying a Bonferroni correction for multiple comparisons. This analysis shows that the success rates are predominantly above chance levels (p < 0.05) for all 5 classifiers.
Significance between classifiers, features sets and subjects. To test the significance between classifier outputs, feature sets and subjects, the same statistical approach is applied. Table 3 shows the significance p values of comparing the same classifier for 4 different conditions: 1) γ low features: healthy against patients; 2) γ high features: healthy against patients; 3) Healthy: γ low against γ high features; 4) Patients: γ low against γ high features. Table 4, on the other hand, shows the significance p values for all the classifier combinations for the 4 graphs of Fig  7. On both tables, the significant p values are marked in bold after applying a Bonferroni correction.

Discussion
On the spatial distribution of the averaged features shown in Fig 5 it is possible to appreciate some differentiation between tasks and frequencies. Focusing on the minimum and maximum spectral power values used to represent each frequency, it is shown that the spectral power of EEG signals gets lower when the frequency evaluated increases. In δ and θ bands, the distributions present huge deviations between tasks and the main activation is located on the frontal lobe usually associated to blinks, head movements and other physiological artifacts [61]. These results fit the statements of Kline et al. in [33] claiming that motion artifacts have higher contributions on lower frequency bands and present huge variations depending on the subject and the walking speed. α band shows a main activation on the motor area which is an expected behavior during walking. According to literature, the motor cortex side activated depends on the real or imagined motor movement performed [62,63]. For that reason, during walking, left  and right areas of motor cortex are alternatively activated disguising any changes related with the attention. In the β band, the spectral power is softer in the motor area and experiences a smooth increase in the occipital area. Finally, for γ low and γ high bands there is a significant reduction on the motor cortex where the power and distribution seems to be related to the attention level classified. Compared with the motor cortex, the occipital area is much more active suggesting that visual information is being processed. These results and those provided in Table 1 about class separability suggest that the measurement of the attention level on the gait is closely related to the selective attention. As it was stated at the beginning of the paper, selective attention is described as the attention paid to an external object, action or stimulus while another object, action or stimulus, simultaneously happening, is ignored. In this case, the gait process and the environmental visual information are the event and action, respectively, confronted, and the gait attention changes depending on which event is more or less ignored or attended.
Classification values on Fig 7 and the test performed show that all classifiers tested provide significant values in the classification of this attention level, confirming the existence of discernible information related to the attention level on the gait obtained from γ band features. On rows 1 and 2 of Table 4, each classifier from healthy users is compared with its analog on patients showing, in most cases, significant differences between them. These results imply that the classification results from patients are significantly different than those obtained from healthy. This significance could be a consequence of the inherent increased attention on gait presented on incomplete SCI patients who are always more focused on the gait than healthy subjects regardless the distraction event presented to them [18,19]. Both healthy subjects and patients performed the same tasks and the classification model was trained according to them, but patients' attention on gait was always higher than healthy's during standard and low  Table 4 show the same classifier comparison but this time between frequency bands, showing no significant values except for classifiers C2 (NB) and C3 (LDA) on healthy users. These values suggest that the information is similar in both bands and their combination in the feature extraction stage could improve the final classification results. Finally, on Table 3, for each subject group (healthy and patients) and frecuency band (γ low and γ high ), the classifiers are compared between them showing in most cases no significant results between them. In that case, the optimal classifiers selection could be performed by the maximum average success rate. Focusing on Fig 7, classifiers SVM, LDA and KNN provide, in general the best results, being SVM performance slightly higher on healthy users and LDA and KNN on patients.

Conclusion
In this paper, different frequency bands extracted from data of 10 healthy users and 3 incomplete SCI patients have been evaluated to find a relationship between cortical signals and the cognitive mechanisms related to the attenton during gait. The study emphasizes the difficulties to decode cortical information on low frequency (δ and θ) bands due to the artifacts generated by movement and also to classify the cortical information on α band related to attention which is covered up by the motor rhythms produced during gait on cortical signals. On the other hand, information related to selective attention mechanisms that relate the attention to the gait and the attention to the external environment have been found in the γ band. The band selection has been performed by computing a separability index between tasks using the Bhattacharyya distance and the features extracted have been tested using 5 classifiers with an 8-fold cross validation. Final results provide an average value of 67% of success rate classification for 4 attention tasks. Success rate values have been compared to chance level, showing clear significance. Healthy users have obtained significantly higher classification results suggesting that there was less class separability on patient's which is related to the higher level of attention they pay on non-attentional tasks due to their motor disabilities.
This work settles a basis for gait attention classification. For future works this offline analysis should be performed with a higher patient population. In addition, the features extraction should be improved combining both low and high γ bands and using mathematical algorithms for feature selection and reduction. The number of tasks used to train the model should be reduced to 3, removing one of the high gait attention tasks and improving this way the final success rate. Current processing algorithms should be changed to fit real time condition in order to perform online tests with healthy users and patients. Also the the reference electrode standardization technique (REST) [64] should be tested on the data to improve the quality of EEG signals in terms of spatial resolution and artifacts' identification. And, as a final goal, an online system that provides the attention level of an incomplete SCI patient under exoskeleton Decoding the Attentional Demands of Gait through EEG Gamma Band Features rehabilitation could be used to change, in real time, the parameters of the rehabilitation to fit better to the cognitive state of the patients and improve the rehabilitation performance by involving them in a deeper level.