BioMed Research International
Volume 2015 (2015), Article ID 825136, 10 pages
http://dx.doi.org/10.1155/2015/825136

Research Article

MIC as an Appropriate Method to Construct the Brain Functional Network

Ziqing Zhang,1,2 Shu Sun,³ Ming Yi,¹ Xia Wu,4,5 and Yiming Ding¹

¹Key Laboratory of Magnetic Resonance in Biological Systems, Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, China
²University of Chinese Academy of Sciences, Beijing 100049, China
³School of Mathematics and Statistics, Wuhan University, Wuhan 430072, China
⁴College of Information Science and Technology, Beijing Normal University, Beijing 100875, China
⁵State Key Laboratories of Transducer Technology, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China

Received 29 July 2014; Revised 12 October 2014; Accepted 14 October 2014

Academic Editor: Yoshiki Kaneoke

Copyright © 2015 Ziqing Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Using an effective method to measure the brain functional connectivity is an important step to study the brain functional network. The main methods for constructing an undirected brain functional network include correlation coefficient (CF), partial correlation coefficient (PCF), mutual information (MI), wavelet correlation coefficient (WCF), and coherence (CH). In this paper we demonstrate that the maximal information coefficient (MIC) proposed by Reshef et al. is relevant to constructing a brain functional network because it performs best in the comprehensive comparisons in consistency and robustness. Our work can be used to validate the possible new functional connection measures.

1. Introduction

Functional connectivity between brain areas is a hotspot in the field of cognitive neuroscience. In brain functional networks, connectivity is often measured using some form of statistical correlation. Brain activity in one region is correlated with activity in another region to quantify the strength and identify statistical dependency. When repeated for every possible pair of regions, the result is a network characterization of the brain’s connectivity, in which brain regions represent network nodes and correlation strengths correspond to connection weights [1, 2].

To date, the main methods of measuring brain functional connectivity include correlation coefficient (CF), partial correlation coefficient (PCF), mutual information (MI), wavelet correlation coefficient (WCF), and coherence (CH). Correlation analysis is the simplest method for analyzing brain functional connectivity; and it is widely used [3–5]. Zalesky et al. [1] recommended using this method to construct a brain functional network. However, contaminations from noise, such as cardiac and blood vessel activities in the brain, could also lead to high correlations [6], and the correlation coefficient only measures linear dependence [7]. Similar to many commonly used statistics, a correlation coefficient is not robust; its value may be misleading if outliers are present [8, 9]. Partial correlation refers to the normalized correlation between two time series after each series has been adjusted by regressing out all other time series of network nodes. One attractive feature of this method is that it attempts to distinguish direct from indirect connections [10]. Marrelec et al. [11] advocated partial correlation as a suitable method to construct a brain functional network. Some researchers have used this method to study the default mode network or the difference between patients and controls [12–15]. Partial correlation is a type of conditional correlation; it still cannot measure a nonlinear association [16, 17]. Coherence is the spectral representation of correlation in the frequency domain and was proposed by Sun et al. [18]. The expression of correlation in the frequency domain enables researchers to study the time course relationship in a natural and intrinsic manner [19]. This method has recently been used by Chang and Glover [20] to investigate nonstationary effects in resting functional magnetic resonance imaging (fMRI) data. However, coherence provides vague information on the actual cortical areas involved because of the complex relationship between the active brain areas and sensor recordings [21]. Wavelet correlation refers to the correlation between wavelet coefficients, which can be obtained from a discrete wavelet transform (DWT) of the time series. Bullmore et al. [22] studied wavelets and the statistical analysis of fMRI of the human brain. Skidmore et al. [23] used wavelet correlation to construct a brain functional network and identified the differences between healthy subjects and subjects with Parkinson’s disease. Mutual information (MI), a method in information theory, quantifies the shared information between two variables and can reflect both linear and nonlinear dependencies [24]. This method is popular for measuring the brain functional connectivity [25–27].

In 2011, Reshef et al. [28] proposed a new measure named maximal information coefficient (MIC), which can capture both linear and nonlinear association between two variables. Because of its outstanding performance in measuring different kinds of dependences, it is considered the correlation for the 21st century [29]. Reshef et al. [30] noted that MIC is more equitable compared with natural alternatives, such as mutual information estimation and distance correlation. Since it was proposed, it has been widely used [31–36]. Su et al. [37] were the first to use MIC to construct a brain functional network. However, they only demonstrated that the brain functional connectivity between healthy subjects and schizophrenia subjects calculated using MIC is suitable for classification. As a novel method applied in cognitive neuroscience, there is no explanation accounting for the rationality or advantages of MIC.

