Application of machine learning and complex network measures to an EEG dataset from DMT experiments
===================================================================================================

* Caroline L. Alves
* Thaise G. L. de O. Toutain
* Joel Augusto Moura Porto
* Aruane M. Pineda
* Eduardo Pondé de Sena
* Francisco A. Rodrigues
* Christiane Thielemann
* Manuel Ciba

## Abstract

There is a growing interest in the medical use of psychedelic substances as preliminary studies using them for psychiatric disorders have shown positive results. In particularly, one of these substances is N,N-dimethyltryptamine (DMT) an agonist serotonergic psychedelic that can induce profound alterations in state of consciousness.

In this work, we propose a computational method based on machine learning as an exploratory tool to reveal DMT-induced changes in brain activity using EEG data and provide new insights into the mechanisms of action of this psychedelic substance. To answer these questions, we propose a two-class classification based on (A) the connectivity matrix or (B) complex network measures derived from it as input to a support vector machine We found that both approaches were able to automatically detect changes in the brain activity, with case (B) showing the highest AUC (89%), indicating that complex network measurements best capture the brain changes that occur due to DMT use. In a second step, we ranked the features that contributed most to this result. For case (A) we found that differences in the high alpha, low beta, and delta frequency band were most important to distinguish between the state before and after DMT inhalation, which is consistent with results described in the literature. Further, the connection between the temporal (TP8) and central cortex (C3) and between the precentral gyrus (FC5) and the lateral occipital cortex (T8) contributed most to the classification result. The connection between regions TP8 and C3 has been found in the literature associated with finger movements that might have occurred during DMT consumption. However, the connection between cortical regions FC5 and P8 has not been found in the literature and is presumably related to emotional, visual, sensory, perceptual, and mystical experiences of the volunteers during DMT consumption. For case (B) closeness centrality was the most important complex network measure. Moreover, we found larger communities and a longer average path length with the use of DMT and the opposite in its absence indicating that the balance between functional segregation and integration was disrupted. This findings supports the idea that cortical brain activity becomes more entropic under psychedelics.

Overall, a robust computational workflow has been developed here with an interpretability of how DMT (or other psychedelics) modify brain networks and insights into their mechanism of action. Finally, the same methodology applied here may be useful in interpreting EEG time series from patients who consumed other psychedelic drugs and can help obtain a detailed understanding of functional changes in the neural network of the brain as a result of drug administration.

## I. INTRODUCTION

N,N-dimethyltryptamine (DMT) is a substance endogenously produced in various mammals [1], including humans [2], and has serotonin agonist properties. Thus, it is able to bind to serotonin receptors, simulating the neurotransmitter [3]. In [4] it was for the first time suggested that DMT is produced by the pineal gland in stress situations such as birth and death. In [5, 6] it seems clear that it is produced in small quantities by this gland [7]. The substance was first synthesized in 1931 [8], while its psychoactive effects have been described for the first time many years later in 1956 by [9]. When administered externally in large quantities, DMT can cause altered states of consciousness [10], hallucinations [11–13] and spiritual experiences such as communication with ‘presences’ or ‘entities’, plus reflections on death [14]. An exogenous ingestion can be done by smoking or injecting. Its effect by oral ingestion depends on the inhibition of monoamine oxidase, an enzyme that degrades the alkaloid DMT in the liver and intestine [12]. This enzyme and DMT are also found in ayahuasca tea that has been used in the Amazon for a couple of hundred years, being part of the traditional medicine of the natives from this region [15].

Recently, the interest in the medical use of psychedelics has increased significantly. Only last year, in [16] it was identified about 100 psychedelic clinical trials currently being conducted worldwide. This shows an increase in the number of clinical trials, compared to 43 assisted psychedelic therapy clinical trials conducted since 1999. One example is the psychedelic 3,4-methylenedioxymethamphetamine (MDMA) which is already in phase 3 clinical trials for the treatment of post-traumatic stress disorder (PTSD) [17] and major depression with positive results [18]. Another notable example is the psychedelic psilocybin whose therapeutic use in the U.S. has come to be considered a revolutionary therapy for treatment-resistant depression and major depressive disorder [19]. These first promising results suggest that other psychedelic substances, such as lysergic acid diethylamide (LSD), ibogaine hydrochloride, salvia divinorum, 5-MeO-DMT, ayahuasca and DMT, which have been less studied so far, should be investigated in more detail [20].

Only a few studies on administration of micro- or low-dose DMT to non-human species (predominantly rats) have been published in the scientific literature [21]. In [22], a low dose of DMT was administered to rats resulting in changes in frequency and amplitude of spontaneous excitatory postsynaptic currents (EPSCs) in the prefrontal cortex (PFC) that lasted long even after the drug was removed from the body. In [23], it was described that chronic, intermittent, low doses of DMT produced an antidepressant effect and increased fear extinction learning in rats without affecting working memory or social interaction. For a high dose of DMT (10 mg/kg), an increase in the density of the dendritic spines in the prefrontal cortex was found in rodents and antidepressant and anxiolytic behavioral effects were observed [24]. In humans, a single dose of 0.1 mg/kg of DMT caused an apparent anxiolytic effect shown in [25]. Other studies using inhaled 5-MeO-DMT also observed complete mystical experiences in 75% of volunteers [26] and improvements in depression and anxiety, which were associated with greater intensity of mystical experiences, with spiritual and personal meaning of the experience, when using this substance [27].

Summarizing, there is some evidence that DMT may be effective for the treatment of depression and post-traumatic stress disorder. However, most studies have been conducted in animals and thus only have a reduced power. Therefore, more in-depth studies on DMT, its mechanisms in the brain, and its potential clinical effects in humans are needed, since there are few studies investigating the use of DMT in humans through EEG [28–30] and fMRI [31].

The application of mathematical methods of graph theory yielded interesting insights into the complex network structure of the human brain. It is known from the literature [32–37] that the topology of brain is a small world network. Networks of this type combine completely random structural characteristics and regular connection topologies. They also preserve a high degree of connectivity between local neighborhoods, while allowing all their nodes 1 which to be connected to surprisingly short paths [32]. Altering this topology is also associated with pathological states [38–41] and the use of substances such as psychedelics [42–44]. Notably, complex networks parameters have been used as biomarker for several diseases [45, 46]. The use of complex networks is widely used in EEG to characterize the functional networks of the Brain [47–50].

In this context, machine learning (ML) has been used for more accurate and automatic medical diagnosis [51–58]. Compared to traditional statistical techniques, this approach has the advantage of not relying on prior assumptions (such as adequate distribution, independence in observations, absence of multicollinearity, and interaction problems) and, moreover, are suited to automatically analyze and capture non-linear complex relationships in data[59, 60]. As brain data are characterized by high complexity and highly correlated brain regions, ML algorithms have been widely used as a important tool capable of detecting acute and permanent abnormalities in the brain [61–63]. On the other hand, ML shows a lack of interpretability and a black box nature that is an especially disadvantageous general limitation when it comes to understanding medical data [64, 65]. In the last years new techniques have emerged to help in the interpretation of machine learning results. Most notable is the SHapley Additive ExPlanations (SHAP)values method [66]. This metric enables the identification and prioritization of features and can be used with any machine learning algorithm [67–69].

The present work aims to investigate EEG data using ML as an exploratory tool to detect temporal changes in the brain functionality of participants after DMT consumption. The study raised the following research questions:

*   Can we automatically detect changes in the functional network structure induced by DMT using ML ?

*   Which new insights into the mechanisms of action of DMT can we draw when we use ML in combination with SHAP values?

To answer these questions, we propose a two-class classification based on (A) the connectivity matrix or (B) complex network measures derived from it as input to a support vector machine (SVM) [70]. SVM has been used with excellent results for the classification of complex network measures before [50, 71, 72]. Furthermore, this ML algorithm can handle problems where the sample size of the data is generally smaller in comparison to the dimensionality of its feature space and is therefore applicable to the study of brain disorders with neuroimaging [73], whose data have these characteristics, and also this is the case of the data in this work (the case A with connectivity matrix).

For a biological interpretation the DMT-induced changes, the SHAP values method was applied to identify the features that contributed most to the classification (feature ranking). A robust workflow has been developed usable for medical professionals who are interested in the interpretation of brain network modifications due to DMT (or other psychedelics).

## II. DATA

The data used for this study has been published in [74] and is public available in a raw format 2. Thirty-five healthy male and female subjects (7 women and 28 men), volunteered to inhale, using pipes, 40 mg of free DMT extracted from the root of Mimosa hostilis. It should be noted that all participants had previous experiences with ayahuasca. Recordings were made with 24 electrodes, following the EEG electrode positions in the standard 10 – 20 location system. These channels are: Fp1, Fp2, Fz, F7, F8, FC1, FC2, Cz, C3, C4, T7, T8, CPz, CP1, CP2, CP5, CP6, TP9, TP10, Pz, P3, P4, O1, and O2. The recordings on the subjects started 10 minutes before DMT inhalation, 5 min with eyes closed and 5 minutes with eyes open. After DMT use, subjects were recorded about 6 min (6 ± 1.4 min).

## III. METHODOLOGY

In an earlier work of the authors [75], ML in combination with complex network measures was successfully applied to EEG data recorded after ayahuasca consumption to detect changes in brain activity. For this purpose, different levels of data abstraction were used as input: (a) the raw EEG time series, (b) the correlation of the EEG time series, and (c) the complex network measures calculated from (b). Several ML algorithms were tested and the best performance was obtained with the SVM at the abstraction levels (b) and (c). Based on this result, we decided to use in the present work connectivity matrices (see subsection III B) and derived complex network measures (see subsection III C as input for a SVM).

More details are displayed in Figure 1 which summarizes the methodology workflow. In short, EEG time series were separated by filtering in eight frequency bands. In a next step, preprocessing of the EEG time series were performed to obtain the connectivity matrices for each frequency band (and the unfiltered signal), see Figure 1-A and B with details of this process described in subsection III A and III B. In a second step, complex networks measures are derived from the connectivity matrices as described in subsection III C, see Figure 1-C and both types of data sets were used as input to a SVM as described in subsection III D. For interpretation of the classification results the feature ranking algorithm SHapley Additive exPlanations (SHAP) is finally applied as described in subsection III E.

