Pith. sign in

REVIEW 3 major objections 5 minor 14 references

The impact of recovery rate heterogeneity in achieving herd immunity

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In heterogeneous SIR populations, herd immunity holds exactly when the average recovery time is finite; the mean recovery rate is not the deciding quantity.

desk verdict The criterion E[1/γ]<∞ is genuine and mostly proved, but the uniform transmission bound βmin>0 is load-bearing, not technical, and the abstract oversells the scope. read the letter →

arxiv 2411.13130 v2 pith:PDNVXFRU submitted 2024-11-20 q-bio.PE math.DS

Gabriel Turinici This is my paper · ORCID
classification q-bio.PEmath.DS MSC 92D30
keywords herdimmunityheterogeneousrecoveryratemeantimeSIRmodelSEIRfinalepidemicsizeheterogeneity
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper asks which single number tells an epidemiologist whether an epidemic in a population with variable recovery speeds will stop before infecting everyone. The answer it defends is the mean recovery time, the average of $1/\gamma(x)$ across people, not the mean recovery rate $\mathbb{E}[\gamma]$. The central result is that, under the model's assumptions, the heterogeneous SIR system has the herd immunity property $S(\infty)>0$ if and only if $\mathbb{E}[1/\gamma(x)]<\infty$, and the same criterion holds for the SEIR extension. The paper also exhibits populations with identical average recovery rate that reach opposite fates, so the mean rate alone is not a reliable metric. If the result is right, forecasts of final epidemic size and vaccination needs should be built from the distribution of recovery times, especially their tail, rather than from the average recovery rate.

What carries the argument

The load-bearing object is the trait-structured SIR system (1)-(3), in which each individual carries a trait $x$ with recovery rate $\gamma(x)$, and the key identity obtained from rewriting the infection term as a recovery-driven integral: $\frac{d}{dt}\ln S(t,x)=-\int_\Omega \frac{\beta(x,y)}{\gamma(y)}\frac{d}{dt}R(t,y)\,dy$. This identity converts the infection history into cumulative recoveries, giving the bound $S(\infty,x)\ge S(0,x)\exp(-\int_\Omega \beta_{\max}/\gamma(y)\,dy)$. The named object deciding the outcome is the mean recovery time $\mathbb{E}[1/\gamma(x)]$: if it is finite the exponential bound is positive, while if it diverges the same structure forces $S(\infty)=0$; the proof handles the case $\gamma=0$ by showing infecteds with that trait grow exponentially and drag all susceptibles to zero.

What would settle it

Run the trait-structured SIR system with $\Omega=[0,1]$, uniform $\beta=1/4$, and $\gamma(x)=x/3$, so $\mathbb{E}[1/\gamma]=\infty$; the theorem predicts $S(\infty)=0$, matching the paper's numerics. To disprove the theorem one would need, under assumptions (8)-(11), any instance with $\mathbb{E}[1/\gamma]=\infty$ and $S(\infty)>0$; a simulation or analytic example producing that would settle the claim negatively.

Watch

Extended reading notes

Core claim

Proposition 3 states the main discovery: assume the transmission function is uniformly positive and bounded, there are no initially recovered individuals, and some initial infectives exist; then the trait-structured SIR equations (1)-(3) have the herd immunity property, meaning a positive fraction of the initially susceptible population remains uninfected, if and only if the population mean of the inverse recovery rate $1/\gamma(x)$ is finite. The direction 'finite mean recovery time implies herd immunity' follows from an exponential lower bound on the susceptible fraction that depends on $\int \beta_{\max}/\gamma(y)\,dy$. The reverse direction shows that if $\mathbb{E}[1/\gamma]=\infty$, the susceptible fraction decays to zero, with the special case $\gamma=0$ on a non-negligible subgroup forcing total infection. Proposition 5 extends the same equivalence to the SEIR model when the exposure rate is bounded away from zero and infinity. The paper's examples with $\gamma=1/6$ versus a two-point or density distribution with the same mean illustrate that equal average recovery rates can produce herd immunity in one population and complete infection in another.

