{"id":"4b03a0ca-32b2-4f46-bf52-3d0643708328","arxiv_id":"2502.04550","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper defines Partial Information Rate Decomposition, a spectral lattice method that decomposes mutual information rate for stationary Gaussian processes into redundant, unique, and synergistic dynamic components.","lead":"This paper introduces a way to split the information flowing between time-varying signals into unique, shared, and synergistic parts. It gives time-series researchers a principled alternative to static information decomposition when data have temporal correlations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Deferred proof of the redundancy axioms for the spectral minimum leaves the N>2 nonnegativity claim unverified; conditional acceptance is appropriate.","rationale":"The reader identified the same weakest assumption: the spectral-pointwise minimum redundancy function must satisfy the PID axioms for all relevant atoms, not only for the N=2 Gaussian examples. I agree that this is the single most load-bearing concern. The manuscript itself flags the gap by deferring the proof to [14], and the provided validations do not exercise N>2, where the Möbius inversion becomes a nontrivial alternating sum over a larger lattice. This is not an accusation of error; the framework is internally consistent for N=2, the frequency-wise construction plausibly reduces to static PID in the i.i.d. case, and the simulations and code are useful. But the central claim that PIRD yields nonnegative unique, redundant, and synergistic information rates for arbitrary N is not self-contained in this paper. Given that the proof is deferred to a companion preprint and no machine-checked or independent N>2 verification is provided, a conditional verdict is the honest outcome. I therefore do not change the reader's verdict, while emphasizing that a concrete numerical test on N=3 Gaussian systems plus a re-derivation of the companion lemma would settle the concern either way.","tokens_in":7931,"tokens_out":31770,"duration_ms":342875,"concrete_test":"Run a numerical stress test: generate many stationary Gaussian VAR(1) systems with N=3 sources plus a target, compute the spectral MIRs via Eq. (7), form redundancy rates via Eqs. (5)-(6) for all 18 antichain nodes of the redundancy lattice, apply Möbius inversion via Eq. (4), and require every atom to be nonnegative and Eqs. (2)-(3) to hold to numerical precision. Simultaneously, independently re-derive from [14] the lemma that the frequency-wise minimum satisfies the lattice axioms for arbitrary N; if any random draw produces a negative atom or if the lemma fails for some N, the PIRD lattice is invalid. If 10^4 random draws all pass, the deferred claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (4), via Möbius inversion of the redundancy rate defined in Eqs. (5)-(6), yields a valid PIRD with nonnegative atoms for any number of sources. The paper asserts after Eq. (7) that the frequency-wise minimum redundancy 'preserves nonnegativity, satisfies the main axioms' and defers the proof to the companion paper [14]. This is the load-bearing step: for N=2, nonnegativity follows from pointwise inequalities such as i_{X1,X2;Y}(ω) ≥ max{i_{X1;Y}(ω), i_{X2;Y}(ω)}, but for N>2 the redundancy lattice has many more nodes and the Möbius transform is an alternating sum. Monotonicity on a non-chain lattice does not by itself guarantee nonnegative Möbius atoms; a monotone function on the Boolean lattice can have a negative Möbius atom. Thus, if the deferred proof fails, Eq. (4) can produce negative or semantically meaningless atoms. The paper provides no N>2 validation: the simulations use N=2 and the climate application is restricted to triplets, so it never exercises the regime where the lattice structure is nontrivial. The central generalization to N>2 therefore rests entirely on an unstated proof in a companion preprint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes Partial Information Rate Decomposition (PIRD), a dynamic extension of Partial Information Decomposition for stationary random processes. It defines the information shared between a scalar target process Y and N source processes through a redundancy rate constructed as the frequency-wise minimum of spectral mutual information rates (Eqs. 5-6), and then obtains unique, redundant, and synergistic information-rate atoms by Möbius inversion over the PID redundancy lattice (Eqs. 2-4). For Gaussian processes, spectral mutual information rates are computed from power spectral densities via Eq. (7). The framework is illustrated on a two-source Gaussian network and on six triplets of climate indices, and it is claimed to reduce to static PID for i.i.d. processes and to transfer-entropy-based PID in strictly causal cases.","tokens_in":8175,"tokens_out":15612,"duration_ms":173637,"significance":"If the construction is valid, PIRD is a useful and computationally tractable dynamic counterpart to PID for stationary Gaussian processes, and the frequency-domain formulation offers the attractive possibility of band-limited analyses. The paper has concrete strengths: a closed-form spectral implementation, explicit reduction to known static and causal special cases, an openly available codebase, and a real-data demonstration on climate indices. The principal weakness is that the axiomatic validity of the spectral-minimum redundancy rate—especially nonnegativity of all Möbius atoms for N>2—is asserted and deferred to a companion paper rather than proved or demonstrated in the present manuscript.","major_comments":[{"comment":"The claim that the frequency-wise minimum in Eq. (6) defines a redundancy rate that is nonnegative and satisfies symmetry, self-redundancy, and monotonicity is the load-bearing premise of the entire construction, because Eq. (4) then produces all atoms by Möbius inversion. This property is not proved in the manuscript; the text says only that the rate is shown in [14] to be bounded by the minimum time-domain MIR and that it preserves nonnegativity and satisfies the main axioms. For N=2 the needed nonnegativity follows from pointwise inequalities, but for N>2 the redundancy lattice is not a chain, and monotonicity of a lattice function does not by itself guarantee nonnegative Möbius atoms—on the Boolean lattice, for example, the monotone function f(S)=1_{|S|>=2} has a negative atom at the full set. Thus the central N>2 generalization is currently unsupported by the present text. Please include a self-contained proof, or state the exact theorem with full hypotheses and a proof outline, so that the Möbius inversion in Eq. (4) is justified for all lattice nodes.","section":"Definition of redundancy rate, Eqs. (5)-(6) and following paragraph"},{"comment":"All numerical evidence in the paper uses N=2 sources: the simulated network has two source processes, and the climate analysis is restricted to triplets (a target plus two sources). The nontrivial region of the redundancy lattice, in which antichains with more than two entries appear and the Möbius inversion is genuinely alternating, is therefore never exercised. I request at least one N=3 or N=4 Gaussian example with known model parameters, reporting all atoms and checking both nonnegativity and the marginal constraints in Eqs. (2)-(3). This is the minimum validation needed to make the claimed generalization beyond N=2 credible.","section":"Simulation and application sections"},{"comment":"The statement that the Gaussian restriction does not invalidate the analysis of nonlinear systems because stationary random processes have a linear (albeit infinite-order) representation is too strong. The Wold representation is linear in the innovations, but the innovations need not be Gaussian, and Eq. (7) computes a spectral mutual information rate that is exactly the MIR only under joint Gaussianity. For non-Gaussian stationary processes, the power spectral density does not determine the mutual information rate. Please restrict the claim to Gaussian processes or provide a separate argument explaining why the spectral MIR formula remains valid beyond Gaussianity.","section":"Summary/limitations paragraph"}],"minor_comments":[{"comment":"The verbal description of A as the collection of subsets of sources such that no subset is a superset of any other is imprecise; what is meant is the lattice of antichains of nonempty subsets of {1,...,N}. Please correct the definition and use notation that distinguishes a source subset such as {i,j} from an antichain such as {{1},{2}}.","section":"Notation in Eqs. (2)-(4)"},{"comment":"Properties that are central to the method, including the inequality I^∩_{Xα;Y} <= min_j I_{Xαj;Y}, the frequency-band decomposition, and the proof of the redundancy axioms, are repeatedly attributed to [14]. Please state explicitly which claims are proven in the present manuscript and which are proven only in the companion paper.","section":"Companion paper [14]"},{"comment":"The notation R(Y(tn); X(tn)) and S(Y(tn); X(tn)) is easy to confuse with the source processes X1 and X2; consider denoting the static redundant and synergistic components with subscripts, for example R_static and S_static, to avoid ambiguity.","section":"Fig. 1 caption"},{"comment":"The statement that the importance of NTA is elicited only using PIRD rests on point estimates; no confidence intervals or statistical tests are reported for the PID values, so this claim should be softened or accompanied by interval estimates for both PID and PIRD.","section":"Climate application"},{"comment":"There is a typo in the ENSO application section: 'athmospheric' should be 'atmospheric'. In addition, Eq. (1) should state the stationarity and regularity assumptions under which the limit defining the mutual information rate exists.","section":"Typo and assumptions"}],"recommendation":"major_revision","confidential_remarks":"The key issue is that the mathematical core of the paper—validity of the spectral-minimum redundancy measure on the full PID lattice, and hence nonnegativity of all Möbius atoms for N>2—is deferred to companion arXiv:2502.04555. If that companion is not part of this submission, the present paper is not self-contained on its central claim. I do not think rejection is warranted: the idea is interesting, the N=2 case is internally consistent, and the computational framework is useful. The authors should be asked to provide a theorem with proof sketch, or at least an explicit statement of exactly which axioms hold and why, and to include an N>2 numerical example. The code availability statement is helpful but should point to an archived, versioned repository rather than a personal webpage."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"PIRD is the first dynamic PID for Gaussian processes that I'd actually consider using. The idea is simple and correct in outline: translate the Williams-Beer redundancy lattice from random variables to stationary processes by defining redundancy at the level of spectral MIR, taking the pointwise minimum across frequencies. This gives a decomposition that reduces to static PID in the i.i.d. case and to joint transfer entropy PID in the strictly causal case. The authors show with a VAR model that zero-lag PID can be badly misled by temporal correlations, and the climate example is a reasonable demonstration.\n\nWhat's genuinely new is the spectral minimum redundancy rate; the lattice translation itself is straightforward once you have that. The paper is honest about deferring the axiom proof to a companion preprint, and it does ship code and a reproducible simulation setup. That counts.\n\nThe soft spots are real but not fatal. The biggest is that the axioms of the redundancy rate -- symmetry, self-redundancy, monotonicity, and the implied nonnegativity of Möbius atoms -- are asserted, not proven here. For N=2 you can verify nonnegativity by hand, but the paper doesn't. For N>2, a monotone function on the lattice is not enough; Möbius inversion of a monotone but not modular function can go negative. So the claimed generalization to N>2 is unsupported in this manuscript. The simulations and the climate application only use two sources, so they never exercise the nontrivial lattice. This is the load-bearing gap, and the stress test is right about it.\n\nThe circularity concern is overblown. The decomposition follows by definition from the redundancy function; that's how all PID works. The validation is qualitative -- expected trends in a designed VAR -- but that's typical for methods papers and not a flaw.\n\nThe surrogate analysis in the climate section is appropriate and partially addresses error-bar concerns; I wouldn't push on that. The missing N>2 validation and the proof discrepancy are the honest reasons for a conditional rather than clean verdict.\n\nWho's it for: anyone doing PID on time series in neuroscience, physiology, or climate. It deserves a serious referee. I'd accept it for review and ask the authors to either include the axiom proof or make the companion paper part of the submission, and to add one N>2 example. If the companion proof is watertight, this is a solid methods contribution.","headline":"A genuine spectral PID for dynamic Gaussian systems, but the N>2 generalization rests on a deferred proof; worth refereeing.","tokens_in":8691,"tokens_out":2995,"would_cite":true,"duration_ms":31784,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["62M10","62B10","94A17"],"pacs":[],"model":"deepseek-v4-flash","headline":"PIRD is a dynamic extension of partial information decomposition that splits the mutual information rate of Gaussian processes into unique, redundant, and synergistic parts via a frequency-wise minimum redundancy rule.","keywords":["partial information decomposition","mutual information rate","spectral redundancy","synergy","redundancy","Gaussian processes","vector autoregressive models","climate time series"],"falsifier":"Fit a three-source stationary Gaussian vector autoregressive model in which two sources carry the same message to the target in two disjoint frequency bands and the third source is a frequency-dependent linear mixture of the other two; compute the redundancy rate of Eq. (6) for every lattice atom, apply the Möbius inversion of Eq. (4), and check whether any atom is negative. A single negative atom, or any violation of monotonicity of $I^\\cap_{X_\\alpha;Y}$ under adding sources, would show that the spectral minimum redundancy function is not valid for $N>2$.","tokens_in":7746,"feed_emoji":"🧩","tokens_out":7429,"duration_ms":67000,"temperature":0.7,"pith_summary":"The paper introduces Partial Information Rate Decomposition (PIRD), a framework that decomposes the mutual information rate between a target random process and a set of source processes into unique, redundant, and synergistic information rates. It argues that the standard static partial information decomposition (PID) misrepresents high-order dynamic interactions whenever the analyzed time series carry temporal correlations, because it treats each time point as an independent observation. PIRD instead applies the PID lattice to information rates and defines redundancy as the frequency-wise minimum of spectral mutual information rates, so that Möbius inversion recovers the individual atoms. For stationary Gaussian processes the decomposition is computable from the joint power spectral density, and it reduces to static PID in the i.i.d. case and to a PID of joint transfer entropy in strictly causal cases. A sympathetic reader should care because PIRD offers a principled dynamic replacement for a method that is widely applied, often incorrectly, to time series from networks.","feed_headline":"Dynamic information rates split into redundancy and synergy","feed_subtitle":"A frequency-wise minimum extends PID to Gaussian time series, where static PID misreads dynamic networks.","key_machinery":"The load-bearing object is the spectral redundancy rate $i^\\cap_{X_\\alpha;Y}(\\omega)=\\min_{j=1,\\dots,J} i_{X_{\\alpha_j};Y}(\\omega)$, the pointwise minimum at each frequency of the spectral mutual information rates between the target and each source in the atom. Integrating this quantity over the normalized angular frequency yields the redundancy rate $I^\\cap_{X_\\alpha;Y}$, and Möbius inversion of Eq. (4) on the PID lattice turns it into the unique, redundant, and synergistic information-rate atoms. The spectral mutual information rate itself is computed from the power spectral density of the joint Gaussian process, which is what makes the framework operational for vector autoregressive models and linear Gaussian systems.","core_discovery":"On its own terms, the paper claims that for a jointly stationary Gaussian process $S=\\{X_1,\\dots,X_N,Y\\}$, the mutual information rate $I_{X;Y}$ can be written as a sum over the PID lattice of information-rate atoms, with each atom obtained by Möbius inversion of a redundancy rate function $I^\\cap_{X_\\alpha;Y}$. The redundancy rate is defined spectrally as the integral over frequency of the pointwise minimum, across the sources in an atom, of the spectral mutual information rates $i_{X_{\\alpha_j};Y}(\\omega)=\\frac{1}{2}\\log\\frac{|P_{X_{\\alpha_j}}(\\omega)|P_Y(\\omega)}{|P_{[X_{\\alpha_j}Y]}(\\omega)|}$. In simulations of a two-source network with only time-lagged effects, the paper shows that static PID reports the lag-zero information as predominantly redundant, while PIRD reveals increasing total information rate and expected net synergy as common-target coupling grows. In the climate application with SOI as target, PIRD brings out redundant and synergistic contributions of the North Tropical Atlantic index that static PID makes negligible, and surrogate data show that temporal correlations cannot be ignored.","pith_inferences":["Beyond the paper, the frequency-wise minimum suggests a band-limited redundancy measure, so one could ask whether redundant information is carried by slow or fast oscillations in a given network.","Beyond the paper, a natural stress test is to apply PIRD to nonlinear stationary processes such as coupled chaotic maps, where the Gaussian spectral formula is an approximation, and check whether all atoms remain nonnegative.","Beyond the paper, the conservative spectral minimum implies that when two sources carry identical information at disjoint frequencies, PIRD will label that information unique rather than redundant, a behavior worth explicit study in frequency-multiplexed systems.","Beyond the paper, one could define a significance test for individual PIRD atoms using block-permutation surrogates that preserve the spectral density, which would make the decomposition hypothesis-testable in practice."],"forward_implications":["Static PID applied to zero-lag variables can mischaracterize temporally correlated dynamics as redundant even when the true dynamic interaction is synergistic; PIRD includes the full temporal structure and corrects this.","In strictly causal systems with no target-to-source feedback, PIRD reduces to the PID of the joint transfer entropy from all sources to the target, linking dynamic information decomposition to established causality measures.","Because redundancy is defined frequency by frequency, PIRD can be restricted to predefined oscillation bands, enabling band-specific analysis of oscillatory networks.","For stationary Gaussian processes the decomposition is estimated from VAR parameters, so it is computationally light and directly applicable to real multivariate time series.","Surrogate tests in the paper show PIRD values can be either larger or smaller than static PID, so time-lagged and instantaneous effects cannot be separated by a simple adjustment."],"supporting_citations":[{"why":"Defines the PID lattice and consistency equations that PIRD adapts from random variables to random processes.","marker":"[1]"},{"why":"Supplies the minimum mutual information principle used to define the pointwise redundancy rate.","marker":"[15]"},{"why":"Provides the spectral mutual information rate formula for Gaussian processes used in Eq. (7).","marker":"[21]"},{"why":"Companion paper that proves the redundancy rate satisfies the PID axioms and derives the inequality against time-domain minimum redundancy; the present paper defers key properties to it.","marker":"[14]"},{"why":"Defines the mutual information rate that is the quantity PIRD decomposes.","marker":"[13]"},{"why":"Connects common-target versus common-drive network structures to dominance of synergy versus redundancy, which the simulations use to interpret PIRD outputs.","marker":"[22]"}],"fun_headline_variants":["PIRD: Decomposing information rates in time series","Spectral PID uncovers synergy static analysis misses","Frequency-wise information decomposition for dynamic networks","Time-aware PID: Redundancy and synergy beyond static snapshots"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes the spectral pointwise minimum in Eq. (6) is a valid redundancy measure across the whole lattice: nonnegative, symmetric, equal to the source's own information rate for a single source, and decreasing when sources are added; the proof is deferred to the companion paper, and if the assumption fails the decomposed atoms can become negative or meaningless.","fun_headline_variants_meta":{"raw":{"variants":["PIRD: Decomposing information rates in time series","Spectral PID uncovers synergy static analysis misses","Frequency-wise information decomposition for dynamic networks","Time-aware PID: Redundancy and synergy beyond static snapshots"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000443,"raw_usage":{"total_tokens":2219,"prompt_tokens":899,"completion_tokens":1320,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":515,"completion_tokens_details":{"reasoning_tokens":1257}},"tokens_in":515,"tokens_out":1320,"duration_ms":12107,"temperature":1.0,"reasoning_tokens":1257,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T22:22:17.089149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit a three-source stationary Gaussian vector autoregressive model in which two sources carry the same message to the target in two disjoint frequency bands and the third source is a frequency-dependent linear mixture of the other two; compute the redundancy rate of Eq. (6) for every lattice atom, apply the Möbius inversion of Eq. (4), and check whether any atom is negative. A single negative atom, or any violation of monotonicity of $I^\\cap_{X_\\alpha;Y}$ under adding sources, would show that the spectral minimum redundancy function is not valid for $N>2$.","supporting_citations":[{"cited_title":"HONEST - High-Order Dynamical Networks in Computational NeuroscienceandPhysiology: anInformation-Theoretic Framework","cited_arxiv_id":null,"evidence_quote":"Defines the PID lattice and consistency equations that PIRD adapts from random variables to random processes."},{"cited_title":"Decomposing Multivariate Information Rates in Networks of Random Processes","cited_arxiv_id":"2502.04555","evidence_quote":"Supplies the minimum mutual information principle used to define the pointwise redundancy rate."},{"cited_title":"Sparacino, Y","cited_arxiv_id":null,"evidence_quote":"Provides the spectral mutual information rate formula for Gaussian processes used in Eq. (7)."}],"review_version":1}