Deep learning and machine learning for EEG signal processing on the example of recognizing the disease of alcoholism ==================================================================================================================== * Rakhmatulin Ildar ## Abstract Alcoholism is one of the most common diseases in the world. This type of substance abuse leads to mental and physical dependence on ethanol-containing drinks. Alcoholism is accompanied by progressive degradation of the personality and damage to the internal organs. Today still not exists a quick diagnosis method to detect this disease. This article presents the method for the quick and anonymous alcoholism diagnosis by neural networks. For this method, don’t need any private information about the subject. For the implementation, we considered various algorithms of machine learning and deep neural networks. In detail analyzed the correlation of the signals from electrodes by neural networks. The wavelet transforms and the fast Fourier transform was considered. The manuscript demonstrates that the deep neural network which operates only with a dataset of EEG correlation signals can anonymously classify the alcoholic and control groups with high accuracy. On the one hand, this method will allow subjects to be tested for alcoholism without any personal data, which will not cause inconvenience or shame in the subject, and on the other hand, the subject will not be able to deceive specialists who diagnose the subject for the presence of the disease. Keywords * EEG alcoholism * EEG machine learning * EEG deep neural networks * machine learning alcoholism * deep neural networks alcoholism * python for EEG * python for BCI ## 1. Introduction According to the World Health Organization, in recent decades the number of patients with alcoholism grew. Research shows that alcohol abuse is associated with behavioral disinhibition, but the neurophysiological mechanisms governing these relationships remain largely unknown. For these reasons, the diagnosis of this disease is difficult. The disease can be identified by the many symptoms. Anuragi et al. [1] and Onarom et al. [2] described the biological processes that occur in the brain during drinking. Ishiguro et al. [3] and Kumar [4] described the physiological consequences of the long-term intake of drinks containing alcohol. These articles showed the complexity of the process operation of the brain during illness and the complexity of diagnosing the presence of this disease. For accurate diagnosis, this disease for specialists needs many private information about patients. But not all patients want to be diagnosed openly. Therefore, the purpose of the research to develop an anonymous method for classifies the alcoholic and control groups by neural networks with the EEG signals dataset. Today, medicine has stepped far enough in this direction. Winterer et al. [5], Patidar et al. [6] and Acharya et al. [7] provided an overview of the EEG signals of patients diagnosed with alcoholism. There is enough information in these works to understand the situation in the field of detection of alcoholism by the EEG signal. But interestingly, despite the seeming knowledge of this issue in the EEG field, many papers on this research have conflicting results. Jeremy et al. [8], showed that compared with men, women are at increased risk of negative physical and neurocognitive correlates of alcohol consumption. In research proved that alcohol abuse has a detrimental effect on the dynamics of EEG suppression of the reaction in the theta range. The opposite conclusion is in the work by Ahmadi et al. [11]. Ahmadi decomposed the EEG signal is into five frequency subbands using the wavelet transform. He showed that there is a lower synchronization in the sub-band of beta frequencies and a loss of lateralization in the sub-band of alpha frequencies in alcoholic subjects. Ahmadi realized classification by machine learning algorithms. But in the research, deep neural networks were not used. But, Paulchamy et al. [15] used all threshold alpha, beta, and theta waves to detect this disease in subject. Ziya et al. [9], implemented software in the Matlab program for classifying by EEG indicator - alcoholic or control subject was presented. The research did not present the result of signal pre-processing. The article does not have enough information about the neural network model. Wajid et al. [10] used EEG data to extract EEG characteristics such as absolute power (AP) and relative power (RP). The classification accuracy of the model is not high. Guohun et al. [12], showed that the areas with electrodes - C1, C3, and FC5 for alcoholic’s groups are significantly different. Mingyue et al. [13] present a new algorithm for analyzing an EEG signal. Mingyue calculated the distinguish non-linear EEG characteristics with alcoholics and controls by the exponential strength Ratio Index (EPRI). But in the research deep neural networks were not used. Joel et al. [15], for the diagnosis of alcoholism by EEG extracted features from four-minute records of EEG of the scalp with eyes closed. In finally the influence of age and gender on the diagnosis of alcoholism was researched. Madhavi et al. [16] noted that increased absolute theta strength in patients with alcohol dependence in all areas of the scalp. Also, Madhavi considered the increase in theta login power in male alcoholics in the central and parietal regions. Anuragi et al. [17] showed that chronic alcoholism is associated with a high frequency of low-voltage recordings. ## 2. Materials and method In this manuscript, the dataset from Henri Begliter (Laboratory of Neurodynamics at the Center for Health at New York State University in Brooklyn, presented publicly, [https://archive.ics.uci.edu/ml/datasets/eeg+database](https://archive.ics.uci.edu/ml/datasets/eeg+database)) was used. This dataset from the research of genetic predisposition to alcoholism. In experiment 64 electrodes placed on the scalp were used (frequency of reading signal of 256 Hz). Two groups of subjects: an alcoholic and a control group were involved in the experiment. Each subject to either one stimulus (S1) or two stimuli (S1 and S2) was subjected. S1 and S2 are a set of images of objects selected from the set of images of 1980 Snodgrass and Vanderwart. The dataset has the following structure, 480 tables in format - csv for training and 480 tables for verification. The table shows 64 electrodes. Each electrode has 256 records with a duration of 1 second, table. 1. View this table: [Table 1.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/T1) Table 1. The structure of the data file for training - Data1, test The table contains the following information: sensor position, sensor value (µV), subject identifier (Alcoholic(a) or Control (c)), matching condition, name(a serial code assigned to each subject), time(inverse of sample num measured in seconds)) The location of the sensors on the head is shown in fig. 1. ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F1) Fig. 1. Arrangement of 64 electrodes during EEG recording Fig. 2. demonstrates the dataset with the library plotly, python3.7, where a - alcoholic group with picture from collection - s1 (alk\_s1), b - alcoholic group with picture from collection - s2 (alk\_s2), c - not alcoholic group with picture from collection-s1 (not\_alk_s1), d - not alcoholic group with picture from collection-s2 (not_alk_s2). ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F2/graphic-3.medium.gif) [](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F2/graphic-3) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F2/graphic-4.medium.gif) [](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F2/graphic-4) Fig. 2. Dataset visualization. X-axis for the number of electrode (from 1 to 64), Y-axis - the record of value for each electrode (256 times with a length of 1 second), Z-axis - the signal amplitude, μV, a -1.alk\_s1, b - 2.1alk\_s2, c - 3.not\_alk\_s1, d - 4.not\_alk\_s2 The ratio of the electrodes on the head and the X-axis (fig.2) is presented in table 2. View this table: [Table 2.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/T2) Table 2. Location of the position of the electrodes on the X-axis for Fig. 1, Fig. 2 This dataset has some artifacts. In fig. 2 - c there is a sudden increase in tension, which can be caused by eye movement or blinking. To exclude this kind of artifact, the principal component method was used. The method of principal components is a multidimensional statistical analysis method used to reduce the dimension of the feature space with minimal loss of useful information. The potential of electrooculography (EOG) is one of the most popular artifacts that occur with eye movement. In this case, the maximum amplitude of artifacts is observed in the frontal leads and decreases towards the occipital leads. In the next researches for [18,19,20] to remove artifacts caused by involuntary eye movements of the subject from a multi-channel EEG, a wide analysis of the main components is used. For these artifacts, it is very difficult to visually find regularity in the presented figures (Fig. 2). Therefore, it is advisable to try neural networks and machine learning for EEG signal recognition. ### 3.2.1. Wavelet transforms We considered the most popular methods of signal preprocessing - wavelet transform and decomposition into a fast Fourier series. In many researches, the fast Fourier transform is used in conjunction with the wavelet transform [25, 26, 27, 28, 29]. The wavelet transform carries a huge amount of information about the signal, but, on the other hand, has a strong redundancy, since each point of the phase plane affects its result. A continuous wavelet transform is defined as the scalar product of the original signal x(t) and the daughter wavelet function ¥*t,a* :(*t*) ![Formula][1] where W(τ,a) - wavelet expansion coefficients; τ, a – п parameters of time shift and scale, respectively; operator * means complex pairing. Child wavelet functions ¥*τ,a*, formed by shear and scale operations of the mother wavelet function ¥*τ* and related to it by the ratio: ![Formula][2] The complex Morlet wavelet, which is the product of a complex sinusoid and a Gaussian, is used as the mother wavelet function. The analytical expression of the Morlet wavelet has the form: ![Formula][3] where ω0 – maternal wavelet center frequency; σ – standard deviation of the envelope of the mother wavelet. Using a stationary wavelet transform, the signal from the electrode to 4 frequency levels delta (0-4 Hz), theta (4-8 Hz), alpha (8-12 Hz), beta (12-20 Hz) were decomposed, fig.3. ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F3) Fig. 3. a - initial signal, b - beta (12-20 Hz), c - alpha (8-12 Hz), d - theta (4-8 Hz), e - delta (0-4 Hz), The received image dataset will be used as input data to the neural network. #### 3.2.2. Fast Fourier Series The continuous Fourier transform, and the discrete Fourier transform have not found wide application in the process of extracting attributes due to their low efficiency, which was explained in the next articles [30,31,32,33,34]. The most popular is the decomposition of the signal into harmonic components using the Fast Fourier transform. For the signal x (n), presented in the form of a sequence of samples, taken with sampling frequency Fs, time moments with numbers n = 0,1, …, N-1, the discrete Fourier transform is defined as: ![Formula][4] where F (k) is the complex amplitude of the sinusoidal signal with a frequency k * Δ f, Δ f = Fs / N resolution (step) in frequency, x (n) are the measured signal values at time instants with numbers n = 0,1, N -1. The result of expanding the signal into a Fast Fourier series with the [scipy.fftpack](http://scipy.fftpack) library is shown in fig.4. ![Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F4.medium.gif) [Fig. 4.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F4) Fig. 4. Expanding the signal into a Fast Fourier series. a -1.alk\_s1, b -2.1alk\_s2, c - 3.not\_alk\_s1, d - 4.not\_alk\_s2 The received image dataset will be used as input data to the neural network. ## 3. Experimental research ### 3.1. Machine Learning The use of machine learning in classification tasks today is becoming less popular due to the development of deep neural networks. But in the next papers, results with high accuracy with machine learning were obtained [21,22,23,24]. For this reason, we tried using machine learning. In our research for machine learning tasks, the available dataset of 480 excel files was converted to a single file, of the following form, table. 3 View this table: [Table 3.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/T3) Table 3. Type of machine learning dataset Algorithms Logical Regression, Naive bayes, k-Nearest Neighbors, Support\_Vector\_Machines shows results of approximately 0.50 accuracy. Algorithms Random Forest Classifier - 0, 75 accuracies. The maximum result was obtained with the use of Classification and Regression Trees (CART) - 0.81 accuracies. ### 3.2. Deep neural networks To increase accuracy, we decided to use deep neural networks. We decided to use a convolutional neural network (CNN) that works with images. Figure 5 shows the two-dimensional (2D) graphics data correlation (Python3.7, matplotlib). ![Fig. 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F5.medium.gif) [Fig. 5.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F5) Fig. 5. Correlation EEG signals. Where, X-axis for the electrode from 1 to 64, Y-axis – value for electrode. Every electrode has 256 records with length 1 sec. a -1.alk\_s1, b -2.1alk\_s2, c - 3.not\_alk\_s1, d - 4.not\_alk\_s2 Image analysis shows a high correlation between regions that are close to each other. Visually we noticed that brain regions show different correlation values between subjects for the following regions PO3-CPZ and F4-C4. Visual observations allow us to conclude that the image data can be used for deep machine learning in the classification task. Today, the following CNN can be used in the EEG signal classification process [35,36,37,38], fig.6. ![Fig. 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F6.medium.gif) [Fig. 6.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F6) Fig. 6. Popular deep neural networks for image classification Initially, for classification tasks, we used the Inception network. The weights originally trained for layers 1 through 205 were used. Further, the network was completed on training images. The result of the classification accuracy was 0.72. The low accuracy is since the scales were initially trained to classify full-fledged images and not graph. Similar results were obtained for networks Resnet, VGG. Therefore, for classification, the new CNN was created. CNN consists of 6 layers of Convolution_2D and MaxPooling2D layers after the second and fourth convolution. On all layers except the output fully connected layer, the ReLU activation function is used, the last layer uses softmax. To regularize our model, after each subsample layer and the first fully connected layer, the Dropout layer was used. The following graphs were used as input images: * - source graphics without conversion; * - correlation graphs; * - Beta charts; * - Alpha charts; * - Theta graphics; * - Delta charts; * - Fast Fourier Series. Correlation graphs showed the highest accuracy. A schematic network when working with correlation graphs is shown in fig. 7. ![Fig. 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/04/2021.06.02.21258251/F7.medium.gif) [Fig. 7.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/F7) Fig. 7. Schematic representation of a neural network when working with correlation graphs The results of the classification accuracy in % when using various input images are presented in table 4. View this table: [Table 4.](http://medrxiv.org/content/early/2021/06/04/2021.06.02.21258251/T4) Table 4. Accuracy of CNN model for different dataset For the developed CNN model, the higher accuracy for classification objects between alcohol and not alcohol objects when we use Correlation EEG images was received. ## 4. Discussion and conclusions The highest accuracy result in classification was obtained using an image with the correlation of signals. The next areas PO3-CPZ and F4-C4 have the highest correlations. From the frequency range, a high classification result when working with the Beta range was obtained. The delta range showed the lowest result, which is associated with the loss of the useful signal in the original frequency. Using images with the Fourier series has accuracy commensurate with accuracy for the machine learning algorithm. In this research, there is no direct pattern between the magnitude of the voltage across the electrode and the group of studies. In 72% of cases, it is observed that in the group of alcoholics the voltage on the electrode is lower than in the control group. But, the location of the electrodes is different and there is no way to establish an exact relationship between the magnitude of the voltage on the electrode and the presence of the disease. In many papers, the beta, alpha, and theta rhythms of the EEG signal were used for classification signals. In this research, it is shown that with the same type of data about the object under research (only electrode voltage), it is preferable to use convolution networks with images of the correlation of EEG signals. For correct research in the field of analysis of EEG data using neural networks, it is necessary to submit as much data as possible from studies: age, gender, medical history, etc. Much research in the field of alcohol recognition by EEG signals has different results. In order to avoid it, it is necessary to develop a standard in the field of using neural networks regulating the number of signs for neural networks for classifying an alcoholism disease in a subject. ## Data Availability approval [https://github.com/Ildaron/3.eeg\_recognation](https://github.com/Ildaron/3.eeg_recognation) ## Dataset and code [https://github.com/Ildaron/3.eeg\_recognation](https://github.com/Ildaron/3.eeg_recognation) ## Conflicts of Interest None ## Funding None ## Ethical Approval Not required ## Footnotes * ildar.o2010{at}yandex.ru * [https://github.com/Ildaron/3.eeg\_recognation](https://github.com/Ildaron/3.eeg_recognation) * Received June 2, 2021. * Revision received June 2, 2021. * Accepted June 4, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory The copyright holder has placed this preprint in the Public Domain. It is no longer restricted by copyright. Anyone can legally share, reuse, remix, or adapt this material for any purpose without crediting the original authors. ## Reference 1. 1.Anuragi, A., Sisodia, D. Empirical wavelet transform based automated alcoholism detecting using EEG signal features. Biomedical Signal Processing and Control Volume 57, March 2020, 101777 2. 2.Oon-arom, A., Likhitsathain, S. and Srisurapanont, M. Efficacy and acceptability of varenicline for alcoholism: A systematic review and meta-analysis of randomized-controlled trials. Drug and Alcohol Dependence Volume 205, 1 December 2019, 107631 3. 3.Ishiguro, H., Higuchi, S., Arinami, T. and Onaivi, E. Association between alcoholism and the gene encoding the endocannabinoid synthesizing enzyme diacylglycerol lipase alpha in the Japanese population. Alcohol Volume 68, May 2018, pp. 59–62 4. 4.Kumar, R., Kumar, K. and Benegal, V. Cognitive and behavioural dispositions in offspring at high risk for alcoholism. Asian Journal of Psychiatry Volume 35, June 2018, pp. 38–44 5. 5.Winterer. GKloppel, B. and Heinza, A. Quantitative EEG (QEEG) predicts relapse in patients with chronic alcoholism and points to a frontally pronounced cerebral disturbance. Psychiatry Research Volume 78, Issues 1–2, 20 March 1998, pp. 101–113 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0165-1781(97)00148-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9579706&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F04%2F2021.06.02.21258251.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072882000010&link_type=ISI) 6. 6.Patidar, S., Pachori, R. and Upadhyay, A. An integrated alcoholic index using tunable-Q wavelet transform based features extracted from EEG signals for diagnosis of alcoholism. Applied Soft Computing Volume 50, January 2017, pp. 71–78 7. 7.Acharya, R., Vidya. S. and Bhat, S. Computer-aided diagnosis of alcoholism-related EEG signals. Epilepsy & Behavior Volume 41, December 2014, pp. 257–263 8. 8.Jeremy Harper Stephen M. Malone William G. Iacono. Impact of alcohol use on EEG dynamics of response inhibition: a cotwin control analysis. addiction biology. 2016 [https://doi.org/10.1111/adb.12481](https://doi.org/10.1111/adb.12481) 9. 9.Ziya Ekşi. Akif Akgül. Mehmet Recep Bozkurt. The Classificaton of EEG Signals Recorded in Drunk and Non-Drunk People. International Journal of Computer Applications (0975 – 8887) Volume 68– No.10, April 2013 10. 10. Wajid Mumtaz. Pham Lam Vuong. Automatic diagnosis of alcohol use disorder using EEG features. Journal Knowledge-Based Systems archive Volume 105 Issue C, August 2016 Pages 48–59 11. 11. N. Ahmadi, Y. Pei, M. Pechenizkiy. Detection of alcoholism based on EEG signals and functional brain network features extraction. 2017 IEEE 30th International Symposium on Computer-Based Medical Systems (CBMS), June 2017, pp. 179–184 12. 12. Guohun Zhu, Yan Li, Peng (Paul) Wen, Shuaifang Wang. Analysis of alcoholic EEG signals based on horizontal visibility graph entropy. Brain Informatics (2014) 1:19–25 DOI 10.1007/s40708-014-0003-x, pp.19–25 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40708-014-0003-x&link_type=DOI) 13. 13. Mingyue Yan. An Exponential Power Ratio Index based Algorithm for Analysis of Alcoholic EEG Signal. ACADEMY PUBLISHER doi:10.4304/jsw.8.9., pp. 2385–2390 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4304/jsw.8.9&link_type=DOI) 14. 14.Paulchamy B and Ila Vennila. A PROFICIENT SYSTEM FOR PREVENTING AND ACKNOWLEDGING ABOUT THE DRUNKEN DRIVE BY ANALYSING THE NEURONAL - ACTIVITIY OF THE BRAIN. ARPN Journal of Engineering and Applied Sciences. VOL. 7, NO. 8, AUGUST 2012 15. 15.Joel Brooks. Discrimination of Alcoholics from Non-Alcoholics using Supervised Learning on Resting EEG. [https://dspace.mit.edu/bitstream/handle/1721.1/87943/880382385-MIT.pdf?sequence=2](https://dspace.mit.edu/bitstream/handle/1721.1/87943/880382385-MIT.pdf?sequence=2) 16. 16. Madhavi Rangaswamy, Bernice Porjesz, David B. Chorlian,Theta Power in the EEG of Alcoholics. 0145-6008/03/2704-0607$03.00/0 ALCOHOLISM: CLINICAL AND EXPERIMENTAL RESEARCH. Vol. 27, No. 4 April 2003 17. 17.Anuragi, A., Sisodia, D. Empirical wavelet transform based automated alcoholism detecting using EEG signal features. Biomedical Signal Processing and Control Volume 57, March 2020, 101777 18. 18.Siuly, S., Li, Y. Designing a robust feature extraction method based on optimum allocation and principal component ana. Computer Methods and Programs in Biomedicine Volume 119, Issue 1, April 2015, pp. 29–42 19. 19.Artoni, F., Delorme, A. and Makeig, S. Applying dimension reduction to EEG data by Principal Component Analysis reduces the quality of its subsequent Independent Component decomposition. NeuroImage Volume 175, 15 July 2018, pp. 176–187 20. 20.Dea, F. and Stecca, M. A Big-Data-Analytics Framework for Supporting Classification of ADHD and Healthy Children via Principal Component Analysis of EEG Sleep Spindles Power Spectra. Procedia Computer Science Volume 159, 2019, pp. 1584–1590 21. 21.Raghu, S., Sriraam, N. Classification of focal and non-focal EEG signals using neighborhood component analysis and machine learning algorithms. Expert Systems with Applications Volume 113, 15 December 2018, pp. 18–32 22. 22.Zhang, Y., Wang, Y., Zhou, G. and Jin, J. Multi-kernel extreme learning machine for EEG classification in brain-computer interfaces. Expert Systems with Applications Volume 96, 15 April 2018, pp. 302–310 23. 23.Fernández-Varela, I., Hernández-Pereira, E. and lvarez-Estévez, A. Combining machine learning models for the automatic detection of EEG arousals. Neurocomputing Volume 268, 13 December 2017, pp. 100–108 24. 24.Hasanzadeh. F., Mohebbi, M. Prediction of rTMS treatment response in major depressive disorder using machine learning techniques and nonlinear features of EEG signal. Journal of Affective Disorders Volume 256, 1 September 2019, pp. 132–142 25. 25.Khan, H., Shanir, P. and Farooq, O. A hybrid Local Binary Pattern and wavelets based approach for EEG classification for diagnosing epilepsy. Expert Systems with Applications Volume 140, February 2020, 112895 26. 26.Balasubramanian, G., Kanagasabai, A. and Mohan, J. Music induced emotion using wavelet packet decomposition—An EEG study. Biomedical Signal Processing and Control Volume 42, April 2018, pp. 115–128 27. 27.Khalaf, A., Sejdic, E. and Akcakaya, M. Common spatial pattern and wavelet decomposition for motor imagery EEG-fTCD brain-computer interface. Journal of Neuroscience Methods Volume 320, 15 May 2019, pp. 98–106 28. 28.Kaleem, M., Guergachi, A. and Krishnan, S. Patient-specific seizure detection in long-term EEG using wavelet decomposition. Biomedical Signal Processing and Control Volume 46, September 2018, pp. 157–165 29. 29.Bajaj, N., Carrión, J. and Berta, R. Automatic and tunable algorithm for EEG artifact removal using wavelet decomposition with applications in predictive modeling during auditory tasks. Biomedical Signal Processing and Control Volume 55, January 2020, 101624 30. 30.Al-Salman, W., Li, Y. and Wen, P. Detecting sleep spindles in EEGs using wavelet fourier analysis and statistical features. Biomedical Signal Processing and Control Volume 48, February 2019, pp. 80–92 31. 31.Hyvärinen, .A, Ramkumar, P., Parkkonen, L. and Hari, R. Independent component analysis of short-time Fourier transforms for spontaneous EEG/MEG analysis. NeuroImage Volume 49, Issue 1, 1 January 2010, pp. 257–271 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2009.08.028&link_type=DOI) 32. 32.Sitnikova, E., Hramov, A., Koronovsky, A. and Luijtelaar, G. Sleep spindles and spike– wave discharges in EEG: Their generic features, similarities and distinctions disclosed with Fourier transform and continuous wavelet analysis. Journal of Neuroscience Methods Volume 180, Issue 2, 15 June 2009, pp. 304–316 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jneumeth.2009.04.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19383511&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F04%2F2021.06.02.21258251.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266953200012&link_type=ISI) 33. 33.Chen, G. Automatic EEG seizure detection using dual-tree complex wavelet-Fourier features. Expert Systems with Applications Volume 41, Issue 5, April 2014, pp. 2391–2394 34. 34.Kato, K., Takahashi, K., Nobuaki, M. and Ushiba, J. Online detection of amplitude modulation of motor-related EEG desynchronization using a lock-in amplifier: Comparison with a fast Fourier transform, a continuous wavelet transform, and an autoregressive algorithm. Journal of Neuroscience Methods Volume 293, 1 January 2018, pp. 289–298 35. 35.Jiao, Z., Gao, X., Wang, Y. and Xu, H. Deep Convolutional Neural Networks for mental load classification based on EEG data. Pattern Recognition Volume 76, April 2018, pp. 582–595 36. 36.Qiao, W., Bi, X. Ternary-task convolutional bidirectional neural turing machine for assessment of EEG-based cognitive workload. Biomedical Signal Processing and Control Volume 57, March 2020, 101745 37. 37.Hua, C., Hong, W. and Chen, J. Novel functional brain network methods based on CNN with an application in proficiency evaluation. Neurocomputing Volume 359, 24 September 2019, pp. 153–162 38. 38.Acharya, U., Oh, S. and Hagiwara, Y. Automated EEG-based screening of depression using deep convolutional neural network. Computer Methods and Programs in Biomedicine Volume 161, July 2018, pp. 103–113 [1]: /embed/graphic-7.gif [2]: /embed/graphic-8.gif [3]: /embed/graphic-9.gif [4]: /embed/graphic-11.gif