![FIG. 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F1.medium.gif)

[FIG. 1:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F1)

FIG. 1: 
Diagram showing the methodology used in the present work. In **(A) Data preprocessing**, described in the subsection III A, the EEG time series are filtered to remove artifacts (in the picture the time series of a subject at the time the DMT was used) and then separated into the frequency bands high alpha, low alpha, high beta, mid-range beta, low beta, gamma, delta and theta (as an example in the picture the topographic map for the frequencies high alpha and high beta for the same subject). For each band the correlation between the channels is calculated using Pearson’s correlation to obtain a 24×24 connectivity matrix. In **(B) Connectivity matrices**, described in the subsection III B, where the connectivity matrices are flattened into a vector that are put into the SVM in order to verify the most important connections with the use of DMT (in the figure the best performing model, using the high alpha, low beta and delta bands, found TP8 and C3 as the main connections). In **(C) complex network measures**, described in the subsection III C, where the connectivity matrices are analyzed as graphs (in the figure for the same subject, the graph for the frequencies high alpha and high beta, where the number of connections in each node varies according to the color bar) and from them are extracted measures of complex networks that are applied in the SVM and the best model found for the delta frequency found the closeness centrality as the main measure.

### A. Data preprocessing

First, a high-pass filter with a cut off frequency of 0.5 Hz was used to remove artifacts such as electrogalvanic signals and motion artifacts [76]. This type of filtering is widely used in the literature [77–80]. To remove eye artifacts, we employed an independent component analysis (ICA) approach in which EEG signals are decomposed to maximize independent components and those with eye activity are identified and eliminated [81]. An example of the ICA analysis for a subject using DMT can be seen in Figure 2. This analysis was done with a python package called MNE [82] using an algorithm based on maximum information (Infomax) perspective [83].

![FIG. 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F2.medium.gif)

[FIG. 2:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F2)

FIG. 2: 
Scalp topography maps generated by ICA analytics recorded from a subject shortly after inhalation of DMT. The EEG signal was decomposed into twenty principal components. The component with ocular activity (see ICA018 with activity in the frontal region near the eyes) is removed and an artifact-free EEG signal is reconstructed.

In the next step, EEG time series were separated by filtering in eight frequency bands: high alpha (10-12 Hz), low alpha (8-10 Hz), low beta waves (12–15 Hz), mid-range beta waves (15–20 Hz), high beta waves (18–40 Hz), gamma (30 - 44 Hz), delta (0.5 – 3 Hz) and theta (4 – 7 Hz).

### B. Connectivity matrices

Connectivity matrices were calculated by the well-known Pearson correlation. It is a widely used and successfully approved measure to capture the correlation of EEG electrodes [84–89]. The Pearson correlation was calculated for all electrode pairs and all frequency bands (including the unfiltered signal). The connectivity matrices serve as input for the following steps as described in subsections III B and III C. To input the data into the ML algorithm, the connectivity matrices were flattened into one vector. Then, all vectors were sequentially merged into a 2D matrix where each column represents a connection between two electrodes and each row represents a subject. Such 2D matrices were generated for each frequency band (and the unfiltered signal).

### C. Complex network measures

For each connectivity matrix, a complex network graph was generated to extract different measures. To input the data into the ML algorithm, the complex network measures were stored in a matrix, where each column represents a complex network measure and each row a subject. Such a 2D matrices were generated for all frequencies bands (and the unfiltered signal).

The following complex network measures were calculated: Assortativity [90, 91], average path length (APL) [92], betweenness centrality (BC) [93], closeness centrality (CC) [94], eigenvector centrality (EC) [95], diameter [96], hub score [97], average degree of nearest neighbors [98] (Knn), mean degree [99], second moment degree (SMD) [100], entropy degree (ED) [101], transitivity [102, 103], complexity, k-core [104, 105], eccentricity [106], density [107], and efficiency [108]. In addition, newly developed metrics (described in detail in [75]) reflecting the number of communities in a complex network are applied. Community detection (also called clustering graph) is one of the fundamental analyses of complex networks aiming to decompose the network in order to find densely connected structures, so-called communities [109–111]. However, the community detection measures need to be transformed into a single scalar value to include them in the matrix. To this aim, we perform the community detection algorithms to find the largest community, then calculate the average path length within this community and receive a single value as the result. The community detection algorithms used were: Fastgreedy community (FC) [112], infomap community (IC) [113], leading eigenvector community (LC) [114], label propagation community (LPC) [115], edge betweenness community (EBC) [116], spinglass (SPC) [117], and multilevel community (MC) [118]. To indicate our approach, we extended the given abbreviations by the letter “A” (for average path length) as follows: AFC, AIC, ALC, ALPC, AEBC, ASPC, and AMC

### D. Machine learning process

In order to classify these two levels of data abstraction, namely the connectivity matrix and the complex network measures, the matrices were sampled by separating them into training (train) and test sets, with 25% of the data composing the test set. Then, for a reliable model, a k-cross validation was used [119], with k = 10 (value widely used in the literature [120–124]). For the training process, the training sets were applied to the SVM. SVM is based on the search for a hyperplane that geometrically divides samples into two classes. Three important hyper-parameters of the SVM have been considered in this work:

*   Kernel function: also known as kernel trick, has the function of projecting the input vectors in higher dimensions, because by increasing the dimension of the problem, the probability of it becoming a linearly separable problem increases, which makes it easier to solve [125, 126].

*   Regularization parameter C: this is the penalty term of the optimization problem and is an added constant that creates flexible margins with respect to the optimal hyperplane found.

*   Gamma: defines how much influence a single training example has. When the gamma value is too small, the model is too restricted and fails to capture the complexity of the data.

To find the best parameters, these hyper-parameters were optimized with the grid search method, widely used in the literature [127–131]. The grid search combines in a comprehensive way all values of the parameters selected for the models using some metrics to evaluate the performance of these combinations, which in the present work was the area under ROC curve (AUC) (for explanation see below). Here, we used the following functions as values for the kernel: gaussian (rbf), polynomial (poly), sigmoid and linear. Optimized values for parameter C and gamma are displayed in Appendix A.

For evaluation, the standard performance metrics accuracy (Acc.) was used as described in [132–136]. As we have a two-class (negative and positive) classification problem, other metrics like Precision and Recall are considered, also common in the literature [137–140]. Precision (also called specificity) corresponds to the hit rate in the negative class (here corresponding no effect induced by DMT). Whereas Recall (also called sensitivity) measures how well a classifier can predict positive examples (hit rate in the positive class), here related with an effect of DMT. Another well-known measure, see [128, 141, 142], is the F1 score which is the harmonic mean of the recall and precision [143]. For visualization of these two latter measures, the receiver operating characteristic (ROC) curve is a common method as it displays the relation between the rate of true positives and false positives. The area below this curve, called area under ROC curve (AUC) has been widely used in classification problems [130, 132, 144, 145]. The value of the AUC varies from 0 to 1, where the value of one corresponds to a classification result free of errors. *AUC* = 0.5 indicates that the classifier is not able to distinguish the two classes equal to the random choice. Furthermore, we consider the micro average of ROC curve, which computes the AUC metric independently for each class (calculate AUC metric for healthy individuals, class zero, and separately calculate for unhealthy subjects, class one) and then the average is computed considering these classes equally. The macro average is also used in our evaluation, which does not consider both classes equally, but aggregates the contributions of the classes separately and then calculates the average.

### E. Feature Ranking

As described in I, most notable technique for interpreting ML results is the SHAP values method based on the Shapley value concept which has its origin in game theory [146, 147], where it aims to assign payoffs to players depending on their contribution to the total payoff in the game. In addition, those who cooperate in a coalition receive a certain profit from this cooperation [148]. Applying this approach to our ML problem, each feature corresponds to a player in a game and the prediction corresponds to the payoff. Thus Shapley’s values tell us how to distribute the payoff fairly among the features [149].

Here, we used this methodology to evaluate which complex network measures and which correlation between electrodes (brain regions) contributed most to the classification result allowing for a biological interpretation of the results obtained with our ML algorithms.

## IV. RESULTS

ML was applied for two different levels of data abstraction: (A) the correlation of EEG time series (connectivity matrix) and (B) the complex network measures calculated from (A). We found that both approaches were able to automatically detect acute changes in the brain activity induced by the inhalation of DMT. The highest classification performance was obtained for the complex network measures with an AUC of 89% (see Table I). The following subsections IV A and IV B describe the results in more detail.

View this table:
[TABLE I:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/T1)

TABLE I: 
Performances of the SVM classifier. The classification of complex network measures capture the changes in the brain due to DMT slightly better than the connectivity matrix. Best performance is highlighted in bold.

### A. Connectivity matrix

EEG data recorded from subjects before DMT inhalation (control with eyes closed) and those after inhalation of DMT were filtered and divided into eight frequency bands as described in III A. Detailed results for each frequency band are given in the appendix B. The best performance was achieved for the low beta frequency band (test sample performance with mean AUC of 0.78, mean precision of 0.78, mean F1 score of 0.78, mean recall of 0.78, and mean Acc. of 0.78) followed by the high alpha and delta frequency bands (test sample performance for the both frequency band was a mean AUC of 0.72, mean precision of 0.72, mean F1 score of 0.72, mean recall of 0.72, and mean Acc. of 0.72). Clearly better results were achieved by combining these frequency bands whose test sample performance a mean AUC of 0.82, mean precision of 0.82, mean F1 score of 0.82, mean recall of 0.82, and mean Acc. of 0.82. Furthermore, see appendix D for similarity of results obtained for each frequency.

In Figure 3, the confusion matrix (3-(a)), the learning curve (Figure 3-(b)), and the ROC curve (3-(c)) are displayed. The learning curve evaluates the predictability of the model by varying the size of the training set [69]. Results show that the entire database is required to achieve the highest validation accuracy.

![FIG. 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F3.medium.gif)

[FIG. 3:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F3)

