Edgewise and subgraph-level tests for brain networks.

Resting-state functional magnetic resonance image is a useful technique for investigating brain functional connectivity at rest. In this work, we develop flexible regression models and methods for determining differences in resting-state functional connectivity as a function of age, gender, drug intervention, or neuropsychiatric disorders. We propose two complementary methods for identifying changes of edges and subgraphs. (i) For detecting changes of edges, we select the optimal model at each edge and then conduct contrast tests to identify the effects of the important variables while controlling the familywise error rate. (ii) We adopt the network-based statistics method to improve power by incorporating the graph topological structure. Both methods have wide applications for low signal-to-noise ratio data. We propose stability criteria for the choice of threshold in the network-based statistics procedure and utilize efficient massive parallel procedure to speed up the estimation and inference procedure. Results from our simulation studies show that the thresholds chosen by the proposed stability criteria outperform the Bonferroni threshold. To demonstrate applicability, we use both methods in the context of the Oxytocin and Aging Study to determine effects of age, gender, and drug treatment on resting-state functional connectivity, as well as in the context of the Autism Brain Imaging Data Exchange Study to determine effects of autism spectrum disorder on functional connectivity at rest.