Use of random matrix theory in the discovery of resting state brain networks

Connectomics identifies brain networks in vivo in resting state functional MRI. However, the presence of noise produces spurious identification of brain networks, which have low test-retest reliability. A Network Based Statistics approach to network identification has been previously proposed that affords much better statistical power relative to Bonferroni method but nevertheless provides a sufficiently conservative, family-wise control for false positives. We propose the use of Random Matrix Theory (RMT) to discover brain networks and to associate those networks with demographic and clinical variables. We parcellated the brain into cortical and subcortical regions using either an anatomical or a functional brain atlas. We applied RMT to study functional connectivity across brain regions by first computing the correlation matrix for time courses in those brain regions and then identifying eigenvalues that deviate from the theoretical random distribution that RMT predicts, on the assumption that real brain networks would produce eigenvalues that differ significantly from the random distribution. We assessed the specificity and test-retest reliability of identified networks through application of this RMT-based approach to (1) synthetic data generated under the null-hypothesis, (2) resting state functional MRI data from 4 real-world cohorts of patients and healthy controls, and (3) synthetic data generated by the addition of increasing amounts of noise to real-world datasets. Our findings showed that RMT method was robust to the atlas used for parcellating the brain and did not discover a brain network in synthetic data when in fact a network was not present (i.e., specificity was high); RMT-identified networks in the real-world dataset had high test-retest reliability; and RMT-based method consistently discovered the same network in the presence of increasing noise in the real-world dataset.