{"id":"dbbe9c72-6e97-4376-a146-35e8481615a1","arxiv_id":"2411.14602","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"An automated, KS-statistic-based way to choose modified diffusion entropy analysis parameters is proposed and used to show complexity synchronization among EEG, ECG, and filtered respiration during a cognitive task.","lead":"This methods paper introduces automated parameter selection for modified diffusion entropy analysis, a technique that measures complexity from brain, heart, and lung signals. It applies the approach to cognitive task data and reports that complexity indices of these signals often move together, especially between brain and heart, after preprocessing choices that are still under debate.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The KS-based stripe-size selection fits an IPL model by construction, so without null-surrogate validation the reported complexity synchronization may reflect parameter-optimization artifacts rather than genuine coupling.","rationale":"The reader's weakest assumption focuses on the post-hoc RESP filter creating an IPL-like artifact. That is a valid and important concern, and the paper's own cautionary text supports it. However, I see a broader and more load-bearing issue: the entire parameter-selection procedure is designed to find an IPL fit and a linear scaling region, so the production of a δ value is guaranteed for almost any input. The paper demonstrates sensitivity on Mittag-Leffler data but never demonstrates specificity on null data. Without a null-surrogate test, the central claim that the observed correlations reflect genuine complexity synchronization is not fully supported. This is not an accusation of circular reasoning in a pejorative sense; it is a standard model-validation gap. The proposed surrogate test would settle whether the method creates spurious CS. If the false-positive rate is low, the method gains real support and the filter concern is largely resolved. If the false-positive rate is high, the headline empirical claim is undermined regardless of whether the RESP filter is physiologically justified. I therefore keep the reader's CONDITIONAL verdict: the paper's contribution is plausible and code is provided, but acceptance should be conditional on null-surrogate validation and, ideally, on a pre-specified filtering decision. My agreement with the reader is partial because I identify the missing negative control as the deeper issue, with the RESP filter as a special case.","tokens_in":15155,"tokens_out":5475,"duration_ms":63385,"concrete_test":"Run the full MDEA+KS pipeline on surrogate data with no true cross-signal coupling: (i) phase-randomized surrogates of the actual EEG, ECG, and filtered RESP that preserve power spectra but destroy nonlinear structure, and (ii) independent realizations of a known non-IPL process (e.g., AR(1) or sine-plus-noise) with the same length and sampling rate. Apply the identical stripe-size search, linear-fit selection, and 30 s overlapping-window correlation analysis, then count how often all three pairwise correlations reach p < .05. If the false-positive rate substantially exceeds the nominal 5%, or if phase-randomized surrogates yield a three-way CS rate close to the observed 24/54, the reported synchronization is an artifact of the estimation pipeline.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the automated MDEA pipeline detects synchronization of multifractal scaling indices across EEG, ECG, and filtered RESP. That claim requires the estimated δ series to track genuine complexity, not just best-fit parameters. But Eq. (5) selects stripe size and μ by minimizing the KS statistic against an IPL CCDF, and the adaptive linear-region selection in §2.3.2 fits the longest linear segment of the entropy plot. Together these steps will produce a scaling exponent for almost any input signal, because there is no goodness-of-fit threshold and no rejection of non-IPL cases. The paper validates the estimator on Mittag-Leffler synthetic data (a positive control), but provides no negative controls: no demonstration that uncoupled or non-IPL signals yield low or absent CS. The post-hoc 2 Hz high-pass filtering of RESP (§2.3.1) is a concrete instance of this broader issue: the filter changes the stripe-size distribution and the scaling estimates, and the authors themselves caution that removing the dominant respiratory oscillation raises interpretability questions. The reader's filter-artifact concern is therefore best understood as one manifestation of a more general missing piece: the pipeline's specificity has not been established. If the KS optimization can make phase-randomized noise or independent AR(1) processes look IPL-like, then the reported correlations among EEG, ECG, and filtered RESP, especially the 24/54 three-way significant datasets, may be artifacts of the estimation procedure rather than evidence of brain-heart-lung complexity synchronization.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an automated pipeline for modified diffusion entropy analysis (MDEA) of neurophysiological time series, with the goal of estimating multifractal scaling indices δ and detecting complexity synchronization (CS) across EEG, ECG, and respiration. Two parameter-selection problems are addressed: automated stripe-size selection via a Kolmogorov–Smirnov (KS) fit to an inverse power law (IPL) complementary cumulative distribution function (Eq. 5), and automated linear-fit-region selection using Matlab's ischange() function. The method is validated on Mittag-Leffler synthetic data and then applied to 54 datasets (27 subjects under low and high time stress) from a Go-NoGo neurofeedback study. The empirical results indicate significant EEG-ECG CS in 46/54 datasets and significant three-way CS with high-pass filtered RESP in 24/54 datasets. The authors provide open-source code and a sample dataset, and they candidly discuss limitations, including the post-hoc decision to filter RESP and unresolved interpretability questions.","tokens_in":15504,"tokens_out":4156,"duration_ms":40988,"significance":"If the central claim is upheld, the paper would make a useful methodological contribution to network physiology by offering a repeatable, automated procedure for measuring complexity synchronization across heterogeneous organ time series. The open-source implementation, Monte Carlo bias/variance evaluation, and honest reporting of limitations are clear strengths. However, the main empirical headline—that brain, heart, and lung complexity 'significantly co-vary'—is only supported in 24 of 54 datasets, and the pipeline's specificity is not established: there are no negative controls or goodness-of-fit rejection criteria, and the post-hoc RESP filter is a potential source of artifact. These issues do not invalidate the method but they do block acceptance as a validated methods paper.","major_comments":[{"comment":"The abstract states that 'the complexity of brain, heart, and lung ONTS significantly co-vary over time during cognitive task performance,' but the three-way CS (EEG-ECG-RESP) was significant in only 24 of 54 datasets, with 30 datasets lacking significance. Under a per-dataset threshold of p<.05, about 2.7 false-positive datasets are expected by chance, so the aggregate result is above chance, but the headline overstates the consistency of the effect. The authors should report effect sizes and confidence intervals for the correlations and ideally perform a multilevel analysis with subjects and conditions as random effects, rather than relying on a count of 'significant' datasets.","section":"Results, after Table 4"},{"comment":"The decision to high-pass filter RESP at 2 Hz with a Kaiser filter of order 8192 was made only after finding that unfiltered RESP did not synchronize with EEG and ECG. The authors themselves raise the interpretability concern ('what does it mean to remove the dominant feature of the RESP time series?'). The paper provides no control analysis to show that this filter preserves an IPL complexity component rather than creating an IPL-like artifact. I request a surrogate validation: apply the identical filter, stripe-size selection, and MDEA pipeline to phase-randomized RESP, to simulated mixtures of periodic and IPL components, and to pure noise, and show that the resulting EEG-RESP and ECG-RESP correlations are not reproduced without genuine coupling.","section":"Section 2.3.1, following Figure 7"},{"comment":"The KS-based stripe-size and μ selection minimizes the KS statistic against an IPL CCDF with two free parameters, and MDEA subsequently estimates δ which is interpreted through the same IPL relation μ = 1 + 1/δ. This couples the measurement to the assumed model: for almost any signal, the grid search will find a stripe size that brings the inter-event intervals close to some IPL curve, unless a goodness-of-fit threshold is enforced. The paper demonstrates positive controls (Mittag-Leffler, Figs. 4-5) but no negative controls. I ask the authors to report the KS statistics and the proportion of windows that would fail an IPL goodness-of-fit test, and to apply the pipeline to non-IPL signals (e.g., AR(1), sinusoids, white noise) to demonstrate that spurious scaling indices and spurious CS are not produced.","section":"Section 2.3.1, Eqs. (3)-(5)"},{"comment":"The use of 30-s sliding windows with 20-s overlap creates strong serial dependence in the δ time series, so the effective sample size for the correlation tests in Tables 3-4 is much smaller than the nominal number of windows. The reported p-values do not appear to account for this dependence and are therefore anti-conservative. The authors should either use non-overlapping windows or correct for autocorrelation, for example with block bootstrap or an effective-degrees-of-freedom adjustment, before claiming statistical significance for the CS correlations.","section":"Section 2.3.1 and Results"}],"minor_comments":[{"comment":"There is a typo in the text before Eq. (4): 'pr oviding' should be 'providing', and the notation is inconsistent because the definition of F_emp uses the indicator 1_{τi ≤ T} while the text describes the probability that 'the inter-arrival time variable τ < T'; please clarify the roles of τ and T.","section":"Section 2.3.1, Eq. (4)"},{"comment":"The reference to Corder and Foreman (2014) lists 'The city' as the publisher location; this placeholder should be replaced with the actual location.","section":"References"},{"comment":"The caption says 'mean (blue) and variance (red)', but the text in Section 2.3.1 refers to bias and variance; please clarify which curve is bias and which is variance, and add axis labels to the figure.","section":"Figure 4"},{"comment":"The supplementary tables are numbered Table 1-6, which duplicates the numbering of the main-text tables; please label them 'Supplementary Table S1' etc. to avoid confusion.","section":"Supplementary Materials"},{"comment":"The description of the adaptive adjustment of the ischange threshold ('adjusted in small increments') does not specify the increment size or the stopping criterion, which reduces reproducibility; please provide the exact algorithm or pseudocode for this procedure.","section":"Section 2.3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a methods-oriented contribution with a plausible core but the empirical validation is incomplete. The circularity concern raised in the reader's report is real: the KS-based stripe selection fits to the IPL model by construction, so the 'discovery' of IPL complexity is partly baked into the estimator. The missing negative controls and the post-hoc filtering of RESP are fixable with additional analyses within the manuscript's scope. I also note that the paper's theory sections largely restate prior work by the same group, so the novel contribution rests on the automated parameter selection and its application; the empirical support for the headline claim is weaker than the abstract suggests. I would be willing to look at a revised version that addresses the four major comments, especially the surrogate/specificity testing and the statistical correction for overlapping windows."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something real: it automates two parameter choices in MDEA that were previously manual — the stripe size via a KS fit to an IPL CCDF, and the linear fit region via Matlab's ischange. That is a genuine step toward reproducibility, and they shipped code and a sample dataset, which deserves credit. The synthetic Mittag-Leffler validation of the stripe-size estimator is also a sensible positive control, and the authors are candid about the RESP filtering dilemma, even asking what it means to remove the main oscillatory component.\n\nThe soft spots are real and mostly concern the empirical claim. First, there are no negative controls: no demonstration that phase-randomized noise, AR(1) processes, or explicitly uncoupled signals produce low CS. Because the KS optimization will fit an IPL-shaped CCDF to almost anything, and the ischange routine then finds a linear segment in the entropy plot, the pipeline will produce scaling exponents for nearly any input. Until you show it does not manufacture synchronization for independent or non-IPL signals, the specificity is unknown. This is the load-bearing gap. Second, the RESP high-pass filter at 2 Hz is post-hoc — chosen after the unfiltered signal failed to synchronize — and the authors themselves flag interpretability concerns. That is honest, but it means the three-way result is conditional on an unvalidated preprocessing choice, not a robust finding. Third, the 30 s windows with 20 s overlap make the p-values anti-conservative; the abstract's statement that brain, heart, and lung complexity \"significantly co-vary\" is stronger than the evidence, given that three-way significance holds in 24 of 54 datasets. The EEG-ECG coupling at 46/54 is more convincing, but the RESP additions are not.\n\nNone of these are fatal — they are standard fixes. What the paper needs is a surrogate analysis (e.g., phase randomization, independent AR(1) or colored noise, and mismatched pairing) to establish that the CS measure is specific, plus a correction for the overlapping-window inflation or a statement that the p-values are descriptive. The authors seem capable of doing this; the limitations section reads as genuinely thoughtful rather than pro forma.\n\nWho should read it: anyone working on network physiology, complexity synchronization, or multifractal methods for physiological time series. The automation is worth citing even if the empirical results are provisional. I would send this to referees. My own take is conditional acceptance: the methods contribution is solid, but the headline claim needs the surrogate controls before it can stand.","headline":"Useful automation of MDEA parameter selection with honest limitations, but the central empirical claim about brain-heart-lung synchronization is not yet supported without null-surrogate controls and corrected statistics.","tokens_in":16044,"tokens_out":1479,"would_cite":true,"duration_ms":17000,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a standardized diffusion-entropy method can detect synchronized complexity across EEG, ECG, and respiration during cognitive tasks.","keywords":["complexity synchronization","modified diffusion entropy analysis","multifractality","crucial events","inverse power law","EEG","ECG","respiration"],"falsifier":"Run the paper's order-8192 2-Hz high-pass filter on a synthetic signal built from a known deterministic oscillator plus a known inverse-power-law component, then apply the KS-stripe-size MDEA pipeline and check whether the recovered δ matches the injected IPL exponent; additionally, phase-shuffle the filtered RESP signal and recompute the EEG–RESP and ECG–RESP correlations, which should vanish if the synchronization reflects genuine coupling.","tokens_in":14963,"feed_emoji":"🧠","tokens_out":9476,"duration_ms":85093,"temperature":0.7,"pith_summary":"This paper tries to establish that complexity synchronization — simultaneous, time-varying changes in the multifractal scaling of different organ signals — is a real and measurable property of brain, heart, and lung time series recorded during cognitive task performance. The vehicle is modified diffusion entropy analysis (MDEA), which converts each signal into a sequence of crucial events (amplitude-threshold crossings) and reads off the slope δ of the diffusion entropy versus log-window plot, a scaling index tied to the inverse power-law exponent of the waiting times between events. The authors add two automatic parameter choices: a Kolmogorov-Smirnov statistic that selects the stripe size and power-law exponent by fitting the event inter-times to an inverse power law, and automatic detection of the linear fit region. Across 54 datasets from a Go-NoGo study, the method finds significant EEG–ECG synchronization in 46 of 54 datasets and significant three-way synchronization among EEG, ECG, and high-pass-filtered respiration in 24 of 54 datasets. A sympathetic reader would care because, if the claim holds, it gives a reproducible, data-driven way to measure coordination among heterogeneous physiological signals without frequency-band decompositions or trial averaging.","feed_headline":"Brain, heart, and lungs sync their complexity during tasks","feed_subtitle":"A standardized entropy-scaling method finds synchronized complexity across EEG, ECG, and filtered respiration.","key_machinery":"The load-bearing machinery is modified diffusion entropy analysis (MDEA). MDEA coarse-grains a time series into stripes by amplitude, marks crucial events when the signal crosses stripe boundaries, builds a binary event series, and accumulates it into a diffusion trajectory; the Shannon-Wiener entropy of segment endpoints, plotted against log window length, has slope δ, related to the inverse power-law exponent µ of the waiting-time distribution by µ = 1 + 1/δ. The paper's main technical additions are (1) a Kolmogorov-Smirnov statistic that scores how well the crucial-event inter-times fit an IPL complementary cumulative distribution, minimized over stripe size and µ by grid search, and (2) an automated change-point-based selection of the longest linear region of the entropy plot, with a special middle-segment rule for ECG. Complexity synchronization is then defined by correlations over sliding windows of the δ values across signals; the resulting time-varying δ tracks are the objects whose co-variation is tested.","core_discovery":"The paper's central claim is that complexity synchronization (CS) — high-order co-variation of the inverse power-law scaling index δ extracted by modified diffusion entropy analysis (MDEA) — is present among brain, heart, and lung organ-network time series during cognitive task performance, and that this can be detected repeatably once two MDEA parameters are chosen automatically. The authors first validate a KS-based estimator of stripe size and of the IPL index µ on synthetic duration-time data with known power-law exponents, showing lower bias and variance as sample size grows. Applied to 54 datasets (27 subjects × 2 stress conditions), the pipeline finds significant EEG–ECG synchronization in 46 of 54 datasets, and significant three-way synchronization among EEG, ECG, and a 2-Hz high-pass-filtered RESP signal in 24 of 54 datasets. A central caveat, stated by the authors, is that the unfiltered RESP signal did not synchronize with EEG and ECG, and meaningful lung synchronization appeared only after the dominant periodic component of respiration was removed, a preprocessing decision whose interpretation they flag as an open question.","pith_inferences":["An explicit test the authors did not run: apply the same order-8192 2-Hz high-pass filter to synthetic periodic-plus-power-law signals to confirm that the filter preserves the true IPL exponent; until then, the RESP results rest on preprocessing.","The correlation between EEG and ECG scaling indices could partly be driven by a shared global state such as arousal, movement, or task engagement rather than organ-network information flow; window-shuffled surrogates or partial correlations controlling for a common driver would discriminate.","If the CS measure generalizes, it could be applied to other multimodal recordings, such as fMRI with peripheral physiology or hyperscanning between individuals, wherever signals are long enough for reliable entropy-slope estimation; the paper's own claims are limited to the neurophysiological data analyzed.","The discussion's 'information force' interpretation goes beyond the data, which establish co-variation of scaling indices but not direction or mechanism of information transfer."],"forward_implications":["With automated parameter selection, MDEA/CS can be run on large-scale datasets: the authors estimate stripe sizes across 190,938 moving windows, and provide code and sample data for independent replication.","The empirical pattern of correlations, with stronger EEG–ECG than EEG–RESP or ECG–RESP coupling, implies that brain–heart complexity coupling is the most consistent organ-network link during cognitive tasks.","For periodic physiological signals such as respiration, CS analysis requires removing the dominant deterministic oscillation; otherwise the scaling index is not informative about complexity synchronization.","Tracking δ over time in sliding windows yields a real-time, single-trial measure of organ-network coordination, applicable to continuous recordings without trial averaging.","Cross-modal synchronization of scaling indices, if replicated, offers a candidate signature of network-level information transfer distinct from linear correlation, coherence, or phase coupling."],"supporting_citations":[{"why":"Introduces diffusion entropy analysis, the basis of MDEA, and the derivation of the entropy slope δ from segment endpoint distributions.","marker":"Scafetta and Grigolini (2002)"},{"why":"Prior proof-of-concept applying MDEA/CS to EEG, ECG, and RESP in two subjects; the approach this paper validates and automates.","marker":"Mahmoodi et al. (2023a)"},{"why":"Sets out the complexity synchronization hypothesis for organ networks on which the present analysis rests.","marker":"West et al. (2023b)"},{"why":"Source of the 30-subject neurofeedback Go-NoGo dataset with simultaneous EEG, ECG, and RESP recordings.","marker":"Kerick et al. (2023)"},{"why":"Provides the Kolmogorov-Smirnov statistic used to score the IPL fit for stripe-size and exponent selection.","marker":"Corder and Foreman (2014)"},{"why":"Defines the Mittag-Leffler map used to generate synthetic duration times for the bias/variance validation of the KS-based estimator.","marker":"Huillet (2016)"},{"why":"Provides the sample dataset and Matlab functions implementing MDEA, stripe-size search, and fit-region detection, making the pipeline reproducible.","marker":"Github (2024)"}],"fun_headline_variants":["Complexity sync links brain, heart, lungs in tasks","Entropy-scaling method shows organ complexity sync","Brain, heart sync complexity; lungs only when filtered","Organ complexity syncs during tasks, but lungs need filtering"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the 2-Hz high-pass filter applied to the respiration signal removes a deterministic oscillatory component while leaving an inverse-power-law complexity component intact; if the filter instead creates or distorts the IPL statistics, the reported EEG–RESP and ECG–RESP synchronization is an artifact of preprocessing.","fun_headline_variants_meta":{"raw":{"variants":["Complexity sync links brain, heart, lungs in tasks","Entropy-scaling method shows organ complexity sync","Brain, heart sync complexity; lungs only when filtered","Organ complexity syncs during tasks, but lungs need filtering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000573,"raw_usage":{"total_tokens":2669,"prompt_tokens":868,"completion_tokens":1801,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":1746}},"tokens_in":484,"tokens_out":1801,"duration_ms":11727,"temperature":1.0,"reasoning_tokens":1746,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:06:40.730203+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the paper's order-8192 2-Hz high-pass filter on a synthetic signal built from a known deterministic oscillator plus a known inverse-power-law component, then apply the KS-stripe-size MDEA pipeline and check whether the recovered δ matches the injected IPL exponent; additionally, phase-shuffle the filtered RESP signal and recompute the EEG–RESP and ECG–RESP correlations, which should vanish if the synchronization reflects genuine coupling.","supporting_citations":[],"review_version":1}