REVIEW 1 major objections 3 minor 2 cited by
Stochastic heterogeneous SIR model with infection-age dependent infectivity on large random graphs
T0 review · 1 major / 3 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read The paper proves a functional law of large numbers for stochastic SIR epidemics on large non-homogeneous random graphs, with a limit described by explicit measure-valued equations.
desk verdict A substantial FLLN for non-Markovian SIR on graphons, with a load-bearing but fixable gap in Lemma 5.6 and a restrictive initial-condition assumption that should be flagged. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The proof constructs an auxiliary or intermediate process in which each susceptible individual is infected at rate $F(t,X^N_i)$ using the limiting force of infection rather than the sum over actual infected neighbours in the graph. Because the auxiliary infection times are conditionally independent, the dependence among individuals is decoupled, and the difference between the true and auxiliary processes is controlled by decomposing it into seven terms: age mismatch, kernel approximation error, conditional-mean and variance fluctuations of the contact rates, and errors in the initial measures. Each term is bounded under Assumptions 2.2, 3.3, and 3.5, with a generic weak-convergence result for kernels that are continuous only almost everywhere, and a Gronwall inequality on $\bar D^N(t)=\mathbb E^N_0[N^{-1}\sum_{i}\sup_{r\le t}|D^N_i(r)-\tilde D^N_i(r)|]$ closes the comparison.
What would settle it
Simulate the finite-$N$ model with a fixed a.e.-continuous kernel such as $\bar\omega(x,x')=\mathbf{1}_{\{|x-x'|\le\delta\}}$ and compare the empirical susceptible measure with $\bar\mu^S_t(dx)=\exp(-\int_0^t F(s,x)\,ds)\,\bar\mu^S_0(dx)$; Theorem 3.6 predicts the difference tends to zero in probability. A reproducible simulation showing a nonzero gap that does not shrink with $N$, or an explicit initial distribution satisfying Assumption 3.1 for which (3.6)-(3.7) admits two distinct solutions, would falsify the claim.
Extended reading notes
Core claim
The central result, Theorem 3.6, states that under Assumptions 2.2, 2.4, 2.6, 3.1, 3.3, and 3.5, the LLN-scaled measure-valued processes $(\bar\mu^{S,N},\bar\mu^{I,N},\bar\mu^{R,N})$ converge in probability to $(\bar\mu^S,\bar\mu^I,\bar\mu^R)$ in $D_1\times D_2\times D_1$. The limit is the unique solution of the system (3.6)-(3.9), where the force of infection $F(t,x)$ is given by (3.7) in terms of the limiting graph kernel $\bar\omega$, the average infectivity $\bar\lambda$, and the current susceptible measure. In particular, the susceptible population satisfies $\bar\mu^S_t(dx)=\exp(-\int_0^t F(s,x)\,ds)\,\bar\mu^S_0(dx)$. When the recovery-time distributions have densities, the infected measure solves a linear transport PDE with boundary condition $\tilde\mu^I_t(dx,0)=F(t,x)\bar\mu^S_t(dx)$, so the entire epidemic reduces to a deterministic evolution on the space of types and infection ages.
Load-bearing premise
The argument needs the initial recovered measure to dominate the remaining infectivity carried by initially infected individuals, a compatibility inequality that fails for a fresh outbreak with no recovered people and infected individuals who already have positive infection age.
Editorial extensions
If this is right
- For large $N$, the random epidemic trajectory is essentially deterministic: fluctuations around the limiting PDE vanish in probability, so a single PDE solve replaces many stochastic simulations.
- The limit supports continuous type spaces and kernels that are continuous only almost everywhere, so spatial or age-structured models with discontinuous contact thresholds are included.
- The boundary condition $\tilde\mu^I_t(dx,0)=F(t,x)\bar\mu^S_t(dx)$ links new infections to the current susceptible measure, making the infection-age distribution computable without tracking the whole graph.
- Non-Markovian infectivity profiles, including random individual infection durations with non-exponential laws, obey the same equations, extending earlier Markovian graphon limits.
Reading between the lines
- Editorial inference: the same intermediate-model comparison suggests a natural route to a central limit theorem for these measure-valued processes, because the variance terms in the seven-term decomposition carry the $N^{-1/2}$ fluctuations; the paper does not pursue that step.
- Editorial inference: the aggregate averaging conditions in Assumption 3.5 leave room for localized superspreading on a negligible fraction of edges, so an extension to unbounded kernels would likely require $L^p$ integral conditions rather than uniform boundedness.
- Editorial inference: the explicit exponential formula for $\bar\mu^S_t(dx)$ gives a directly testable prediction for the fraction of each type remaining susceptible, which could be compared with simulations or contact-tracing data.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies an individual-based stochastic SIR model on a large random graph with individual characteristics, random edge weights, and infection-age-dependent infectivity. The main theorem (Theorem 3.6) asserts that, under Assumptions 2.2, 2.4, 2.6, 3.1, 3.3, and 3.5, the normalized measure-valued processes (μ̄S,N, μ̄I,N, μ̄R,N) converge in probability to a deterministic limit characterized by the measure-valued equations (3.6)-(3.9), with an additional PDE-on-graphon representation. The proof constructs an auxiliary mean-field process driven by the limiting force of infection, decomposes the discrepancy between the original and auxiliary processes into seven terms, and controls them through a series of lemmas and a Gronwall argument. Appendix A supplies technical supporting results, including a generic convergence lemma (Proposition 5.9) of independent interest.
Significance. If the result holds, this is a substantial extension of functional law of large numbers results for SIR epidemics: it combines non-Markovian infectivity, individual heterogeneity, and non-homogeneous random graph connectivity in a single graphon limit. The paper is methodical and self-contained, and the generic convergence result Proposition 5.9 is of independent interest. There is no parameter fitting or circular prediction; the deterministic limit is derived from explicit probabilistic assumptions. The main proof has a localized algebraic gap in the error bounds of Lemmas 5.6-5.7 that should be corrected before publication.
major comments (1)
- [Section 5.1, Lemma 5.6 and Lemma 5.7, Eq. (5.27) vs. (5.34)] The proof of Lemma 5.6 establishes the bound U^N ≤ λ* sqrt(Υ_N + ω* γ̄_N), but the displayed V^N in (5.27) is λ* sqrt(Υ_N) + ω* γ̄_N. The inequality λ* sqrt(Υ_N + ω* γ̄_N) ≤ λ* sqrt(Υ_N) + ω* γ̄_N is not true in general; for instance, with Υ_N = 0 and γ̄_N = 1/N, the left side is λ* sqrt(ω*/N) while the right side is ω*/N. Lemma 5.7 contains the same mismatch, and Proposition 5.11 and Proposition 5.14 rely on V^N as a vanishing error term in the Gronwall estimate. The proof as written therefore has a genuine algebraic gap. The gap is local and repairable: replacing V^N by λ* sqrt(Υ_N + ω* γ̄_N), or by λ*(sqrt(Υ_N) + sqrt(ω* γ̄_N)), preserves the required convergence to zero under Assumption 3.5.
minor comments (3)
- [Assumption 3.1 and Remark 3.2] The compatibility inequality (3.2) is a substantive restriction on admissible initial conditions; it fails when μ̄R_0 = 0 and μ̄I_0 has positive mass on ages a with F_x(a) > 0. Since this inequality is used in Proposition 4.1 to bound the force of infection and close the uniqueness argument, the authors should state this restriction prominently as a modeling limitation or self-consistency condition, rather than only providing the heuristic in Remark 3.2.
- [Proof of Proposition 5.16] In the paragraph following (5.60), the claim that the right-hand side converges by 'Proposition 4.1' should refer to Proposition 5.11, since it is Proposition 5.11 that establishes the convergence of D^N(t).
- [Appendix A, Lemma A.1 and references] The citation 'Pólya and Szegő [?]' contains a placeholder and should be completed. Additionally, reference [23] appears poorly formatted, and the text contains several OCR-like artifacts (e.g., '/BD' in place of indicator functions, doubled letters, and missing accents); these should be cleaned up in the final version.
Circularity Check
No significant circularity: the FLLN is derived from explicit probabilistic assumptions and a self-contained coupling argument.
full rationale
The paper's derivation chain is a standard FLLN proof: define the individual-based SIR dynamics (Section 2), state LLN-scaled measure-valued processes, and prove convergence to the unique solution of the limit system (3.6)-(3.9). The limit F is not fitted from the finite-N process; it is defined as the solution of an integral equation (3.7) from the limiting initial measures and kernels. The auxiliary model in (5.1)-(5.2) uses this limiting F, but only as a proof device: Lemma 5.1 verifies convergence of the auxiliary process, Lemmas 5.6-5.10 bound the seven error terms in Lemma 5.4, and Proposition 5.11 uses Gronwall's inequality to show the difference D^N(t) between original and auxiliary processes vanishes. No fitted parameter is renamed as a prediction, and no quantity is defined in terms of the result it is supposed to establish. Uniqueness is proved directly in Proposition 4.1 rather than imported from prior work. Self-citations (e.g., [20], [27], [46]) are used for model context and as background for the intermediate-model technique, but the convergence argument is self-contained and does not rely on a self-citation as its load-bearing step. The compatibility inequality in Assumption 3.1 is an explicit hypothesis on the initial condition that bounds the force of infection; it restricts the class of admissible epidemics but does not make the theorem circular. A possible algebraic gap in Lemma 5.6's bound V^N is a correctness issue, not a circularity issue. Overall the central claim has independent content and is not forced by definition.
Assumptions & free parameters
assumptions (6)
- domain assumption Assumption 2.2: uniform bound λ* on all infectivity functions λ_i^N(a).
- domain assumption Assumption 2.4 and 2.6: existence of duration distribution F_x and mean infectivity λ̄ with χ = λ̄/F^c bounded by λ*.
- domain assumption Assumption 3.1: weak convergence of initial measures plus the compatibility inequality (3.2).
- domain assumption Assumption 3.3: ω̄_N = N κ_N γ_N → ω̄ uniformly, with ω̄ bounded and μ̄X^2-a.e. continuous, plus continuity of F_x and λ̄.
- domain assumption Assumption 3.5: γ̄_N and Υ_N converge to zero.
- standard math Standard weak convergence, Skorokhod, Portmanteau, and tightness criteria, including Perkins' Theorem II.4.1 (Proposition 5.12).
Cite this review
Pith. "Pith review of Stochastic heterogeneous SIR model with infection-age dependent infectivity on large random graphs." pith.science (2026). https://pith.science/paper/JSSPT5L7
@misc{pith2026250204225,
author = {Pith},
title = {Pith review of: Stochastic heterogeneous SIR model with infection-age dependent infectivity on large random graphs},
year = {2026},
howpublished = {\url{https://pith.science/paper/JSSPT5L7}},
note = {Machine review of arXiv:2502.04225}
}
read the original abstract
We study an individual-based stochastic SIR epidemic model with infection-age dependent infectivity on a large random graph, capturing individual heterogeneity and non-homogeneous connectivity. Each individual is associated with particular characteristics (for example, spatial location and age structure), which may not be i.i.d., and is represented by a particular node. The connectivities among the individuals are given by a non-homogeneous random graph, whose connecting probabilities may depend on the individual characteristics of the edge. To each individual is associated a random infectivity function of its infection age, which is allowed to depend upon the individual characteristics. We use measure-valued processes to describe the epidemic evolution dynamics, tracking the infection age of all individuals, and their associated characteristics. We consider the epidemic dynamics as the population size grows to infinity under a specific scaling of the connectivity graph related to the convergence to a graphon. In the limit, we obtain a system of measure-valued equations, which can be also represented as a PDE model on graphon, and reflects the heterogeneities in individual characteristics and social connectivity.
Forward citations
Cited by 2 Pith papers
-
Functional Central Limit Theorem and SPDE for epidemic model with memory of the last infection and waning immunity
A functional central limit theorem and SPDE are derived for fluctuations of a structured epidemic model with memory of the last infection and waning immunity.
-
Information Spreading in Random Graphs Evoving by Norros-Reittu Model
The paper attempts to derive finite-time message-spreading distributions for Norros-Reittu graphs, but the Poisson-binomial derivation is invalid because success probabilities depend on the random set of informed nodes.
Reference graph
Works this paper leans on
- [20]
-
[37]
D. Keliger, I. Horv´ ath, and B. Tak´ acs. Local-densitydependent Markov processes on graphons with epidemiological applications. Stoch. Proc. Appl. , 148:324–352, 2022
work page 2022
- [1]
-
[2]
F. Ball and T. Britton. Epidemics on networks with preven tive rewiring. Random Structures & Algo- rithms, 61(2):250–297, 2022
work page 2022
-
[3]
F. Ball, T. Britton, K. Y. Leung, and D. Sirl. A stochastic SIR network epidemic model with preventive dropping of edges. Journal of Mathematical Biology , 78(6):1875–1951, 2019
work page 1951
-
[4]
F. Ball and D. Sirl. An SIR epidemic model on a population w ith random network and household structure, and several types of individuals. Advances in Applied Probability , 44(1):63–86, 2012
work page 2012
-
[5]
F. Ball and D. Sirl. Acquaintance vaccination in an epide mic on a random graph with specified degree distribution. Journal of Applied Probability , 50(4):1147–1168, 2013
work page 2013
-
[6]
F. Ball and D. Sirl. Evaluation of vaccination strategie s for SIR epidemics on random networks incor- porating household structure. Journal of Mathematical Biology , 76(1):483–530, 2018
work page 2018
Show all 55 references
-
[7]
F. Ball, D. Sirl, and P. Trapman. Threshold behaviour and final outcome of an epidemic on a random network with household structure. Advances in Applied Probability , 41(3):765–796, 2009
2009
-
[8]
Billingsley
P. Billingsley. Convergence of Probability Measures . Wiley Series in Probability and Statistics. John Wiley & Sons, Inc., New York, second edition edition, 1999. 42 GUODONG PANG, ´ETIENNE PARDOUX, AND AUR ´ELIEN VELLERET
1999
-
[9]
Bohman and M
T. Bohman and M. Picollelli. SIR epidemics on random grap hs with a fixed degree sequence. Random Structures & Algorithms , 41(2):179–214, 2012
2012
-
[10]
Borgs, K
C. Borgs, K. Huang, and C. Ikeokwu. A law of large numbers for SIR on the stochastic block model: A proof via herd immunity. arXiv preprint arXiv:2410.07097 , 2024
2024 arXiv
-
[11]
Britton, F
T. Britton, F. Ball, and P. Trapman. A mathematical mode l reveals the influence of population hetero- geneity on herd immunity to SARS-CoV-2. Science, 369(6505):846–849, 2020
2020
-
[12]
Britton, M
T. Britton, M. Deijfen, and F. Liljeros. A weighted confi guration model and inhomogeneous epidemics. Journal of Statistical Physics , 145:1368–1384, 2011
2011
-
[13]
Britton and D
T. Britton and D. Lindenstrand. Inhomogeneous epidemi cs on weighted networks. Mathematical Bio- sciences, 240(2):124–131, 2012
2012
-
[14]
W. Chen, Y. Hou, and D. Yao. SIR epidemics on evolving Erd ¨ os R´ enyi graphs.arXiv preprint arXiv:2208.11923, 2022
2022 arXiv
-
[15]
Clancy Jr
D. Clancy Jr. Epidemics on critical random graphs with h eavy-tailed degree distribution. Stochastic Processes and their Applications , 179:104510, 2025
2025
-
[16]
D. L. Cohn. Measure Theory: Second Edition . Birkh¨ auser Advanced Texts Basler Lehrb¨ ucher. Springer, New York, NY, 2013
2013
-
[17]
Decreusefond, J.-S
L. Decreusefond, J.-S. Dhersin, P. Moyal, and V. C. Tran . Large graph limit for an SIR process in random network with heterogeneous connectivity. Ann. Appl. Probab. , 22(2):541–575, 2012
2012
-
[18]
M. Deijfen. Epidemics and vaccination on weighted grap hs. Mathematical Biosciences , 232(1):57–65, 2011
2011
-
[19]
Delmas, D
J.-F. Delmas, D. Dronnier, and P.-A. Zitt. An infinite-d imensional metapopulation SIS model. J. Differ. Equ., 313:1–53, 2022
2022
-
[21]
Duchamps, F
J.-J. Duchamps, F. Foutel-Rodier, and E. Schertzer. Ge neral epidemiological models: Law of large numbers and contact tracing. Electronic Journal of Probability , 28:1–37, 2023
2023
-
[22]
Durrett and D
R. Durrett and D. Yao. Susceptible–infected epidemics on evolving graphs. Electronic Journal of Prob- ability, 27:1–66, 2022
2022
-
[23]
A. S. Elie B, Selinger C. The source of individual hetero geneity shapes infectious disease outbreaks. Proc Biol Sci. , 289: 20220232(1974), May 2022
1974
-
[24]
Fernley, P
J. Fernley, P. M¨ orters, and M. Ortgiese. The contact process on a graph adapting to the infection. arXiv preprint arXiv:2312.06251, 2023
2023 arXiv
-
[25]
E. J. Ferreyra, M. Jonckheere, and J. P. Pinasco. SIR dynamics with vaccination in a large configuration model. Applied Mathematics & Optimization , 84:1769–1818, 2021
2021
-
[26]
Forien, G
R. Forien, G. Pang, and ´E. Pardoux. Epidemic models with varying infectivity. SIAM Journal on Applied Mathematics, 81(5):1893–1930, 2021
1930
-
[27]
Forien, G
R. Forien, G. Pang, and ´E. Pardoux. Multi-patch multi-group epidemic model with va rying infectivity. Probability, Uncertainty and Quantitative Risk , 7(4):333–364, 2022
2022
-
[28]
Forien, G
R. Forien, G. Pang, and ´E. Pardoux. Recent advances in epidemic modeling: Non-Mark ov stochastic models and their scaling limits. The Graduate Journal of Mathematics , 7(2):19–75, 2022
2022
-
[29]
Huang and A
Y. Huang and A. R¨ ollin. The SIR epidemic on a dynamic Erd ¨ os R´ enyi random graph.arXiv preprint arXiv:2404.12566, 2024
2024 arXiv
-
[30]
K. A. Jacobsen, M. G. Burch, J. H. Tien, and G. A. Rempa/suppress la. The large graph limit of a stochastic epidemic model on a dynamic multilayer network. Journal of Biological Dynamics , 12(1):746–788, 2018
2018
-
[31]
Jacod and A
J. Jacod and A. Shiryaev. Limit Theorems for Stochastic Processes . Springer-Verlag, Berlin, 2 edition, 2003
2003
-
[32]
Janson, M
S. Janson, M. Luczak, and P. Windridge. Law of large numbers for the SIR epidemic on a random graph with given degrees. Random Structures & Algorithms , 45(4):726–763, 2014
2014
-
[33]
Janson, M
S. Janson, M. Luczak, P. Windridge, and T. House. Near-c ritical SIR epidemic on a random graph with given degrees. Journal of Mathematical Biology , 74:843–886, 2017
2017
-
[34]
Jiang, R
Y. Jiang, R. Kassem, G. York, M. Junge, and R. Durrett. SI R epidemics on evolving graphs. arXiv preprint arXiv:1901.06568, 2019
1901 arXiv
-
[35]
Kallenberg
O. Kallenberg. Foundations of Modern Probability. Probability and Its Applications. Springer New York, New York, NY, 2002. 43
2002
-
[36]
Kanga and ´E
A. Kanga and ´E. Pardoux. Spatial SIR epidemic model with varying infecti vity without movement of individuals: Law of Large Numbers. arXiv preprint arXiv:2412.01673 , 2024
2024
-
[38]
W. R. Khudabukhsh, C. Woroszylo, G. A. Rempa/suppress la, and H. Koeppl. A functional central limit theorem for SI processes on configuration model graphs. Advances in Applied Probability , 54(3):880–912, 2022
2022
-
[39]
M. Kubasch. Large population limit for a multilayer SIR model including households and workplaces. arXiv preprint arXiv:2305.17064 , 2023
2023 arXiv
-
[40]
A. A. Lashari, A. Serafimovi´ c, and P. Trapman. The durat ion of a supercritical SIR epidemic on a configuration model. Electronic Journal of Probability , 26:1–49, 2021
2021
-
[41]
Lloyd-Smith, S
J. Lloyd-Smith, S. Schreiber, P. Kopp, and W. Getz. Supe rspreading and the effect of individual varia- tion on disease emergence. Nature, 438(7066):355—359, November 2005
2005
-
[42]
M. Luczak. A note on the Markovian SIR epidemic on a rando m graph with given degrees. Journal of Theoretical Probability, 37(2):1039–1051, 2024
2024
-
[43]
Mossong, N
J. Mossong, N. Hens, M. Jit, P. Beutels, K. Auranen, R. Mi kolajczyk, M. Massari, S. Salmaso, G. S. Tomba, J. Wallinga, J. Heijne, M. Sadkowska-Todys, M. Rosinska, and W. J. Edmunds. Social contacts and mixing patterns relevant to the spread of infectious dis eases. PLOS Medic...
2008
-
[44]
Pang and ´E
G. Pang and ´E. Pardoux. Functional limit theorems for non-Markovian ep idemic models. The Annals of Applied Probability , 32(3):1615–1665, 2022
2022
-
[45]
Pang and E
G. Pang and E. Pardoux. Functional central limit theorems for epidemic models with varying infectivity. Stochastics, 95(5):819–866, 2023
2023
-
[46]
Pang and ´E
G. Pang and ´E. Pardoux. Functional law of large numbers and PDEs for epidemic models with infection- age dependent infectivity. Applied Mathematics and Optimization , 87:Article 50, 2023
2023
-
[47]
Pang and ´E
G. Pang and ´E. Pardoux. SPDE for stochastic SIR epidemic models with inf ection-age dependent infectivity. arXiv:2410.17200, 2024
2024 arXiv
-
[48]
E. A. Perkins. Dawson-Watanabe superprocesses and mes ure-valued diffusions. In Springer, editor, Ecole d’Et´ e de Probabilit´ es de Saint-Flour, volume 1781 of Lectures Notes in Math. , pages 125–329, New York, 1999
1999
-
[49]
G. Reinert. The asymptotic evolution of the general sto chastic epidemic. The Annals of Applied Proba- bility, 5(4):1061–1086, 1995
1995
-
[50]
Y. Shang. Mixed SI (R) epidemic dynamics in random graph s with general degree distributions. Applied Mathematics and Computation , 219(10):5042–5048, 2013
2013
-
[51]
Spricer and T
K. Spricer and T. Britton. An SIR epidemic on a weighted n etwork. Network Science , 7(4):556–580, 2019
2019
-
[52]
Sznitman
A.-S. Sznitman. Topics in propagation of chaos. In Ecole d’´ et´ e de probabilit´ es de Saint-Flour XIX—1989, pages 165–251. Springer, 1991
1989
-
[53]
E. Volz. SIR dynamics in random networks with heterogen eous connectivity. Journal of mathematical biology, 56:293–310, 2008
2008
-
[54]
F. J. Wang. Limit theorems for age and density dependent stochastic population models. Journal of Mathematical Biology, 2(4):373–400, 1975
1975
-
[55]
F. J. Wang. Gaussian approximation of some closed stoch astic epidemic models. Journal of Applied Probability, 14(2):221–231, 1977. Department of Computational Applied Mathematics and Opera tions Research, George R. Brown School of Engineering and Computing, Rice University , ...
1977
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.