Skip to content
DraconiticResearch
Search Sources API
197,768 studies31 sourcesUpdated 58 min ago
PreprintarXiv licencearXiv

Graphon Hawkes Processes: Sparse Estimation and Network Recovery

Abstract

Learning interaction structures from spatiotemporal event data is a key statistical problem in fields such as neuroscience, epidemiology, and seismology. We study this problem for Hawkes processes on a continuous spatial domain, where events may trigger future events at different locations. Following the framework developed in ``Spatiotemporal Hawkes processes with a graphon-induced connectivity structure'' by Baars, Laeven, and Mandjes (2026), we model the spatial interaction mechanism using an underlying graphon that characterizes how the intensities' excitation strength varies across locations. By expanding the baseline intensity, temporal excitation kernel and graphon in orthogonal bases, we represent the problem as a high-dimensional parametric inference problem. We propose a sparse penalized least-squares estimator that recovers the effective interaction structure while controlling model complexity. Our analysis establishes non-asymptotic error bounds in a high-dimensional regime where the number of parameters may grow with the observation horizon, without requiring the conditional intensity to be almost surely bounded. We further establish sparsistency results, showing that the latent interaction network can be consistently recovered. Simulation studies demonstrate the good finite sample performance of the proposed method and an application to large-scale neuronal spike data, comprising 1.81 million events, showcases its practical effectiveness.

Related

  1. Parametric estimation of Hawkes processes based on ordinary least squares
    Benjamin Poignard, Yoann Potiron
    arXivPreprintarXiv licenceFull textStatisticsarXiv
  2. nethist: An R package for Nonparametric Graphon Estimation via Network Histograms
    Youngseok Song, Sofia C. Olhede
    arXivPreprintCC BY 4.0Full textStatisticsarXiv
  3. Gaussian approximation for the integrated intensity of counting processes: application to Hawkes processes
    Benjamin Massat
    HALPreprintCC BY 4.0Full textProbabilityHALarXiv
  4. Local stability and rates of convergence to equilibrium for the Nonlinear Renewal Equation; applications to Hawkes processes
    Céline Duval, Eric Luçon
    NonlinearityPeer-reviewedCC BY 4.0Full textDynamical systemsProbabilityHAL
  5. Estimation and Recovery of a Planted Dense Subgraph from a Single Network Cascade
    Maximilien Dreveton, Paula Mürmann, Patrick Thiran
    arXivPreprintCC BY 4.0Full textProbabilityStatisticsarXiv
BibTeX
@misc{baars2026graphon,
  title = {{Graphon Hawkes Processes: Sparse Estimation and Network Recovery}},
  author = {Justin R. Baars and Roger J. A. Laeven and Michel R. H. Mandjes},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2610.16677},
  eprint = {2610.16677},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2610.16677},
}
RIS
TY  - UNPB
TI  - Graphon Hawkes Processes: Sparse Estimation and Network Recovery
AU  - Justin R. Baars
AU  - Roger J. A. Laeven
AU  - Michel R. H. Mandjes
PY  - 2026
DA  - 2026/10/12
PB  - arXiv
DO  - 10.48550/arxiv.2610.16677
UR  - https://arxiv.org/abs/2610.16677
AB  - Learning interaction structures from spatiotemporal event data is a key statistical problem in fields such as neuroscience, epidemiology, and seismology. We study this problem for Hawkes processes on a continuous spatial domain, where events may trigger future events at different locations. Following the framework developed in ``Spatiotemporal Hawkes processes with a graphon-induced connectivity structure'' by Baars, Laeven, and Mandjes (2026), we model the spatial interaction mechanism using an underlying graphon that characterizes how the intensities' excitation strength varies across locations. By expanding the baseline intensity, temporal excitation kernel and graphon in orthogonal bases, we represent the problem as a high-dimensional parametric inference problem. We propose a sparse penalized least-squares estimator that recovers the effective interaction structure while controlling model complexity. Our analysis establishes non-asymptotic error bounds in a high-dimensional regime where the number of parameters may grow with the observation horizon, without requiring the conditional intensity to be almost surely bounded. We further establish sparsistency results, showing that the latent interaction network can be consistently recovered. Simulation studies demonstrate the good finite sample performance of the proposed method and an application to large-scale neuronal spike data, comprising 1.81 million events, showcases its practical effectiveness.
ER  -