FIG. 3: 
ML results using connectivity matrices. (a) Confusion matrix indicating a true negative rate of 80.8% (purple according to the color bar) and a true positive rate of 84.6% (blue according to the color bar). (b) Learning curve for the training accuracy (blue) and for test accuracy (green). (c) ROC curve with class 0 (control) and class 1 (after inhalation of DMT).

In order to reveal the importance of the connections between electrode pairs (brain connections) by considering the combination of the best EEG frequency bands (high alpha, low beta and delta), the SHAP values were calculated. The results are shown in Figure 4. Clearly, the most important connection was between electrodes TP8 (temporal and parietal region) and C3 (central region). In addition, the presentation of the data in Figure 4 shows that for the connection between TP8 and C3, low values of correlation (blue dots) were important for detecting the presence of DMT (positive SHAP values), and high values of correlation (red dots) were important for detecting the absence of DMT (negative SHAP values). The second most important connection was between the electrodes FC5 and P8. The corresponding brain regions are depicted in Figure 5.

![FIG. 4:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F4.medium.gif)

[FIG. 4:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F4)

FIG. 4: 
Feature importance ranking for the SVM classifier with electrode correlation (brain regions) ranked in descending order of importance by considering the combination of the best EEG frequency bands (high alpha, low beta and delta). The connection between the regions TP8 and C3 is the most important to classify the effect of DMT. (a) Feature ranking based on the average of absolute SHAP values over all subjects considering both classes (gray: control, cyan: after inhalation of DMT). (b) Same as (a), but additionally showing details of the impact of each feature on the model.

![FIG. 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F5.medium.gif)

[FIG. 5:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F5)

FIG. 5: 
Brain plot displaying these electrode pairs which contributed most to the classification result based on connectivity matrices. The brain plot was made using Braph tool [150], based on the coordinates in [151, 152].

### B. Complex network measures

We received the best performance considering complex network measures for the delta frequency band (test sample performance with mean AUC of 0.89, mean precision of 0.91, mean F1 score of 0.88, mean recall of 0.88 and mean Acc. of 0.89), see Table I. The precision measure is related to the positive class (with DMT). Thus, since the precision was higher than the recall, we conclude that the model slightly better detects the presence of DMT than its absence. Furthermore, see appendix D for similarity of results obtained for each frequency.

In Figure 6, the confusion matrix (Figure 6 (a)), the learning curve (Figure 6 (b)), and the ROC curve (Figure 6 (c)) are plotted. Again, the entire database is necessary in order to get the highest accuracy. All the other results can be found in the Appendix C

![FIG. 6:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F6.medium.gif)

[FIG. 6:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F6)

FIG. 6: 
ML results using complex network measures. (a) Confusion matrix indicating a true negative rate of 77.8% (purple according to the color bar) and a true positive rate of 100.0% (blue according to the color bar). (b) Learning curve for the training accuracy (blue) and for test accuracy (green). (c) ROC curve with class 0 (control) and class 1 (after inhalation of DMT).

Based on the SHAP values in Figure 7 it can be seen that the most important measure for the model was the CCy, followed by the ALPC measure and the APL. In addition, high values of the CC measure (pink dots) indicate its importance for the detection of the absence of DMT (negative SHAP values), see Figure 7 (b).

![FIG. 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F7.medium.gif)

[FIG. 7:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F7)

FIG. 7: 
Feature importance ranking for SVM classifier with features ranked in descending order of importance. The CC measure is the most important to classify the effect of DMT. (a) Feature ranking based on the average of absolute SHAP values over all subjects considering both classes (gray: control, cyan: after inhalation of DMT). (b) Same as (a) but additionally showing details of the impact of each feature on the model.

## V. DISCUSSION

In the previous sections, we presented a computational workflow including data preprocessing and ML algorithm revealing acute differences in the brain activity before and after the consumption of the psychedelic drug DMT. As a result we achieved a classification accuracy of at least 82%. We further showed that the classification accuracy based on complex network measures (89%) was higher than that based on the connectivity matrix alone (see Table I), suggesting that these measures are important to capture differences in brain activity.

Further, we searched for descriptive parameters related to changes in the functional network structure by ranking the importance of features which contributed to the classification result. The results are discussed in this section with the aim to get insights into the effects of DMT consumption on the brain in terms of EEG frequency band (subsection V A), connection of the most activated brain regions (subsection V B), and measures of complex networks (subsection V C).

### A. Frequency bands

With our workflow, we were able to identify these frequency bands which were mostly modified after the intake of DMT. We found that classification results received with the connectivity matrices and the complex network measures are strongly based on acute changes in the delta band. Thus, changes in the delta band were most robust for both input data types.

This observation corresponds to what was also found in the literature [28, 74]. Delta band activity is usually associated with states where there is no wakefulness, such as sleep [153, 154] and coma [155]. However, some studies such in [156] observed that the delta frequency is present even when there are behavioral responses, such as in propofol anesthesia, postoperative delirium, and in powerful psychedelic states. Moreover, delta band activity has also been detected in studies involving spiritual experiences [157–159] and meditation states [160, 161].

Although the increase of delta band activity points clearly to an altered state of consciousness after the inhalation of DMT, also other frequency bands were affected. We found for the connectivity matrices that in addition to delta, high alpha and low beta bands were important features. This finding is supported by [74] who describes that inhalation of DMT reduces the alpha band activity while increasing the delta and gamma band at the same time [162]. According to the authors, the increase in gamma is associated with subjective perceptions typical of mystical experiences. In our data, we observed no changes in the gamma band.

### B. Connection between brain regions

With the connectivity matrices, we found that classification results are strongly based on a decreasing correlation between the temporal/parietal (TP8) and the central brain (C3) region after DMT uptake. These brain areas correspond to occipitotemporal (Right BA37), primary somatosensory cortex and the motor cortex (Left BA01/02), via Brodmann’s map [163]. The temporal lobe is associated with perception and production of speech, hearing, memory and emotional processes, because it is connected to the amygdala and the limbic system [145]. The right temporal region, TP8 found here, is associated with the recognition of familiar faces, with participation from the frontal cortex [164]. Furthermore, in humans, TP8 region contributes to the global processing of visual information [165]. The connection between this two regions, TP8 and C3, is involved in visual and tactile perceptions, and finger movements [166, 167].

In addition, the correlation between the electrodes FC5 (frontal region) and P8 (parietal region) contributed significantly to the classification result. These regions correspond to Left BA6 and Right BA19 of Brodmann’s map, respectively. The frontal region is involved in cognitive processing, planning behavior, and has connections to the somatosensory cortex, motor, and auditory areas [168, 169] and limbic system, and is also involved in emotions. The placement of FC5 electrode encompasses the region of the precentral gyrus in the pre motor region, which is responsible for controlling voluntary motor movement of the body. This region also includes a portion of the supplementary motor cortex, responsible for planning the voluntary movement of the limbs [170]. The P8 region, on the other hand, is located in the lateral occipital cortex, responsible for integrating different types of information so that our interaction with the environment is efficient, forming representational spaces through perception, semantics, through perception, semantics, and motor functions [171].

Studies using other psychedelics, such as LSD, have found reduced functional connectivity in the anterior medial prefrontal cortex, and time-specific effects were correlated with different aspects of subjective experiences under the effect of psychedelics [172]. Psilocybin consumption, on the other hand, was related to decreased functional connectivity between the medial temporal lobe and high-level cortical regions. The changes found in the cortical regions reported above, are related to the visual, sensory, perceptual, and motor type experiences experienced by the volunteers during the during the use of these two psychedelics (LSD and Psilocybin). Correlating our findings with previous studies, the FC5 and P8 regions also found here are part of the cortical region and a possible inference is that they are related to the participants’ experience with DMT in [74] in which it was reported that 13 of 35 participants (equivalent to 37%) accessed a complete mystical experience 3. Although these two cortical regions, FC5 and P8, may be related to the participants’ sensory and visual experiences when using DMT, no interpretation of the connectivity between these two regions has been obtained, since there is no information in the literature.

### C. Complex network measures

Concerning the measures of complex networks, the most important was the CC. CC is a centrality measure which is defined as the inverse of the average length of the shortest path from one node to all other nodes in the network [173]. The idea is that important nodes participate in many shortest paths within a network and, therefore, play an important role in the flow of information in the brain [94]. ALPC was the second important measure, which is associated to the size of the largest community found by the label propagation community (LPC) detection algorithm. Increased values (compared to controls) of this metric are associated with the effect of DMT (see Figure 7 (b)) indicating communities with increased average path lengths after the use of DMT, in other words, larger communities. The third important metric was the APL which is the average of all shortest paths. The shortest path *d**ij* (also known as the geodesic path) between two nodes *i* and *j*, is defined as the shortest of all possible paths between these vertices. Increased values for APL were associated with the presence of DMT (Figure 7(b)).

In the brain of large vertebrates there are two contrasting concepts: functional segregation (or specialization) and integration (or distributed processes) [174]. Anatomical and functional segregation refers to the existence of specialized neurons and brain areas organized in modules [175] which correspond to communities where their members have high connectivity among themselves and few connections with members of other modules [176]. As opposed to segregation, neuron units don’t operate in isolation [175], there are regions of the brain (distributed system of the cerebral cortex) capable of combining specialized information, characterizing the concept of integration [173]. These regions have an executing function, benefiting from a high global efficiency of information transfer throughout the entire network [177]. The fact that we found larger communities and a longer average path with the use of DMT and the opposite in its absence, indicates a decrease of the brain integration, which might slow down the distribution of information. Larger brain communities were also found in [44] after the use of ayahuasca, a mixture containing DMT.

Furthermore, when looking at the transitivity, which is a measure of the propensity of nodes to be grouped together, and efficiency measure, which is a measure of how effective the exchange of information within a network is, both also presented in the rank of the most important measures for the model in Figure 7 (b), the presence of DMT decrease the values of these two measures. The transitivity is a measure of the efficiency of information transfer between all pairs of nodes in the graph [178] and a higher value of these measures indicates greater segregation [179]. On other hand, higher values of efficiency indicates greater integration of networks. Thus, we can infer that the integration and segregation decreased with the use of the DMT considering the delta frequency. A decrease in brain segregation has been found in studies using other psychedelics such as LSD [42, 172, 180]. Specifically in [180], the authors concluded that the use of LSD caused a decrease in the integration and segregation of brain networks, supporting the idea that cortical brain activity becomes more “entropic” under psychedelics [181]. However, as pointed out in [42], psychedelics not only render the brain more random, but with normal organization disruption, they also produce strong functional and topologically far-reaching connections not seen in the normal state. Thus, even though our results show that integration and segregation have been disrupted, further experiments should be conducted to verify if there have been new long-distance connections as shown in the literature.

