We consider exploratory methods for the discovery of cortical functional
connectivity. Typically, data for the i-th subject (i=1...NS) is represented as
an NVxNT matrix Xi, corresponding to brain activity sampled at NT moments in
time from NV cortical voxels. A widely used method of analysis first
concatenates all subjects along the temporal dimension, and then performs an
independent component analysis (ICA) for estimating the common cortical
patterns of functional connectivity. There exist many other interesting
variations of this technique, as reviewed in [Calhoun et al.
Understanding of normal and pathological brain function requires the
identification and localization of functional connections between specialized
regions. The availability of high time resolution signals of electric neuronal
activity at several regions offers information for quantifying the connections
in terms of information flow. When the signals cover the whole cortex, the
number of connections is very large, making visualization and interpretation
very difficult.