{"id":"089567f2-4b0e-46e1-a3d4-ae22989d8114","arxiv_id":"1908.01286","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Recurrence networks from short ECG recordings show bimodal degree distributions and two scaling regions in link density, with disease-specific variations.","lead":"Researchers built recurrence networks from one-minute ECG recordings of 96 patients and found that the networks' degree distributions are bimodal, with two peaks linked to separate spatial scales in the heart's reconstructed dynamics. The study also reports two scaling regimes in network link density versus threshold and suggests these measures could help distinguish cardiac disease classes from healthy controls.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'bimodality in all ECG types' claim is contradicted by the paper's own text: BB is said to be unimodal at the chosen ε=0.1, making the headline feature threshold-dependent.","rationale":"The reader's weakest assumption identifies threshold choice as the critical vulnerability, and the paper's text confirms rather than merely suggests the problem: BB is described as unimodal at ε=0.1 and bimodal only at ε=0.2, while the abstract and conclusion claim bimodality is common to all ECG types. This is a direct, load-bearing contradiction in the central claim. The scaling claim is similarly fragile because the authors report classes where the two regions are not discernible, and Table 2's large error bars make the disease-specific distinctions statistically unsupported. I am not raising an outside-consensus objection; the concern follows from the manuscript's own reported observations. The appropriate remedy is not rejection, because the study is exploratory and the proposed features could still be useful if properly quantified, but the conditional verdict is correct: the claims need a systematic, threshold-swept, statistically tested reanalysis before they can be accepted. My stress-test therefore leaves the reader's verdict unchanged, since the reader already required exactly this kind of additional validation.","tokens_in":7296,"tokens_out":3564,"duration_ms":41473,"concrete_test":"Re-analyze all 96 PTB records, not just typical cases, at ε = 0.1, 0.2, and 0.3, and apply a quantitative bimodality test (e.g., Hartigan dip test or Gaussian-mixture BIC) to each degree distribution, with phase-randomized surrogates to set the null. Record the fraction of bimodal records per class and per threshold. If BB records are mostly unimodal at ε=0.1, or if no single ε value yields bimodality for the large majority of records in every class, then the 'common feature of all ECG types' claim should be revised. Separately, fit the LD-versus-ε curves with piecewise linear regression and report breakpoint uncertainties to test whether the two scaling regions are statistically distinguishable per class.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that bimodality of the recurrence-network degree distribution is a stable, common feature of RNs from all ECG classes. The paper's own results undermine this. In the section 'Degree distribution', after stating that 'the distribution appears to be distorted for data from unhealthy cases', the authors write that 'For BB, the distribution is unimodal at this threshold', and later that for MI 'the first peak disappears as we increase ε to 0.5'. At the fixed construction threshold ε=0.1, one of the five classes is explicitly unimodal; the second peak appears only at ε=0.2. The abstract and conclusion nevertheless assert that bimodality is common to all types of ECG data. This is an internal inconsistency, not a matter of outside consensus. Because ε=0.1 is justified only by 'most of the networks become just connected', not by any bimodality criterion, the claimed two-scale structure is not shown to be a property of cardiac dynamics independent of threshold choice. The same problem weakens the scaling claim: for MI the authors say it is 'difficult to discern the distinct scaling regions', and for BB 'there is almost no change of slope'. Table 2's scaling indices also overlap heavily across classes (e.g., healthy γ2 = 0.387 ± 0.296 vs MI γ2 = 0.372 ± 0.276), so the reported numbers do not support disease-specific quantification. The load-bearing assumption is therefore not merely that ε=0.1 is convenient; it is that the observed bimodality and two-regime scaling survive a fixed, pre-specified threshold protocol across all classes. The text supplies direct counterexamples.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript constructs recurrence networks (RNs) from one-minute clinical ECG recordings taken from the PTB database, for healthy subjects and four disease classes (Bundle Branch Block, Cardiomyopathy, Dysrhythmia, Myocardial Infarction), and compares them with RNs from chaotic, hyperchaotic, and noise data. The two central claims are, first, that the RN degree distribution of ECG data is bimodal, reflecting two spatial scales in the reconstructed attractor, and second, that the link density as a function of the recurrence threshold displays two scaling regions, quantified by indices gamma1 and gamma2. The authors further suggest that these features and their averaged network measures can distinguish disease classes from healthy cases. The analysis uses standard time-delay embedding with fixed embedding dimension m=4 and a fixed recurrence threshold epsilon=0.1, and reports averaged network measures in Table 1 and scaling indices in Table 2.","tokens_in":7712,"tokens_out":3564,"duration_ms":38506,"significance":"If the bimodality and two-regime scaling claims were rigorously established, this would be a useful contribution to the analysis of short clinical ECG recordings, since recurrence-network measures are relatively simple to compute and do not require long data sets. The manuscript has strengths: it applies a well-defined, standard recurrence-network construction to public data; it includes comparisons with standard dynamical systems and noise; and it reports average values over a reasonably large set of 96 recordings. However, the central claims are not currently supported by the evidence presented. The paper's own text states that at the chosen threshold epsilon=0.1 the degree distribution for Bundle Branch Block is unimodal, directly contradicting the abstract's claim that bimodality is a common feature of RNs from all types of ECG data. The scaling claim is similarly weakened by statements that for Myocardial Infarction the scaling regions are difficult to discern and for Bundle Branch Block there is almost no change of slope. No quantitative bimodality test, no threshold-invariance analysis, and no surrogate-data comparison are given.","major_comments":[{"comment":"The central claim that bimodality is a common feature of all types of ECG data is internally inconsistent with the reported results. In the section 'Degree distribution', the authors state for BB: 'the distribution is unimodal at this threshold', and in the threshold-variation section they state that for MI 'the first peak disappears as we increase epsilon to 0.5'. The abstract and conclusion nevertheless assert that both bimodality and two-region scaling are 'common features of RNs from all types of ECG data'. Since the fixed construction threshold epsilon=0.1 is the only threshold used for the headline degree distributions in Fig. 3, the universal claim is contradicted by the manuscript's own observations. This must be resolved, either by revising the claims or by providing a systematic threshold-dependent definition of bimodality that is tested across all 96 datasets.","section":"Degree distribution and Abstract/Conclusion"},{"comment":"The choice of a single recurrence threshold epsilon=0.1 and a single embedding dimension m=4 for all classes is not justified for the central claims. The text justifies epsilon=0.1 only by saying 'most of the networks become just connected', which is a connectivity criterion, not a criterion for preserving bimodality or scaling behavior. The authors themselves note that BB is unimodal at epsilon=0.1 and becomes bimodal only at epsilon=0.2, while MI changes qualitatively as epsilon increases. Because the claimed two-scale attractor structure and the scaling indices depend on threshold, the manuscript needs to show that the qualitative features are stable over a meaningful range of epsilon and m, or to explain explicitly why a single fixed threshold is appropriate for the universal claims made in the abstract.","section":"Construction of Recurrence Networks from embedded phase space attractors of ECG data"},{"comment":"The two-scaling-region claim is not supported as a common feature. The text states that for MI 'it is difficult to discern the distinct scaling regions' and for BB 'there is almost no change of slope'. These statements contradict the general assertion that two scaling regions are observed for ECG data. In addition, the scaling indices in Table 2 show large standard deviations and heavy overlap across classes; for example, healthy gamma2 = 0.387 +/- 0.296 and MI gamma2 = 0.372 +/- 0.276 are statistically indistinguishable. The manuscript also does not state how the boundaries of the scaling regions are chosen. The authors need to provide a quantitative, reproducible procedure for identifying scaling regions, and should test whether gamma1 and gamma2 actually separate disease classes using a statistical test rather than visual inspection of Fig. 5.","section":"Scaling of link density with recurrence threshold, Table 2"},{"comment":"Fig. 3 shows degree distributions for single typical cases without error bars, uncertainty bands, or a quantitative bimodality test. The claim that the bimodal shape is a stable property of the recurrence-network ensemble requires at minimum a measure such as the dip test, Hartigan's dip test, or a Gaussian-mixture fit applied to all 96 datasets, together with a comparison against surrogate data or randomized thresholds. As presented, the evidence for bimodality is visual inspection of a few examples, which is not sufficient for a novelty claim of the strength made in the abstract.","section":"Degree distribution, Fig. 3"}],"minor_comments":[{"comment":"There is a typo in the abstract: 'we also show that that there is relevant information' should read 'we also show that there is relevant information'.","section":"Abstract"},{"comment":"The Fig. 3 caption states 'The degree distribution is bimodal in nature in all cases of ECG data sets', but the text in the same section says that BB is unimodal at the chosen threshold. The caption and the text should be brought into agreement.","section":"Fig. 3 caption and Degree distribution"},{"comment":"The downsampling procedure from 60000 to 6000 points by binning is described only in one sentence. The authors should specify the bin width, whether the binning uses averaging or another rule, and show that the downsampling does not remove the small-scale structures that are central to the bimodality interpretation.","section":"Construction of Recurrence Networks from embedded phase space attractors of ECG data"},{"comment":"The manuscript states that epsilon is varied from 0.1 to 1.0 in steps of 0.01 but results are shown only up to 0.5 because 'there is no significant change after 0.5'. The criterion for 'no significant change' is not given, and the scaling regions in Fig. 5 appear to be selected by eye; a reproducible method for locating the two regions is needed.","section":"Scaling of link density with recurrence threshold"},{"comment":"The conclusion says that the analysis 'is to be applied to larger number of data sets so that disease specific ranges of measures and scaling indices can be derived' and that this work is 'already in progress'. This statement is in tension with the earlier claims that disease-specific variations can be quantified from the present 96 datasets; the manuscript should clarify which claims are established now and which are prospective.","section":"Conclusion"},{"comment":"No data availability or code availability statement is provided, although all data come from a public database. A statement with the exact PTB record identifiers and, where possible, the analysis code would substantially increase reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope as an application of nonlinear time-series and network methods to physiological data, but the novelty claim ('for the first time') should be checked carefully against the existing recurrence-network and ECG literature, which is not discussed in the introduction. The main issue is not the method itself but the gap between the universal statements in the abstract and the class-by-class exceptions documented in the body; this gap can be addressed with additional quantitative analysis, but the revision will need to be substantial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, this paper reports something I haven't seen in the recurrence-network literature: bimodal degree distributions in RNs built from short ECG recordings, plus a two-regime scaling of link density with recurrence threshold. Second, the paper's own results contradict its headline. In Section 4 it says that for BB the distribution is unimodal at ε=0.1 and only becomes bimodal at 0.2, and for MI the first peak disappears as ε approaches 0.5. So \"bimodality common to all types\" is an overstatement.\n\nWhat's actually new and useful: the observation itself is new relative to the cited RN reviews and multifractal ECG work, and the comparison with Rössler, Chen, and white noise gives a point of reference. The idea of varying ε and extracting scaling indices is a reasonable way to get more from short data. The authors are upfront about secondary diagnoses and about needing larger samples. They use public PTB data, standard preprocessing, and there is no circularity in the method.\n\nSoft spots, in proportion. The central claim rests on visual inspection of Fig. 3; there is no quantitative bimodality test, no error bars on the distributions, and no surrogate comparison. That matters because the feature is threshold-dependent for at least one class. The scaling indices in Table 2 overlap heavily—healthy γ2 = 0.387 ± 0.296 vs MI γ2 = 0.372 ± 0.276—so disease-specific quantification is not actually demonstrated. The paper also acknowledges that for BB there is almost no slope change and for MI the regions are hard to discern. Those are honest statements, but they undercut the \"two scaling regions\" as a universal signature.\n\nSo the core findings are plausible and worth pursuing, but the claims as written are stronger than the evidence. This is a proof-of-concept, not a classification tool.\n\nWho should read it: people working on network-based physiological time-series analysis who want a short-data alternative to multifractal methods. It deserves a serious referee. I'd recommend sending it to peer review with a clear request: replace visual inspection with quantitative bimodality tests, add surrogate data, fit the scaling regions with an explicit algorithm, and temper the \"all types\" language. If the authors do that, the observation could become a solid contribution. As it stands, I wouldn't cite it as evidence for disease-specific indices, but I'd cite it as an early report of a novel empirical feature.","headline":"A promising but overclaimed proof-of-concept: the bimodality and scaling observations are new, yet the paper's own text shows the headline features are threshold-dependent and not class-separating as stated.","tokens_in":8152,"tokens_out":2469,"would_cite":false,"duration_ms":23571,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["05.45.Tp","64.60.aq","05.45.-a"],"model":"deepseek-v4-flash","headline":"This paper claims that one-minute ECG recordings, turned into recurrence networks, show a two-peaked degree distribution and two distinct scaling regimes in link density, with disease-specific values.","keywords":["recurrence networks","ECG","degree distribution","bimodality","link density","scaling exponents","phase space embedding","cardiac dynamics"],"falsifier":"Take the 96 recordings and, for each, compute the degree distribution over a grid of thresholds from 0.05 to 0.5 and embedding dimensions from 2 to 6. If a disease class is bimodal at one threshold and unimodal at neighboring thresholds, the bimodality is not robust. Separately, fit the link-density-versus-threshold curve on log-log axes without pre-chosen breakpoints; if two distinct linear segments cannot be identified consistently, the claimed two scaling regions do not exist.","tokens_in":7090,"feed_emoji":"❤️","tokens_out":10814,"duration_ms":98856,"temperature":0.7,"pith_summary":"This paper sets out to show that a one-minute clinical ECG recording, when converted into a recurrence network, carries two stable structural signatures of cardiac dynamics: a degree distribution with two separated peaks, and a link density that scales with the recurrence threshold in two distinct regimes. The authors report that both features appear across all ECG classes they study, while networks from standard chaotic, hyperchaotic, and noise signals are unimodal and scale differently. They interpret the two peaks as evidence of two spatial scales in the heart's reconstructed phase-space attractor, and the two scaling exponents as finer quantifiers of that structure. If the claim holds, recurrence networks give a way to quantify cardiac complexity from short, nonstationary clinical data using the full ECG waveform, not just heart-rate variability.","feed_headline":"Heart signals form networks with two connection peaks, unlike chaos","feed_subtitle":"The two peaks and two scaling regimes separate healthy from diseased hearts in one-minute recordings.","key_machinery":"The central object is the recurrence network: each point of the time-delay-embedded phase-space trajectory is a node, and two nodes are linked when their distance falls below a chosen threshold $\\varepsilon$, giving an adjacency matrix $A_{ij} = \\Theta(\\varepsilon - \\|\\vec{v}_i - \\vec{v}_j\\|) - \\delta_{ij}$. The embedding uses dimension $m=4$ with time delay taken from the autocorrelation falling to $1/e$, and $\\varepsilon=0.1$ is chosen as the value where most networks become just connected. The machinery works by translating geometric structure into network statistics: the small dense loop and larger ring of the cardiac attractor produce the two peaks in degree, and increasing $\\varepsilon$ from 0.1 to 0.5 gradually connects the network, exposing two scaling regimes in link density whose log-log slopes define $\\gamma_1$ and $\\gamma_2$.","core_discovery":"The paper's central discovery is that recurrence networks built from short clinical ECG time series have a bimodal degree distribution, with two well-separated peaks in the probability of a node having a given number of links, and the authors state this is the first time bimodality has been reported for recurrence networks derived from time series. Underlying the claim is the reconstructed phase-space attractor: healthy ECGs show a small dense loop-like region and a larger ring, and the dense region produces the second peak at higher degrees. The paper further finds that link density $LD$ varies with the recurrence threshold $\\varepsilon$ as $LD \\sim \\varepsilon^{\\gamma}$ in two distinct scaling regions, with exponents $\\gamma_1$ and $\\gamma_2$ that differ across diagnostic classes. Disease-specific average values of clustering coefficient, average path length, link density, and the two scaling exponents are presented as quantifiable signatures, with healthy cases showing the least variability.","pith_inferences":["If the bimodality is robust, a simple classifier could use the positions and heights of the two degree peaks alone; the paper does not test classification accuracy, so this is an extension, not a claim.","The two scaling regimes might correspond to fast local recurrences within the dense loop and slower returns around the larger ring, a dynamical interpretation the paper does not explicitly make.","A direct test of the mechanism would be to generate synthetic ECG-like signals with and without the small-scale loop; if the second peak disappears when the loop is removed, the geometric interpretation would be confirmed.","Applying the same threshold-scanning analysis to heart-rate-variability series could connect this network signature to the larger HRV classification literature."],"forward_implications":["Healthy ECGs show the least variability in clustering coefficient and average path length, so those measures provide a stable baseline for detecting abnormality.","Because the two scaling exponents differ across diagnostic classes, $\\gamma_1$ and $\\gamma_2$ can serve as compact quantitative features for separating cardiac conditions.","The method works on one-minute recordings, allowing complexity quantification for short and nonstationary clinical data where multifractal analysis would need longer series.","The bimodality separates ECG recurrence networks from networks of standard chaotic, hyperchaotic, and noise signals, placing cardiac dynamics in a distinct structural class.","In bundle branch block cases the second peak appears only when the threshold is increased, showing that the underlying structure exists but is less dense than in healthy cases."],"supporting_citations":[{"why":"Introduces the recurrence network construction and the adjacency-matrix definition used throughout the paper.","marker":"[2]"},{"why":"Surveys recurrence network measures and the rationale for choosing the recurrence threshold $\\varepsilon$.","marker":"[3]"},{"why":"Supplies the clinical ECG recordings from 96 cases used for all network calculations.","marker":"[12]"},{"why":"Provides the data repository through which the ECG records are accessed and validated.","marker":"[13]"},{"why":"Gives the time-delay embedding method used to reconstruct the phase-space attractor from the ECG time series.","marker":"[14]"},{"why":"Provides the false nearest neighbours algorithm used to select the embedding dimension $m=4$.","marker":"[17]"},{"why":"Defines the standard network measures, including degree distribution, clustering coefficient, and average path length, that the study computes.","marker":"[19]"},{"why":"Supplies the chaotic and hyperchaotic comparison systems whose recurrence networks are contrasted with ECG networks.","marker":"[20]"}],"fun_headline_variants":["ECG recurrence networks show two-scale degree bimodality","Heart data networks reveal two connection peaks","First glimpse of bimodal degree in ECG networks","Two scaling regimes in ECG recurrence networks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes that one fixed embedding dimension and one fixed recurrence threshold are appropriate for every ECG class, so if that choice creates or hides the second peak in some classes, the bimodality is an artifact of parameter selection rather than a property of cardiac dynamics.","fun_headline_variants_meta":{"raw":{"variants":["ECG recurrence networks show two-scale degree bimodality","Heart data networks reveal two connection peaks","First glimpse of bimodal degree in ECG networks","Two scaling regimes in ECG recurrence networks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000469,"raw_usage":{"total_tokens":2320,"prompt_tokens":911,"completion_tokens":1409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":1352}},"tokens_in":527,"tokens_out":1409,"duration_ms":9741,"temperature":1.0,"reasoning_tokens":1352,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:17:11.511050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the 96 recordings and, for each, compute the degree distribution over a grid of thresholds from 0.05 to 0.5 and embedding dimensions from 2 to 6. If a disease class is bimodal at one threshold and unimodal at neighboring thresholds, the bimodality is not robust. Separately, fit the link-density-versus-threshold curve on log-log axes without pre-chosen breakpoints; if two distinct linear segments cannot be identified consistently, the claimed two scaling regions do not exist.","supporting_citations":[{"cited_title":"V., Zou Y., Donges J","cited_arxiv_id":null,"evidence_quote":"Introduces the recurrence network construction and the adjacency-matrix definition used throughout the paper."},{"cited_title":"V., Marwan N., Donges J","cited_arxiv_id":null,"evidence_quote":"Surveys recurrence network measures and the rationale for choosing the recurrence threshold $\\varepsilon$."},{"cited_title":"and Schnabel A","cited_arxiv_id":null,"evidence_quote":"Supplies the clinical ECG recordings from 96 cases used for all network calculations."},{"cited_title":"and Kantz H., Chaos: An Interdisciplinary Journal of Nonlinear Science , 25 (2015) 097610","cited_arxiv_id":null,"evidence_quote":"Gives the time-delay embedding method used to reconstruct the phase-space attractor from the ECG time series."},{"cited_title":"and Schreiber T","cited_arxiv_id":null,"evidence_quote":"Provides the false nearest neighbours algorithm used to select the embedding dimension $m=4$."},{"cited_title":", Networks: an introduction (Oxford univer- sity press) 2010","cited_arxiv_id":null,"evidence_quote":"Defines the standard network measures, including degree distribution, clustering coefficient, and average path length, that the study computes."},{"cited_title":"and Yuan Z","cited_arxiv_id":null,"evidence_quote":"Supplies the chaotic and hyperchaotic comparison systems whose recurrence networks are contrasted with ECG networks."}],"review_version":1}