The degree centrality (DC) was subsequently computed, followed by the calculation of Pearson correlation coefficients between all pairs of voxels within the gray matter mask to establish a time series correlation among them. This matrix represents the total number of voxels in the whole-brain mask. The graph for each subject was constructed by creating a binary undirected adjacency matrix, with correlations thresholded at r > 0.25. This threshold excludes voxels with low temporal correlation due to signal noise and serves as the default setting for calculating the DC map. Subsequently, The weighted DC was calculated by summing up positive Pearson correlation coefficients between a given voxel and other voxels, as negative correlations were not interpreted in this study. Finally, the individual-level voxel-wise DC was transformed into a z-score map.