## VI. CONCLUSION AND FUTURE WORK

In summary, our results demonstrated that the application of ML methods was able to automatically reveal changes in brain functional connectivity induced by DMT consumption considering a two-class classification based on (A) the connectivity matrix or (B) complex network measures. The workflow developed here was indeed powerful for detecting the brain changes caused due to the psychedelic substance, with case (B) showing the highest AUC (89%), indicating that complex network measurements best capture the brain changes that occur due to DMT use. In terms of frequency, the workflow employed here detected that the delta frequency was most associated with DMT. Although DMT induces an altered state of consciousness with the presence of delta, other frequencies were important for recognizing the pattern of brain activity with the use of this substance, such as high alpha and low beta, through the connectivity matrix. This may suggest that the combination between the brain frequencies may represent an important point to be investigated, to further define the altered state of consciousness induced by DMT.

Furthermore, by using the SHAP value it was possible to interpret the results of the ML algorithms with a biological interpretation associated with the use of DMT on EEG data. The most important connections found with the use of DMT were between the temporal (TP8) and central cortex (C3) regions, followed by the connection between the precentral gyrus (FC5) and the lateral occipital cortex (T8). The connection between regions TP8 and C3 has been found in the literature associated with finger movements that might have occurred during DMT consumption. However, the connection between cortical regions FC5 and P8 has not been found in the literature and is presumably related to emotional, visual, sensory, perceptual, and mystical experiences of the volunteers during DMT consumption.

Concerning the measures of complex networks similar to what was found with the use of ayahuasca in [75], the most important was the centrality measure CC. Also, the fact that we found larger communities and a longer average path with the use of DMT and the opposite in its absence, indicates that this balance between functional segregation and integration was disrupted. This suggests that the distribution of information is slower. This findings supports the idea that cortical brain activity becomes more entropic under psychedelics. However, from the literature, psychedelics don’t simply make the brain more random, but after the normal organization is disrupted, strong and topologically far-reaching functional connections emerge that are not present in the normal state. Therefore, we would like to investigate in the long term how psychedelics change the functional connectivity of the brain using our workflow. Overall, a robust computational workflow has been developed here with an interpretability of how DMT (or other psychedelics) modify brain networks and insights into their mechanism of action. Finally, the same methodology applied here may be useful in interpreting EEG time series from patients who consumed other psychedelic drugs and can help obtain a detailed understanding of functional changes in the neural network of the brain as a result of drug administration. Thus, in future work we intend to use this methodology on the psychedelic drug called ketamine [182].

## Data Availability

All data are available online at [https://zenodo.org/record/3992359](https://zenodo.org/record/3992359)

## VII. ACKNOWLEDGEMENTS

F.A.R. acknowledges CNPq (grant 309266/2019-0) and FAPESP (grant 19/23293-0) for the financial support given for this research. T.G.L.O.T acknowledges FAPESB (grant number 307/2020 – Cota 2020; BOL0202/2020) for the financial support given this research. A.M.P. acknowledges FAPESP (grant 2019/22277-0) for the financial support given this research.

## Appendix A:

### Grid search hyperparameter tuning

The Figure 8 contains the values used in the present work where for one of the models (considering all frequencies and comparation of dmt and open eyes control) the combination of hyperparameter values and the grid search was plotted in relation to the AUC metric.

![FIG. 8:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F8.medium.gif)

[FIG. 8:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F8)

FIG. 8: 
Figure containing the values of each hyperparameter of the SVM varied with the grid search. In (A) for the model considering all frequencies and comparing the subject DMT and without DMT with the eye closed, the two-dimensional plot with the x-axis being the values of the parameter C and the y-axis being the values of the kernel and gamma function. For each combination of values and hyperparameters, AUC performance was obtained (whose obtained values are illustrated in the color table). In (B), the three-dimensional plot of (A) in which each hyperparameter corresponds to an axis.

## Appendix B:

### Results comparing different band frequencies

View this table:
[Table2](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/T2)

## Appendix C:

### Results considering complex network measures and different frequencies band

View this table:
[Table3](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/T3)

## Appendix D:

### Similarity of results obtained for each frequency

The SHAP value calculated for each frequency band was also considered. For each band a vector of connection between electrode pairs and its respective SHAP value found by the model is generated. For each of these vectors the Euclidean distance between them is then calculated generating a distance matrix of these vectors. This aims to quantify how close resulting vectors are. The distance matrix is displayed in form of a cluster map, see Figure 9, where vectors with a distance less than 0.2 are connected hierarchically in a dendrogram indicating clusters. Here, the cluster is most prominent between the low beta and the delta frequency band vector.

![FIG. 9:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F9.medium.gif)

[FIG. 9:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F9)

FIG. 9: 
Cluster map with the Euclidean distance of each band SHAP value vector. The delta bands and low beta frequency are the closest frequencies, forming in the cluster map a connection.

The same was made for the complex network measure and, then, a cluster map with the Euclidean distance between vectors containing SHAP values for each complex network measure is generated (see Figure 10). All vectors, except low beta, are very close to each other. This proximity indicates that the results obtained were similar, with other words the connections between the electrodes and their respective SHAP value were similar for all frequency bands. Thus, it can be seen that the results of the SHAP value vectors of each frequency, with the exception of the low beta frequency, were very close, which means that they show a great similarity between them.

![FIG. 10:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/06/16/2022.06.14.22276410/F10.medium.gif)

[FIG. 10:](http://medrxiv.org/content/early/2022/06/16/2022.06.14.22276410/F10)

FIG. 10: 
Cluster map showing the Euclidean distance of the SHAP value vectors for different frequency bands. Frequency band vectors, except low beta, are very close to each other, indicating great similarity.

## Footnotes

*   1 The network nodes can be representations from of neurons (¡1*µm*, microscale) to brain regions (*≈* 10 cm, macro scale)

*   2 Avaiable on Zenodo. [https://doi.org/10.5281/zenodo.3992359](https://doi.org/10.5281/zenodo.3992359)

*   3 According to [74], there was a significant positive correlation between items of the following scales, used with the volunteers in this study: the 5D altered states of consciousness scale (5D-ASC), mystical experience scales such as the affective component of the NDE scale and MEQ-30, and the social interactions on the psychedelic experience (post-experience questionnaire). Among the experiences accessed by the volunteers during the use of DMT, the scores with the highest percentages (mean; SD) for the 5D-ASC scale were: elementary imagery (85.27%; 20.72%), blissful (61.77%; 25.9%), complex imagery (50.21%;20.4%), spiritual (49.61%; SD 29%) and disembodiment (47.58%; 31.3%). For the NDE scale, related to positive mood, it was 60.94% (26.39%) for affect experience. For MEQ-30 scale, was found awe, with 46.29% (14.57) and for post-experience questionnaire scale, setting (76.48%; 25.38) and social (61.79; 3.91%)

*   Received June 14, 2022.
*   Revision received June 14, 2022.
*   Accepted June 16, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  [1]. S. T. Christian,  R. Harrison,  E. Quayle,  J. Pagel, and  J. Monti, The in vitro identification of dimethyltryptamine (dmt) in mammalian brain and its characterization as a possible endogenous neuroregulatory agent, Biochemical medicine 18, 164 (1977).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0006-2944(77)90088-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20877&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1977DV98700005&link_type=ISI) 

2.  [2]. J. R. Smythies,  R. Morin, and  G. Brown, Identification of dimethyltryptamine and o-methylbufotenin in human cerebrospinal fluid by combined gas chromatography/mass spectrometry, Biological psychiatry 14, 549 (1979).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=289421&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GZ84400012&link_type=ISI) 

3.  [3]. R. L. Smith,  H. Canton,  R. J. Barrett, and  E. Sanders-Bush, Agonist properties of n, n-dimethyltryptamine at serotonin 5-ht2a and 5-ht2c receptors, Pharmacology Biochemistry and Behavior 61, 323 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0091-3057(98)00110-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9768567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

4.  [4]. R. Strassman, Dmt: the spirit molecule. rochester, vt at¡ [http://scroungehound.com/sellsheets](http://scroungehound.com/sellsheets), B0145 DMT. pdf (2001).
    
    

5.  [5]. S. A. Barker,  E. H. McIlhenny, and  R. Strassman, A critical review of reports of endogenous psychedelic n, n-dimethyltryptamines in humans: 1955–2010, Drug testing and analysis 4, 617 (2012).
    
    

6.  [6]. S. A. Barker,  J. Borjigin,  I. Lomnicka, and  R. Strassman, Lc/ms/ms analysis of the endogenous dimethyltryptamine hallucinogens, their precursors, and major metabolites in rat pineal gland microdialysate, Biomedical Chromatography 27, 1690 (2013).
    
    

7.  [7]. D. E. Nichols, N, n-dimethyltryptamine and the pineal gland: Separating fact from myth, Journal of Psychopharmacology 32, 30 (2018).
    
    

8.  [8]. R. H. Manske, A synthesis of the methyltryptamines and some derivatives, Canadian Journal of Research 5, 592 (1931).
    
    

9.  [9]. S. Szara, Dimethyltryptamin: Its metabolism in man; the relation of its psychotic effect to the serotonin metabolism, Experientia 12, 441 (1956).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF02157378&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=13384414&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1956WS96400021&link_type=ISI) 

10. [10]. S. A. Barker, N, n-dimethyltryptamine (dmt), an endogenous hallucinogen: Past, present, and future research to determine its role and function, Frontiers in neuroscience, 536 (2018).
    
    

11. [11]. A. Gaujac,  J. L. Ford,  N. M. Dempster,  J. B. de Andrade, and  S. D. Brandt, Investigations into the polymorphic properties of n, n-dimethyltryptamine by x-ray diffraction and differential scanning calorimetry, Microchemical Journal 110, 146 (2013).
    
    

12. [12]. J. Ott, Pharmahuasca: human pharmacology of oral dmt plus harmine, Journal of psychoactive drugs 31, 171 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02791072.1999.10471741&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10438001&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

13. [13]. M. M. Schartner and  C. Timmermann, Neural network models for dmt-induced visual hallucinations, Neuroscience of Consciousness 2020, iaa024 (2020).
    
    

14. [14]. C. Timmermann,  L. Roseman,  L. Williams,  D. Erritzoe,  C. Martial,  H. Cassol,  S. Laureys,  D. Nutt, and  R. Carhart-Harris, Dmt models the near-death experience, Frontiers in psychology, 1424 (2018).
    
    

15. [15]. D. B. de Arauújo, Evidence for the therapeutic effects of ayahuasca, Advances in psychedelic medicine. State-of-the-art therapeutic applications, 103 (2019).
    
    

16. [16]. D. Perkins,  J. Sarris,  S. Rossell,  Y. Bonomo,  D. Forbes,  C. Davey,  D. Hoyer,  C. Loo,  G. Murray,  S. Hood, et al., Medicinal psychedelics for mental health and addiction: Advancing research of an emerging paradigm, Australian & New Zealand Journal of Psychiatry 55, 1127 (2021).
    
    

17. [17]. J. M. Mitchell,  M. Bogenschutz,  A. Lilienstein,  C. Harrison,  S. Kleiman,  K. Parker-Guilbert,  M. Ot’alora G,  W. Garas,  C. Paleos,  I. Gorman, et al., Mdma-assisted therapy for severe ptsd: a randomized, double-blind, placebo-controlled phase 3 study, Nature Medicine 27, 1025 (2021).
    
    

18. [18]. T. Andrews and  K. Wright, The frontiers of new psychedelic therapies: A survey of sociological themes and issues, Sociology Compass 16, e12959 (2022).
    
    

19. [19]. D. E. Nichols, Psilocybin: from ancient magic to modern medicine, The Journal of antibiotics 73, 679 (2020).
    
    

20. [20]. A. N. Siegel,  S. Meshkat,  K. Benitah,  O. Lipsitz,  H. Gill,  L. M. Lui,  K. M. Teopiz,  R. S. McIntyre, and  J. D. Rosenblat, Registered clinical studies investigating psychedelic drugs for psychiatric disorders, Journal of Psychiatric Research 139, 71 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpsychires.2021.05.019&link_type=DOI) 