The paper demonstrates that MIC is a suitable method to measure brain functional connectivity and illustrates its advantages in constructing a brain functional connectivity network. Thirteen healthy subjects with minimal differences were selected from 75 healthy subjects, and the brain functional network was constructed using MIC, as well as other methods (CF, PCF, MI, WAF, and CH). Based on the node importance, we compared the consistency and robustness of the methods from different aspects of the network. Compared with other methods, MIC performed better in terms of consistency and robustness. Although there are many measures provided to capture the functional connections between brain areas, there is no work to compare the performance of the measures. Our work can be used to validate the possible new functional connection measures.

2. Materials and Methods

2.1. Subject Information and Data Preprocessing

The data utilized for this paper is available for download at http://fcon_1000.projects.nitrc.org/indi/retro/cobre.html. The study comprised 75 healthy samples (ages ranged from 18 to 65 years in each group). To control external variables for our research, we carefully selected 13 samples with the smallest differences. The ID of samples are 40020, 40051, 40056, 40065, 40091, 40093, 40104, 40114, 40115, 40120, 40127, 40128, and 40129. All subjects were male and right handed according to the Edinburgh Handedness Inventory. Mean age was for the samples. Detailed subject information is provided in Table 1.

Table 1: Sample information.

Age range (mean age)
Sex
Right handed

Healthy samples
13
22–27 (24.6)
All male
All

