Analysis of Decision-Making Process Using Methods of Quantitative Electroencephalography and Machine Learning Tools

The electroencephalographic activity of particular brain areas during the decision making process is still little known. This paper presents results of experiments on the group of 30 patients with a wide range of psychiatric disorders and 41 members of the control group. All subjects were performing the Iowa Gambling Task that is often used for decision process investigations. The electroencephalographical activity of participants was recorded using the dense array amplifier. The most frequently active Brodmann Areas were estimated by means of the photogrammetry techniques and source localization algorithms. The analysis was conducted in the full frequency as well as in alpha, beta, gamma, delta, and theta bands. Next the mean electric charge flowing through each of the most frequently active areas and for each frequency band was calculated. The comparison of the results obtained for the subjects and the control groups is presented. The difference in activity of the selected Brodmann Areas can be observed in all variants of the task. The hyperactivity of amygdala is found in both the patients and the control group. It is noted that the somatosensory association cortex, dorsolateral prefrontal cortex, and primary visual cortex play an important role in the decision-making process as well. Some of our results confirm the previous findings in the fMRI experiments. In addition, the results of the electroencephalographic analysis in the broadband as well as in specific frequency bands were used as inputs to several machine learning classifiers built in Azure Machine Learning environment. Comparison of classifiers' efficiency is presented to some extent and finding the most effective classifier may be important for planning research strategy toward finding decision-making biomarkers in cortical activity for both healthy people and those suffering from psychiatric disorders.

The electroencephalographic activity of particular brain areas during the decision making process is still little known. This paper presents results of experiments on the group of 30 patients with a wide range of psychiatric disorders and 41 members of the control group. All subjects were performing the Iowa Gambling Task that is often used for decision process investigations. The electroencephalographical activity of participants was recorded using the dense array amplifier. The most frequently active Brodmann Areas were estimated by means of the photogrammetry techniques and source localization algorithms. The analysis was conducted in the full frequency as well as in alpha, beta, gamma, delta, and theta bands. Next the mean electric charge flowing through each of the most frequently active areas and for each frequency band was calculated. The comparison of the results obtained for the subjects and the control groups is presented. The difference in activity of the selected Brodmann Areas can be observed in all variants of the task. The hyperactivity of amygdala is found in both the patients and the control group. It is noted that the somatosensory association cortex, dorsolateral prefrontal cortex, and primary visual cortex play an important role in the decision-making process as well. Some of our results confirm the previous findings in the fMRI experiments. In addition, the results of the electroencephalographic analysis in the broadband as well as in specific frequency bands were used as inputs to several machine learning classifiers built in Azure Machine Learning environment. Comparison of classifiers' efficiency is presented to some extent and finding the most effective classifier may be important for planning research strategy toward finding decision-making biomarkers in cortical activity for both healthy people and those suffering from psychiatric disorders.

