REVIEW 2 major objections 5 minor 1 cited by
A stochastic epidemic model with memory of the last infection and waning immunity
T0 review · 2 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read An epidemic with waning immunity stays endemic precisely when a single computable threshold built from infectivity, susceptibility, and memory of the last infection falls below 1.
desk verdict A careful extension of the FPPZ framework that earns a serious referee; the endemicity threshold is real, and the main scope limit (no reinfection while infectious) is explicit and acknowledged. 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 machinery is a piecewise deterministic Markov process on (age, trait), with infection rate $F(t)\gamma(a,\theta)$, age reset to zero at infection, and trait jumping according to the kernel $K$; its large-population limit is a nonlinear transport PDE whose boundary condition feeds reinfections back in. The threshold is carried by the memory operator $T(B)(\theta)=\int_{\Theta} K(\tilde{\theta},\theta)B(\tilde{\theta})\nu(d\tilde{\theta})$: its Perron eigenfunction $S^*$ is the stationary trait distribution generated by reinfections, and the function $H(x)=x\int_{\mathbb{R}_+\times\Theta} \exp(-x\int_0^a \gamma(s,\theta)ds)\, S^*(\theta)\, da\, \nu(d\theta)$ satisfies $H(0)=\int_\Theta S^*/\gamma^* \, d\nu$, so solving $H(x)=1$ yields the equilibrium force of infection.
What would settle it
Drop the separation assumption $\lambda\gamma\equiv 0$ while keeping everything else, and numerically solve the stationary equations (5.1) for a two-trait example: if an endemic equilibrium appears while the computed threshold (3.9) is above 1, or fails to appear while the threshold is below 1, then the separation assumption is load-bearing and the threshold claim does not extend beyond it.
Extended reading notes
Core claim
The central claim is Theorem 3.6: under assumptions that force susceptibility to vanish while an individual is infectious and to have a well-defined long-run Cesàro limit $\gamma^*(\theta)$, the PDE limit has an endemic equilibrium when $\int_{\Theta} \frac{1}{\gamma^*(\theta)} S^*(\theta) \nu(d\theta) < 1$, where $S^*$ is the unique positive eigenfunction of the memory operator $T$ with eigenvalue 1, normalized so that $\int_{\mathbb{R}_+ \times \Theta} \lambda(a,\theta) S^*(\theta) \, da\, \nu(d\theta) = 1$. Equivalently, after reweighting the trait distribution by $S^*$, the average number of infections produced in a fully susceptible population exceeds $E^*_\nu[1/\gamma^*]$. Under a monotonicity condition on susceptibility, the same inequality is necessary and sufficient; when it fails, only the disease-free equilibrium exists. The paper also proves local stability of the endemic equilibrium for a memory-free SIS-type model with bounded infectivity duration, a step that earlier approaches could not reach.
Load-bearing premise
The load-bearing premise is that no one can be reinfected while still infectious, so an infected person's susceptibility is zero during the entire infectious period; if infectivity and susceptibility overlap, the stationary state constructed in the proof no longer satisfies the equilibrium equations and the threshold is not proved.
Editorial extensions
If this is right
- Under the theorem's assumptions, the classical $R_0>1$ condition is replaced by $E^*_\nu[1/\gamma^*] < R^*_0$, with the memory kernel entering only through the reweighted trait distribution $S^*\nu$.
- With susceptibility non-decreasing in age, the threshold is necessary and sufficient: the disease-free state is the only equilibrium above it.
- Vaccination enters through susceptibility curves; for a renewal vaccination scheme the paper recovers the previously known threshold without assuming infectivity and susceptibility are independent.
- For the one-shot-vaccine model there are parameter regions with two endemic equilibria, so the same model can admit both a mild and a severe persistent state.
- Local stability is proved for a memory-free SIS-type model with bounded infectivity duration, under the explicit parameter condition $\lambda_* \le 2\rho e^{\rho a_*}$.
Reading between the lines
- Editorial inference: the threshold formula suggests that memory affects endemicity only through the stationary trait distribution $S^*$, so transient correlations between consecutive infections should not move the endemicity boundary, only the shape of $S^*$.
- Editorial inference: if infectivity and susceptibility overlap in real pathogens, the separation assumption fails; a natural numerical test is to solve the stationary equations without $\lambda\gamma\equiv 0$ and see whether a modified threshold emerges.
- Editorial inference: when the auxiliary function $H$ is non-monotone, the model can have two endemic equilibria, which hints at hysteresis—temporarily lowering the force of infection could push a population from the higher persistent state to the disease-free state.
- Editorial inference: the stability proof is limited to bounded infectivity duration, so a plausible extension is to analyze the eigenvalue equation (5.29) for unbounded durations using Laplace-transform sign arguments.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces a stochastic individual-based epidemic model with age since last infection and a trait parameter that carries memory of the last infection. At each reinfection, the trait is redrawn according to a Markov kernel K, and the infectivity and susceptibility curves depend on age and trait. The authors prove a functional law of large numbers (Theorem 3.1) showing that the empirical measure of ages and traits converges to the unique solution of a nonlinear age-and-trait structured PDE. The main result is Theorem 3.6: under Assumptions A1 and A2, an endemic equilibrium exists exactly when the threshold condition ∫Θ (1/γ*(θ)) S*(θ) ν(dθ) < 1 holds, where S* is the positive eigenfunction of the memory kernel's integral operator normalized by ∫ λ S* = 1, and γ* is the long-run Cesàro limit of susceptibility. Uniqueness is obtained under a monotonicity condition (3.10). The paper also analyzes local stability in the memory-free case, conditional on a spectral condition (5.29), and verifies it for a step-susceptibility SIS-type model. Applications to vaccination policies are given, including a renewal-process vaccination model for which the threshold of [12] is recovered.
Significance. If correct, the paper provides a genuine extension of the Forien--Pang--Pardoux--Zotsa framework by allowing the new infectivity and susceptibility curves to depend on the previous infection via a trait Markov chain. The threshold is parameter-free and explicit: it is expressed in terms of the eigenfunction of the memory kernel, the susceptibility long-run limit, and the infectivity curves. The FLLN proof follows a standard tightness-and-martingale argument and is coherent. The equilibrium construction is algebraically consistent under A1 and A2, and the examples, including the vaccination models, are concrete and give interpretable thresholds. The paper also delivers, to my knowledge, a nontrivial local stability result for a memory-free SIS-type model, albeit only under an explicit no-eigenvalue condition.
major comments (2)
- [Theorem 3.6, Eq. (5.11)] The endemicity threshold is established under Assumption A2(1)(c), which imposes that the support of λ precedes the support of γ and that λγ ≡ 0. The proof uses Eq. (5.11), where the exponential waning factor is dropped on the support of λ. If the infectivity and susceptibility supports overlap, the exponential factor is strictly less than 1 on part of the λ-support, and the proposed stationary state does not satisfy the force-of-infection equation. The assumption is explicit and is consistent with [11,12], so this is a scoping limitation rather than an internal contradiction; however, the abstract and introduction should state prominently that the threshold is proved only in the no-reinfection-while-infectious regime.
- [Theorem 3.8, Eq. (5.29)] The local stability result is conditional on the non-existence of eigenvalues of the linearized operator with real part ≥ 0, as stated in Eq. (5.29). For the general memory-free model this condition is not verified; it is checked only for the special SIS-type model in Proposition 5.9. The text acknowledges this, but the statement of Theorem 3.8 should be more explicit that it is a conditional stability criterion for general curves, with a concrete verification in one model class, rather than a fully general stability theorem.
minor comments (5)
- [Section 2 and Eq. (4.1)] The notation Qk is used both for the original Poisson measure and for its compensated version; please distinguish the two, for example by writing Q̃k for the compensated measure.
- [Proposition 5.3, Eq. (5.7)] The normalization condition ∫ λ S* = 1 requires that ∫ λ S* > 0; the case λ ≡ 0, which would make the threshold ill-defined, should be excluded or explicitly discussed.
- [Section 6.2] The class of infectivity curves λ used in the renewal vaccination model is not specified precisely; please state the assumptions on λ needed both for Assumption A2 and for the definition of R0.
- [Throughout] The manuscript contains numerous typographical and OCR-related artifacts, for example 'memor y', 'W ANING', 'd/greaterorequalslant1', and broken mathematical spacing. The final version should be carefully proofread.
- [Remark 3.7] The probability measure E*ν is introduced only after the condition (3.9); it may be clearer to define E*ν before rewriting the threshold as E*ν[1/γ*] < R*0.
Circularity Check
No significant circularity: the endemicity threshold is derived from the model's Perron eigenfunction and Cesàro susceptibility limit, and the authors' earlier papers are used only as baseline comparisons.
full rationale
The paper's central claim, Theorem 3.6, is not circular. The equilibrium is constructed from the stationary equations (5.1), which are reduced to the system (5.3)-(5.5). The function S* is defined as the normalized positive eigenfunction of the integral operator T from the memory kernel K, using the external spectral theorem of Schaefer [30], not from the desired threshold. The function H in (5.8)-(5.10) is computed from S* and the susceptibility curves, and its value at zero is obtained from the Cesàro limit in Assumption A2-(2), giving H(0)=∫(1/γ*)S*dν. The threshold (3.9) is therefore an output of the equilibrium equations, not an input. No parameter is fitted to the conclusion, and the change of measure in Remark 3.7 is merely a rewriting of (3.9). Equation (5.11), which drops the exponential waning factor on the support of λ, is a direct consequence of the explicitly stated disjoint-support assumption A2-(1)(c), not of a prior result by the same authors. The FLLN, Theorem 3.1, is proved in Section 4 by tightness, uniqueness, and martingale estimates, not by citation. The local stability result, Theorem 3.8, is stated conditionally on the spectral condition (5.29), and the paper explicitly acknowledges in Section 5.2 that a complete stability proof for the memory model remains open; this is an honest scope limitation rather than a circular reduction. The self-citations [11] and [12] are used as baselines: Remark 3.4 notes that the memory-free system coincides with [11], Example 5.2 recovers the [11] threshold after independent computation, and Section 6.2 recovers the [12] threshold by direct calculation. These prior-author results are not load-bearing assumptions in the proofs; they are confirmatory comparisons. Thus no circularity step can be exhibited from the paper's own equations.
Assumptions & free parameters
assumptions (9)
- domain assumption Assumption A1: K is a measurable kernel with ∫Θ K(θ, dθ~)ν(dθ~)=1 for each θ.
- domain assumption Assumption A2-(1c): λγ≡0 and the support of λ lies before the support of γ.
- domain assumption Assumption A2-(2): the Cesàro average of γ has a limit γ*(θ) ν-a.e.
- domain assumption Assumption A2-(3): ess supΘ ∫ exp(-x∫_0^a γ) da < ∞ for every x>0.
- domain assumption Assumption A2-(4): K>0 on Θ² and θ ↦ sup K ∈ L1(ν).
- domain assumption Assumption A3: γ(a,θ) ≥ σ 1_{a>a*} uniformly in θ.
- domain assumption The population is homogeneous, demographic turnover is neglected, and at t=0 all individuals have been infected at least once.
- standard math Standard tools: tightness criteria for measure-valued processes, Gronwall, semigroup perturbation theory, spectral theory of positive integral operators.
- domain assumption For the stability section, Θ is open in Rd and ν is absolutely continuous with respect to Lebesgue measure.
Cite this review
Pith. "Pith review of A stochastic epidemic model with memory of the last infection and waning immunity." pith.science (2026). https://pith.science/paper/GAGQW3OS
@misc{pith2026250500601,
author = {Pith},
title = {Pith review of: A stochastic epidemic model with memory of the last infection and waning immunity},
year = {2026},
howpublished = {\url{https://pith.science/paper/GAGQW3OS}},
note = {Machine review of arXiv:2505.00601}
}
read the original abstract
We adapt the article of Forien, Pang, Pardoux and Zotsa: Arxiv preprint Arxiv2210.04667(2022), on epidemic models with varying infectivity and waning immunity, to incorporate the memory of the last infection. To this end, we introduce a parametric approach and consider a piecewise deterministic Markov process modeling both the evolution of the parameter, also called the trait, and the age of infection of individuals over time. At each new infection, a new trait is randomly chosen for the infected individual according to a Markov kernel, and their age is reset to zero. In the large population limit, we derive a partial differential equation (PDE) that describes the density of traits and ages. The main goal is to study the conditions under which endemic equilibria exist for the deterministic PDE model and to establish an endemicity threshold that depends on the model parameters. The local stability of these equilibria is also analyzed. The endemicity threshold is computed for several examples, including models that incorporate a vaccination policy, and a local stability result is obtained for a memory-free SIS-type model.
Forward citations
Cited by 1 Pith paper
-
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.
Reference graph
Works this paper leans on
- [11]
-
[12]
F. Foutel-Rodier, A. Charpentier, and H. Gu´ erin. Optimal vaccination policy to prevent endemicity: A stochastic model. J. Math. Biol. , 90(10), 2025
work page 2025
- [1]
-
[2]
H. Brezis. Functional analysis, Sobolev spaces and partial differentia l equations. Universitext. Springer, New York, 2011
work page 2011
-
[3]
T. Britton and E. Pardoux eds. Stochastic epidemic models with inference . Springer, 2019
work page 2019
-
[4]
A. Calsina and J. M. Palmada. Steady states of a selection -mutation model for an age structured population. J. Math. Anal. Appl. , 400(2):386–395, 2013
work page 2013
-
[5]
D. D. Chaplin. Overview of the immune response. Journal of allergy and clinical immunology , 125(2):S3–S23, 2010
work page 2010
-
[6]
M. A. Chowdhury, N. Hossain, M. A. Kashem, M. A. Shahid, an d A. Alam. Immune response in covid-19: A review. Journal of infection and public health , 13(11):1619–1629, 2020
work page 2020
Show all 37 references
-
[7]
Crump and C
K. Crump and C. J. Mode. A general age-dependent branchin g process. ii. Journal of mathematical analysis and applications, 25(1):8–17, 1969
1969
-
[8]
K. S. Crump and C. J. Mode. A general age-dependent branch ing process. i. Journal of Mathematical Analysis and Applications , 24(3):494–508, 1968
1968
-
[9]
M. H. A. Davis. Piecewise-deterministic Markov process es: a general class of nondiffusion stochastic models. J. Roy. Statist. Soc. Ser. B , 46(3):353–388, 1984. With discussion
1984
-
[10]
El Khalifi and T
M. El Khalifi and T. Britton. Extending susceptible-inf ectious-recovered-susceptible epidemics to allow for gra d- ual waning of immunity. Journal of The Royal Society Interface , 20(206):20230042, 2023
2023
-
[13]
P. Gabriel. Measure solutions to the conservative rene wal equation. ESAIM: Proceedings and Surveys , 62:68–78, 2018
2018
-
[14]
D. Gray. Immunological memory. Annual review of immunology , 11(1):49–77, 1993
1993
-
[15]
Hamza, P
K. Hamza, P. Jagers, and F. C. Klebaner. The age structur e of population-dependent general branching processes in environments with a high carrying capacity. Proceedings of the Steklov Institute of Mathematics , 282(1):90–105, 2013
2013
-
[16]
J. M. Heffernan and M. J. Keeling. Implications of vaccin ation and waning immunity. Proceedings of the Royal Society B: Biological Sciences , 276(1664):2071–2080, 2009
2009
-
[17]
H. Inaba. Kermack and McKendrick revisited: The variab le susceptibility model for infectious diseases. Japan journal of industrial and applied mathematics , 18(2):273, 2001
2001
-
[18]
H. Inaba. Variable Susceptibility, Reinfection, and I mmunity. In Age-Structured Population Dynamics in Demog- raphy and Epidemiology , pages 379–442. Springer, 2017
2017
-
[19]
J¨ orgens
K. J¨ orgens. Linear integral operators , volume 7 of Surveys and Reference Works in Mathematics . Pitman (Ad- vanced Publishing Program), Boston, Mass.-London, 1982. T ranslated from the German by G. F. Roach
1982
-
[20]
W. O. Kermack and A. G. McKendrick. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London. Series A. , 115(772):700–721, 1927
1927
-
[21]
W. O. Kermack and A. G. McKendrick. Contributions to the mathematical theory of epidemics. II.—The problem of endemicity. Proceedings of the Royal Society of London. Series A, contai ning papers of a mathematical and physical character, 138(834):55–83, 1932
1932
-
[22]
W. O. Kermack and A. G. McKendrick. Contributions to the mathematical theory of epidemics–III. Further studies of the problem of endemicity. 1933. Proceedings of the Royal Society of London. Series A, contai ning papers of a mathematical and physical character , 141(843):89–118, 1933
1933
-
[23]
Magal and S
P. Magal and S. Ruan. Theory and applications of abstract semilinear Cauchy prob lems, volume 201 of Applied Mathematical Sciences. Springer, Cham, 2018. With a foreword by Glenn Webb. 30 H. GU ´ERIN AND A.B. ZOTSA–NGOUF ACK
2018
-
[24]
M´ el´ eard and V
S. M´ el´ eard and V. C. Tran. Trait substitution sequence process and canonical equation for age-structured pop- ulations. Journal of mathematical biology , 58:881–921, 2009
2009
-
[25]
Oelschlager
K. Oelschlager. Limit theorems for age-structured pop ulations. The Annals of Probability , pages 290–318, 1990
1990
-
[26]
Pakdaman, B
K. Pakdaman, B. Perthame, and D. Salort. Relaxation and self-sustained oscillations in the time elapsed neuron network model. SIAM Journal on Applied Mathematics , 73(3):1260–1279, 2013
2013
-
[27]
Pang and ´E
G. Pang and ´E. Pardoux. Functional central limit theorems for epidemic models with varying infectivity. Stochas- tics, pages 1–48, 2022
2022
-
[28]
Perthame
B. Perthame. Transport equations in biology . Springer Science & Business Media, 2006
2006
-
[29]
Roelly-Coppoletta
S. Roelly-Coppoletta. A criterion of convergence of me asure-valued processes: application to measure branching processes. Stochastics, 17(1-2):43–65, 1986
1986
-
[30]
H. H. Schaefer. Banach lattices and positive operators , volume Band 215 of Die Grundlehren der mathematischen Wissenschaften. Springer-Verlag, New York-Heidelberg, 1974
1974
-
[31]
H. R. Thieme. Semiflows generated by Lipschitz perturba tions of non-densely defined operators. Differential Integral Equations, 3(6):1035–1066, 1990
1990
-
[32]
H. R. Thieme and J. Yang. An endemic model with variable r e-infection rate and applications to influenza. Mathematical Biosciences, 180(1-2):207–235, 2002
2002
-
[33]
Torres, B
N. Torres, B. Perthame, and D. Salort. A multiple time re newal equation for neural assemblies with elapsed time model. Nonlinearity, 35(10):5051, 2022
2022
-
[34]
V. C. Tran. Mod` eles particulaires stochastiques pour des probl` emes d’´ evolution adaptative et pour l’approximation de solutions statistiques . PhD thesis, Universit´ e de Nanterre - Paris X, 2006
2006
-
[35]
G. F. Webb. Theory of nonlinear age-dependent population dynamics , volume 89 of Monographs and Textbooks in Pure and Applied Mathematics . Marcel Dekker, Inc., New York, 1985
1985
-
[36]
A. B. Zotsa Ngoufack. Stochastic Epidemic model with varying infectivity and wan ing immunity: Law of large numbers and Central limit theorem. Weighted norm inequalit y in the variable Lebesgue spaces for Bergman Projector on the unit ball of Cn. PhD thesis, Aix-Marseille Univ...
2024
-
[37]
A. B. Zotsa-Ngoufack. Functional central limit theore ms for epidemic models with varying infectivity and waning immunity. ESAIM: P&S , 29:45 – 112, 2025. D´epartement de math ´ematiques, Universit ´e du Qu ´ebec `a Montr ´eal Email address : guerin.helene@uqam.ca D´epartement...
2025
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.