21. [21]. S. A. Barker, Administration of n, n-dimethyltryptamine (dmt) in psychedelic therapeutics and research and the study of endogenous dmt, Psychopharmacology, 1 (2022).
    
    

22. [22]. C. Ly,  A. C. Greb,  L. P. Cameron,  J. M. Wong,  E. V. Barragan,  P. C. Wilson,  K. F. Burbach,  S. S. Zarandi,  A. Sood,  M. R. Paddy, et al., Psychedelics promote structural and functional neural plasticity, Cell reports 23, 3170 (2018).
    
    

23. [23]. L. P. Cameron,  C. J. Benson,  B. C. DeFelice,  O. Fiehn, and  D. E. Olson, Chronic, intermittent microdoses of the psychedelic n, n-dimethyltryptamine (dmt) produce positive effects on mood and anxiety in rodents, ACS chemical neuroscience 10, 3261 (2019).
    
    

24. [24]. L. P. Cameron,  C. J. Benson,  L. E. Dunlap, and  D. E. Olson, Effects of n, n-dimethyltryptamine on rat behaviors relevant to anxiety and depression, ACS chemical neuroscience 9, 1582 (2018).
    
    

25. [25]. R. J. Strassman,  C. R. Qualls,  E. H. Uhlenhuth, and  R. Kellner, Dose-response study of n, n-dimethyltryptamine in humans: Ii. subjective effects and preliminary results of a new rating scale, Archives of general psychiatry 51, 98 (1994).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archpsyc.1994.03950020022002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8297217&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994MW11800002&link_type=ISI) 

26. [26]. J. Barsuglia,  A. K. Davis,  R. Palmer,  R. Lancelotta,  A.-M. Windham-Herman,  K. Peterson,  M. Polanco,  R. Grant, and  R. R. Griffiths, Intensity of mystical experiences occasioned by 5-meo-dmt and comparison with a prior psilocybin study, Frontiers in Psychology, 2459 (2018).
    
    

27. [27]. A. K. Davis,  S. So,  R. Lancelotta,  J. P. Barsuglia, and  R. R. Griffiths, 5-methoxy-n, n-dimethyltryptamine (5-meo-dmt) used in a naturalistic group setting is associated with unintended improvements in depression and anxiety, The American journal of drug and alcohol abuse 45, 161 (2019).
    
    

28. [28]. C. Timmermann,  L. Roseman,  M. Schartner,  R. Milliere,  L. T. Williams,  D. Erritzoe,  S. Muthukumaraswamy,  M. Ashton,  A. Bendrioua,  O. Kaur, et al., Neural correlates of the dmt experience assessed with multivariate eeg, Scientific reports 9, 1 (2019).
    
    

29. [29]. E. Tagliazucchi,  F. Zamberlan,  F. Cavanna,  L. De La Fuente,  C. Romero,  Y. S. Perl, and  C. Pallavicini, Baseline power of theta oscillations predicts mystical-type experiences induced by dmt in a natural setting, Frontiers in Psychiatry 12 (2021).
    
    

30. [30]. A. Alamia,  C. Timmermann,  D. J. Nutt,  R. VanRullen, and  R. L. Carhart-Harris, Dmt alters cortical travelling waves, Elife 9, e59784 (2020).
    
    

31. [31]. J. Daumann,  K. Heekeren,  A. Neukirch,  C. M. Thiel,  W. Möller-Hartmann, and  E. Gouzoulis-Mayfrank, Pharmacological modulation of the neural basis underlying inhibition of return (ior) in the human 5-ht2a agonist and nmda antagonist model of psychosis, Psychopharmacology 200, 573 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00213-008-1237-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18649072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

32. [32]. O. Sporns and  J. D. Zwi, The small world of the cerebral cortex, Neuroinformatics 2, 145 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1385/NI:2:2:145&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15319512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223582200003&link_type=ISI) 

33. [33]. K. R. Van Dijk,  T. Hedden,  A. Venkataraman,  K. C. Evans,  S. W. Lazar, and  R. L. Buckner, Intrinsic functional connectivity as a tool for human connectomics: theory, properties, and optimization, Journal of neurophysiology 103, 297 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.00783.2009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19889849&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273526000029&link_type=ISI) 

34. [34]. J. Scannell,  G. Burns,  C. Hilgetag,  M. O’Neil, and  M. P. Young, The connectional organization of the corticothalamic system of the cat, Cerebral Cortex 9, 277 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/9.3.277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10355908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000080255100009&link_type=ISI) 

35. [35]. C.-C. Hilgetag,  G. A. Burns,  M. A. O’Neill,  J. W. Scannell, and  M. P. Young, Anatomical connectivity defines the organization of clusters of cortical areas in the macaque and the cat, Philosophical Transactions of the Royal Society of London B: Biological Sciences 355, 91 (2000).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2000.0551&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10703046&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085153300007&link_type=ISI) 

36. [36]. A. Fornito,  A. Zalesky, and  E. Bullmore, Fundamentals of brain network analysis (Academic Press, 2016).
    
    

37. [37]. J. G. White,  E. Southgate,  J. N. Thomson, and  S. Brenner, The structure of the nervous system of the nematode caenorhabditis elegans: the mind of a worm, Philosophical Transactions of the Royal Society of London 314, 1 (1986).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.1986.0056&link_type=DOI) 

38. [38]. D. S. Bassett and  M. S. Gazzaniga, Understanding complexity in the human brain, Trends in cognitive sciences 15, 200 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2011.03.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21497128&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291135800006&link_type=ISI) 

39. [39]. A. M. Pineda,  F. M. Ramos,  L. E. Betting, and  A. S. Campanharo, Quantile graphs for eeg-based diagnosis of alzheimer’s disease, Plos one 15, e0231169 (2020).
    
    

40. [40]. O. Sporns, Graph theory methods: applications in brain networks, Dialogues in clinical neuroscience 20, 111 (2018).
    
    

41. [41]. D. S. Bassett,  P. Zurn, and  J. I. Gold, On the nature and use of models in network neuroscience, Nature Reviews Neuroscience 19, 566 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41583-018-0038-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30002509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

42. [42]. D. E. Nichols,  M. W. Johnson, and  C. D. Nichols, Psychedelics as medicines: an emerging new paradigm, Clinical Pharmacology & Therapeutics 101, 209 (2017).
    
    

43. [43]. M. Girn,  C. Mills,  L. Roseman,  R. L. Carhart-Harris, and  K. Christoff, Updating the dynamic framework of thought: creativity and psychedelics, Neuroimage 213, 116726 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2020.116726&link_type=DOI) 

44. [44]. A. Viol,  F. Palhano-Fontes,  H. Onias,  D. B. de Araujo, and  G. Viswanathan, Shannon entropy of brain functional complex networks under the influence of the psychedelic ayahuasca, Scientific reports 7, 1 (2017).
    
    

45. [45]. M. Hayashida and  T. Akutsu, Complex network-based approaches to biomarker discovery, Biomarkers in medicine 10, 621 (2016).
    
    

46. [46]. T. Fekete,  N. Zach,  L. R. Mujica-Parodi, and  M. R. Turner, Multiple kernel learning captures a systems-level functional connectivity biomarker signature in amyotrophic lateral sclerosis, PloS one 8, e85190 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0085190&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24391997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

47. [47]. E. Bullmore and  O. Sporns, Complex brain networks: graph theoretical analysis of structural and functional systems, Nature reviews neuroscience 10, 186 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn2575&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19190637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263556500012&link_type=ISI) 

48. [48]. R. Baravalle,  N. Guisande,  M. Granado,  O. A. Rosso, and  F. Montani, Characterization of visuomotor/imaginary movements in eeg: An information theory and complex network approach, Frontiers in Physics, 115 (2019).
    
    

