Spectral learning of predictive dynamics AI

How past–future relationships become a compact dynamical representation and spatial network patterns.
Published

October 21, 2026

A predictive representation describes the recent history of a signal by what it helps predict next. In fMRI, this changes the modeling target: instead of finding only the largest fluctuations, we seek a compact set of temporal relationships linking regional activity across successive windows.

What the decomposition learns

Let ptp_t stack a short history of ROI measurements and ftf_t stack the following measurements. Their empirical cross-moment, Cfp=⟨ftpt𝖳⟩C_{fp}=\langle f_t p_t^\mathsf{T}\rangle, summarizes relationships between those windows. The past moment Cpp=⟨ptpt𝖳⟩C_{pp}=\langle p_t p_t^\mathsf{T}\rangle describes the scale and redundancy of the histories. These are uncentered moments unless the inputs have been centered.

In the implementation used here, a regularized whitening of the past precedes a singular value decomposition of the cross-moment. Keeping kk singular directions defines a reduced predictive space. This resembles PCA in its use of an SVD, but the decomposed matrix relates past and future rather than simply measuring variance within one observation vector. The retained singular values therefore are not percentages of BOLD variance explained.

History and future windows feed a whitened cross-moment decomposition; retained coordinates support spatial and temporal readouts.

A schematic of the continuous-signal implementation. Rank sets the dimension of the predictive representation.

Where networks enter

A history is mapped into coordinates xk(t)x_k(t). Regression then associates these coordinates with observation moments and their evolution. For the spatial readout,

M(t)=W0+∑k=1kmaxxk(t)Wk.M(t)=W_0+\sum_{k=1}^{k_{\max}}x_k(t)W_k.

Each WkW_k is an ROI-by-ROI coefficient pattern, and several patterns contribute simultaneously. Opposing signed weights are part of this linear representation: a negative coordinate reverses the contribution of its entire pattern. They do not independently establish physiological inhibition or anticorrelation.

Predictive-state methods provide a broader framework for representing dynamics through observable predictions; see Boots, Siddiqi and Gordon. The continuous moment model here is an adaptation, so probability guarantees from a particular discrete-state formulation cannot simply be transferred to it. Its predictive value requires evaluation on recordings withheld from fitting.