Since the conventional mass-univariate testing might be underpowered with a low contrast-to-noise ratio of the synchronization matrix (Zalesky et al., 2010), network-based statistics (NBS; Zalesky et al., 2010, 2012) were performed to identify the node pairs between which the SL value was significantly changed in BD participants compared to controls. For this purpose, two-sample t-tests were independently performed at each synchronization value, and t-statistics larger than an uncorrected threshold of t = 2.61 (p = 0.005) were extracted into a set of supra-threshold connections. Then we identified all connected components in the adjacency matrix of supra-threshold links and saved the number of links. Finally, a permutation test was performed 10,000 times to estimate the null distribution of maximal component size, and the corrected p-value was calculated as the proportion of permutations for which the most connected components consists of two or more links. Methodological details for NBS can be found in Zalesky et al. (2010).
