A Mixed Self-Exciting Process to Model Epileptic Seizures
Pith reviewed 2026-05-22 04:15 UTC · model grok-4.3
The pith
A Bayesian mixed Hawkes process quantifies seizure clustering and shows bias from ignoring patient heterogeneity in epilepsy data.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We propose a Bayesian mixed Hawkes process model that addresses seizure clustering and heterogeneity between individuals. In the Hawkes process, the intensity is accelerated each time an event occurs, through the composition of background and excitation intensity functions. The proposed model incorporates a Weibull baseline intensity to model a trend in background seizure rates over time, while the excitation process accounts for seizure clustering within individuals. We model heterogeneity among individuals by including covariates and random effects in both the background and excitation intensities. In the HEP study, the average time between primary and secondary seizures within an ividuals
What carries the argument
Mixed Hawkes process intensity that adds a Weibull baseline trend to an excitation kernel, with random effects and covariates for individual heterogeneity in background and excitation.
If this is right
- The model yields specific estimates for the average time between primary and secondary seizures of 1.57 days and 2.20 seizures per cluster.
- Including random effects corrects for underestimation of background intensity when heterogeneity is present.
- Including random effects corrects for overestimation of excitation rates when heterogeneity is present.
- The approach allows individualized modeling of seizure patterns through covariates and random effects.
Where Pith is reading between the lines
- The clustering estimates could support clinical decisions on adjusting treatment or monitoring after an observed seizure.
- Similar mixed self-exciting models could be tested on other recurrent medical events such as migraine attacks to measure their clustering properties.
- Replacing the Weibull baseline with other flexible time trends would test whether the self-excitation findings remain stable.
Load-bearing premise
The observed seizures arise from a self-exciting intensity that combines a smooth Weibull trend in the background rate with a decaying boost after each seizure, plus patient-specific random adjustments to those components.
What would settle it
Collect an independent set of seizure diary records from a comparable group of patients and verify whether the mean number of seizures per cluster falls within the credible interval 1.96 to 2.47.
Figures
read the original abstract
Epilepsy is a neurological disorder characterized by recurrent seizures affecting more than 70 million people worldwide. Often, an individual with epilepsy is more likely to experience subsequent seizures following an initial seizure, a process we call seizure clustering. Motivated by seizure diary data collected over three years from 407 individuals newly diagnosed with focal epilepsy in the Human Epilepsy Project (HEP), we propose a Bayesian mixed Hawkes process model that addresses seizure clustering and heterogeneity between individuals. In the Hawkes process, the intensity is accelerated each time an event occurs, through the composition of background and excitation intensity functions. The proposed model incorporates a Weibull baseline intensity to model a trend in background seizure rates over time, while the excitation process accounts for seizure clustering within individuals. We model heterogeneity among individuals by including covariates and random effects in both the background and excitation intensities. In the HEP study, the average time between primary and secondary seizures within an individual is 1.57 (95\% CrI: 1.43, 1.70) days, with an average of 2.20 (1.96, 2.47) seizures per cluster. We demonstrate that omitting random effects in the presence of heterogeneity leads to underestimation of the background intensity and overestimation of excitation rates.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a Bayesian mixed Hawkes process for modeling seizure clustering and individual heterogeneity in diary data from 407 newly diagnosed focal epilepsy patients in the Human Epilepsy Project. The intensity combines a Weibull baseline trend, an excitation kernel, covariates, and random effects on both background and excitation components. Posterior summaries include an average of 1.57 days (95% CrI 1.43–1.70) between primary and secondary seizures and 2.20 seizures per cluster (1.96–2.47); a simulation study is used to show that omitting random effects underestimates background intensity and overestimates excitation rates.
Significance. If the simulation results are robust, the work supplies a practical Bayesian framework for recurrent-event data with self-excitation and unobserved heterogeneity. Explicit demonstration of bias induced by ignoring random effects is a useful methodological warning for applied Hawkes-process analyses in biostatistics and neurology.
major comments (2)
- [Simulation Study] Simulation study: the reported direction of bias (underestimation of background intensity, overestimation of excitation) is shown only when data are generated from the exact random-effects distribution and variance components used in the simulation. It is not demonstrated that the same bias pattern would arise under plausible alternative heterogeneity structures (heavier tails, correlation between background and excitation random effects, or time-varying heterogeneity) that might better match the HEP diaries.
- [Model Specification] Model definition: the precise functional form of the intensity (Weibull baseline plus excitation kernel plus random effects) and the identifiability constraints imposed on the random-effects variances are not fully specified in a single equation block. Without this, it is difficult to verify that the posterior estimates of cluster size and inter-seizure time are invariant to reparameterization.
minor comments (2)
- [Results] Table or figure presenting the posterior means and credible intervals for the Weibull shape/scale and excitation parameters would improve readability of the main results.
- [Data Description] The description of the HEP data (number of seizures per patient, observation lengths, covariate distributions) is summarized only briefly; a supplementary table with basic descriptive statistics would aid reproducibility.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive comments on our manuscript. We address each major comment below with honest assessments of the current work and clear plans for revision. We believe these changes will strengthen the presentation and robustness discussion without altering the core contributions.
read point-by-point responses
-
Referee: [Simulation Study] Simulation study: the reported direction of bias (underestimation of background intensity, overestimation of excitation) is shown only when data are generated from the exact random-effects distribution and variance components used in the simulation. It is not demonstrated that the same bias pattern would arise under plausible alternative heterogeneity structures (heavier tails, correlation between background and excitation random effects, or time-varying heterogeneity) that might better match the HEP diaries.
Authors: We agree that the simulation demonstrates the bias specifically under the random-effects distribution and variance components used in our model. This choice was deliberate to illustrate the consequences of omitting random effects when they are present in the data-generating process, providing a direct methodological warning for applied Hawkes process analyses. While we acknowledge that additional simulations under alternative structures (e.g., heavier tails, correlated random effects, or time-varying heterogeneity) would offer broader robustness evidence, such extensions involve substantial computational cost and are unlikely to reverse the qualitative finding that ignoring heterogeneity biases background and excitation estimates. In revision, we will add an explicit discussion of this scope limitation and note that the demonstrated bias pattern is representative for the model class considered. revision: partial
-
Referee: [Model Specification] Model definition: the precise functional form of the intensity (Weibull baseline plus excitation kernel plus random effects) and the identifiability constraints imposed on the random-effects variances are not fully specified in a single equation block. Without this, it is difficult to verify that the posterior estimates of cluster size and inter-seizure time are invariant to reparameterization.
Authors: We thank the referee for identifying this presentational gap. The intensity function is described across sections but not consolidated. In the revised manuscript, we will add a single, self-contained equation block that explicitly defines the full intensity as the sum of the Weibull baseline, the excitation kernel, covariate terms, and the random effects on both background and excitation components. We will also state the identifiability constraints (e.g., any centering or fixed-scale assumptions on the random-effects variances) in the same block. This will make it straightforward to verify that the reported posterior summaries for cluster size and inter-seizure times remain invariant under reparameterization. revision: yes
Circularity Check
No significant circularity detected in derivation or demonstration
full rationale
The paper defines a Bayesian mixed Hawkes process with Weibull baseline intensity, excitation kernel, covariates, and random effects for individual heterogeneity, then fits this model to the external HEP diary data to report posterior summaries such as mean inter-seizure time and cluster size. The demonstration that omitting random effects biases background and excitation estimates is performed via a standard simulation study that generates data under the fitted mixed model and refits the misspecified version; this illustrates consequences of misspecification under the paper's own assumptions rather than reducing any reported prediction or first-principles result to an input quantity by construction. No self-definitional equations, fitted-input predictions, or load-bearing self-citations that collapse the central claims are present. The derivation chain remains self-contained against the external HEP data and conventional Hawkes-process literature.
Axiom & Free-Parameter Ledger
free parameters (3)
- Weibull baseline parameters
- Excitation kernel parameters
- Random effects variances
axioms (2)
- domain assumption Seizure times follow a self-exciting point process whose intensity is the sum of a baseline and excitation terms.
- domain assumption Bayesian hierarchical modeling with random effects and priors is appropriate for the clustered, heterogeneous data.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
λ*_i(t) = ν_i μ0(t) exp(∑ β_p x_ip) + ω_i exp(∑ ζ_q z_iq) ∑ exp(−δ(t−t_ij)) with μ0(t)=α t^{α−1}, ν_i~Gamma(ϕ,ϕ), ω_i~Gamma(ξ,ξ)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
branching ratio r_i = ω_i κ_i / δ; average cluster size 1/(1−r_i)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Handling missing data in self-exciting point process models , volume =. Spatial Statistics , author =. 2019 , pages =. doi:10.1016/j.spasta.2018.12.004 , abstract =
-
[2]
Intelligent Data Analysis , author =
Bayesian estimation of decay parameters in. Intelligent Data Analysis , author =. 2023 , pages =. doi:10.3233/IDA-216283 , abstract =
-
[3]
The Annals of Statistics , author =
Nonparametric. The Annals of Statistics , author =. 2020 , file =. doi:10.1214/19-AOS1903 , language =
-
[4]
Deutsch, Isabella and Ross, Gordon J. , month = jun, year =. arXiv preprint , publisher =
-
[5]
Proceedings of the AAAI Conference on Artificial Intelligence , author =
Hawkes. Proceedings of the AAAI Conference on Artificial Intelligence , author =. 2018 , file =. doi:10.1609/aaai.v32i1.12116 , abstract =
-
[6]
Journal of Computational and Graphical Statistics , author =
A. Journal of Computational and Graphical Statistics , author =. 2020 , pages =. doi:10.1080/10618600.2019.1617159 , abstract =
-
[7]
Statistical Science , author =
Bayesian. Statistical Science , author =. 2018 , file =. doi:10.1214/17-STS630 , abstract =
-
[8]
Journal of the American Statistical Association , author =
A. Journal of the American Statistical Association , author =. 2015 , pages =. doi:10.1080/01621459.2014.969424 , language =
-
[9]
Journal of Computational and Graphical Statistics , author =
Parameter. Journal of Computational and Graphical Statistics , author =. 2022 , pages =. doi:10.1080/10618600.2022.2050247 , abstract =
-
[10]
Journal of Futures Markets , author =
Hawkes. Journal of Futures Markets , author =. 2014 , pages =. doi:10.1002/fut.21644 , abstract =
-
[11]
and Lee, Young and Pollett, Philip K
Laub, Patrick J. and Lee, Young and Pollett, Philip K. and Taimre, Thomas , month = may, year =. Hawkes
-
[12]
Prospective validation study of an epilepsy seizure risk system for outpatient evaluation , volume =. Epilepsia , author =. 2020 , pages =. doi:10.1111/epi.16397 , abstract =
-
[13]
Does my posterior look big in this?—. Epilepsia Open , author =. 2019 , pages =. doi:10.1002/epi4.12328 , language =
-
[14]
A big data approach to the development of mixed‐effects models for seizure count data , volume =. Epilepsia , author =. 2017 , pages =. doi:10.1111/epi.13727 , abstract =
-
[15]
Epilepsy as a dynamic disease:. Epilepsia Open , author =. 2018 , pages =. doi:10.1002/epi4.12112 , abstract =
-
[16]
Annals of Clinical and Translational Neurology , author =
Is seizure frequency variance a predictable quantity? , volume =. Annals of Clinical and Translational Neurology , author =. 2018 , pages =. doi:10.1002/acn3.519 , abstract =
-
[17]
The Annals of Applied Statistics , author =
Bayesian non-homogeneous hidden. The Annals of Applied Statistics , author =. 2023 , file =. doi:10.1214/22-AOAS1630 , language =
-
[18]
Auvin, Stéphane and French, Jacqueline and Dlugos, Denis and Knupp, Kelly G and Perucca, Emilio and Arzimanoglou, Alexis and Whalen, Ed and Shellhaas, Renée A , file =. Novel study design to assess the efficacy and tolerability of antiseizure medications for focal‐onset seizures in infants and young children:
-
[19]
Simple discrete-time self-exciting models can describe complex dynamic processes:. PLOS ONE , author =. 2021 , pages =. doi:10.1371/journal.pone.0250015 , abstract =
-
[20]
A flexible, random histogram kernel for discrete-time
Browning, Raiha and Rousseau, Judith and Mengersen, Kerrie , month = aug, year =. A flexible, random histogram kernel for discrete-time
-
[21]
Miscouridou, Xenia and Caron, Francois and Teh, Yee Whye , file =. Modelling sparsity, heterogeneity, reciprocity and community structure in temporal interaction data , abstract =
-
[22]
Miscouridou, Xenia and Bhatt, Samir and Mohler, George and Flaxman, Seth and Mishra, Swapnil , month = oct, year =. Cox-
- [23]
-
[24]
Engelhard, Matthew and Xu, Hongteng and Carin, Lawrence and Oliver, Jason A and Hallyburton, Matthew and McClernon, F Joseph , year =. Predicting
-
[25]
Journal of the Royal Statistical Society Series C: Applied Statistics , author =
Network. Journal of the Royal Statistical Society Series C: Applied Statistics , author =. 2022 , pages =. doi:10.1111/rssc.12581 , abstract =
-
[26]
and Curley, James and Zheng, Tian , month = aug, year =
Wu, Jing and Ward, Owen G. and Curley, James and Zheng, Tian , month = aug, year =. Markov-
-
[27]
Probability and random processes , isbn =
Grimmett, Geoffrey and Stirzaker, David , year =. Probability and random processes , isbn =
-
[28]
Journal of Computational and Graphical Statistics , author =
Programming with models: writing statistical algorithms for general model structures with. Journal of Computational and Graphical Statistics , author =. 2017 , pages =. doi:10.1080/10618600.2016.1172487 , number =
-
[29]
and Vere-Jones, David , year =
Daley, Daryl J. and Vere-Jones, David , year =. An introduction to the theory of point processes. 1:
-
[30]
Non-uniform random variate generation , isbn =
Devroye, Luc , year =. Non-uniform random variate generation , isbn =
-
[31]
Hawkes, Alan G. , year =. Spectra of. Biometrika , publisher =. doi:10.2307/2334319 , abstract =
-
[32]
Hawkes, Alan G. and Oakes, David , year =. A. Journal of Applied Probability , publisher =. doi:10.2307/3212693 , abstract =
-
[33]
Long-term seizure diary tracking habits in clinical studies:. Epilepsy Research , author =. 2024 , pages =. doi:10.1016/j.eplepsyres.2024.107379 , abstract =
-
[34]
Epileptic Disorders , author =
Electronic seizure diaries for clinical care and research , volume =. Epileptic Disorders , author =. 2022 , pages =. doi:10.1684/epd.2022.1451 , abstract =
-
[35]
Pénichoux, Juliette and Moreau, Thierry and Latouche, Aurélien , month = mar, year =. Simulating recurrent events that mimic actual data: a review of the literature with emphasis on event-dependence , shorttitle =
-
[36]
Pediatric Neurology , author =
Responsive. Pediatric Neurology , author =. 2023 , pages =. doi:10.1016/j.pediatrneurol.2023.03.001 , language =
-
[37]
Seizure: European Journal of Epilepsy , author =
Epilepsy surgery in children with genetic etiologies:. Seizure: European Journal of Epilepsy , author =. 2023 , pages =. doi:10.1016/j.seizure.2023.10.011 , abstract =
-
[38]
Leveraging expertise and optimizing clinical research:. Epilepsia , author =. 2023 , pages =. doi:10.1111/epi.17579 , abstract =
-
[39]
Surgical evaluation in children. Epilepsia , author =. 2022 , pages =. doi:10.1111/epi.17124 , abstract =
-
[40]
Statistics in Medicine , author =
Generating survival times to simulate. Statistics in Medicine , author =. 2005 , pages =. doi:10.1002/sim.2059 , abstract =
-
[41]
Lim, Kar Wai and Lee, Young and Hanlen, Leif and Zhao, Hongbiao , month = mar, year =. Simulation and
-
[42]
and Lee, Young and Taimre, Thomas , year =
Laub, Patrick J. and Lee, Young and Taimre, Thomas , year =. The. doi:10.1007/978-3-030-84639-8 , language =
-
[43]
Juarez-Colunga, E and Dean, C B , file =
-
[44]
Cook, Richard J. and Lawless, Jerald F. , year =. The statistical analysis of recurrent events , isbn =
-
[45]
Post-traumatic epilepsy: an overview , volume =. Therapy , author =. 2010 , pages =. doi:10.2217/thy.10.57 , abstract =
-
[46]
Standards for epidemiologic studies and surveillance of epilepsy , volume =. Epilepsia , author =. 2011 , pages =. doi:10.1111/j.1528-1167.2011.03121.x , abstract =
-
[47]
Methodology and Computing in Applied Probability , author =
Bayesian. Methodology and Computing in Applied Probability , author =. 2013 , pages =. doi:10.1007/s11009-011-9272-5 , abstract =
-
[48]
Management of post‐traumatic epilepsy:. Epilepsia Open , author =. 2017 , pages =. doi:10.1002/epi4.12049 , abstract =
-
[49]
Obral, Katarzyna , file =. A
-
[50]
Methodology and Computing in Applied Probability , author =
Approximate. Methodology and Computing in Applied Probability , author =. 2006 , pages =. doi:10.1007/s11009-006-7288-z , abstract =
-
[51]
Advances in Applied Probability , author =
Perfect simulation of. Advances in Applied Probability , author =. 2005 , pages =. doi:10.1239/aap/1127483739 , abstract =
- [52]
-
[53]
Methodology and Computing in Applied Probability , author =
Hawkes. Methodology and Computing in Applied Probability , author =. 2022 , pages =. doi:10.1007/s11009-022-09938-1 , abstract =
-
[54]
and Wong, Victoria and Nurse, Ewan S
Hannon, Timothy and Fernandes, Kiran M. and Wong, Victoria and Nurse, Ewan S. and Cook, Mark J. , month = mar, year =. Over- and underreporting of seizures:. Epilepsia , publisher =. doi:10.1111/epi.17930 , abstract =
-
[55]
Rasmussen, Jakob Gulddahl , month = jun, year =. Lecture
-
[56]
Electronic Communications in Probability , author =
Exact simulation of. Electronic Communications in Probability , author =. 2013 , file =. doi:10.1214/ECP.v18-2717 , abstract =
-
[57]
Tracking epilepsy with an electronic diary , volume =. Acta Paediatrica , author =. 2010 , pages =. doi:10.1111/j.1651-2227.2010.01694.x , language =
-
[58]
Focal nonmotor versus motor seizures:. Epilepsia , author =. 2020 , pages =. doi:10.1111/epi.16707 , abstract =
-
[59]
Fifty years later: new directions in. SORT. Statistics and Operations Research Transactions , author =. 2022 , pages =. doi:10.2436/20.8080.02.116 , abstract =
-
[60]
Choi, Edward and Du, Nan and Chen, Robert and Song, Le and Sun, Jimeng , month = nov, year =. Constructing. 2015. doi:10.1109/ICDM.2015.144 , abstract =
-
[61]
Journal of the American Statistical Association , author =
Self-. Journal of the American Statistical Association , author =. 2011 , pages =. doi:10.1198/jasa.2011.ap09546 , language =
-
[62]
International Journal of Forecasting , author =
Hawkes process modeling of. International Journal of Forecasting , author =. 2022 , pages =. doi:10.1016/j.ijforecast.2021.07.001 , abstract =
-
[63]
Epilepsy in adults , volume =. The Lancet , author =. 2019 , pages =. doi:10.1016/S0140-6736(18)32596-0 , language =
-
[64]
The burden of epilepsy and unmet need in people with focal seizures , volume =. Brain and Behavior , author =. 2022 , pages =. doi:10.1002/brb3.2589 , abstract =
-
[65]
The prevalence of epilepsy and pharmacoresistant epilepsy in adults:. Epilepsia , author =. 2008 , pages =. doi:10.1111/j.1528-1167.2008.01579.x , abstract =
-
[66]
Neurobiology of Disease , author =
Long-term seizure dynamics are determined by the nature of seizures and the mutual interactions between them , volume =. Neurobiology of Disease , author =. 2021 , pages =. doi:10.1016/j.nbd.2021.105347 , abstract =
-
[67]
Nature Communications , author =
Critical slowing down as a biomarker for seizure susceptibility , volume =. Nature Communications , author =. 2020 , pages =. doi:10.1038/s41467-020-15908-3 , abstract =
-
[68]
Nature Communications , author =
Multi-day rhythms modulate seizure risk in epilepsy , volume =. Nature Communications , author =. 2018 , pages =. doi:10.1038/s41467-017-02577-y , abstract =
-
[69]
Disease Models & Mechanisms , author =
Semiology, clustering, periodicity and natural history of seizures in an experimental occipital cortical epilepsy model , copyright =. Disease Models & Mechanisms , author =. 2018 , pages =. doi:10.1242/dmm.036194 , abstract =
-
[70]
Interictal spikes and epileptic seizures: their relationship and underlying rhythmicity , volume =. Brain , author =. 2016 , pages =. doi:10.1093/brain/aww019 , language =
-
[71]
Seizure. Epilepsia , author =. 2005 , pages =. doi:10.1111/j.0013-9580.2005.29004.x , abstract =
-
[72]
Frontiers in Neurology , author =
The. Frontiers in Neurology , author =. 2014 , file =. doi:10.3389/fneur.2014.00217 , abstract =
-
[73]
Experimental Neurology , author =
Endogenous multidien rhythm of epilepsy in rats , volume =. Experimental Neurology , author =. 2019 , pages =. doi:10.1016/j.expneurol.2019.02.006 , abstract =
-
[74]
Seizure. Epilepsia , author =. 1991 , pages =. doi:10.1111/j.1528-1157.1991.tb04703.x , abstract =
-
[75]
Ogata, Yosihiko , year =. Statistical. Journal of the American Statistical Association , publisher =. doi:10.2307/2288914 , abstract =
-
[76]
Hawkes process model with a time-dependent background rate and its application to high-frequency financial data , volume =. Physical Review E , author =. 2017 , keywords =. doi:10.1103/PhysRevE.96.012303 , abstract =
-
[77]
Journal of Quantitative Criminology , author =
Terrorism. Journal of Quantitative Criminology , author =. 2013 , pages =. doi:10.1007/s10940-012-9181-y , abstract =
-
[78]
SN Computer Science , author =
Sports. SN Computer Science , author =. 2025 , pages =. doi:10.1007/s42979-025-03939-w , abstract =
-
[79]
Real-time analysis of neuronal firing patterns via
Gerencser, Laszlo and Perczel, Gyorgy , month = oct, year =. Real-time analysis of neuronal firing patterns via. 2021 55th. doi:10.1109/IEEECONF53345.2021.9723129 , abstract =
-
[80]
Mohler, George and Raje, Rajeev and Carter, Jeremy and Valasik, Matthew and Brantingham, Jeffrey , month = oct, year =. A. 2018. doi:10.1109/SMC.2018.00421 , abstract =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.