49. [49]. S. Das and  S. D. Puthankattil, Complex network analysis of mci-ad eeg signals under cognitive and resting state, Brain Research 1735, 146743 (2020).
    
    

50. [50]. M. Diykh,  Y. Li, and  P. Wen, Classify epileptic eeg signals using weighted complex networks based community structure detection, Expert Systems with Applications 90, 87 (2017).
    
    

51. [51]. C.-W. Song,  H. Jung, and  K. Chung, Development of a medical big-data mining process using topic modeling, Cluster Computing 22, 1949 (2019).
    
    

52. [52]. M. Mozaffarinya,  A. R. Shahriyari,  M. K. Bahadori,  A. Ghazvini,  S. S. Athari, and  G. Vahedi, A data-mining algorithm to assess key factors in asthma diagnosis, Revue FranÇaise d’Allergologie (2019).
    
    

53. [53]. N. Ilyasova,  A. Kupriyanov,  R. Paringer, and  D. Kirsh, Particular use of big data in medical diagnostic tasks, Pattern Recognition and Image Analysis 28, 114 (2018).
    
    

54. [54]. J. G. Richens,  C. M. Lee, and  S. Johri, Improving the accuracy of medical diagnosis with causal machine learning, Nature communications 11, 1 (2020).
    
    

55. [55]. C. J. Lynch and  C. Liston, New machine-learning technologies for computer-aided diagnosis, Nature medicine 24, 1304 (2018).
    
    

56. [56]. R. Alizadehsani,  M. Roshanzamir,  M. Abdar,  A. Beykikhoshk,  A. Khosravi,  M. Panahiazar,  A. Koohestani,  F. Khozeimeh,  S. Nahavandi, and  N. Sarrafzadegan, A database for using machine learning and data mining techniques for coronary artery disease diagnosis, Scientific data 6, 1 (2019).
    
    

57. [57]. P. A. Keane and  E. J. Topol, With an eye to ai and autonomous diagnosis (2018).
    
    

58. [58]. C. Bhatt,  I. Kumar,  V. Vijayakumar,  K. U. Singh, and  A. Kumar, The state of the art of deep learning models in medical science and their challenges, Multimedia Systems 27, 599 (2021).
    
    

59. [59]. R. Li,  A. Shinde,  A. Liu,  S. Glaser,  Y. Lyou,  B. Yuh,  J. Wong, and  A. Amini, Machine learning–based interpretation and visualization of nonlinear interactions in prostate cancer survival, JCO Clinical Cancer Informatics 4, 637 (2020).
    
    

60. [60]. H. S. R. Rajula,  G. Verlato,  M. Manchia,  N. Antonucci, and  V. Fanos, Comparison of conventional statistical methods with machine learning in medicine: diagnosis, drug development, and treatment, Medicina 56, 455 (2020).
    
    

61. [61]. R. C. Fong,  W. J. Scheirer, and  D. D. Cox, Using human brain activity to guide machine learning, Scientific reports 8, 1 (2018).
    
    

62. [62]. P. A. Kragel and  K. S. LaBar, Decoding the nature of emotion in the brain, Trends in cognitive sciences 20, 444 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2016.03.011&link_type=DOI) 

63. [63]. A. Boutet,  R. Madhavan,  G. J. Elias,  S. E. Joel,  R. Gramer,  M. Ranjan,  V. Paramanandam,  D. Xu,  J. Germann,  A. Loh, et al., Predicting optimal deep brain stimulation parameters for parkinson’s disease using functional mri and machine learning, Nature communications 12, 1 (2021).
    
    

64. [64]. C. Rudin, Stop explaining black box machine learning models for high stakes decisions and use interpretable models instead, Nature Machine Intelligence 1, 206 (2019).
    
    

65. [65]. I. Ekanayake,  D. Meddage, and  U. Rathnayake, A novel approach to explain the black-box nature of machine learning in compressive strength predictions of concrete using shapley additive explanations (shap), Case Studies in Construction Materials, e01059 (2022).
    
    

66. [66]. S. M. Lundberg and  S.-I. Lee, A unified approach to interpreting model predictions, in Proceedings of the 31st international conference on neural information processing systems (2017) pp. 4768–4777.
    
    

67. [67]. D. Bowen and  L. Ungar, Generalized shap: Generating multiple types of explanations in machine learning, arXiv preprint arxiv:2006.07155 (2020).
    
    

68. [68]. R. Rodríguez-Pérez and  J. Bajorath, Interpretation of compound activity predictions from complex machine learning models using local approximations and shapley values, Journal of Medicinal Chemistry 63, 8761 (2019).
    
    

69. [69]. G. Spadon,  A. C. de Carvalho,  J. F. Rodrigues-Jr, and  L. G. Alves, Reconstructing commuters network using machine learning and urban indicators, Scientific reports 9, 1 (2019).
    
    

70. [70]. L. Bottou and  C.-J. Lin, Support vector machine solvers, Large scale kernel machines 3, 301 (2007).
    
    

71. [71]. M. Mazrooyisebdani,  V. A. Nair,  C. Garcia-Ramos,  R. Mohanty,  E. Meyerand,  B. Hermann,  V. Prabhakaran, and  R. Ahmed, Graph theory analysis of functional connectivity combined with machine learning approaches demonstrates widespread network differences and predicts clinical variables in temporal lobe epilepsy, Brain connectivity 10, 39 (2020).
    
    

72. [72]. S. Dey,  A. R. Rao, and  M. Shah, Attributed graph distance measure for automatic detection of attention deficit hyperactive disordered subjects, Frontiers in neural circuits 8, 64 (2014).
    
    

73. [73]. D. A. Pisner and  D. M. Schnyer, Support vector machine, in Machine learning (Elsevier, 2020) pp. 101–121.
    
    

74. [74]. C. Pallavicini,  F. Cavanna,  F. Zamberlan,  L. A. de la Fuente,  Y. Ilksoy,  Y. S. Perl,  M. Arias,  C. Romero,  R. Carhart-Harris,  C. Timmermann, et al., Neural and subjective effects of inhaled n, n-dimethyltryptamine in natural settings, Journal of Psychopharmacology 35, 406 (2021).
    
    

75. [75]. C. L. Alves,  R. G. Cury,  K. Roster,  A. M. Pineda,  F. A. Rodrigues,  C. Thielemann, and  M. Ciba, Application of machine learning and complex network measures to an eeg dataset from ayahuasca experiments, medRxiv (2022).
    
    

76. [76]. P. Kumar,  R. Saini,  P. P. Roy, and  D. P. Dogra, A bio-signal based framework to secure mobile devices, Journal of Network and Computer Applications 89, 62 (2017).
    
    

77. [77]. N. J. Stevenson,  K. Tapani,  L. Lauronen, and  S. Vanhatalo, A dataset of neonatal eeg recordings with seizure annotations, Scientific data 6, 1 (2019).
    
    

78. [78]. M. Louwerse and  S. Hutchinson, Neurological evidence linguistic processes precede perceptual simulation in conceptual processing, Frontiers in psychology 3, 385 (2012).
    
    

79. [79]. C. Daftari,  J. Shah, and  M. Shah, Detection of epileptic seizure disorder using eeg signals, in Artificial Intelligence-Based Brain-Computer Interface (Elsevier, 2022) pp. 163–188.
    
    

80. [80]. M. Klug and  K. Gramann, Identifying key factors for improving ica-based decomposition of eeg data in mobile and stationary experiments, European Journal of Neuroscience 54, 8406 (2021).
    
    

81. [81]. T.-P. Jung,  S. Makeig,  C. Humphries,  T.-W. Lee,  M. J. Mckeown,  V. Iragui, and  T. J. Sejnowski, Removing electroencephalographic artifacts by blind source separation, Psychophysiology 37, 163 (2000).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0048577200980259&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10731767&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000085709100003&link_type=ISI) 

82. [82]. A. Gramfort,  M. Luessi,  E. Larson,  D. A. Engemann,  D. Strohmeier,  C. Brodbeck,  R. Goj,  M. Jas,  T. Brooks,  L. Parkkonen, et al., Meg and eeg data analysis with mne-python, Frontiers in neuroscience 7, 267 (2013).
    
    

83. [83]. S.-i. Amari,  A. Cichocki, and  H. Yang, A new learning algorithm for blind signal separation, Advances in neural information processing systems 8 (1995).
    
    

84. [84]. G. M. Rojas,  C. Alvarez,  C. E. Montoya,  M. de la Iglesia-Vayá,  J. E. Cisternas, and  M. Gálvez, Study of resting-state functional connectivity networks using eeg electrodes position as seed, Frontiers in neuroscience 12, 235 (2018).
    
    

85. [85]. L. Wang,  W. Wang,  T. Yan,  J. Song,  W. Yang,  B. Wang,  R. Go,  Q. Huang, and  J. Wu, Beta-band functional connectivity influences audiovisual integration in older age: an eeg study, Frontiers in aging neuroscience 9, 239 (2017).
    
    

86. [86]. M. Jalili, Functional brain networks: does the choice of dependency estimator and binarization method matter?, Scientific reports 6, 1 (2016).
    
    

87. [87]. C. Han,  X. Sun,  Y. Yang,  Y. Che, and  Y. Qin, Brain complex network characteristic analysis of fatigue during simulated driving based on electroencephalogram signals, Entropy 21, 353 (2019).
    
    

88. [88]. A. Tokariev,  J. A. Roberts,  A. Zalesky,  X. Zhao,  S. Vanhatalo,  M. Breakspear, and  L. Cocchi, Large-scale brain modes reorganize between infant sleep states and carry prognostic information for preterms, Nature communications 10, 1 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-09238-2&link_type=DOI) 

89. [89]. F. De Vico Fallani,  F. A. Rodrigues,  L. da Fontoura Costa,  L. Astolfi,  F. Cincotti,  D. Mattia,  S. Salinari, and  F. Babiloni, Multiple pathways analysis of brain functional networks from eeg signals: an application to real data, Brain topography 23, 344 (2011).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20614232&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