Load-bearing premise

The argument that finite mean recovery time is also necessary assumes every trait group can both infect and be infected by every other group, through the lower bound $\beta_{\min}>0$ on transmission rates; if some subgroup is epidemiologically isolated, it can remain susceptible no matter how long other groups stay infectious.

Editorial extensions

If this is right

  • In any heterogeneous SIR population meeting the assumptions, a positive fraction of susceptibles escapes infection exactly when the population-averaged infectious period $\mathbb{E}[1/\gamma]$ is finite.
  • The mean recovery rate $\mathbb{E}[\gamma]$ is not a reliable predictor: the paper exhibits three settings with the same mean rate $1/6$ whose final susceptible fractions are approximately $0.42$, $0$, and $0$.
  • Any non-negligible subgroup that never recovers ($\gamma=0$) eliminates herd immunity for the whole population, no matter how quickly everyone else recovers.
  • Herd immunity can fail even when every individual has finite recovery time, because distributions like $\gamma(x)=x/3$ on $[0,1]$ give $\mathbb{E}[1/\gamma]=\infty$ and empty the susceptible class.
  • The identical finite-mean-recovery-time criterion governs the SEIR model with an exposed class, provided exposure rates are bounded away from $0$ and $\infty$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A practical consequence the paper leaves implicit: surveillance should track the tail of long infectious periods, since an infinite-mean tail alone destroys herd immunity even when most patients recover quickly.
  • In contact networks with disconnected components the criterion likely localizes: each epidemiologically closed component needs finite mean recovery time for its own susceptibles to be protected.
  • The theorem yields a testable ordering: among populations matched for mean $\gamma$, final susceptible fraction should increase as $\mathbb{E}[1/\gamma]$ decreases; stratified recovery data could check this.
  • Vaccination strategy could exploit the result by prioritizing subgroups with the longest infectious periods, since shortening their infectious period directly restores herd immunity even if it barely moves the population mean $\gamma$.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper studies a heterogeneous SIR model (and then an SEIR extension) in which the recovery rate is a trait-dependent function γ(x). Herd immunity is defined as S(∞)>0, the aggregate susceptible fraction remaining positive at infinite time. The paper first gives numerical examples showing that the mean recovery rate E[γ] does not predict herd immunity, then states its central result (Proposition 3): under assumptions (8)–(11), the system has the herd immunity property if and only if the mean recovery time E[1/γ] is finite. Proposition 5 extends this characterization to the SEIR model. The paper concludes that the mean recovery time, rather than the mean recovery rate, is the relevant metric for herd immunity in heterogeneous populations.

Significance. The proposed criterion is elegant and, if established, practically relevant: it gives a sharp, parameter-free condition for whether a positive susceptible fraction survives an epidemic in a heterogeneous population. The forward direction is clean, the proof is self-contained, the examples are informative, and the code is made available. The main caveat is that the equivalence relies on the full-contact assumption 0<βmin≤β(x,y), and the paper's abstract and discussion state the result without this qualifier; as written, the theorem is correct under its hypotheses, but the scope claimed outside the assumptions is broader than what is proved.

