THOI: An efficient and accessible library for computing higher-order interactions enhanced by batch-processing

TL;DR

THOI is a Python library that efficiently computes higher-order interactions in complex systems, leveraging Gaussian copula and batch processing to improve speed, scalability, and applicability to large datasets.

Contribution

We introduce THOI, a novel library that enhances the estimation of high-order interactions using Gaussian copula and parallel processing, addressing computational challenges in complex system analysis.

Findings

THOI outperforms existing tools in speed and scalability.

It accurately estimates high-order interactions in synthetic and real datasets.

THOI enables analysis of large-scale systems like fMRI data.

Abstract

Complex systems are characterized by nonlinear dynamics, multi-level interactions, and emergent collective behaviors. Traditional analyses that focus solely on pairwise interactions often oversimplify these systems, neglecting the higher-order interactions critical for understanding their full collective dynamics. Recent advances in multivariate information theory provide a principled framework for quantifying these higher-order interactions, capturing key properties such as redundancy, synergy, shared randomness, and collective constraints. However, two major challenges persist: accurately estimating joint entropies and addressing the combinatorial explosion of interacting terms. To overcome these challenges, we introduce THOI (Torch-based High-Order Interactions), a novel, accessible, and efficient Python library for computing high-order interactions in continuous-valued systems. THOI leverages the well-established Gaussian copula method for joint entropy estimation, combined with state-of-the-art batch and parallel processing techniques to optimize performance across CPU, GPU, and TPU environments. Our results demonstrate that THOI significantly outperforms existing tools in terms of speed and scalability. For larger systems, where exhaustive analysis is computationally impractical, THOI integrates optimization strategies that make higher-order interaction analysis feasible. We validate THOI accuracy using synthetic datasets with parametrically controlled interactions and further illustrate its utility by analyzing fMRI data from human subjects in wakeful resting states and under deep anesthesia. Finally, we analyzed over 900 real-world and synthetic datasets, establishing a comprehensive framework for applying higher-order interaction (HOI) analysis in complex systems.