90. [90]. M. E. Newman, The structure and function of complex networks, SIAM review 45, 167 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/s003614450342480&link_type=DOI) 

91. [91]. M. E. Newman, Assortative mixing in networks, Physical review letters 89, 208701 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevLett.89.208701&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12443515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

92. [92]. R. Albert and  A.-L. Barabási, Statistical mechanics of complex networks, Reviews of modern physics 74, 47 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/RevModPhys.74.47&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174548700003&link_type=ISI) 

93. [93]. L. C. Freeman, A set of measures of centrality based on betweenness, Sociometry, 35 (1977).
    
    

94. [94]. L. C. Freeman, Centrality in social networks conceptual clarification, Social networks 1, 215 (1978).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8733(78)90021-7&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GL19400002&link_type=ISI) 

95. [95]. P. Bonacich, Power and centrality: A family of measures, American journal of sociology 92, 1170 (1987).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2780000&link_type=DOI) 

96. [96]. R. Albert,  H. Jeong, and  A.-L. Barabási, Diameter of the world-wide web, nature 401, 130 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/43601&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000082458800041&link_type=ISI) 

97. [97]. J. M. Kleinberg, Hubs, authorities, and communities, ACM computing surveys (CSUR) 31, 5 (1999).
    
    

98. [98]. D. Eppstein,  M. S. Paterson, and  F. F. Yao, On nearestneighbor graphs, Discrete & Computational Geometry 17, 263 (1997).
    
    

99. [99]. J. Doyle and  J. Graver, Mean distance in a graph, Discrete Mathematics 17, 147 (1977).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0012-365X(77)90144-3&link_type=DOI) 

100.[100]. T. A. Snijders, The degree variance: an index of graph heterogeneity, Social networks 3, 163 (1981).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8733(81)90014-9&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1981MV68200001&link_type=ISI) 

101.[101]. M. Dehmer and  A. Mowshowitz, A history of graph entropy measures, Information Sciences 181, 57 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ins.2010.08.041&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000284511100004&link_type=ISI) 

102.[102]. D. J. Watts and  S. H. Strogatz, Collective dynamics of ‘small-world’networks, Nature 393, 440 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/30918&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9623998&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074020000035&link_type=ISI) 

103.[103]. M. E. Newman,  D. J. Watts, and  S. H. Strogatz, Random graph models of social networks, Proceedings of the National Academy of Sciences 99, 2566 (2002).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxNToiOTkvc3VwcGxfMS8yNTY2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMTYvMjAyMi4wNi4xNC4yMjI3NjQxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

104.[104]. S. B. Seidman, Network structure and minimum degree, Social networks 5, 269 (1983).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8733(83)90028-X&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1983RS06200002&link_type=ISI) 

105.[105]. M. Newman, Networks: an introduction (Oxford university press, 2010).
    
    

106.[106]. P. Hage and  F. Harary, Eccentricity and centrality in networks, Social networks 17, 57 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0378-8733(94)00248-9&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QM57300003&link_type=ISI) 

107.[107]. B. S. Anderson,  C. Butts, and  K. Carley, The interaction of size and density with graph-level indices, Social networks 21, 239 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0378-8733(99)00011-8&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084385700002&link_type=ISI) 

108.[108]. V. Latora and  M. Marchiori, Economic small-world behavior in weighted networks, The European Physical Journal B-Condensed Matter and Complex Systems 32, 249 (2003).
    
    

109.[109]. M. E. Newman, Communities, modules and large-scale structure in networks, Nature physics 8, 25 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nphys2162&link_type=DOI) 

110.[110]. J. Kim and  J.-G. Lee, Community detection in multilayer graphs: A survey, ACM SIGMOD Record 44, 37 (2015).
    
    

111.[111]. X. Zhao,  J. Liang, and  J. Wang, A community detection algorithm based on graph compression for large-scale social networks, Information Sciences 551, 358 (2021).
    
    

112.[112]. A. Clauset,  M. E. Newman, and  C. Moore, Finding community structure in very large networks, Physical review E 70, 066111 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.70.066111&link_type=DOI) 

113.[113]. M. Rosvall,  D. Axelsson, and  C. T. Bergstrom, The map equation, The European Physical Journal Special Topics 178, 13 (2009).
    
    

114.[114]. M. E. Newman, Finding community structure in networks using the eigenvectors of matrices, Physical review E 74, 036104 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.74.036104&link_type=DOI) 

115.[115]. U. N. Raghavan,  R. Albert, and  S. Kumara, Near linear time algorithm to detect community structures in largescale networks, Physical review E 76, 036106 (2007).
    
    

116.[116]. M. Girvan and  M. E. Newman, Community structure in social and biological networks, Proceedings of the national academy of sciences 99, 7821 (2002).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiOTkvMTIvNzgyMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA2LzE2LzIwMjIuMDYuMTQuMjIyNzY0MTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

117.[117]. J. Reichardt and  S. Bornholdt, Statistical mechanics of community detection, Physical review E 74, 016110 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1103/PhysRevE.74.016110&link_type=DOI) 

118.[118]. V. D. Blondel,  J.-L. Guillaume,  R. Lambiotte, and  E. Lefebvre, Fast unfolding of communities in large networks, Journal of statistical mechanics: theory and experiment 2008, P10008 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1742-5468/2008/10/P10008&link_type=DOI) 

119.[119]. P. Refaeilzadeh,  L. Tang, and  H. Liu, Cross-validation., Encyclopedia of database systems 5, 532 (2009).
    
    

120.[120]. D. Berrar, Cross-validation. (2019).
    
    

121.[121]. Y. Bengio and  Y. Grandvalet, No unbiased estimator of the variance of k-fold cross-validation, Journal of machine learning research 5, 1089 (2004).
    
    

122.[122]. A. A. Shah and  Y. D. Khan, Identification of 4-carboxyglutamate residue sites based on position based statistical feature and multiple classification, Scientific Reports 10, 1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-73501-6&link_type=DOI) 

123.[123]. T. Kawamoto and  Y. Kabashima, Cross-validation estimate of the number of clusters in a network, Scientific reports 7, 1 (2017).
    
    

124.[124]. J. Chan,  T. Rea,  S. Gollakota, and  J. E. Sunshine, Contactless cardiac arrest detection using smart devices, NPJ digital medicine 2, 1 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41746-019-0103-3&link_type=DOI) 

125.[125]. Y. Xu,  S. Zomer, and  R. G. Brereton, Support vector machines: a recent method for classification in chemometrics, Critical Reviews in Analytical Chemistry 36, 177 (2006).
    
    

126.[126]. M. Awad and  R. Khanna, Support vector machines for classification, in Efficient learning machines (Springer, 2015) pp. 39–66.
    
    

127.[127]. M. Sato,  K. Morimoto,  S. Kajihara,  R. Tateishi,  S. Shiina,  K. Koike, and  Y. Yatomi, Machine-learning approach for the development of a novel predictive model for the diagnosis of hepatocellular carcinoma, Scientific reports 9, 1 (2019).
    
    

128.[128]. Z. Zhong,  X. Yuan,  S. Liu,  Y. Yang, and  F. Liu, Machine learning prediction models for prognosis of critically ill patients after open-heart surgery, Scientific Reports 11, 1 (2021).
    
    

129.[129]. F. Arcadu,  F. Benmansour,  A. Maunz,  J. Willis,  Z. Haskova, and  M. Prunotto, Author correction: Deep learning algorithm predicts diabetic retinopathy progression in individual patients, NPJ digital medicine 3, 1 (2020).
    
    

130.[130]. C. Krittanawong,  H. U. H. Virk,  A. Kumar,  M. Aydar,  Z. Wang,  M. P. Stewart, and  J. L. Halperin, Machine learning and deep learning to predict mortality in patients with spontaneous coronary artery dissection, Scientific reports 11, 1 (2021).
    
    

131.[131]. H. H. Rashidi,  S. Sen,  T. L. Palmieri,  T. Blackmon,  J. Wajda, and  N. K. Tran, Early recognition of burn-and trauma-related acute kidney injury: a pilot comparison of machine learning techniques, Scientific reports 10, 1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-73501-6&link_type=DOI) 

132.[132]. A. Mincholé and  B. Rodriguez, Artificial intelligence for the electrocardiogram, Nature medicine 25, 22 (2019).
    
    

133.[133]. Y. Tolkach,  T. Dohmgörgen,  M. Toma, and  G. Kristiansen, High-accuracy prostate cancer pathology using deep learning, Nature Machine Intelligence 2, 411 (2020).
    
    

134.[134]. J. Dukart,  S. Weis,  S. Genon, and  S. B. Eickhoff, Towards increasing the clinical applicability of machine learning biomarkers in psychiatry, Nature Human Behaviour 5, 431 (2021).
    
    

135.[135]. R. C. Li,  S. M. Asch, and  N. H. Shah, Developing a delivery science for artificial intelligence in healthcare, NPJ digital medicine 3, 1 (2020).
    
    

136.[136]. Y. Park and  M. Kellis, Deep learning for regulatory genomics, Nature biotechnology 33, 825 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nbt.3313&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26252139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

137.[137]. Y. Ito,  M. Unagami,  F. Yamabe,  Y. Mitsui,  K. Nakajima,  K. Nagao, and  H. Kobayashi, A method for utilizing automated machine learning for histopathological classification of testis based on johnsen scores, Scientific reports 11, 1 (2021).
    
    

138.[138]. J. Kim,  J. Lee,  E. Park, and  J. Han, A deep learning model for detecting mental illness from user content on social media, Scientific reports 10, 1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-73501-6&link_type=DOI) 

139.[139]. Y. Li,  C. M. Nowak,  U. Pham,  K. Nguyen, and  L. Bleris, Cell morphology-based machine learning models for human cell state classification, NPJ systems biology and applications 7, 1 (2021).
    
    

140.[140]. X. Yu,  W. Pang,  Q. Xu, and  M. Liang, Mammographic image classification with deep fusion learning, Scientific Reports 10, 1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-73501-6&link_type=DOI) 

141.[141]. S. Berryman,  K. Matthews,  J. H. Lee,  S. P. Duffy, and  H. Ma, Image-based phenotyping of disaggregated cells using deep learning, Communications Biology 3, 1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s42003-020-01468-1&link_type=DOI) 