major comments (3)
  1. [§2 (Assumption (9)) and Proposition 3] The assumption 0<βmin≤β(x,y) is load-bearing for the 'only if' direction and is not merely technical. Without it the equivalence is false: take Ω={1,2}, β(1,1)=β(2,2)=β>0, β(1,2)=β(2,1)=0, γ(1)=0, γ(2)=1, S0(1)=S0(2)=1/2, I0(1)>0, I0(2)=0. Group 2 is never exposed, so S(∞)≥1/2>0 (herd immunity holds), yet E[1/γ]=∞ because P(γ=0)=1/2. The paper's statement in §2 that (9) mainly prevents 'totally immune individuals' mischaracterizes the issue; zero cross-group transmission is a structural feature of many realistic contact matrices. Please restrict the abstract and discussion to models satisfying the full-contact condition, or state clearly which parts of the theorem survive when β has zero entries.
  2. [§3, proof of Proposition 3, Eq. (25)] The chain after (23)–(24) is not written correctly: from R(∞,y)≥1−C3 one obtains S(∞,x)≤S(0,x)exp(−∫ βmin γ(y)^{-1}(1−C3)dy), not the bound with C3 in the exponent. The displayed inequality with C3 is not justified unless C3≤1/2, which is not established. The argument is easily repaired by replacing C3 with any positive constant bounded away from zero, but the displayed equation should be fixed.
  3. [§3, proof of Proposition 3, Part II] The statement 'it can be easily checked that lim_{t→∞}I(t,x)=0' is used to pass from S(∞,x)≤C3 to R(∞,x)≥1−C3. Please include the short argument: S(t,x) and R(t,x) are monotone, so I(t,x) has a pointwise limit; since ∫0∞ γ(x)I(t,x)dt=R(∞,x)≤1 and γ(x)>0 a.e., the limit must be zero. Similarly, in the SEIR proof the conclusion E_Z(∞)=0, I_Z(∞)=1 from (35) needs a few lines of justification rather than being asserted.
minor comments (5)
  1. [§2.1, Lemma 2 proof, Eq. (14)] The formula for I(t) appears to contain a typo: I(t)=(I(t,1)+S(t,2))/2 should presumably be I(t)=(I(t,1)+I(t,2))/2.
  2. [§2.1, Lemma 2 proof] The sentence 'since I(0)>0 we have I(0,1)>0' is not automatic; if I(0,1)=0, positivity follows from I′(0,1)=βS(0,1)I(0)>0.
  3. [§3, after Eq. (23)] Notation such as I(0) and I(t) in the proof of Proposition 3 refers to the aggregate variables defined in (4), but this is not restated at the point of use; please clarify to avoid confusion with the trait-dependent I(t,x).
  4. [§2.1, Figure 1 and surrounding text] In the middle panel, the caption says S(∞)=0.0, but for the γ=0 group the susceptible and infected fractions do not vanish individually; the text should specify that this is the aggregate S(∞).
  5. [Data availability] The GitHub link contains spaces: 'epidemiology heterogeneous gamma' should be a proper URL or percent-encoded, otherwise it may not be clickable.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the herd-immunity criterion is derived from the model equations via direct inequalities.

full rationale

The paper's central result, Proposition 3, is an if-and-only-if statement about the heterogeneous SIR system (1)-(3). Both directions are proved from the equations themselves: the forward direction uses the identity d/dt ln(S(t,x)) = -∫ β(x,y) I(t,y) dy together with the upper bound β ≤ βmax to obtain S(∞) ≥ S(0) exp(-∫ βmax/γ(y) dy) > 0 when E[1/γ] < ∞; the reverse direction uses the lower bound βmin > 0 to show that a non-negligible set with γ = 0 forces the infected fraction in that set to converge to 1, hence S(∞) = 0, and then uses the same logarithmic identity to conclude that S(∞) > 0 implies E[1/γ] < ∞. No parameter is fitted to data and then renamed a prediction; the numerical examples in Section 2.1 are illustrations rather than statistical inference. The assumption βmin > 0 in (9) restricts the theorem's scope and is flagged by the paper as technical, but this is a matter of hypothesis validity, not circularity: the equivalence is derived from that hypothesis, not assumed as the conclusion. The paper does not rely on self-citations for load-bearing content; the cited prior works [1,14] are used only to motivate the model formulation. The concern that real contact matrices may violate (9) is a correctness or scoping issue, not a circularity issue. Consequently, the derivation is self-contained given its assumptions, and no circular step can be exhibited.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

No parameters are fitted to data; the examples use hand-chosen distributions for illustration. The axioms are the model assumptions stated in the paper.

