Extracting Phase Coupling Functions between Collectively Oscillating Networks from Time-Series Data

Abstract

Many real-world systems are often regarded as weakly coupled limit-cycle oscillators, in which each oscillator corresponds to a dynamical system with many degrees of freedom that have collective oscillations. One of the most practical methods for investigating the synchronization properties of such a rhythmic system is to statistically extract phase coupling functions between limit-cycle oscillators directly from observed time-series data. In Particular, using a method that combines phase reduction theory and Bayesian inference, the phase coupling functions can be extracted from the time-series data of even just one variable in each oscillatory dynamical system with many degrees of freedom. However, it remains unclear how the choice of the observed variables affects the statistical inference for the phase coupling functions. In this study, we examine the influence of observed variable types on the extraction of phase coupling functions using some typical dynamical elements under various conditions. We demonstrate that our method can consistently extract the macroscopic phase coupling functions between two phases representing collective oscillations in a fully locked state, regardless of the observed variable types; for example, even using one variable of any element in one system and the mean-field value over all the elements in another system. We also study the case of globally coupled phase oscillators in a partially locked state. Our results reveal directional asymmetry in the robustness of extracting the macroscopic phase coupling function between two networks. For instance, when an asynchronous oscillator in network $A$ and the macroscopic collective oscillation of network $B$ is observed, the macroscopic phase coupling function from network $A$ to network $B$ can be extracted more robustly than in the opposite direction.