142.[142]. S. Yang,  J. Kweon,  J.-H. Roh,  J.-H. Lee,  H. Kang,  L.-J. Park,  D. J. Kim,  H. Yang,  J. Hur,  D.-Y. Kang, et al., Deep learning segmentation of major vessels in x-ray coronary angiography, Scientific reports 9, 1 (2019).
    
    

143.[143]. A. Y. Hannun,  P. Rajpurkar,  M. Haghpanahi,  G. H. Tison,  C. Bourn,  M. P. Turakhia, and  A. Y. Ng, Cardiologist-level arrhythmia detection and classification in ambulatory electrocardiograms using a deep neural network, Nature medicine 25, 65 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-018-0268-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

144.[144]. M. Bracher-Smith,  K. Crawford, and  V. Escott-Price, Machine learning for genetic prediction of psychiatric disorders: a systematic review, Molecular Psychiatry 26, 70 (2021).
    
    

145.[145]. D. Patel,  V. Kher,  B. Desai,  X. Lei,  S. Cen,  N. Nanda,  A. Gholamrezanezhad,  V. Duddalwar,  B. Varghese, and  A. A. Oberai, Machine learning based predictors for covid-19 disease severity, Scientific Reports 11, 1 (2021).
    
    

146.[146]. R. Rodríguez-Pérez and  J. Bajorath, Interpretation of machine learning models using shapley values: application to compound potency and multi-target activity predictions, Journal of computer-aided molecular design 34, 1013 (2020).
    
    

147.[147]. L. S. Shapley, A value for n-person games, contributions to the theory of games, 2, 307–317 (1953).
    
    

148.[148]. I. Parrachino, Cooperative game theory and its application to natural, environmental and water resource issues, Vol. 4072 (World Bank Publications, 2012).
    
    

149.[149]. C. Molnar, Interpretable machine learning (Lulu. com, 2020).
    
    

150.[150]. M. Mijalkov,  E. Kakaei,  J. B. Pereira,  E. Westman,  G. Volpe, and  A. D. N. Initiative, Braph: a graph theory software for the analysis of brain connectivity, PloS one 12, e0178798 (2017).
    
    

151.[151]. C. M. Michel and  D. Brunet, Eeg source imaging: a practical review of the analysis steps, Frontiers in neurology 10, 325 (2019).
    
    

152.[152]. E. E. Asher,  M. Plotnik,  M. Günther,  S. Moshel,  O. Levy,  S. Havlin,  J. W. Kantelhardt, and  R. P. Bartsch, Connectivity of eeg synchronization networks increases for parkinson’s disease patients with freezing of gait, Communications biology 4, 1 (2021).
    
    

153.[153]. F. Amzica and  M. Steriade, Electrophysiological correlates of sleep delta waves, Electroencephalography and clinical neurophysiology 107, 69 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0013-4694(98)00051-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9751278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000075760900001&link_type=ISI) 

154.[154]. G. Bernardi,  M. Betta,  E. Ricciardi,  P. Pietrini,  G. Tononi, and  F. Siclari, Regional delta waves in human rapid eye movement sleep, Journal of Neuroscience 39, 2686 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjEwOiIzOS8xNC8yNjg2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMTYvMjAyMi4wNi4xNC4yMjI3NjQxMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

155.[155]. N. I. Ardeshna, Eeg and coma, The Neurodiagnostic Journal 56, 1 (2016).
    
    

156.[156]. J. Frohlich,  D. Toker, and  M. M. Monti, Consciousness among delta waves: a paradox?, Brain 144, 2257 (2021).
    
    

157.[157]. D. Biello, Searching for god in the brain, Scientific American Mind 18, 38 (2007).
    
    

158.[158]. A. Sadeghi Habibabad,  J.-e.-D. MahdiNejad,  H. Azemati, and  P. Matracchi, Using neurology sciences to investigate the color component and its effect on promoting the sense of spirituality in the interior space of the vakil mosque of shiraz (using quantitative electroencephalography wave recording), Journal of Religion and Health, 1 (2019).
    
    

159.[159]. M. Beauregard and  V. Paquette, Eeg activity in carmelite nuns during a mystical experience, Neuroscience letters 444, 1 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neulet.2008.08.028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18721862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

160.[160]. J.-P. Banquet, Spectral analysis of the eeg in meditation, Electroencephalography and clinical neurophysiology 35, 143 (1973).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0013-4694(73)90170-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=4124606&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1973Q506200004&link_type=ISI) 

161.[161]. P. Kora,  K. Meenakshi,  K. Swaraja,  A. Rajani, and  M. S. Raju, Eeg based interpretation of human brain activity during yoga and meditation using machine learning: A systematic review, Complementary Therapies in Clinical Practice 43, 101329 (2021).
    
    

162.[162]. C. Pallavicini,  F. Cavanna,  F. Zamberlan,  L. A. de la Fuente,  M. Arias,  M. C. Romero,  R. Carhart-Harris,  C. Timmermann, and  E. Tagliazucchi, Neural and subjective effects of inhaled dmt in natural settings, bioRxiv (2020).
    
    

163.[163]. C. L. Scrivener and  A. T. Reader, Variability of eeg electrode positions and their underlying brain regions: visualizing gel artifacts from a simultaneous eeg-fmri dataset, Brain and Behavior, e2476 (2022).
    
    

164.[164]. G. Gainotti, Face familiarity feelings, the right temporal lobe and the possible underlying neural mechanisms, Brain research reviews 56, 214 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.brainresrev.2007.07.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17822771&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000251871200011&link_type=ISI) 

165.[165]. J. Doyon and  B. Milner, Right temporal-lobe contribution to global visual processing, Neuropsychologia 29, 343 (1991).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0028-3932(91)90024-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1886679&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

166.[166]. J. Brand,  M. Piccirelli,  M.-C. Hepp-Reymond,  K. Eng, and  L. Michels, Brain activation during visually guided finger movements, Frontiers in Human Neuroscience, 309 (2020).
    
    

167.[167]. A. C. Pierno,  F. Tubaldi,  L. Turella,  P. Grossi,  L. Barachino,  P. Gallo, and  U. Castiello, Neurofunctional modulation of brain regions by the observation of pointing and grasping actions, Cerebral Cortex 19, 367 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/bhn089&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18534989&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000262518800011&link_type=ISI) 

168.[168]. D. J. Cobia,  M. J. Smith,  L. Wang, and  J. G. Csernansky, Longitudinal progression of frontal and temporal lobe changes in schizophrenia, Schizophrenia research 139, 1 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2012.05.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22647883&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

169.[169]. M. Catani, The anatomy of the human frontal lobe, Handbook of clinical neurology 163, 95 (2019).
    
    

170.[170]. L. Banker and  P. Tadi, Neuroanatomy, precentral gyrus.[updated 2020 jul 31], StatPearls [Internet]. StatPearls Publishing: Treasure Island (FL) (2021).
    
    

171.[171]. A. Lingnau and  P. E. Downing, The lateral occipitotemporal cortex in action, Trends in cognitive sciences 19, 268 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2015.03.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25843544&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

172.[172]. A. I. Luppi,  R. L. Carhart-Harris,  L. Roseman,  I. Pappas,  D. K. Menon, and  E. A. Stamatakis, Lsd alters dynamic integration and segregation in the human brain, NeuroImage 227, 117653 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2020.117653&link_type=DOI) 

173.[173]. M. Rubinov and  O. Sporns, Complex network measures of brain connectivity: uses and interpretations, Neuroimage 52, 1059 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2009.10.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19819337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280181800027&link_type=ISI) 

174.[174]. G. Tononi,  O. Sporns, and  G. M. Edelman, A measure for brain complexity: relating functional segregation and integration in the nervous system, Proceedings of the National Academy of Sciences 91, 5033 (1994).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiOTEvMTEvNTAzMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA2LzE2LzIwMjIuMDYuMTQuMjIyNzY0MTAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

175.[175]. O. Sporns, Network analysis, complexity, and brain function, Complexity 8, 56 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cplx.10047&link_type=DOI) 

176.[176]. O. Sporns, Network attributes for segregation and integration in the human brain, Current opinion in neurobiology 23, 162 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.conb.2012.11.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23294553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

177.[177]. E. Bullmore and  O. Sporns, The economy of brain network organization, Nature Reviews Neuroscience 13, 336 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn3214&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22498897&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F16%2F2022.06.14.22276410.atom) 

178.[178]. D. Rangaprakash,  M. N. Dretsch,  J. S. Katz,  T. S. Denney Jr., and  G. Deshpande, Dynamics of segregation and integration in directional brain networks: illustration in soldiers with ptsd and neurotrauma, Frontiers in neuroscience, 803 (2019).
    
    

179.[179]. W. Luo,  A. S. Greene, and  R. T. Constable, Within node connectivity changes, not simply edge changes, influence graph theory measures in functional connectivity studies of the brain, NeuroImage 240, 118332 (2021).
    
    

180.[180]. R. L. Carhart-Harris,  S. Muthukumaraswamy,  L. Roseman,  M. Kaelen,  W. Droog,  K. Murphy,  E. Tagliazucchi,  E. E. Schenberg,  T. Nest,  C. Orban, et al., Neural correlates of the lsd experience revealed by multimodal neuroimaging, Proceedings of the National Academy of Sciences 113, 4853 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEzLzE3LzQ4NTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNi8xNi8yMDIyLjA2LjE0LjIyMjc2NDEwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

181.[181]. R. L. Carhart-Harris,  R. Leech,  P. J. Hellyer,  M. Shanahan,  A. Feilding,  E. Tagliazucchi,  D. R. Chialvo, and  D. Nutt, The entropic brain: a theory of conscious states informed by neuroimaging research with psychedelic drugs, Frontiers in human neuroscience, 20 (2014).
    
    

182.[182]. N. Farnes,  B. E. Juel,  A. S. Nilsen,  L. G. Romundstad, and  J. F. Storm, Increased signal diversity/complexity of spontaneous eeg, but not evoked eeg responses, in ketamine-induced psychedelic state in humans, Plos one 15, e0242056 (2020).