assumptions (6)
  • domain assumption Standard SIR/SEIR compartment dynamics with trait-dependent parameters (equations 1-3 and 26-29).
    The model is inherited from cited prior work [14,1] and is the object of study, not derived.
  • domain assumption Nonnegative recovery rate γ(x) ≥ 0 (hypothesis 8).
    Needed to define 1/γ with the convention 1/0=∞.
  • domain assumption Uniform positive bounds on transmission: 0<βmin≤β(x,y)≤βmax<∞ (hypothesis 9).
    Load-bearing for both directions of Proposition 3; without βmin>0 the iff fails.
  • domain assumption R(0)=0 and I(0)>0 (hypotheses 10-11).
    Conventions that make the epidemic start meaningful and simplify the lower bound.
  • domain assumption No coupling between recovery rate and epidemic state.
    Stated in Section 2; allows γ to be a fixed function of trait x.
  • domain assumption For SEIR, bounded transition α between exposed and infectious (hypothesis 30).
    Used in the SEIR extension, Proposition 5.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The impact of recovery rate heterogeneity in achieving herd immunity." pith.science (2026). https://pith.science/paper/PDNVXFRU

@misc{pith2026241113130,
  author       = {Pith},
  title        = {Pith review of: The impact of recovery rate heterogeneity in achieving herd immunity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PDNVXFRU}},
  note         = {Machine review of arXiv:2411.13130}
}
abstract

Herd immunity is a critical concept in epidemiology, describing a threshold at which a sufficient proportion of a population is immune, either through infection or vaccination, thereby preventing sustained transmission of a pathogen. In the classic Susceptible-Infectious-Recovered (SIR) model, which has been widely used to study infectious disease dynamics, the achievement of herd immunity depends on key parameters, including the transmission rate ($\beta$) and the recovery rate ($\gamma$), where $\gamma$ represents the inverse of the mean infectious period. While the transmission rate has received substantial attention, recent studies have underscored the significant role of $\gamma$ in determining the timing and sustainability of herd immunity. Additionally, it is becoming increasingly evident that assuming $\gamma$ as a constant parameter might oversimplify the dynamics, as variations in recovery times can reflect diverse biological, social, and healthcare-related factors. In this paper, we investigate how heterogeneity in the recovery rate affects herd immunity. We show empirically that the mean of the recovery rate is not a reliable metric for determining the achievement of herd immunity. Furthermore, we provide a theoretical result demonstrating that it is instead the mean recovery time, which is the mean of the inverse $1/\gamma$ of the recovery rate that is critical in deciding whether herd immunity is achievable within the SIR framework. A similar result is proved for the SEIR model. These insights have significant implications for public health interventions and theoretical modeling of epidemic dynamics.

Figures

Figures reproduced from arXiv: 2411.13130 by the authors.

