Multiple‐region directed functional connectivity based on phase delays

Network analysis is increasingly advancing the field of neuroimaging. Neural networks are generally constructed from pairwise interactions with an assumption of linear relations between them. Here, a high‐order statistical framework to calculate directed functional connectivity among multiple regions, using wavelet analysis and spectral coherence has been presented. The mathematical expression for 4 regions was derived and used to characterize a quartet of regions as a linear, combined (nonlinear), or disconnected network. Phase delays between regions were used to obtain network's temporal hierarchy and directionality. The validity of the mathematical derivation along with the effects of coupling strength and noise on its outcomes were studied by computer simulations of the Kuramoto model. The simulations demonstrated correct directionality for a large range of coupling strength and low sensitivity to Gaussian noise compared with pairwise coherences. The analysis was applied to resting‐state fMRI data of 40 healthy young subjects to characterize the ventral visual system, motor system and default mode network (DMN). It was shown that the ventral visual system was predominantly composed of linear networks while the motor system and the DMN were composed of combined (nonlinear) networks. The ventral visual system exhibits its known temporal hierarchy, the motor system exhibits center ↔ out hierarchy and the DMN has dorsal ↔ ventral and anterior ↔ posterior organizations. The analysis can be applied in different disciplines such as seismology, or economy and in a variety of brain data including stimulus‐driven fMRI, electrophysiology, EEG, and MEG, thus open new horizons in brain research. Hum Brain Mapp 38:1374–1386, 2016. © 2016 Wiley Periodicals, Inc.