Maximum likelihood estimation of burst-merging kernels for bursty time series
Pith reviewed 2026-05-23 00:46 UTC · model grok-4.3
The pith
Maximum likelihood estimation recovers the burst-merging kernel from observed time series.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The burst-tree structure can be simply characterized by a burst-merging kernel that dictates which bursts are merged together as the timescale increases. We develop the maximum likelihood estimation method of the burst-merging kernel from time series, which is successfully tested against the time series generated using several model kernels. We also apply our method to some empirical time series from various backgrounds.
What carries the argument
The burst-merging kernel, a function that determines which bursts merge as timescale coarsens, recovered by maximum likelihood estimation from the observed event times.
If this is right
- The estimated kernel supplies a precise quantitative characterization of any given time series.
- Comparison of kernels across datasets becomes possible once the estimation procedure is applied.
- Underlying generative mechanisms can be studied more accurately once the kernel is known.
- The method works for both synthetic series drawn from known kernels and real empirical series.
Where Pith is reading between the lines
- Kernels recovered from different domains could be compared to test whether bursty processes share common merging rules.
- The estimation procedure could be extended to time series with slowly varying statistics by allowing the kernel to depend on external covariates.
- If the kernel turns out to be low-dimensional in many systems, it would reduce the effective degrees of freedom needed to model burst hierarchies.
Load-bearing premise
The hierarchical merging pattern of bursts is fully captured by a single kernel depending on burst properties.
What would settle it
Forward simulation of new event series from the estimated kernel that produces burst trees whose merging statistics deviate from those measured in the original data.
Figures
read the original abstract
Various time series in natural and social processes have been found to be bursty. Events in the time series rapidly occur within short time periods, forming bursts, which are alternated with long inactive periods. As the timescale defining bursts increases, individual events are sequentially merged to become small bursts and then bigger ones, eventually leading to the single burst containing all events. Such a merging pattern has been depicted by a tree that fully reveals the hierarchical structure of bursts, thus called a burst tree. The burst-tree structure can be simply characterized by a burst-merging kernel that dictates which bursts are merged together as the timescale increases. In this work, we develop the maximum likelihood estimation method of the burst-merging kernel from time series, which is successfully tested against the time series generated using several model kernels. We also apply our method to some empirical time series from various backgrounds. Our method provides a useful tool to precisely characterize the time series data, hence enabling to study their underlying mechanisms more accurately.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a maximum likelihood estimation (MLE) method to recover the burst-merging kernel that governs the hierarchical merging of events into bursts as the observation timescale increases. The kernel is estimated directly from observed event times. The procedure is validated by generating synthetic time series from several chosen model kernels and recovering the input kernels; the method is then applied to empirical bursty time series from multiple domains.
Significance. If the estimator is correctly derived and implemented, the approach supplies a compact, data-driven characterization of burst-tree structure that goes beyond conventional inter-event time statistics. Successful synthetic recovery under the modeling assumptions demonstrates that the numerical procedure is feasible and can serve as a quantitative tool for comparing bursty processes across systems.
major comments (1)
- [synthetic tests and empirical application] The validation consists exclusively of recovering kernels from data generated by the same family of model kernels used to define the likelihood (abstract and synthetic-test section). This confirms numerical correctness of the MLE under the modeling assumption but leaves untested whether any single kernel suffices to describe merging statistics in data generated by other processes or in empirical series; if the single-kernel assumption fails, the estimated kernel is not a sufficient statistic for the observed burst tree.
minor comments (2)
- [abstract] The abstract states that the method is 'successfully tested' but supplies no explicit form of the likelihood, no optimization details, and no quantitative recovery metrics (e.g., parameter error or likelihood ratio); including these would allow immediate assessment of the central technical contribution.
- [empirical results] When the estimated kernel is applied to empirical series, report at least the maximized log-likelihood value and a comparison against a baseline (e.g., uniform or exponential kernel) to quantify the improvement.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive feedback on our manuscript. We respond to the major comment below.
read point-by-point responses
-
Referee: [synthetic tests and empirical application] The validation consists exclusively of recovering kernels from data generated by the same family of model kernels used to define the likelihood (abstract and synthetic-test section). This confirms numerical correctness of the MLE under the modeling assumption but leaves untested whether any single kernel suffices to describe merging statistics in data generated by other processes or in empirical series; if the single-kernel assumption fails, the estimated kernel is not a sufficient statistic for the observed burst tree.
Authors: We agree that the synthetic tests verify numerical correctness and unbiased recovery of the kernel when data are generated exactly from the assumed model family; this is the conventional validation for an MLE procedure. The empirical applications demonstrate practical use of the estimator on real bursty series, yielding a compact characterization under the model. The manuscript does not claim that a single kernel is always sufficient for arbitrary generative processes or that the estimate is invariably a sufficient statistic; it is presented as the maximum-likelihood description within the model class. When the assumption is violated, the fitted kernel remains the best approximation under the model (with possible diagnostics via likelihood or residuals), but we do not test model misspecification here. Broadening the validation to other generative mechanisms lies beyond the present scope. revision: no
Circularity Check
No circularity: MLE estimator derives kernel from observed event times without self-definition or fitted-input reduction
full rationale
The paper introduces a maximum likelihood procedure that takes raw event timestamps as input and outputs an estimated burst-merging kernel. Validation consists of forward simulation from chosen kernels followed by recovery; this tests numerical correctness of the estimator under the modeling assumption but does not make the target kernel a function of itself or rename fitted parameters as predictions. No self-citation chain, ansatz smuggling, or uniqueness theorem is invoked to close the derivation. The central claim therefore remains an independent statistical construction rather than a tautology.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The burst-tree structure can be simply characterized by a burst-merging kernel that dictates which bursts are merged together as the timescale increases.
Reference graph
Works this paper leans on
-
[1]
At the initial time step s = 0, each event makes its own burst of size one
-
[2]
At each time step s, two bursts of sizes b and b′ are chosen at random with a probability proportional to Kbb′ [Eq. (5)]
-
[3]
Two merged bursts are assigned to left and right children nodes at random
Those two bursts are merged to make another burst of size b + b′, which becomes a parent node. Two merged bursts are assigned to left and right children nodes at random
-
[4]
Steps 2–3 are repeated until all events belong to the single burst of size n. Once the ordinal burst tree is generated, we randomly draw n − 1 IETs from the power-law IET distribution as P (τ) = τ −α Pτc τ =1 τ −α for τ = 1, . . . , τc, (15) where α is the power-law exponent and τc is the upper bound. These IETs are sorted in a descending order to obtain ...
-
[5]
A.-L. Barab´ asi and M. P´ osfai,Network Science (Cam- bridge University Press, Cambridge, 2016)
work page 2016
-
[6]
M. E. J. Newman, Networks, second edition ed. (Oxford University Press, Oxford, United Kingdom; New York, NY, United States of America, 2018)
work page 2018
-
[7]
F. Menczer, S. Fortunato, and C. A. Davis, A First Course in Network Science (Cambridge University Press, Cambridge, 2020)
work page 2020
-
[8]
S. N. Dorogovtsev and J. F. F. Mendes, The Nature of Complex Networks (Oxford University Press, Oxford, UK; New York, NY, 2022)
work page 2022
-
[9]
D. Easley and J. Kleinberg, Networks, Crowds, and Mar- kets: Reasoning About a Highly Connected World (Cam- bridge University Press, 2010)
work page 2010
- [10]
-
[11]
P. Holme and J. Saram¨ aki, Temporal networks, Physics Reports 519, 97 (2012)
work page 2012
-
[12]
P. Holme and J. Saram¨ aki, eds.,Temporal Network The- ory, Computational Social Sciences (Springer Interna- tional Publishing, Cham, 2019)
work page 2019
-
[13]
N. Masuda and R. Lambiotte, A Guide to Temporal Net- works, Series on Complexity Science (World Scientific, New Jersey, 2016)
work page 2016
-
[14]
Barab´ asi, The origin of bursts and heavy tails in human dynamics, Nature 435, 207 (2005)
A.-L. Barab´ asi, The origin of bursts and heavy tails in human dynamics, Nature 435, 207 (2005)
work page 2005
-
[15]
M. Karsai, H.-H. Jo, and K. Kaski, Bursty Human Dy- namics (Springer International Publishing, Cham, 2018)
work page 2018
- [16]
-
[17]
P. Bak, C. Tang, and K. Wiesenfeld, Self-organized crit- icality: An explanation of the 1/f noise, Physical Review Letters 59, 381 (1987)
work page 1987
-
[18]
H. J. Jensen, Self-Organized Criticality: Emergent Com- plex Behavior in Physical and Biological Systems , 1st ed. (Cambridge University Press, 1998)
work page 1998
-
[19]
M. S. Wheatland, P. A. Sturrock, and J. M. McTiernan, The Waiting-Time Distribution of Solar Flare Hard X- Ray Bursts, The Astrophysical Journal 509, 448 (1998)
work page 1998
-
[20]
´A. Corral, Long-term clustering, scaling, and universal- ity in the temporal occurrence of earthquakes, Physical Review Letters 92, 108501 (2004)
work page 2004
-
[21]
J. M. Beggs and D. Plenz, Neuronal Avalanches in Neo- cortical Circuits, The Journal of Neuroscience 23, 11167 (2003)
work page 2003
-
[22]
T. Petermann, T. C. Thiagarajan, M. A. Lebedev, M. A. L. Nicolelis, D. R. Chialvo, and D. Plenz, Spon- taneous cortical activity in awake monkeys composed of neuronal avalanches, Proceedings of the National Academy of Sciences 106, 15921 (2009)
work page 2009
-
[23]
T. Kemuriyama, H. Ohta, Y. Sato, S. Maruyama, M. Tandai-Hiruma, K. Kato, and Y. Nishida, A power- law distribution of inter-spike intervals in renal sym- pathetic nerve activity in salt-sensitive hypertension- induced chronic heart failure, BioSystems 101, 144 (2010)
work page 2010
-
[24]
K.-I. Goh and A.-L. Barab´ asi, Burstiness and memory in complex systems, EPL (Europhysics Letters) 81, 48002 (2008)
work page 2008
-
[25]
R. Crane and D. Sornette, Robust dynamic classes re- vealed by measuring the response function of a social system, Proceedings of the National Academy of Sciences 105, 15649 (2008)
work page 2008
- [26]
-
[27]
H.-H. Jo, T. Hiraoka, and M. Kivel¨ a, Burst-tree decom- position of time series reveals the structure of temporal correlations, Scientific Reports 10, 12202 (2020)
work page 2020
-
[28]
J. Choi, T. Hiraoka, and H.-H. Jo, Individual-driven versus interaction-driven burstiness in human dynamics: The case of Wikipedia edit history, Physical Review E 104, 014312 (2021). 8
work page 2021
- [29]
-
[30]
T. Birhanu and H.-H. Jo, Burst-tree structure and higher-order temporal correlations, Physical Review E 111, 014308 (2025)
work page 2025
-
[31]
T. Pham, P. Sheridan, and H. Shimodaira, PAFit: A Sta- tistical Method for Measuring Preferential Attachment in Temporal Complex Networks, PLoS ONE 10, e0137796 (2015)
work page 2015
-
[32]
A.-L. Barab´ asi and R. Albert, Emergence of Scaling in Random Networks, Science 286, 509 (1999)
work page 1999
-
[33]
P. L. Krapivsky, S. Redner, and F. Leyvraz, Connectivity of Growing Random Networks, Physical Review Letters 85, 4629 (2000)
work page 2000
-
[34]
M. E. J. Newman, Clustering and preferential attach- ment in growing networks, Physical Review E64, 025102 (2001)
work page 2001
- [35]
-
[36]
P. Sheridan, Y. Yagahara, and H. Shimodaira, Measur- ing preferential attachment in growing networks with missing-timelines using Markov chain Monte Carlo, Phys- ica A: Statistical Mechanics and its Applications 391, 5031 (2012)
work page 2012
-
[37]
J. Kunegis, M. Blattner, and C. Moser, Preferential at- tachment in online networks: Measurement and explana- tions, in Proceedings of the 5th Annual ACM Web Science Conference (ACM, Paris France, 2013) pp. 205–214
work page 2013
-
[38]
W. H. Stockmayer, Theory of Molecular Size Distribu- tion and Gel Formation in Branched-Chain Polymers, The Journal of Chemical Physics 11, 45 (1943)
work page 1943
-
[39]
Lushnikov, Evolution of coagulating systems, Journal of Colloid and Interface Science 45, 549 (1973)
A. Lushnikov, Evolution of coagulating systems, Journal of Colloid and Interface Science 45, 549 (1973)
work page 1973
-
[40]
W. H. White, On the form of steady-state solutions to the coagulation equations, Journal of Colloid and Interface Science 87, 204 (1982)
work page 1982
-
[41]
E. M. Hendriks, M. H. Ernst, and R. M. Ziff, Coagulation equations with gelation, Journal of Statistical Physics31, 519 (1983)
work page 1983
-
[42]
D. J. Aldous, Deterministic and Stochastic Models for Coalescence (Aggregation and Coagulation): A Review of the Mean-Field Theory for Probabilists, Bernoulli 5, 3 (1999)
work page 1999
-
[43]
M. H. Lee, A survey of numerical solutions to the coagu- lation equation, Journal of Physics A: Mathematical and General 34, 10219 (2001)
work page 2001
-
[44]
F. Leyvraz, Scaling theory and exactly solved models in the kinetics of irreversible aggregation, Physics Reports 383, 95 (2003)
work page 2003
-
[45]
F. Leyvraz, Rigorous Results in the Scaling Theory of Irreversible Aggregation Kinetics:, Journal of Nonlinear Mathematical Physics 12, 449 (2005)
work page 2005
-
[46]
J. A. Wattis, An introduction to mathematical models of coagulation–fragmentation processes: A discrete deter- ministic mean-field approach, Physica D: Nonlinear Phe- nomena 222, 1 (2006)
work page 2006
-
[47]
T. Birhanu and H.-H. Jo, Codes for the maximum likelihood estimation method of burst-merging kernels, https://github.com/tibebe22/Kernel-estimate (2025)
work page 2025
-
[48]
English Wikipedia, https://dumps.wikimedia.org/
-
[49]
J. Yang and J. Leskovec, Patterns of temporal variation in online media, in Proceedings of the Fourth ACM In- ternational Conference on Web Search and Data Mining (ACM Press, Hong Kong, China, 2011) pp. 177–186
work page 2011
-
[50]
PhysioBank, https://physionet.org/physiobank/
-
[51]
Japan University Network Earthquake Catalog, http://wwweic.eri.u-tokyo.ac.jp/CATALOG/junec/
-
[52]
D. R. Hunter and K. Lange, A Tutorial on MM Algo- rithms, The American Statistician 58, 30 (2004)
work page 2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.