Figure 1
Figure 1. Left: Solution of the standard SIR model for β = 1/4, γ = 1/6. Here S(∞) ≃ 0.42. Middle: solution of a two group SIR model for β = 1/4, γ = 0 or 1/3 (each with probability 1/2). Here S(∞) = 0.0. Right: solution of a SIR model with β = 1/4, Ω = [0, 1], γ(x) = x/3. Again S(∞) = 0.0. Consider first the simple example where β = 1/4, γ = 1/6 (both are con￾stants) that corresponds to a R0 = 3/2. The numerical result is gi… view at source ↗
Figure 2
Figure 2. 3.1 Herd immunity : extension to the SEIR model Similar results can be proved for the SEIR model. Recall that in this case the “Infected and infectious” class is split in two parts, the “Exposed” class that 8 [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 2
Figure 2. Solution of a SIR model with β = 1/4, Ω = [0, 1], γ(x) = 3 2 √ x. Herd immunity (S(∞) > 0) is observed numerically, coherent with the theoretical results. is infected but not yet infectious and the “Infectious” class that is infected and infectious. The equations analogue to (1)-(3) read now: ∂tS(t, x) = −S(t, x) Z Ω β(x, y)I(t, y)dy, S(0, x) = S0(x) (26) ∂tE(t, x) = S(t, x) Z Ω β(x, y)I(t, y)dy − α(x)E(t, x), E(0, … view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 14 canonical work pages

  1. [1]

    Final size and convergence rate for an epidemic in heterogeneous populations

    Luis Almeida, Pierre-Alexandre Bliman, Gr´ egoire Nadin, Beno ˆ ıt Perthame, and Nicolas Vauchelet. Final size and convergence rate for an epidemic in heterogeneous populations. Mathematical Models and Methods in Applied Sciences, 31(05):1021–1051, 2021

    work page 2021

  2. [2]

    Anderson and Robert M

    Roy M. Anderson and Robert M. May. Infectious Diseases of Humans: Dynamics and Control . Oxford University Press, 1991

    work page 1991

  3. [3]

    Weitz, and Jean-Marc Oury

    Henri Berestycki, Beno ˆ ıt Desjardins, Joshua S. Weitz, and Jean-Marc Oury. Epidemic modeling with heterogeneity and social diffusion. Journal of Mathematical Biology, 86(4):60, March 2023

    work page 2023

  4. [4]

    A mathematical model re- veals the influence of population heterogeneity on herd immunity to SARS- CoV-2

    Tom Britton, Frank Ball, and Pieter Trapman. A mathematical model re- veals the influence of population heterogeneity on herd immunity to SARS- CoV-2. Science, 369(6505):846–849, 2020

    work page 2020

  5. [5]

    Rodrigues, and Yamir Moreno

    Guilherme Ferraz de Arruda, Giovanni Petri, Francisco A. Rodrigues, and Yamir Moreno. Impact of the distribution of recovery rates on disease spreading in complex networks. Phys. Rev. Res. , 2:013046, Jan 2020

    work page 2020

  6. [6]

    How heterogeneous susceptibility and recovery rates affect the spread of epidemics on networks

    Wei Gou and Zhen Jin. How heterogeneous susceptibility and recovery rates affect the spread of epidemics on networks. Infectious Disease Modelling , 2(3):353–367, 2017

    work page 2017

  7. [7]

    SIR Model with Stochastic Transmission

    Christian Gourieroux and Yang Lu. SIR model with stochastic transmis- sion, 2020. arXiv:2011.07816

    work page Pith review arXiv 2020

  8. [8]

    Heterogene- ity matters: Contact structure and individual variation shape epidemic dynamics

    Gerrit Grossmann, Michael Backenkohler, and Verena Wolf. Heterogene- ity matters: Contact structure and individual variation shape epidemic dynamics. PLOS ONE , 16(7):1–19, 07 2021

    work page 2021

Show all 14 references
  1. [9]

    Keeling and Pejman Rohani

    Matt J. Keeling and Pejman Rohani. Modeling Infectious Diseases in Hu- mans and Animals . Princeton University Press, 2008

    2008

  2. [10]

    Kermack and Anderson G

    William O. Kermack and Anderson 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

  3. [11]

    The dynamics of epidemiological systems with nonautonomous and random coefficients

    Peter E Kloeden and V Kozyakin. The dynamics of epidemiological systems with nonautonomous and random coefficients. MESA, 2(2):159–172, 2011

    2011

  4. [12]

    A. L. Lloyd. Destabilization of epidemic models with the inclusion of real- istic distributions of infectious periods. Proceedings of the Royal Society of London. Series B: Biological Sciences , 268(1470):985–993, 2001

    2001

  5. [13]

    Alun L. Lloyd. Realistic Distributions of Infectious Periods in Epidemic Models: Changing Patterns of Persistence and Dynamics. Theoretical Pop- ulation Biology, 60(1):59–71, 2001. 11

    2001

  6. [14]

    Novozhilov

    Artem S. Novozhilov. On the spread of epidemics in a closed heterogeneous population. Mathematical Biosciences, 215(2):177–185, 2008. 12

    2008

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.