- 1Chair of Neuroinformatics and Biomedical Engineering, Faculty of Mathematics, Physics and Computer Science, Institute of Computer Science, Maria Curie-Sklodowska University, Lublin, Poland
- 2Neurophysiological Independent Unit of the Department of Psychiatry, Medical University of Lublin, Lublin, Poland
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, 2013, 2014; Teruel 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 et al., 2007; 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 and Kaminski, 2007; 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., 2008, 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., 1994, 1999; Pascual-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 software1, 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 averaging2 and source analysis and imaging3 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 EGI4. 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 environment5.
We investigated 30 patients, 9 females and 21 males (avg. age 28.1, s.d. 12.4). They have been diagnosed with a wide range of psychiatric disorders. The disorders are classified in ICD-10 as: 12 × F41 (Panic disorder), 5 × F32.1 (Major depressive episode), 5 × F84.5 (Asperger syndrome), 3 × F40 (Social anxiety disorders), 2 × F31 (Bipolar affective disorder), 2 × F42 (Obsessive-compulsive disorder co-occurrent with the patients with F84.5), 2 × F51.1 (Non-organic hypersomnia), and 1 × F20 (Schizophrenia). The control group of 30 healthy volunteers were also examined. The control group were only males (avg. age 22.4, s.d. 1.7). It is worth noting that about 30% more subjects were investigated from both control and patients' groups as the signal of all those for whom the recordings were too noisy or incomplete had to be eliminated. All participants were right-handed and measured by a handedness questionnaire (Chapman and Chapman, 1987).
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 reward-punishment 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.
Figure 1. Typical screens shown to participants during the IGT experiment. Card decks—in the left and reward/punishment with account state in the right. The screens are generated by the PST e-Prime 2.0.8.90 which is synchronized with Net Station 4.5.4 recordings.
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 decision-making 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 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).
Figure 2. Typical results of GeoSource BA activity visualization on the brain cortex so-called Flat Map. The increase of activity in BA36 (Uncus Lobe, Limbic Lobe) for good choice and BA37 (Fusiform Gyrus, Temporal Lobe) for bad choice are indicated.
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.
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.
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.
Table 1. Most active BA in particular subjects of the patients' group while receiving a reward during the IGT experiment in the full and in the alpha, beta, gamma, delta, and theta EEG bands.
Table 2. Most active BA in particular subjects of the patients' group while receiving a punishment during the IGT experiment in the full and in the alpha, beta, gamma, delta, and theta EEG bands.
Table 3. The ι for the most active BA in particular patients receiving a reward during the IGT experiment obtained from the sLORETA quantitative analysis.
Table 4. The ι for the most active BA in particular patients receiving a punishment during the IGT experiment obtained from the sLORETA quantitative analysis.
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.
Table 5. The most active BA in particular subjects of the control group while receiving a reward during the IGT experiment in the full and in the alpha, beta, gamma, delta, and theta EEG bands.
Table 6. The most active BA in particular subjects of the control group while receiving a punishment during the IGT experiment in the full and alpha, beta, gamma, delta, and theta EEG bands.
Table 7. The ι for the most active BA in particular subjects of the control group receiving a reward during the IGT experiment obtained from the sLORETA quantitative analysis.
Table 8. The ι for the most active BA in particular subjects of the control group receiving a punishment during the IGT experiment obtained from the sLORETA quantitative analysis.
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.
Figure 4. The comparison of the most frequently active BAs in the patients and the control groups in the “reward” variant of IGT response for (A)—full, (B)—alpha, (C)—beta, (D)—gamma, (E)—delta, (F)—theta frequency bands. The y-axis refers to the N—number of subjects with particular BA maximum activity noted.
Figure 5. The comparison of the most frequently active BAs in the patients and the control groups in the “punishment” variant of IGT response for (A)—full, (B)—alpha, (C)—beta, (D)—gamma, (E)—delta, (F)—theta frequency bands. The y-axis refers to the N—number of subjects with particular BA maximum activity noted.
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:
• 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 decision-making 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 locally-deep 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.
Table 9. Comparison of the discussed classifiers efficiency for all frequency bands, including the broadband in the STD (punishment) and TGT (reward) classification tasks.
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 (Bolanowski and Paszkiewicz, 2015; 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). The protocol was approved by the Medical University of Lublin Bioethical Commission. All subjects gave written informed consent in accordance with the GCP. Permission No. KE-0254/138/2015 and No. KE-0254/140/2015 given by the Medical University of Lublin Bioethical Commission on May 28th, 2015.
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.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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.
Footnotes
1. ^GeoSource 3 electrical source imaging packages: https://www.egi.com/research-division/electrical-source-imaging/geosource
2. ^BESA: ERP analysis and averaging, http://www.besa.de/products/besa-research/features/erp-analysis-and-averaging/
3. ^BESA: Source analysis and imaging, http://www.besa.de/products/besa-research/features/source-analysis-and-imaging/
4. ^Electrical Geodesic Systems, Inc., 500 East 4th Ave. Suite 200, Eugene, OR, 97401, USA.
5. ^Psychology Software Tools, Inc. PST, Sharpsburg Business Park, 311 23rd Street Ext., Suite 200, Sharpsburg, PA, 15215-2821, USA.
References
Barbey, A. K., Koenigs, M., and Grafman, J. (2013). Dorsolateral prefrontal contributions to human working memory. Cortex 49, 1195–1205. doi: 10.1016/j.cortex.2012.05.022
Bechara, A., Damasio, A. R., Damasio, H., and Anderson, S. W. (1994). Insensitivity to future consequences following damage to human prefrontal cortex. Cognition 50, 7–15.
Bechara, A., Damasio, H., and Damasio, A. R. (2000). Emotion, decision making and the orbitofrontal cortex. Cereb. Cortex 10, 295–307. doi: 10.1093/cercor/10.3.295
Bechara, A., Damasio, H., and Damasio, A. R. (2003). Role of the amygdala in decision-making. Ann. N. Y. Acad. Sci. 985, 356–369. doi: 10.1111/j.1749-6632.2003.tb07094.x
Bechara, A., Damasio, H., Damasio, A. R., and Lee, G. P. (1999). Different contributions of the human amygdala and ventromedial prefrontal cortex to decision-making. J. Neurosci. 19, 5473–5481.
Bigdely-Shamlo, N., Makeig, S., and Robbins, K. A. (2016). Preparing laboratory and real-world EEG data for large-scale analysis: a containerized approach. Front. Neuroinform. 10:7. doi: 10.3389/fninf.2016.00007
Bolanowski, M., and Paszkiewicz, A. (2015). “The use of statistical signatures to detect anomalies in computer network,” in L. Gołȩbiowski and M. Mazur (eds.) Analysis and Simulation of Electrical and Computer Systems (Springer International), 251–260.
Brevers, D., Bechara, A., Cleeremans, A., and Noël, X. (2013). Iowa gambling task (IGT): twenty years after–gambling disorder and igt. Front. Psychol. 4:665. doi: 10.3389/fpsyg.2013.00665
Campanella, S. (2013). Why it is time to develop the use of cognitive event-related potentials in the treatment of psychiatric diseases. Neuropsychiatr. Dis. Treat. 9:1835. doi: 10.2147/NDT.S53687
Cavanagh, J. F., Napolitano, A., Wu, C., and Mueen, A. (2017). The patient repository for EEG data+ computational tools (PRED+CT). Front. Neuroinform. 11:67. doi: 10.3389/fninf.2017.00067
Cella, M., Dymond, S., and Cooper, A. (2010). Impaired flexible decision-making in major depressive disorder. J. Affect. Disord. 124, 207–210. doi: 10.1016/j.jad.2009.11.013
Chan, R. C., Shum, D., Toulopoulou, T., and Chen, E. Y. (2008). Assessment of executive functions: review of instruments and identification of critical issues. Archiv. Clin. Neuropsychol. 23, 201–216. doi: 10.1016/j.acn.2007.08.010
Chaouachi, M., Chalfoun, P., Jraidi, I., and Frasson, C. (2010). “Affect and mental engagement: towards adaptability for intelligent systems,” in Proceedings of the 23rd International FLAIRS Conference (Daytona Beach, FL: Citeseer).
Chapman, R. M., and Bragdon, H. R. (1964). Evoked responses to numerical and non-numerical visual stimuli while problem solving. Nature 203:1155.
Chow, Y., Masiak, J., Mikołajewska, E., Mikołajewski, D., Wójcik, G. M., Wallace, B., et al. (2018). Limbic brain structures and burnout–a systematic review. Adv. Med. Sci. 63, 192–198. doi: 10.1016/j.advms.2017.11.004
Cui, J.-F., Chen, Y.-h., Wang, Y., Shum, D. H., and Chan, R. C. (2013). Neural correlates of uncertain decision making: ERP evidence from the iowa gambling task. Front. Hum. Neurosci. 7:776. doi: 10.3389/fnhum.2013.00776
Donner, T. H., Siegel, M., Fries, P., and Engel, A. K. (2009). Buildup of choice-predictive activity in human motor cortex during perceptual decision making. Curr. Biol. 19, 1581–1585. doi: 10.1016/j.cub.2009.07.066
Elliott, R. (2003). Executive functions and their disorders: imaging in clinical neuroscience. Br. Med. Bull. 65, 49–59. doi: 10.1093/bmb/65.1.49
Fellows, L. K., and Farah, M. J. (2004). Different underlying impairments in decision-making following ventromedial and dorsolateral frontal lobe damage in humans. Cereb. Cortex 15, 58–63. doi: 10.1093/cercor/bhh108
Goldenholz, D. M., Ahlfors, S. P., Hämäläinen, M. S., Sharon, D., Ishitobi, M., Vaina, L. M., et al. (2009). Mapping the signal-to-noise-ratios of cortical sources in magnetoencephalography and electroencephalography. Hum. Brain Mapp. 30, 1077–1086. doi: 10.1002/hbm.20571
Golonka, K., Mojsa-Kaja, J., Popiel, K., Marek, T., and Gawlowska, M. (2017). Neurophysiological markers of emotion processing in burnout syndrome. Front. Psychol. 8:2155. doi: 10.3389/fpsyg.2017.02155
Grabowski, F., Paszkiewicz, A., and Bolanowski, M. (2015). “Wireless networks environment and complex networks,” in Analysis and Simulation of Electrical and Computer Systems (Springer International Publishing), 261–270.
Grech, R., Cassar, T., Muscat, J., Camilleri, K. P., Fabri, S. G., Zervakis, M., et al. (2008). Review on solving the inverse problem in EEG source analysis. J. Neuroeng. Rehabil. 5:25. doi: 10.1186/1743-0003-5-25
Gupta, R., Koscik, T. R., Bechara, A., and Tranel, D. (2011). The amygdala and decision-making. Neuropsychologia 49, 760–766. doi: 10.1016/j.neuropsychologia.2010.09.029
Haaland, V. Ø., and Landrø, N. I. (2007). Decision making as measured with the iowa gambling task in patients with borderline personality disorder. J. Int. Neuropsychol. Soc. 13, 699–703. doi: 10.1017/S1355617707070890
Heekeren, H. R., Marrett, S., Bandettini, P. A., and Ungerleider, L. G. (2004). A general mechanism for perceptual decision-making in the human brain. Nature 431:859. doi: 10.1038/nature02966
Heekeren, H. R., Marrett, S., Ruff, D. A., Bandettini, P., and Ungerleider, L. G. (2006). Involvement of human left dorsolateral prefrontal cortex in perceptual decision making is independent of response modality. Proc. Natl. Acad. Sci. U.S.A. 103, 10023–10028. doi: 10.1073/pnas.0603949103
Hsu, M., Bhatt, M., Adolphs, R., Tranel, D., and Camerer, C. F. (2005). Neural systems responding to degrees of uncertainty in human decision-making. Science 310, 1680–1683. doi: 10.1126/science.1115327
Huk, A. C., and Shadlen, M. N. (2005). Neural activity in macaque parietal cortex reflects temporal integration of visual motion signals during perceptual decision making. J. Neurosci. 25, 10420–10436. doi: 10.1523/JNEUROSCI.4684-04.2005
John, E. R., Prichep, L., Fridman, J., and Easton, P. (1988). Neurometrics: computer-assisted differential diagnosis of brain dysfunctions. Science 239, 162–169.
Kakiashvili, T., Koczkodaj, W. W., and Woodbury-Smith, M. (2012). Improving the medical scale predictability by the pairwise comparisons method: evidence from a clinical data study. Comput. Methods Prog. Biomed. 105, 210–216. doi: 10.1016/j.cmpb.2011.09.011
Kamarajan, C., and Porjesz, B. (2015). Advances in electrophysiological research. Alcohol Res. Curr. Rev. 37:53.
Kaminski, W. A., and Wojcik, G. M. (2004). Liquid state machine built of Hodgkin–Huxley neurons. Informatica 15, 39–44. Available online at: https://content.iospress.com/articles/informatica/inf15-1-03
Koczkodaj, W. W., and Szybowski, J. (2015). Pairwise comparisons simplified. Appl. Math. Comput. 253, 387–394. doi: 10.1016/j.amc.2014.12.069
Kotyra, S., and Wojcik, G. M. (2017a). “The station for neurofeedback phenomenon research,” in Polish Conference on Biocybernetics and Biomedical Engineering (Cham: Springer), 32–43.
Kotyra, S., and Wojcik, G. M. (2017b). “Steady state visually evoked potentials and their analysis with graphical and acoustic transformation,” in Polish Conference on Biocybernetics and Biomedical Engineering (Cham: Springer), 22–31.
Krishnaswamy, P., Obregon-Henao, G., Ahveninen, J., Khan, S., Babadi, B., Iglesias, J. E., et al. (2017). Sparsity enables estimation of both subcortical and cortical activity from MEG and EEG. Proc. Natl. Acad. Sci. U.S.A. 114, E10465–E10474. doi: 10.1073/pnas.1705414114
Kufel, D. S., and Wojcik, G. M. (2018). Analytical modelling of temperature effects on an ampa-type synapse. J. Comput. Neurosci. 44, 379–391. doi: 10.1007/s10827-018-0684-x
Lubar, J. F., Swartwood, M. O., Swartwood, J. N., and O'Donnell, P. H. (1995). Evaluation of the effectiveness of eeg neurofeedback training for adhd in a clinical setting as measured by changes in tova scores, behavioral ratings, and WISC-R performance. Biofeedback Selfregul. 20, 83–99. doi: 10.1007/BF01712768
Mann, C. A., Lubar, J. F., Zimmerman, A. W., Miller, C. A., and Muenchen, R. A. (1992). Quantitative analysis of eeg in boys with attention-deficit-hyperactivity disorder: controlled study with clinical implications. Pediatr. Neurol. 8, 30–36.
Mapelli, D., Di Rosa, E., Cavalletti, M., Schiff, S., and Tamburin, S. (2014). Decision and dopaminergic system: an ERPS study of Iowa gambling task in Parkinson's disease. Front. Psychol. 5:684. doi: 10.3389/fpsyg.2014.00684
Martínez-Rodrigo, A., Fernández-Sotos, A., Latorre, J. M., Moncho-Bogani, J., and Fernández-Caballero, A. (2017). Neural correlates of phrase rhythm: an EEG study of bipartite vs. rondo sonata form. Front. Neuroinform. 11:29. doi: 10.3389/fninf.2017.00029
Mikołajewska, E., and Mikołajewski, D. (2012). Neuroprostheses for increasing disabled patients' mobility and control. Adv. Clin. Exp. Med. 21, 263–272. Available online at: http://www.advances.umed.wroc.pl/pdf/2012/1410.pdf
Mikołajewska, E., and Mikołajewski, D. (2013). Ethical considerations in the use of brain-computer interfaces. Central Eur. J. Med. 8, 720–724. doi: 10.2478/s11536-013-0210-5
Mikołajewska, E., and Mikołajewski, D. (2014). The prospects of brain computer interface applications in children. Open Med. 9, 74–79. doi: 10.2478/s11536-013-0249-3
Monsell, S. (2003). Task switching. Trends Cogn. Sci. 7, 134–140. doi: 10.1016/S1364-6613(03)00028-7
Moretti, D. V., Babiloni, C., Binetti, G., Cassetta, E., Dal Forno, G., Ferreric, F., et al. (2004). Individual analysis of eeg frequency and band power in mild Alzheimer's disease. Clin. Neurophysiol. 115, 299–308. doi: 10.1016/S1388-2457(03)00345-6
Nidal, K., and Malik, A. S. (2014). EEG/ERP Analysis: Methods and Applications. Boca Raton, FL; London; New York, NY: CRC Press.
Niedermeyer, E., and da Silva, F. L. (2005). Electroencephalography: Basic Principles, Clinical Applications, and Related Fields. New York, NY: Lippincott Williams & Wilkins.
Ogiela, L., Tadeusiewicz, R., and Ogiela, M. R. (2008). Cognitive techniques in medical information systems. Comput. Biol. Med. 38, 501–507. doi: 10.1016/j.compbiomed.2008.01.017
Ozga, W. K., Zapała, D., Wierzgała, P., Augustynowicz, P., Porzak, R., and Wójcik, G. M. (2018). Acoustic neurofeedback increases beta ERD during mental rotation task. Appl. Psychophysiol. Biofeedback 44, 103–115. doi: 10.1007/s10484-018-9426-0
Pascual-Marqui, R. D. (2002). Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. Methods Find Exp. Clin. Pharmacol. 24, 5–12. Available online at: http://www.institutpsychoneuro.com/wp-content/uploads/2015/10/sLORETA2002.pdf
Pascual-Marqui, R. D., Lehmann, D., Koenig, T., Kochi, K., Merlo, M. C., Hell, D., et al. (1999). Low resolution brain electromagnetic tomography (LORETA) functional imaging in acute, neuroleptic-naive, first-episode, productive schizophrenia. Psychiatry Res. Neuroimaging 90, 169–179.
Pascual-Marqui, R. D., Michel, C. M., and Lehmann, D. (1994). Low resolution electromagnetic tomography: a new method for localizing electrical activity in the brain. Int. J. Psychophysiol. 18, 49–65.
Paulus, M. P. (2007). Decision-making dysfunctions in psychiatry altered homeostatic processing? Science 318, 602–606. doi: 10.1126/science.1142997
Pope, A. T., Bogart, E. H., and Bartolome, D. S. (1995). Biocybernetic system evaluates indices of operator engagement in automated task. Biol. Psychol. 40, 187–195.
Saletu, B., Anderer, P., and Saletu-Zyhlarz, G. (2010). EEG topography and tomography (LORETA) in diagnosis and pharmacotherapy of depression. Clin. EEG Neurosci. 41, 203–210. doi: 10.1177/155005941004100407
Sand, T., Bjørk, M., and Vaaler, A. (2013). Is EEG a useful test in adult psychiatry? Tidsskr. Norske Laegeforen. 133, 1200–1204. doi: 10.4045/tidsskr.12.1253
Sutton, S., Braren, M., Zubin, J., and John, E. (1965). Evoked-potential correlates of stimulus uncertainty. Science 150, 1187–1188.
Szaleniec, J., Wiatr, M., Szaleniec, M., SkłAdzień, J., Tomik, J., Oleś, K., et al. (2013). Artificial neural network modelling of the results of tympanoplasty in chronic suppurative otitis media patients. Comput. Biol. Med. 43, 16–22. doi: 10.1016/j.compbiomed.2012.10.003
Szaleniec, M., Tadeusiewicz, R., and Witko, M. (2008). How to select an optimal neural model of chemical reactivity? Neurocomputing 72, 241–256. doi: 10.1016/j.neucom.2008.01.003
Tamburin, S., Maier, A., Schiff, S., Lauriola, M. F., Di Rosa, E., Zanette, G., et al. (2014). Cognition and emotional decision-making in chronic low back pain: an erps study during iowa gambling task. Front. Psychol. 5:1350. doi: 10.3389/fpsyg.2014.01350
Teruel, M. A., Navarro, E., Romero, D., García, M., Fernández-Caballero, A., and González, P. (2017). “An innovative tool to create neurofeedback games for ADHD treatment,” in International Work-Conference on the Interplay Between Natural and Artificial Computation (Cham: Springer), 183–192.
Tohka, J., and Ruotsalainen, U. (2012). Imaging brain change across different time scales. Front. Neuroinform. 6:29. doi: 10.3389/fninf.2012.00029
Trivedi, J. (2006). Cognitive deficits in psychiatric disorders: current status. Indian J. Psychiatry 48:10. doi: 10.4103/0019-5545.31613
Ważny, M., and Wojcik, G. M. (2014). Shifting spatial attention numerical model of posner experiment. Neurocomputing 135, 139–144. doi: 10.1016/j.neucom.2013.12.043
Wierzgała, P., Zapala, D., Wójcik, G. M., and Masiak, J. (2018). Most popular signal processing methods in motor-imagery bci: a review and meta-analysis. Front. Neuroinform. 12:78. doi: 10.3389/fninf.2018.00078
Wojcik, G. M. (2012). Electrical parameters influence on the dynamics of the Hodgkin–Huxley liquid state machine. Neurocomputing 79, 68–74. doi: 10.1016/j.neucom.2011.10.007
Wojcik, G. M., and Garcia-Lazaro, J. A. (2010). Analysis of the neural hypercolumn in parallel pcsim simulations. Proc. Comput. Sci. 1, 845–854. doi: 10.1016/j.procs.2010.04.092
Wojcik, G. M., and Kaminski, W. A. (2007). Liquid state machine and its separation ability as function of electrical parameters of cell. Neurocomputing 70, 2593–2597. doi: 10.1016/j.neucom.2006.12.015
Wojcik, G. M., and Kaminski, W. A. (2008). “Self-organised criticality as a function of connections' number in the model of the rat somatosensory cortex,” in International Conference on Computational Science (Berlin; Heidelberg: Springer), 620–629.
Wojcik, G. M., Kaminski, W. A., and Matejanka, P. (2007). “Self-organised criticality in a model of the rat somatosensory cortex,” in International Conference on Parallel Computing Technologies (Berlin; Heidelberg: Springer), 468–476.
Wojcik, G. M., Masiak, J., Kawiak, A., Kwasniewicz, L., Schneider, P., Polak, N., et al. (2018b). Mapping the human brain in frequency band analysis of brain cortex electroencephalographic activity for selected psychiatric disorders. Front. Neuroinform. 12:73. doi: 10.3389/fninf.2018.00073
Wojcik, G. M., Masiak, J., Kawiak, A., Schneider, P., Kwasniewicz, L., Polak, N., et al. (2018a). New protocol for quantitative analysis of brain cortex electroencephalographic activity in patients with psychiatric disorders. Front. Neuroinform 12:27. doi: 10.3389/fninf.2018.00027
Wojcik, G. M., and Ważny, M. (2015). Bray-curtis metrics as measure of liquid state machine separation ability in function of connections density: [Procs 51C (2015) 2948–2951]. Proc. Comput. Sci. 51:2978. doi: 10.1016/j.procs.2015.07.446
Wu, L., Wu, L., Chen, Y., and Zhou, J. (2014). A promising method to distinguish vascular dementia from Alzheimer's disease with standardized low-resolution brain electromagnetic tomography and quantitative EEG. Clin. EEG Neurosci. 45, 152–157. doi: 10.1177/1550059413496779
Keywords: electroencephalography, sLORETA, psychiatric disorders, frequency band analysis, machine learning, Iowa Gambling Task, decision-making
Citation: Wojcik GM, Masiak J, Kawiak A, Kwasniewicz L, Schneider P, Postepski F and Gajos-Balinska A (2019) Analysis of Decision-Making Process Using Methods of Quantitative Electroencephalography and Machine Learning Tools. Front. Neuroinform. 13:73. doi: 10.3389/fninf.2019.00073
Received: 25 April 2019; Accepted: 14 November 2019;
Published: 27 November 2019.
Edited by:
Ludovico Minati, Tokyo Institute of Technology, JapanReviewed by:
Stavros I. Dimitriadis, Cardiff University, United KingdomRyszard Tadeusiewicz, AGH University of Science and Technology, Poland
Copyright © 2019 Wojcik, Masiak, Kawiak, Kwasniewicz, Schneider, Postepski and Gajos-Balinska. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Grzegorz M. Wojcik, gmwojcik@live.umcs.edu.pl