INTRODUCTION
Decision-making is an essential skill in everyday life but currently there is little systematic knowledge about how decision-making is affected in people with a diagnosis of psychiatric disorders. Decision-making is a process in which many cognitive functions are engaged. Probably that is why the IGT was often chosen as a task for investigating the behavior of the people with psychiatric disorders, however, there is relatively not much known about the cortical activity of individuals while making decisions in both healthy people and those with psychiatric disorders diagnosis. Some research has been done on the patients with major depressive disorder (Cella et al., 2010;Brevers et al., 2013). Similarly, the cohort of subjects with borderline personality disorder was investigated using IGT (Haaland and Landrø, 2007) as well as bipolar disorder (Paulus, 2007). IGT applications for a variety of research and different disorders are presented to some extent in a review by Brevers et al. (2013) and originally in Bechara (2007). With many applications in psychiatry, we decided to choose the IGT out of many other tasks for this stage of our research.
Quantitative electroencephalography is at its Renaissance stage in last decades (Sand et al., 2013) and has developed toward some forms of research in modern psychiatry (Kamarajan and Porjesz, 2015;Martínez-Rodrigo et al., 2017).
The rapid increase in the number of publications concerning Brain-Computer Interfaces (BCI) is observed (Mikołajewska and Mikołajewski, 2012, 2013Teruel et al., 2017;Ozga et al., 2018;Wierzgała et al., 2018) and the EEG activity can be recognized as one of possible solutions in BCI engineering (Kotyra and Wojcik, 2017a,b). In addition, any ideas for finding biomarkers of psychiatric disorders (Chapman and Bragdon, 1964;Sutton et al., 1965;Campanella, 2013;Golonka et al., 2017) are in demand as the interview is still the most often used tool in psychiatry to make the diagnosis.
The expansion of computational modeling techniques applied to neuroscience makes it possible to simulate selected parts of the brain tissues which we are familiar with Wojcik and Kaminski, 2008;Wojcik and Garcia-Lazaro, 2010) or even investigate the influence of electrophysiological parameters of single cells on the dynamics of the whole simulated system Wojcik, 2012). However, we are still very far from explaining complex phenomena like psychiatric disorders or syndromes e.g., burn-out (Chow et al., 2018). Higher cognitive functions are sometimes a source of inspiration in biomedical engineering (Kaminski and Wojcik, 2004;Ważny and Wojcik, 2014;Wojcik and Ważny, 2015;Kufel and Wojcik, 2018) and artificial intelligence (Ogiela et al., 2008;Szaleniec et al., 2008Szaleniec et al., , 2013 mixed with cognitive science methodology provides some explanation or leads to the construction of classification tools. Nevertheless, we are still in demand for verification theory in the experiment. There are different electroencephalographic methods that allow visualization of recorded activity on the brain model. One of them is the standardized low-resolution brain electromagnetic tomography algorithm (sLORETA) (Pascual-Marqui et al., 1994Pascual-Marqui, 2002). This method advantages come from the high temporal resolution of modern electroencephalographs (Tohka and Ruotsalainen, 2012) and makes possible to compute the subjects brain activity distributed in time and put it on brain topography with the tomography-like quality of detail. Applications of sLORETA were reported e.g., for the attention-deficit-hyperactivity disorder (ADHD) (Mann et al., 1992) and neurodegenerative diseases (Wu et al., 2014). The sLORETA can be also applied in the frequency band analysis (Moretti et al., 2004;Saletu et al., 2010).
Using EEG based source localization techniques for the measurement of subcortical activity can be controversial. We are aware of the fact that for example in Krishnaswamy et al. (2017) authors state that subcortical structures produce smaller scalp EEG signals. This happens because they are farther from the head surface than cortical structures. To make matters worse, subcortical neurons can have a closed-field geometry that further weakens the observed distant fields and subcortical structures are surrounded by the cortical mantle. So measurements of activity in deep brain structures can potentially be explained by a surrogate distribution of currents on the cortex. That is why it can be very difficult to measure subcortical activity when cortical activity is occurring at the same time (Krishnaswamy et al., 2017). However, there are various mathematical models (Grech et al., 2008) that allow us to make some estimation of such kind of activity. Our lab is equipped with the very sophisticated and developed for 25 years GeoSource software 1 , where such models are implemented and based on the results given by it, having access to the photogrammetry station which generates the head model with high accuracy, we are able to draw some conclusions that are some extrapolated indicators for subcortical areas increased activity. The GeoSource is not the only software with subcortical areas activity algorithmic detectors. We have done some comparative analysis with BESA and its: ERP analysis and averaging 2 and source analysis and imaging 3 packages getting the same quality of results.
The investigations of Event-Related Potentials are often chosen by experimental psychologists as well as clinicians and biomedical engineers. One of the best-recognized ERP experiments in which decision-making is investigated was proposed by Bechara et al. (1994). It is known as the Iowa Gambling Task (IGT) and is described in detail in the Materials and Methods section of this contribution.
IGT was used in many clinical experiments (Cui et al., 2013;Mapelli et al., 2014;Tamburin et al., 2014). In Tamburin et al. (2014) the patients with chronic low back pain were investigated and the authors tried to find correlations between the ERP responses and the cognitive measures taken on them. On the other hand, in Cui et al. (2013) the students were investigated during IGT and the amplitudes of P3 potential were observed and discussed. Similar research is reported (Mapelli et al., 2014) but in this case it was focused on the people with Parkinson's disease making decisions and after that their ERP potentials were analyzed. The research mentioned above is concentrated on the analysis of the shape of statistically averaged potential and there are no source localization procedures applied to the analysis.
The aim of the research presented herein was to apply the methodology proposed in Wojcik et al. (2018a) and Wojcik et al. (2018b) to the quantitative electroencephalographic analysis of cortical activity from the patients in different frequency bands as well as in the full spectrum of the EEG signal. We used source localization techniques and having measured the average amperage in time for particular Brodmann Areas (BA) the mean electric charge flowing through them during the experiment was conducted for each patient and member of the control group. For this contribution, the brain activity of a group of patients with selected psychiatric disorders was measured using dense array EEG. These results were compared with those obtained from the participants of the control group. Both groups performed IGT.
Additionally, the results gathered for both healthy and disordered people in the broad and particular frequency EEG bands were taken as inputs to seven different machine learning classifiers in order to distinguish two types of responses in IGT, basing only on BA activity. The efficiency of these classifiers was compared and is presented to some extent.

MATERIALS AND METHODS
The Department of Neuroinformatics is equipped with the dense array amplifier recording the cortical activity with up to 500 Hz frequency through 256 channels HydroCel GSN 130 Geodesic Sensor Nets provided by EGI 4 . In addition, in the EEG Laboratory the Geodesic Photogrammetry System (GPS) was used. Eleven cameras placed in the corners of GPS take a set of subject's photos and then it is possible to make a model of the particular subject brain based on its calculated size, proportion and shape. Next the software imposes all computed activity results on this model with a very good accuracy. The amplifier operates on the Net Station 4.5.4 software, GPS is under control of the Net Local 1.00.00 and GeoSource 2.0. The eye blinks and saccades elimination as well as gaze calibration are obtained owing to the application of dedicated eye-tracker operated by SmartEye 5.9.7. The Event-Related Potentials (ERP) experiments are conducted in the PST e-Prime 2.0.8.90 environment 5 .
The IGT was introduced by Bechara et al. (1994) and since then it has become one of the favorite tasks given to the subjects participating in a wide range of experimental psychology experiments. Originating from the research first carried out at the University of Iowa the IGT was intended to get hold of mechanisms of decision-making process during the rewardpunishment oriented card game. The aim of the task is to choose one card deck symbol out of four in each of 100 trials. The participants are told to earn as much of virtual money as possible starting with 0 dollars. In each set of four cards (or symbols) there is a couple of so-called good cards for which there is a reward and a couple of so-called bad cards for which there is a punishment. The participants do not know which card is good and which is bad but they can conclude it from the game behavior. However, at the beginning all cards seem to be good, but for two of them they make impression to be better as the reward for choosing them is remarkably higher than for choosing the others. After several choices of the better cards, the punishments for choosing the next are extremely high. On the other hand the punishment for choosing cards after the initial selection of those worse at the beginning is very low which finally gives the better financial results when compared to the other case. The typical screens shown on the computer on which our participants make decisions is shown in Figure 1.
The electroencephalographic signal was recorded. After the test the photo of each participant was taken using the GPS. Such a technique allows obtaining spatial resolutions comparable to 1.5 T MRI without the necessity of MRI brain scanning for each participant. Thus, the anatomical models of participants' brains are generated using the GeoSource software and GPS photogrammetry which allows us to estimate the activity of particular BAs with satisfactory precision. Note, that in our approach we do not use the default model of the head, which is also possible. We make use of GPS to achieve the best possible accuracy of source localization.
Such methodology allowed us to apply the source localization algorithm with a satisfactory accuracy and estimate the most active Brodmann Areas in each participant during the decisionmaking process. The Net Station software along with the GeoSource tool has implemented the most popular version of the sLORETA algorithm which is described in the chapter titled Brain Source Localization Using EEG Signals in Nidal and Malik (2014). It is based on standardization of the current density assumption. That means that both the variance of the noise in the signal and the biological variance in the actual signal are taken into account (Goldenholz et al., 2009;Nidal and Malik, 2014). Independent and uniform distribution of the biological signal variance across the brain cortex is taken into consideration and this results in a linear imaging localization technique having exact, zero-localization error (Goldenholz et al., 2009;Nidal and   Malik, 2014). For more details see the comparison of different types of LORETA in Nidal and Malik (2014).
The literature reports a few bands that cover typical rhythmical activity of the brain (Niedermeyer and da Silva, 2005) described as follows: δ-delta band (<4 Hz), θ -theta (4-7 Hz), α-alpha (8-15 Hz), β-beta (16-31 Hz), γ -gamma (more than 31 Hz), and sometimes µ-mu (8-12 Hz) bands. Sometimes the frequency ranges that define each band are slightly different. In Frontiers in Neuroinformatics | www.frontiersin.org FIGURE 3 | Diagram of the IGT research protocol used in this paper. All scripts used for the preprocessing data in Net Station and post-processing in GeoSource are listed. Participation of the subject in the experiment begins when the Sensor Net is put on and ends when it is taken off. All data is collected by the Mac Pro workstation which is the central computational unit of the lab. Statistical analysis, finding the most active BAs in full or each of α, β, γ , δ, and θ frequency bands as well as ι estimations can be conducted on other machines. our lab by default the frequency bands are set as follow: δdelta band (0.1-3 Hz), θ -theta (4-7 Hz), α-alpha (7-12 Hz), β-beta (12-30 Hz), γ -gamma (more than 32 Hz).
One of the most useful functions of the GeoSource software is the possibility of estimation of the amperage of the most active areas (Figure 2) varying in time using source localization. The most active BA is indicated by the GeoSource as the BA with the highest electric current flowing through it in time. The activity of a particular BA could last at its maximum value for a longer or shorter period and it could appear more than once during each epoch. The signal was divided into epochs, as usual in ERP, in this case, IGT experiments, then averaged giving amperage in the function of time. Based on the electrical current measured by the EEG amplifier, most active BAs precisely indicated by the photogrammetry station and having precisely estimated time intervals owing to the perfect EEG time resolution, one of many numerical methods for integration can be applied to calculate the mean electric charge ι with good precision (Wojcik et al., 2018a) by integrating the electrical current in time. The details of mean electric charge estimation are described in detail in Wojcik et al. (2018a). The sLORETA can be run for the full EEG frequency band above 0.1 Hz including the γ spectrum and for the selected frequency band analysis. Besides the full band there were taken into consideration each of the following: alpha, beta, gamma, delta, and theta. For each band including the full band, the varying in time value of amperage of particular BAs was obtained from the GeoSource. Having the amperage in the function of time one can calculate the mean electric charge ι flowing through the given BA as described in Wojcik et al. (2018a). The typical visualization of the GeoSource application to the signal is shown in the flat maps in Figure 2. The time interval in which the BA activity was calculated was set to 5 ms and there was chosen the 800 ms segmentation (each segment starting with the stimuli) for signal averaging.
The BA1, BA2, and BA3 were eliminated from our analysis as they are part of Primary Somatosensory Cortex (S1) which was hyperactive owing to the subject's fingertips contact with the response pad during the experiment.
The scheme of the methodology and research protocol are presented in Figure 3. The full band analysis protocol in the case of P300 experiments was presented in Wojcik et al. (2018a) and the frequency band analysis protocol was described in detail in DIGITS related paper in Wojcik et al. (2018b). For this contribution the mixture of both methods proposed in Wojcik et al. (2018a,b) is applied. The so-called Waveform Tools package from the Net Station ecosystem was used and all scripts shown in Figure 3 originate from it. The description of algorithms used in the preprocessing and post-processing stages of the research is given in Electrical Geodesics (2003) and the procedures used in the photogrammetry Net Local are described in the EGI Lab documentation (Electrical Geodesics, 2009, 2011. There were 100 trials for each participant, duration of the experiment was around 12 min. For the preprocessing we used the following and suggested by EGI engineers rules: As an average reference the average of all electrodes was taken. The artifact correction parameters were set as follow: bad channels filtering-Max-Min > 200 µV; eye-blinks-Max-Min > 140 µV, eye movements-Max-Min > 140 µV-all performed in moving average of 80 ms. Filtration settings were set to passband gain 99%, stopband gain 1% and rollof 2 Hz, The segmentation was performed from 100 ms before stimulus to 700 after stimulus with offset 13 ms. The baseline correction was set to 100 ms from portion of segment. The rejection of trials took place when there were more than 10 bad channels. The Holy Grail for the quantitative EEG based psychiatry is finding the biomarkers of particular psychiatric disorders based on the measured electrical activity of the brain. We proposed some idea to find such biomarkers in Wojcik et al. (2018b) by using the frequency band analysis and estimating the most active BAs in an above-mentioned way. Some research was also reported in Zolubak et al. (2019) where authors were investigating low frequency markers in neurofeedback therapy. But indicating the most active BAs can be not enough. In the last decades, we can, however, observe the rapid growth of data science methods applied to big datasets. One of the most important of them are machine learning tools and our idea was to check whether applying different classifiers to our results, both in broadband and specific frequency band analysis, could shed some light on solving diagnoses problems. If there is a secret code of particular disorders to be found in EEG activity-the application of machine learning tools, like classifiers, seems to be the best way to decode this. As the input to classifiers, the activity of BAs in the spectrum of the mean electric charge flowing through them should be considered. Because our patients' group consisted of only 30 subjects and with a wide range of disorders it was impossible to design classifiers that could distinguish the particular disorder from the another. However, the possibility of distinguishing the reward cortical states from the punishment activity was investigated and the efficiency of selected classifiers will be discussed in the following sections to some extent.

RESULTS
For each participant from both the patient and control groups, we estimated the most active BAs in each EEG frequency band during the reward and punishment phases of the IGT experiment. Thus the most active BAs for the reward variant in the patient group are presented in Table 1 and for the punishment in Table 2. The mean electric charge ι flowing through particular BA when receiving a reward by a patient is shown in Table 3 and for punishment in Table 4.
In analogy for the control group the reward associated most active BAs are presented in Table 5 and the punishment responses in Table 6. Tables 7, 8 present the values of ι calculated for each member of the control group in the reward and punishment variants of response, respectively.
As follows from Tables 1-8 in all cases the Amygdala is hyperactive and the order of the value of ι tends to be similar for reward and punishment in both the subject and control groups in each frequency band. However, the comparison shown in Figure 4 for the rewards and in Figure 5 for the punishment can shed some light on the main differences in cortical responses of people with psychiatric disorders and members of the control group.
As far as the frequency of BA appearances in the IGT reward cortical responses are concerned (Figure 4) one can note that: • The one of the amygdala is significantly higher in the patients compared to the control group when observed in the full, beta and delta frequency bands. • The one of BA17 is significantly higher in the control group than in the patients when observed in the gamma and theta bands. • The one of BA09 is higher in the control compared to the patients' group in the full and gamma frequency bands.
When considering the frequency of BA appearance in the IGT punishment responses one can observe that: For detail see the Discussion section in the text.
• The one of the amygdala is significantly higher in the patients than in the control group members in the full, alpha, beta, delta and gamma frequency bands. • The one of BA17 is significantly higher in the control than in the patients' group in the full and gamma bands. • The BA09 is significantly higher in the control group members than in patients when looking at the full and theta frequency bands. • The BA05 in the control group is higher than in the patients in the full and beta frequency bands while in the alpha band it is lower.
The role of the amygdala during the decision-making process was discussed before even by the authors of IGT (Bechara et al., 2003). It is known that both the amygdala and orbitofrontal cortex are parts of a neural circuit critical for judgement and decision-making being under influence of "primary inducers" defined as stimuli that unconditionally, or through learning (e.g., conditioning and semantic knowledge), can (perceptually or subliminally) produce states that are pleasurable or aversive (Bechara et al., 2003). In order to verify the somatic marker hypothesis which proposes that decision-making is a process depending on emotion, some research of the destroyed amygdala was carried out (Bechara et al., 1999;Gupta et al., 2011). During some fMRI studies it was shown that amygdala influences the decisionmaking process in the risk-taking experiments involving information ambiguity (Hsu et al., 2005).
Our experiments show that the people with psychiatric disorders have the amygdala more frequently hyperactive when compared to healthy participants from the control group.
The Azure Machine Learning Studio was used to construct seven different classifiers and next to compare their efficiency in the reward/punishment characteristic cortical activity detection and classification tasks. Our classifiers were designed in order to separate two classes (reward/punishment) and were as follow: logistic regression, decision jungle, support vector machine, boosted decision tree, averaged perceptron, Bayes point machine, classic neural network and locallydeep support vector algorithms. Each classifier had its own characteristics which are shown in Table 10. The registered activity of particular BAs manifesting itself by the mean electric charge ι in all discussed EEG frequency bands, including the broadband was taken as inputs to the classifier. Under ideal conditions, it would be expected to construct effective classifiers for particular psychiatric disorders, but having only 30 diagnosed subjects in our cohort with so many different diagnoses is far from being enough to perform any statistics. For the machine learning tasks, the control group was extended by an additional 11 healthy males finally reaching 41 subjects. Thus, the joined cohort consisted of 41 healthy males and 30 males with some disorders, in a total of 71 participants. Note that in typical ERP experiments there are standard (STD) and target (TGT) stimuli. In the case of IGT, the punishment is treated as STD and the reward is TGT as practically everyone dares to win. For each participant registered the reward and punishment states, finally giving 2 × 71 = 142 averaged responses to the investigated set.
The 5-fold Cross-Validation method was used to validate the efficiency of all classifiers.
Then the values of classification accuracy, recall and precision were calculated and results are presented in Table 9.
As one can see in Table 9 there is no ideal classifier that could be applied to all of the EEG frequency bands, including the broadband.
For the broadband, the best results were achieved by the Locally-Deep Support Vector (acc. 0.698) and the Average Perceptron Classifier (acc. 0.684) methods.
In the α band the Logistic Regression (acc. 0.690) and Neural Network (acc. 0.669) turned out to be the best classifiers  When one looks at classifiers' results in the β band he notes the Logistic Regression (acc. 0.697), Decision Jungle (acc. 0.690) and Support Vector Machine (0.690) as the best, however the Logistic Regression has the highest recall value of 0.676, while the highest precision of 0.714 is achieved by Support Vector Machine.
If one studies the activity in the γ band he finds the highest efficiency for the Support Vector Machine (acc. 0.662) and again Logistic Regression (acc. 0.634).
For the δ band the highest accuracy 0.704 was achieved by the Support Vector Machine and Average Perceptron Classifier.
In case of the θ band, the best three ones were Logistic Regression (acc. 0.704), Support Vector Machine (0.669), and Locally Deep Support Vector (0.669), the Logistic Regression with the highest precision 0.723.
Note that the Bayes Point Machine did not perform well in any of EEG frequency bands.

DISCUSSION
In our experiment the relations between the decision-making process and the emotional responses given by the soma under such experimental conditions are also visible. Somatosensory association cortex (SAC) is mentioned in some papers on decisions making where it is even stated that somatosensory pattern marks the scenario as good or bad (Bechara et al., 2000;Donner et al., 2009). Our results also show that BA05 is one of the few most frequently active BAs in the patients and the control groups members., Moreover, the activity is qualitatively different in different frequency bands.
As well the dorsolateral prefrontal cortex (BA09) is often reported as engaged in decision-making processes. It was even found that damage of this structure results in poor performance in IGT (Fellows and Farah, 2004) and the fMRI studies have shown that the dorsolateral prefrontal cortex plays a role of negotiator establishing the link among sensory evidence, decision, and action during the decision making (Heekeren et al., 2006). Hyperactive BA09 is also reported to be found in other cognitive processes (Elliott, 2003), like working memory (Barbey et al., 2013), cognitive flexibility (Monsell, 2003), and planning (Chan et al., 2008). In our experiments the BA09 seems to be much more active in the control group when compared to the patients.
The visual processing areas provide the sensory evidence for a decision (Heekeren et al., 2004) and our results confirmed that the primary visual cortex is one of the most engaged areas in such processes, much more active in the control than in the patients' group. Some experiments involve the visual motion detection in the decision-making process among macaques (Huk and Shadlen, 2005) and this can be an inspiration for our future research.
The research protocol has been proposed to record the electroencephalographical cortical activity of the human brain during the decision making process. We chose the IGT as one of the tasks that are most often used to investigate people making decisions. The sLORETA was then applied to find the most frequently active BA in the brain cortex both in the patients and the control group.
Some attempts to find biomarkers in the quantitative EEG signals were made for example by John et al. (1988). The frequency band analysis is often used in real-time computing of the engagement index (Lubar et al., 1995;Pope et al., 1995;Chaouachi et al., 2010). Moreover, some cognitive functions in patients with psychiatric disorders are different from those in healthy members of control groups (Trivedi, 2006).
Even though the cohorts were not large we could prove some findings reported after performing such experiments by means of much more sophisticated techniques including fMRI. We had 30 subjects with several different diagnoses. That is why it is hard to apply any more sophisticated statistical analysis. Collecting neurophysiological data is a real challenge for neuroinformatics (Bigdely-Shamlo et al., 2016;Cavanagh et al., 2017). In future it would require building separate cohorts for each particular disorder, for all genders and age ranges. Then it would be possible to make quantitative comparisons of cortical activity which hopefully could even lead to building psychiatric disorders classifiers.
The additional aim of this paper was to check whether it is possible to assess without looking into logs the subject's response in the IGT experiments using only the brain cortical electric activity as the input to the algorithm. The effectiveness of seven different tools from the Azure Machine Learning environment was investigated. The summary of the results is presented in Table 9.
It was shown that there is no universal classifier for each frequency band. However, for the future analysis the Logistic Regression in the α, β, and θ bands should be considered as well as the Support Vector Machine in the β, γ , and δ. Very interesting behavior can be observed for the Averaged Perceptron Classifier in the δ band which together with the Support Vector Machine has one the best recall and precision characteristics in the discussed research.
It is expected that for the larger dataset the efficiency would be much higher. This is the initial stage of our research but one can take it for granted that tuning-up the parameters would also improve the method performance. Now it is hard to predict which methods would be best for additional improving such classifiers. Probably the applications of sophisticated tools offered by applied mathematics (Kakiashvili et al., 2012;Koczkodaj and Szybowski, 2015) or even solutions found for engineering applications in computer science Grabowski et al., 2015) along with big data analysis in such case could add some value.
At this stage, we had access to a limited number of patients. In our methodology, we decided to choose only those who had not taken any psychotropic medicines before. The aim of the research presented in this paper was to show the way in which the biomarkers can be searched. Putting the representatives of several disorders into one group by some readers can be recognized as controversial. On the other hand, we did not intend to give the final answer and to satisfy the definition of a biomarker in the full range of its properties. This would require at least 30 cases for each disorder and if one takes into consideration males and females, different age ranges, handedness-we get the number of about 400 patients for one problem, not saying about the control group. Consideration only the one disorder based on several patients does not make much sense as it would be hard to do the serious statistical analysis. But our results show that there can be quantitative methods to start the hunt for psychiatric disorders biomarkers.
Remembering that the interview is still the most important tool used in current psychiatry we are aware of the fact that developing tools and methods able to support the psychiatrist in the process of diagnosing are in a great demand and would improve the comfort of patients' life in the future.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of Guidelines for Good Clinical Practice (GCP).

AUTHOR CONTRIBUTIONS
GW: project idea and coordination, experiment design, subjects' recruitment, interpretation of results. JM: project idea, experiment design, subjects' recruitment, psychiatric diagnosis, interpretation of results. AK: work in laboratory, cleaning signal, computations, statistical analysis. PS, FP, and LK: statistical analysis, writing scripts, work in laboratory, cleaning signal. AG-B: work in laboratory.

ACKNOWLEDGMENTS
Special thanks to Mr. Slawomir Kotyra, M.Sc. from the Institute of Computer Science, Maria Curie-Sklodowska University in Lublin for solving the problem of electromagnetic noise in our laboratory.