All fMRI data were processed using SPM8 (http://www.fil.ion.ucl.ac.uk/spm/) and DPARSF-V2.0 (http://www.restfmri.net/forum/index.php) [38]. For each subject, we removed the first 10 volume images from the RS-fMRI data for scanner stabilization and subject adaptation to the environment, which left 140 volumes for further analysis. Then, we performed slice timing to correct for the acquisition time delay between slices within the same TR; realignment to the first volume to correct inter-TR head motions was performed, followed by spatial normalization to a standard MNI template and resampling to a voxel size of 3 × 3 × 3 mm³. No spatial smoothing was applied based on methods from previous studies [39–41]. Finally, we performed bandpass filtering for each voxel in the frequency of 0:01–0:08 Hz to reduce low-frequency drift and high-frequency physiological noise. The RS-fMRI data for each subject were checked for head motion. No subject was excluded according to the criteria that the translation and rotation of head motion in any direction were not more than 1.5 mm or 1.5°. To obtain signals for each region, we applied an automated anatomical labeling (AAL) atlas [42] to parcellate the brain into 90 regions of interest (ROIs) (45 per hemisphere). The names of the ROIs and their corresponding abbreviations are listed in Table 2. The time series for each ROI was calculated by averaging the signals of all voxels within that region.

Table 2: Names and abbreviations of the ROIs [42].

Regions
Abbr.
Regions
Abbr.

Superior frontal gyrus, dorsolateral
SFGdor
Superior parietal gyrus
SPG

Superior frontal gyrus, orbital
SFGorb
Paracentral lobule
PCL

Superior frontal gyrus, medial
SFGmed
Postcentral gyrus
PoCG

Superior frontal gyrus, medial orbital
SFGmorb
Inferior parietal gyrus
IPG

Middle frontal gyrus
MFG
Supramarginal gyrus
SMG

Middle frontal gyrus, orbital
MFGorb
Angular gyrus
ANG

Inferior frontal gyrus, opercular
IFGoper
Precuneus
PCNU

Inferior frontal gyrus, triangular
IFGtri
Posterior cingulate gyrus
PCC

Inferior frontal gyrus, orbital
IFGorb

Gyrus rectus
REG
Insula
INS

Anterior cingulate gyrus
ACC
Thalamus
THA

Olfactory cortex
OLF

Superior temporal gyrus
STG

Precentral gyrus
PreCG
Superior temporal gyrus, temporal pole
STGp

Supplementary motor area
SMA
Middle temporal gyrus
MTG

Rolandic operculum
ROL
Middle temporal gyrus, temporal pole
MTGp

Median and paracingulate gyrus
MCC
Inferior temporal gyrus
ITG

Calcarine fissure and surrounding cortex
CAL

Heschl gyrus
HES
Hippocampus
HIP

Cuneus
CUN
Parahippocampal gyrus
PHIP

Lingual gyrus
LING
Amygdala
AMYG

Superior occipital gyrus
SOG

Middle occipital gyrus
MOG
Caudate nucleus
CAU

Inferior occipital gyrus
IOG
Lenticular nucleus, putamen
PUT

Fusiform gyrus
FG
Lenticular nucleus, pallidum
PAL

2.2. Maximal Information Coefficient

The MIC, introduced by Reshef et al. [28] in 2011, was used as a measure of association between two random variables and . The MIC can capture wide range of relationships. The is the mutual information between random variables and normalized by the minimum entropy of and , which can be written as

where is the conditional entropy, which is the amount of information needed to describe the outcome of given that the value of is known. For a pair of variables and , and if and only if and are independent. The MIC is robust to outliers because the estimations of Shannon entropy and conditional entropy are robust [36, 43].

2.3. Construct Network

For each subject, we obtained a 90 × 90 dependence matrix by calculating the connection strength using one of the 6 previously described methods (CF, PCF, MIC, MI, WCF, and CH) between all ROI pairs. One can check the details about the definitions and computations of the methods in the references. Figure 1 is the visualization of the six dependence matrices of a randomly selected sample. The correlation matrix was thresholded into a binary matrix. By taking each ROI as a node and the functional connectivity as an edge, we obtained a 90 90 adjacency matrix for each subject. The adjacency matrix can be defined as

where is the connection matrix, is the threshold value, and is the adjacency matrix under threshold . In other words, if the absolute value of between a pair of brain regions and exceeds a given threshold , an edge is constructed to connect the brain regions; otherwise, there is no edge between them. We know that is an external variable which determines the size of the network. If is too large, we will obtain a network with fewer edges, which may lead to a disconnected network. If is too small, some of the connection strength may be too weak to be significant. To balance these two aspects, 1200 edges were selected as a representative network size. We also have obtained results on different network sizes, and the results were nonsensitive to network size.

Figure 1: The connection matrices (90 × 90) of a randomly chosen subject from 13 samples were calculated by the 6 previously described methods (CF, PCF, CH, WCF, MIC, MI). The horizontal and vertical axes are the 90 ROIs. The color intensity represents the connection strength. One can get an intuitive understanding of the brain functional connectivity. These connection matrices are the basis for our subsequent research.

2.4. Comparison Scheme

We try to conduct both intermethod and intersample comparisons to validate the consistency and robustness. The important results, say key nodes, obtained by a good method should not conflict with that obtained by known methods. So we conduct intermethod comparisons to see which method is more consistent with other methods. Such an aspect we compare is called consistency. Furthermore, a good method should not be sensitive to its operation objects. With 13 subjects of similar physical condition, a good method should obtain similar network properties. This approach represents a type of repeated trials; that is, a good method should always be good, regardless of the experimental targets. Such an aspect we compare is called robustness.

2.4.1. Node Importance

We use two criterions to evaluate the relative importance of a node in a network [43]. One criterion is the degree centrality (DC), which is proportional to the degree of the node. Another criterion is Shannon-Parry centrality (SPC). The former is a popular method. The latter is based on the Shannon-Parry measure of a network and the relative importance of the th node is proportional to , where and are the left and right eigenvectors of the adjacency matrix of the network. The SPC can effectively illustrate the node importance by synthesizing the node properties and network topology structure [43, 44].

2.4.2. Comparisons on Consistency and Robustness

We conduct the comparisons of the consistency and robustness of the six methods by the flowchart in Figure 2.

Figure 2: The flowchart summarizes the process of our comparison scheme on two aspects: consistency and robustness. Looking from the top to the bottom following the arrow and the branches, we can have a clear vision of our steps and the steps can be easily realized. We will conduct ten comparisons, four for consistency and six for robustness. To make it clear, “voting” is a metaphor. For instance, when we use 6 methods to extract the top 9 important nodes, we think that the 6 methods are “holding” a vote for their top 9 important nodes from the 90 nodes. The same goes for our samples. In the left branch of consistency part, the six methods have their own important nodes sets after “voting,” but the sets are different; that is to say, the six methods have different “opinions” on important nodes. We give each method a score, respectively, to decide which method’s “opinion” is the best. One method’s score is the total votes it receives from other methods, including itself. The method which acquires the most votes pools other methods’ “opinions” together and is considered more reliable. In the first left branch of the robustness part, the 13 samples “vote”    important nodes from the 90 nodes. After “voting,” the 90 nodes have their voting numbers and voting rates (voting number/sum of vote). We can calculate the entropy according to the probability distribution induced by the voting rates. If the entropy is small, it means that the 13 samples have consensus on important nodes, and the corresponding method is more robust.

3. Results

3.1. Consistency

We compared the consistency from top 10% important nodes (A) and the total Euclidean distance between importance vectors (B). Since there were two importance criterions DC and SPC, there were four comparisons. The computational results are listed in Tables 3–5.

Table 3: In Table 3(a), importance of nodes was defined by degree centrality (DC). In Table 3(b), importance of nodes was defined by the Shannon-Parry centrality (SPC). Each column included the important nodes and the score information for a method. In this table, the top 9 (10%) nodes were chosen as the important nodes. If a method obtained a higher score, it was more consistent with all other methods. In Table 3(a), CF, MI, and MIC had the relatively higher score. Although MIC was ranked the third, its score was very close to CF which is ranked the first. In general, the network constructed by MIC had the relatively better consistency compared with other methods in terms of important nodes. In Table 3(b), MIC had the highest score. Thus, the network constructed by MIC had the highest consistency compared with other methods in terms of important nodes.

(a)

CF
PCF
MIC
MI
WCF
CH

MFG. L
ORBsup. L
STG. L
MFG. L
SMG. R
PreCG. L

MFG. R
MFG. R
MFG. R
ROL. L
ORBmid. L
SFGdor. L

ROL. L
ORBmid. L
ROL. L
SFGmed. L
ROL. L
MFG. L

SFGmed. L
CUN. R
SFGmed. L
SFGmed. R
STG. R
MFG. R

SFGmed. R
PCUN. L
SFGmed. R
INS. L
SFGmed. L
ROL. L

INS. L
TPOsup. R
PoCG. R
INS. R
SFGmed. R
SFGmed. R

INS. R
MTG. L
STG. R
PoCG. R
INS. R
PCG. L

STG. L
TPOmid. R
ANG. L
STG. L
PoCG. R
FFG. R

STG. R
ITG. R
MFG. L/SMG. R
STG. R
ROL. R/PreCG. R
STG. L

Score
34.5
13
33
33.5
28.5
25.5

(b)

CF
PCF
MIC
MI
WCF
CH

MFG. L
ORBsup. L
MFG. L
MFG. L
PreCG. R
PreCG. L

ANG. L
MFG. R
MFG. R
ROL. L
ORBmid. L
PreCG. R

ROL. L
ORBmid. L
ROL. L
SFGmed. L
ROL. L
MFG. L

SFGmed. L
CUN. R
SFGmed. L
SFGmed. R
ROL. R
MFG. R

SFGmed. R
PCUN. L
SFGmed. R
SFGdor. L
SFGmed. L
SFGdor. R

INS. L
TPOsup. R
PoCG. R
INS. R
SFGmed. R
PoCG. L

INS. R
MTG. L
STG. R
PoCG. R
INS. R
PCG. L

STG. L
ORBsup. R
ANG. L
STG. L
PoCG. R
FFG. R

STG. R
ORBinf. L
STG. L
STG. R
STG. R
STG. L

Score
30
12
32
31
27
18

Table 4: In Table 4(a), we used degree centrality (DC) to measure the importance of node in network. In Table 4(b), we used Shannon-Parry centrality (SPC) to measure the importance of node in network. So each method got a 1 90 vector to measure the importance of nodes by calculating nodes’ average importance from 13 sample networks. The value in the middle of table was the Euclidean distance between vectors from different method. The value of last row was the sum of Euclidean distance from one method to the others. If this sum of distance was smaller, corresponding method had better consistency. From Table 4(a), CF and MIC had relatively small and very similar sum of distance, so MIC had better consistency than other methods. From Table 4(b), MIC had the smallest sum of distance, so MIC had the best consistency with other methods.

(a)

Method
CF
PCF
MIC
MI
WCF
CH

CF
0
48.77
13.90
11.10
22.00
30.70

PCF
48.77
0
47.48
50.36
44.32
40.34

MIC
13.90
47.48
0
14.85
23.80
27.81

MI
11.10
50.36
14.85
0
22.61
31.73

WCF
22.00
44.32
23.80
22.61
0
34.30

CH
30.70
40.34
27.81
31.73
34.30
0

Total
126.4
231.27
127.84
130.65
147.03
164.88

(b)

Method
CF
PCF
MIC
MI
WCF
CH

CF
0
0.0454
0.0117
0.0094
0.0191
0.0299

PCF
0.0454
0
0.0425
0.0442
0.0413
0.0332

MIC
0.0117
0.0425
0
0.0115
0.0192
0.0267

MI
0.0094
0.0442
0.0115
0
0.0180
0.0288

WCF
0.0191
0.0413
0.0192
0.0180
0
0.0328

CH
0.0299
0.0332
0.0267
0.0288
0.0328
0

Total
0.1155
0.2066
0.1116
0.1119
0.1304
0.1514

Table 5: A summary of the comparison of consistency. A: consistency in the important nodes; B: consistency in the importance of all nodes. The DC column indicated the importance of the node as defined by degree centrality (DC). The SPC column indicated the importance of the node as defined by the Shannon-Parry centrality (SPC). For the consistency comparison in four aspects, MIC had two times of the first ranking, one time of the second ranking, and one time of the third ranking. Overall, MIC had the best performance. Actually, if we pay attention to Tables 3 and 4, we will find that MIC, CF, and MI have a very similar conclusion.

Aspects
DC
SPC

CF
PCF
MIC
MI
WCF
CH
CF
PCF
MIC
MI
WCF
CH

A
1
6
3
2
4
5
3
6
1
2
4
5

B
1
6
2
3
4
5
3
6
1
2
4
5

Table 3 shows the comparison of important nodes; for MIC, the important nodes defined by DC included MEG. L, MEG. R, ROL. L, SFGmed. L, SFGmed. R, INS. L, INS. R, PoCG. R, STG. L, and STG. R. When the node importance was defined by the SPC, nearly the same important nodes were obtained (STG. R instead of SMG. R). MIC had the best rank if the node importance criterion was SPC and ranked the third if SPC was replaced by DC. In the latter case, the score of MIC (33) was very close to the highest score 34.5.

The comparison of the total Euclidean distances of all nodes is shown in Table 4; MIC had the best rank if the node importance criterion is SPC and ranked the second if SPC was replaced by DC. In the latter case, the total distance of MIC was 127.84, which was almost the same as the smallest total distance 126.4.

Table 5 collects the ranking information of the comparisons of consistency.

3.2. Robustness

Since robustness is an important feature of a method, we conducted six comparisons to examine the robustness of the methods. We calculated the total Euclidean distances between the importance vectors (C) and the importance ranking vectors (D) of the 13 samples. We had four comparisons since there were two importance criterions. On the other hand, we calculated the VC (variation coefficient) of the degrees of the 90 nodes, as well as the voting entropy of all nodes. The results are listed in Tables 6–10 and Figures 3 and 4.

Table 6: C: robustness in the importance of all nodes. We used degree centrality (DC) and Shannon-Parry centrality (SPC) to measure the importance of node in network and put the importance vector calculated in the same method from different sample as the row to construct matrix, so each method obtained a 13 90 matrix. The value in the table was the sum of Euclidean distance between different row vectors for one method. D: robustness in the ranking of all nodes in importance. We used degree and Parry measure to measure the importance of node in network and put the importance vector calculated in the same method from different sample as the row to construct matrix, so each method obtained a 13 90 matrix. Then we got the rank of node in each row, so each method got a 13 90 rank matrix. The value in the table was the sum of Euclidean distance between different row vectors from rank matrix of one method. If this sum of distance was smaller, corresponding method had better robustness. In the aspect C, MIC is ranked the second and is only bigger than PCF. In the aspect D, MIC had the smallest sum of distance, so MIC had better robustness than other methods.

Aspects

CF
PCF
MIC
MI
WCF
CH

C
DC
9148.4
7282.8
8344.8
9536.6
9750.4
10119

SPC
8.30
6.54
7.59
8.08
8.90
9.45

D
DC
24910
26483
24588
24889
25705
26544

SPC
24669
26574
24435
24677
25600
26680

Table 7: C: robustness in importance of all nodes. D: robustness in the ranking of all nodes in importance. The DC column indicated the importance of the node defined by degree centrality (DC). The SPC column indicated the importance of the node defined by the Shannon-Parry centrality (SPC). For the robustness comparison, MIC had the best ranking in D and was ranked second in C. Actually, PCF was very unstable. In aspect D, PCF had poor performance. In the part of discussion, we will proceed to explain this phenomenon. Overall, compared with other methods, MIC had the best performance in robustness.

Aspects
DC
SPC

CF
PCF
MIC
MI
WCF
CH
CF
PCF
MIC
MI
WCF
CH

Robustness
C
3
1
2
4
5
6
4
1
2
3
5
6

D
3
5
1
2
4
6
3
5
1
2
4
6

Table 8: We conducted one side two-sample -test between each two methods under the confidence level of 95%. In this chart we showed the value of the test in which the method’s mean in the row was smaller than that in the column. We arranged the methods according to their rank. The events in the symmetric position of the table were complementary, so we just showed the bold font data. We could see PCF passed 5 tests and was ranked first. MIC passed 4 tests and was ranked second.

PCF
MIC
CF
MI
WCF
CH

PCF

>0.9999
>0.9999
>0.9999
>0.9999
>0.9999

MIC

0.9966
0.9999
>0.9999
>0.9999

CF

0.8758
0.9415
0.9936

MI

0.5986
0.8714

WCF

0.8447

CH

Table 9: We conducted one side two-sample -test between each two methods under the confidence level of 95%. In this chart we showed the value of the test in which the method’s mean in the row was smaller than that in the column. We arranged the methods according to their rank. The events in the symmetric position of the table were complementary, so we just showed the bold font data. We could see MIC passed 5 tests, which performed best.

MIC
CF
MI
WCF
PCF
CH

MIC

0.7668
0.9556
0.9757
>0.9999
>0.9999

CF

0.8277
0.9011
0.9999
>0.9999

MI

0.6796
0.9998
>0.9999

WCF

0.9984
>0.9999

PCF

0.9870

CH

Table 10: In the robustness comparison, MIC was ranked the second in the aspect of VC of degree centrality (DC); MIC had the best ranking in terms of entropy of important nodes. Overall, compared with other methods, MIC performed best in terms of robustness.

Approach
CF
PCF
MIC
MI
WCF
CH

VC of DC
3
1
2
4
5
6

Entropy of important nodes
2
5
1
3
4
6

Figure 3: We used degree centrality (DC) to measure node importance and calculated VC of 90 nodes’ degree as to 13 samples. Thus, a 1 × 90 vector was obtained for each method. Box plots of these vectors were constructed to determine which method had a smaller average value. When the average VC was smaller, the method was considered to perform better in terms of robustness. The box plot shows the robustness rankings of PCF, MIC, CF, MI, WCF, and CH. Comparisons regarding this aspect indicated that MIC had a good performance. For statistical tests, see Table 8.

Figure 4: The horizontal axis is the number of important nodes for every subject according to the Shannon-Parry centrality (SPC). The vertical axis is the entropy as a metric to measure results; smaller entropy indicates better performance. Performance results were compared, and MIC performed best. For statistical tests, see Table 9.

As shown in Tables 6 and 7, in the comparisons of the total distance of all nodes, MIC was ranked the second. In the comparison of the distance between the importance ranking vectors, MIC was ranked first. The box plot (Figure 3) shows the average VC of all nodes, MIC performed the second, and PCF was ranked the first. The results were verified by the two-sample t-tests under the confidence level of 95%. In the comparison of the voting entropy of nodes between the 13 samples, MIC performed the best when the top 9–18 (10–20%) nodes were regarded as the important nodes, as shown in Figure 4. MIC passed all the two-sample t-tests under the confidence level of 95%. The ranking information is listed in Table 10. Overall, for the robustness comparisons, MIC was consistently ranked in the top two methods.

4. Discussions

We compared MIC to five existing methods comprehensively from consistency and robustness. According to the results in Tables 3–5, MIC, CF, and MI are more consistent than WCF, CH, and PCF. The consistency scores of MIC, CF, and MI are very close. Combine the two comparisons in important nodes with criterions DC and SPC in Table 3, the total scores of MIC, CF, and MI are 65, 64.5, and 64.5, respectively. So MIC is more consistent than CF and MI. For the comparisons based on total distance, one can see from Table 4 that MIC also is more consistent than CF and MI. From the four comparisons in consistency, MIC was ranked 3, 1, 2, and 1, CF was ranked 1, 3, 1, and 3, and MI was ranked 2, 2, 3, and 2. The sum of the ranks is 7, 8, and 9 for MIC, CF, and MI, respectively. We conclude that MIC was more consistent than CF and MI because the sum of the ranks is smaller.

For the robustness comparison, we compare six aspects. MIC was ranked the first and PCF was ranked the fifth in the comparisons of the total Euclidean distance between the importance (DC, SPC) ranking vectors and the voting entropy. In the remaining three comparisons, MIC was ranked the second, while PCF was ranked the first. PCF is very robust in the comparisons of total Euclidean distances between the importance (DC, SPC) vectors and variation coefficients of degrees. But it does not mean it is a good method because it performs the worst in all of the consistency comparisons. In fact, it regresses out all the other 88 nodes’ influence when calculating the correlation coefficient between a pair of nodes. This leads to relatively uniform results that lack discriminability. We refer to Figure 1 for a typical visualization of a PCF matrix.

5. Conclusions

In this paper, we demonstrate that MIC can be used to construct a brain functional network. We compared MIC with five other methods (CF, PCF, MI, WCF, and CH) and ensured that it is suitable for brain functional network construction. In the comprehensive comparisons in consistency and robustness, MIC performs the best, and the results are convincing.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Ziqing Zhang and Shu Sun contributed equally to this work.

Acknowledgments

This work was partially supported by the National Basic Research Program of China (973 program, 2011CB707802, 2013CB910200), the National Natural Science Foundation of China under 11275259 and 91330113, and the General Program of the National Natural Science Foundation of China (61222113), the Program for New Century Excellent Talents in University (NCET-12-0056), Open Foundation of State Key Laboratory of Networking and Switching Technology (Beijing University of Posts and Telecommunications) (SKLNST-2013-1-03), and State Key Laboratories of Transducer Technology (SKT1303). The authors are grateful to Professor Chaogan Yan for the use of his software DPARSF.

References

A. Zalesky, A. Fornito, and E. Bullmore, “On the use of correlation as a measure of network connectivity,” NeuroImage, vol. 60, no. 4, pp. 2096–2106, 2012.

K. J. Friston, “Functional and effective connectivity: a review,” Brain Connectivity, vol. 1, no. 1, pp. 13–36, 2011.

J. Cao and K. Worsley, “The geometry of correlation fields with an application to functional connectivity of the brain,” The Annals of Applied Probability, vol. 9, no. 4, pp. 1021–1057, 1999.

A. G. He, L. H. Tan, Y. Tang et al., “Modulation of neural connectivity during tongue movement and reading,” Human Brain Mapping, vol. 18, no. 3, pp. 222–232, 2003.

M. J. Lowe, M. Dzemidzic, J. T. Lurito, V. P. Mathews, and M. D. Phillips, “Correlations in low-frequency BOLD fluctuations reflect cortico-cortical connections,” NeuroImage, vol. 12, no. 5, pp. 582–587, 2000.

K. J. Friston, A. P. Holmes, K. J. Worsley, J.-P. Poline, C. D. Frith, and R. S. J. Frackowiak, “Statistical parametric maps in functional imaging: a general linear approach,” Human Brain Mapping, vol. 2, no. 4, pp. 189–210, 1994.

S. M. Stigler, “Francis Galton's account of the invention of correlation,” Statistical Science, vol. 4, no. 2, pp. 73–86, 1989.

S. J. Devlin, R. Gnanadesikan, and J. R. Kettenring, “Robust estimation and outlier detection with correlation coefficients,” Biometrika, vol. 62, no. 3, pp. 531–545, 1975.

P. J. Huber, Robust Statistics, John Wiley & Sons, New York, NY, USA, 2004.

S. M. Smith, K. L. Miller, G. Salimi-Khorshidi et al., “Network modelling methods for FMRI,” NeuroImage, vol. 54, no. 2, pp. 875–891, 2011.

G. Marrelec, A. Krainik, H. Duffau et al., “Partial correlation for functional brain interactivity investigation in functional MRI,” NeuroImage, vol. 32, no. 1, pp. 228–237, 2006.

P. Fransson and G. Marrelec, “The precuneus/posterior cingulate cortex plays a pivotal role in the default mode network: evidence from a partial correlation network analysis,” NeuroImage, vol. 42, no. 3, pp. 1178–1184, 2008.

P. Liu, Y. Zhang, G. Zhou et al., “Partial correlation investigation on the default mode network involved in acupuncture: an fMRI study,” Neuroscience Letters, vol. 462, no. 3, pp. 183–187, 2009.

J. Zhang, J. Wang, Q. Wu et al., “Disrupted brain connectivity networks in drug-naive, first-episode major depressive disorder,” Biological Psychiatry, vol. 70, no. 4, pp. 334–342, 2011.

Y. Liu, M. Liang, Y. Zhou et al., “Disrupted small-world networks in schizophrenia,” Brain, vol. 131, no. 4, pp. 945–961, 2008.

K. Baba, R. Shibata, and M. Sibuya, “Partial correlation and conditional correlation as measures of conditional independence,” Australian and New Zealand Journal of Statistics, vol. 46, no. 4, pp. 657–664, 2004.

X. Liang, H. J. Wang, and Y. He, “Human connectome: structural and functional brain networks,” Chinese Science Bulletin, vol. 55, pp. 1565–1583, 2010 (Chinese).

F. T. Sun, L. M. Miller, and M. D'Esposito, “Measuring interregional functional connectivity using coherence and partial coherence analyses of fMRI data,” NeuroImage, vol. 21, no. 2, pp. 647–658, 2004.

K. Li, L. Guo, J. Nie, G. Li, and T. Liu, “Review of methods for functional brain connectivity detection using fMRI,” Computerized Medical Imaging and Graphics, vol. 33, no. 2, pp. 131–139, 2009.

C. Chang and G. H. Glover, “Time-frequency dynamics of resting-state brain connectivity measured with fMRI,” NeuroImage, vol. 50, no. 1, pp. 81–98, 2010.

J. Gross, J. Kujala, M. Hämäläinen, L. Timmermann, A. Schnitzler, and R. Salmelin, “Dynamic imaging of coherent sources: studying neural interactions in the human brain,” Proceedings of the National Academy of Sciences of the United States of America, vol. 98, no. 2, pp. 694–699, 2001.

E. Bullmore, J. Fadili, M. Breakspear, R. Salvador, J. Suckling, and M. Brammer, “Wavelets and statistical analysis of functional magnetic resonance images of the human brain,” Statistical Methods in Medical Research, vol. 12, no. 5, pp. 375–399, 2003.

F. Skidmore, D. Korenkevych, Y. Liu, G. He, E. Bullmore, and P. M. Pardalos, “Connectivity brain networks based on wavelet correlation analysis in Parkinson fMRI data,” Neuroscience Letters, vol. 499, no. 1, pp. 47–51, 2011.

C. E. Shannon, “A mathematical theory of communication,” Bibliometrics, vol. 5, no. 1, pp. 3–55, 2001.

J. Jeong, J. C. Gore, and B. S. Peterson, “Mutual information analysis of the EEG in patients with Alzheimer's disease,” Clinical Neurophysiology, vol. 112, no. 5, pp. 827–835, 2001.

C.-C. Chen, J.-C. Hsieh, Y.-Z. Wu et al., “Mutual-information-based approach for neural connectivity during self-paced finger lifting task,” Human Brain Mapping, vol. 29, no. 3, pp. 265–280, 2008.

M.-E. Lynall, D. S. Bassett, R. Kerwin et al., “Functional connectivity and brain networks in schizophrenia,” The Journal of Neuroscience, vol. 30, no. 28, pp. 9477–9487, 2010.

D. N. Reshef, Y. A. Reshef, H. K. Finucane et al., “Detecting novel associations in large data sets,” Science, vol. 334, no. 6062, pp. 1518–1524, 2011.

T. Speed, A Correlation for the 21st Century, Science, 2011.

D. Reshef, Y. Reshef, M. Mitzenmacher, and P. Sabeti, “Equitability analysis of the maximal information coefficient, with comparisons,” In press, http://arxiv.org/abs/1301.6314

D. Albanese, M. Filosi, R. Visintainer, S. Riccadonna, G. Jurman, and C. Furlanello, “Minerva and minepy: a C engine for the MINE suite and its R, Python and MATLAB wrappers,” Bioinformatics, vol. 29, no. 3, pp. 407–408, 2013.

V. N. Kovacheva, A. M. Khan, M. Khan, D. B. A. Epstein, and N. M. Rajpoot, “DiSWOP: a novel measure for cell-level protein network analysis in localized proteomics image data,” Bioinformatics, vol. 30, no. 3, pp. 420–427, 2014.

Y. Zhang, Q. Hu, W. Zhang, and J. Liu, “A novel Bayesian network structure learning algorithm based on maximal information coefficient,” in Proceedings of the IEEE 5th International Conference on Advanced Computational Intelligence (ICACI '12), pp. 862–867, IEEE, Nanjing, China, October 2012.

Y. Zhang, W. Zhang, and Y. Xie, “Improved heuristic equivalent search algorithm based on maximal information coefficient for bayesian network structure learning,” Neurocomputing, vol. 117, pp. 186–195, 2013.

C. Lin, H. Canhao, T. Miller et al., “Maximal information coefficient for feature selection for clinical document classification,” in Proceedings of the ICML Workshop on Machine Learning for Clinical Data, 2012.

Q. Tan, H. Jiang, and Y. Ding, “Model selection method based on maximal information coefficient of residuals,” Acta Mathematica Scientia, Series B: English Edition, vol. 34, no. 2, pp. 579–592, 2014.

L. Su, L. Wang, H. Shen, G. Feng, and D. Hu, “Discriminative analysis of non-linear brain connectivity in schizophrenia: an fMRI Study,” Frontiers in Human Neuroscience, vol. 7, article 702, 2013.

Y. Chao-Gan and Z. Yu-Feng, “DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI,” Frontiers in Systems Neuroscience, vol. 4, article 13, 2010.

U. Braun, M. M. Plichta, C. Esslinger et al., “Test-retest reliability of resting-state connectivity network characteristics using fMRI and graph theoretical measures,” NeuroImage, vol. 59, no. 2, pp. 1404–1412, 2012.

J. Wang, L. Wang, Y. Zang et al., “Parcellation-dependent small-world brain functional networks: a resting-state fMRi study,” Human Brain Mapping, vol. 30, no. 5, pp. 1511–1523, 2009.

S. Achard and E. Bullmore, “Efficiency and cost of economical brain functional networks,” PLoS Computational Biology, vol. 3, article e17, 2007.

N. Tzourio-Mazoyer, B. Landeau, D. Papathanassiou et al., “Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain,” NeuroImage, vol. 15, no. 1, pp. 273–289, 2002.

S. Huang, H. Cui F, and Y. Ding M, “Evaluation of node imporyance in complex networks,” http://arxiv.org/abs/1402.5743

C. E. Shannon and W. Weaver, A Mathematical Theory of Communication, 1948.