This paper proposes a scalable method for identifying interactions in higher-order networks from observations of nodal processes. Finding such dependencies is important in many disciplines, including neuroscience, social influence modeling, and beyond. However, current approaches are either limited to extracting pairwise dependencies or struggle with scalability, as estimating higher-order dependencies becomes computationally prohibitive. To overcome these challenges, we introduce a tensorbased graph Volterra model that leverages low-rank decomposition techniques to estimate higher-order interactions efficiently. Our approach not only reduces computational and storage complexity but also acts as an implicit regularizer, improving network estimation in ill-posed settings. We validate our method through simulations and real data experiments, demonstrating competitive performance and enhanced scalability compared to existing techniques.

Identifying overlapping communities from data is crucial for grasping the complex structure and dynamics of networks, amongst others in fields such as computational neuroscience. Research using fMRI has demonstrated that brain regions can change their functional network membership over time using temporal independent component analysis (tICA). However, reproducibility of such overlapping communities remains a challenge. Recently, several alternative approaches have been proposed to identify such overlapping communities. While results are promising, less is known about the model and assumptions that underlie these approaches. This paper shows that the bilinear model, combined with the assumption of quasi-stationary and uncorrelated sources, underlies novel methods for identifying overlapping brain networks. Furthermore, we propose a new algorithm, and through simulations, we investigate the robustness of our algorithm and several existing methods to solve the problem in noisy conditions with few available data samples. We conclude that quasi-stationary blind source separation-based techniques can have a promising advantage over tICA in terms of identifiability of overlapping communities and thus have the potential to contribute towards greater reproducibility of results.

Functional ultrasound (fUS) is a novel neuroimaging technique that measures brain hemodynamics through a time series of Doppler images. The measured spatiotemporal hemodynamic changes reflect changes in neural activity through the neurovascular coupling (NVC). Often, such image time series is used to analyze dynamic functional connectivity (dFC) by directly computing a connectivity metric between the measured hemodynamic signals, ignoring the functional connectomics of underlying neural populations. This work proposes a novel fUS signal model, consisting of a hidden Markov model (HMM) cascaded with a convolutive model, that captures how fUS signals arise from a generative perspective while incorporating high-level biological functioning of neural populations. Consequently, the developed model enables inference of functional connectivity networks, being co-activation patterns (CAPs) of neural populations. Our results show that our methods can identify biologically plausible networks of functional connectivity. Furthermore, this method captures a difference in brain dynamics between wild-type and ${Shank2}^{-/-}$ mouse mutants.