---
title: Stochastic Kermack-McKendrick Model
url: https://www.emergentmind.com/topics/stochastic-kermack-mckendrick-model
type: topic
---

# Stochastic Kermack-McKendrick Model

Searching arXiv for the cited topic and papers to ground the article.
The stochastic Kermack–McKendrick model denotes, in the cited literature, a family of epidemic constructions that preserve the Kermack–McKendrick concern with susceptible depletion, infection age, and renewal structure, while introducing randomness at different levels of description. These levels include finite-population Markov jump epidemics, stochastic network processes, individual-based models with random infectivity and waning-immunity trajectories, and stochastic microscopic derivations whose deterministic limits are Kermack–McKendrick-type PDE or Volterra systems [1602.01730], [1605.03555], [1511.02545], [2210.04667]. The same literature also makes clear that some works adjacent to this theme are deterministic age-of-infection models supplemented by stochastic simulations or early-phase stochastic surrogates rather than literal stochastic analogues of the full Kermack–McKendrick infection-age theory [2205.15634], [2006.01586], [2104.07793].

## 1. Scope of the concept

At its narrowest, a stochastic Kermack–McKendrick model is a probabilistic epidemic process whose transmission mechanism depends on time since infection and whose deterministic closure or limit recovers a Kermack–McKendrick renewal or infection-age equation. This is the role played by the directed continuous-time random-walk constructions in which an infected individual contributes force of infection according to an infection-age kernel, and by individual-based reinfection models in which each infection episode draws a random infectivity path and a random susceptibility path [1511.02545], [2210.04667].

At a broader level, the term also includes stochastic SIR models that retain the Kermack–McKendrick compartmental backbone but randomize infection and recovery events directly. The classical finite-population “general stochastic epidemic” on the complete graph is of this type: it is Markovian, homogeneous-mixing, and does not carry a general infection-age density, yet it is explicitly compared with the deterministic Kermack–McKendrick SIR system and shown to lie on one side of it in expectation [1602.01730].

The literature also identifies a third layer, in which Kermack–McKendrick equations arise as limits or reductions of richer stochastic mechanisms. Message-passing and pairwise systems derived from stochastic network epidemics converge to deterministic Kermack–McKendrick equations under dense symmetric-network scaling; gathering-based jump Markov epidemics yield generalized-incidence SIR ODEs as mean-field limits; kinetic particle systems reduce to classical SIR under spatial homogenization and Maxwell-type closure [1605.03555], [2203.08260], [2005.09967].

A recurrent terminological issue is that “stochastic” is not used uniformly. Some papers are deterministic in their main analytical object and use stochasticity only in simulation-based robustness checks or as hidden probabilistic interpretation of residence-time assumptions. This is explicit in the age-of-infection identifiability paper, the SARS-CoV-2 modeling review, and the reinfection-age Volterra paper [2205.15634], [2104.07793], [2512.05889].

## 2. The classical finite-population stochastic SIR formulation

A central rigorous benchmark is the finite-population general stochastic SIR epidemic on a complete graph. The population is a finite set $\mathcal V$ of size $N$, with deterministic initial susceptible set $\mathcal S_0$ and deterministic initial infected set $\mathcal I_0=\mathcal V\setminus\mathcal S_0$. While infected, each individual makes infectious contacts to any specified other individual according to an independent Poisson process of rate $\beta$, infection occurs immediately upon infectious contact, and the infectious period is an independent exponential random variable with parameter $\gamma$ [1602.01730].

If $X(t)$ and $Y(t)$ are the numbers susceptible and infected at time $t$, then $(X(t),Y(t))$ is a continuous-time Markov chain with transitions
$$
(x,y)\to(x-1,y+1)\quad\text{at rate }\beta xy,
$$
$$
(x,y)\to(x,y-1)\quad\text{at rate }\gamma y,
$$
and $Z(t)=N-X(t)-Y(t)$ is the recovered count. The deterministic comparator is the classical Markovian Kermack–McKendrick system
$$
\dot S(t)=-\beta S(t)I(t),\qquad
\dot I(t)=\beta S(t)I(t)-\gamma I(t),\qquad
\dot R(t)=\gamma I(t),
$$
with $S(0)=|\mathcal S_0|$, $I(0)=|\mathcal I_0|$, and $R(0)=0$ [1602.01730].

The principal result is one-sided and strict: for all $t>0$,
$$
\mathbb E[X(t)]>S(t),\qquad \mathbb E[Z(t)]<R(t).
$$
Thus the deterministic Kermack–McKendrick model is a strict lower bound on the expected number of susceptibles and a strict upper bound on the expected number of recoveries in the Markovian stochastic epidemic. At the end of the epidemic, if $T=\inf\{t>0:Y(t)=0\}$, then $\mathbb E[X(T)]>S(\infty)$, so the expected stochastic final size is strictly smaller than the deterministic final size [1602.01730].

The proof inserts a message-passing approximation between the stochastic process and the deterministic ODE. On the complete graph,
$$
\mathbb E[X(t)]\ge S_{\mathrm{mes}}(t)>S(t),\qquad t>0.
$$
Here $S_{\mathrm{mes}}(t)$ is built from cavity-state non-transmission probabilities $F^{i\leftarrow j}(t)$, which are exact on trees and lower bounds on general graphs. In this setting the deterministic Kermack–McKendrick trajectory is not merely a heuristic large-$N$ average; it is an exact finite-$N$ lower bound for susceptibility and an exact finite-$N$ upper bound for recoveries under Poisson transmission and exponential recovery [1602.01730].

## 3. Network epidemics, message passing, and dense-network limits

A more general stochastic Kermack–McKendrick framework is obtained on contact networks. In the network SIR model, each individual $i$ has an initial state $Y_i$, an infectious period $\mu_i$, and neighbour-specific contact times $W_{ji}$; infectious periods and post-infection contact times may be correlated within an individual, while different individuals remain independent. For this model, the message-passing variable $F^{i\leftarrow j}(t)$ denotes the probability that $i$, placed in cavity state, receives no infectious contact from $j$ by time $t$, and the susceptible approximation is
$$
S_{\mathrm{mes}}^{(i)}(t)=z_i\prod_{j\in N_i}F^{i\leftarrow j}(t).
$$
The message-passing system has a unique feasible solution, and it provides a rigorous upper bound for expected epidemic size at any fixed time $t>0$ because it underestimates susceptibility and overestimates recovery; equality holds on trees and forests [1605.03555].

In homogeneous $n$-regular symmetric networks the system reduces to four equations. If $F_{\mathrm{sym}}(t)$ is the common message, then
$$
S_{\mathrm{mes}}(t)=z\,F_{\mathrm{sym}}(t)^n.
$$
This reduction supports several structural conclusions. Cycles inhibit spread relative to the Bethe lattice, final size satisfies an explicit message-passing equation in terms of $p=\int_0^\infty f(\tau)\,d\tau$, and the threshold parameter is $R_0=(n-1)p$ in the small-initial-infection regime [1605.03555].

For Poisson contact processes, the same message-passing system is equivalent to a non-Markovian pairwise model, with standard pairwise ODE and delay systems emerging as special cases for exponential and fixed infectious periods. The deterministic Kermack–McKendrick equations arise as the limit of a sequence of homogeneous message-passing systems when $n\to\infty$ and $n f_n(\tau)\to f^*(\tau)$ uniformly on compact time intervals. In this sense, deterministic Kermack–McKendrick dynamics appear as a dense-network asymptotic description of an underlying stochastic epidemic rather than as an isolated deterministic postulate [1605.03555].

The effect of infectious-period variability has also been analyzed in a broad stochastic SEIR network model that contains classic Kermack–McKendrick, pairwise, message passing, and spatial models as special cases. Under convexity conditions on contact-time survival functions, decreasing infectious-period variability while keeping its mean fixed increases the probability that infection reaches any given subset of the population by any given time; at fixed per-edge transmission probability, greater variability in posterior transmissibility likewise tends to dampen spread. For Poisson contact processes and arbitrarily distributed infectious periods, delay-differential and ordinary-differential systems can provide lower and upper bounds, respectively, for the probability that any given individual has been infected by any given time [1710.05730].

## 4. Infection age, memory, and fractional stochastic derivations

A specifically Kermack–McKendrick notion of stochasticity appears when infection age is modeled explicitly and infectivity depends on time since infection. In the continuous-time random-walk construction, an individual infected at time $t'$ infects a particular susceptible in $(t,t+\delta t)$ with probability
$$
\sigma(t,t')\,\delta t+o(\delta t),
$$
with the factorization
$$
\sigma(t,t')=\omega(t)\rho(t-t').
$$
If $\Phi(t,t')=\phi(t,t')\theta(t,t')$ is the probability that the individual remains infectious and alive from $t'$ to $t$, then the incidence flux satisfies the renewal equation
$$
q^+(I,t)=\int_{-\infty}^t \sigma(t,t')\,S(t)\,\Phi(t,t')\,q^+(I,t')\,dt'.
$$
This is already a Kermack–McKendrick Volterra structure. The infected population can be written as
$$
I(t)=I_0(t)+\int_0^t \Phi(t,t')q^+(I,t')\,dt',
$$
and the corresponding infection-age density is
$$
i(a,t)=\Phi(t,t-a)q^+(I,t-a).
$$
The model is therefore equivalent to an infection-age-structured Kermack–McKendrick system with age-specific infectivity kernel $\nu(t,a)=\omega(t)\rho(a)$ [1511.02545].

The same construction identifies a memory kernel
$$
\kappa(t)=\mathcal L^{-1}\{s\mathcal L[\rho](t)\},
$$
so that the force of infection becomes a nonlocal Volterra operator acting on $I(t)$. When
$$
\rho(t)=\frac{t^{\alpha-1}}{\Gamma(\alpha)},\qquad 0<\alpha\le 1,
$$
the Laplace multiplier of the memory kernel is $s^{1-\alpha}$, and the convolution becomes a Riemann–Liouville fractional derivative of order $1-\alpha$. The resulting fractional infectivity SIR system is not an ad hoc fractional replacement of $d/dt$; the fractional operator is induced by an explicit stochastic infectivity-age law [1511.02545].

A complementary stochastic derivation places the heavy tail in recovery rather than transmission. In the directed CTRW model with births and deaths, infectious survival is
$$
\Phi(t,t')=\phi(t-t')\theta(t,t'),
$$
where $\phi$ is the probability of not yet recovering and $\theta$ is survival against death. The memory kernel
$$
K(t)=\mathcal L^{-1}\left\{\frac{\mathcal L\{\psi(t)\}}{\mathcal L\{\phi(t)\}}\right\}
$$
enters the generalized master equation for $I(t)$. If the recovery waiting time is Mittag–Leffler with power-law tail, then the recovery term becomes fractional and the model is consistent with both the age-structured Kermack–McKendrick formulation and the Hethcote–Tudor integral-equation SIR model. Here the non-Markovianity is in recovery, whereas in the infectivity-fractional model it is in transmission [1505.02492].

These two constructions show that stochastic Kermack–McKendrick theory need not be confined to Markovian Poisson infection and exponential recovery. It can also arise from explicit stochastic processes with infection-age-dependent infectivity, heavy-tailed residence times, and memory kernels that remain exactly interpretable in renewal terms [1511.02545], [1505.02492].

## 5. Individual-based stochastic reinfection and waning immunity

A more general individual-based extension introduces repeated infection, random infectivity profiles, and random waning-immunity trajectories. In a population of size $N$, individual $k$ and infection episode $i$ draw an i.i.d. pair
$$
(\lambda_{k,i},\gamma_{k,i}),
$$
where $\lambda_{k,i}(t)$ is infectivity $t$ units after the $i$-th infection and $\gamma_{k,i}(t)$ is susceptibility at that same elapsed time. If $A_k^N(t)$ counts infections of individual $k$ on $(0,t]$ and $\varsigma_k^N(t)$ is the age since the last infection, then the current infection intensity of $k$ is
$$
\Upsilon_k^N(t)=\gamma_{k,A_k^N(t)}(\varsigma_k^N(t))\,\overline{\mathfrak F}^N(t),
$$
where
$$
\overline{\mathfrak F}^N(t)=\frac1N\sum_{k=1}^N \lambda_{k,A_k^N(t)}(\varsigma_k^N(t))
$$
is average infectivity. The counting processes are constructed from independent Poisson random measures via
$$
A_k^N(t)=\int_{[0,t]\times\mathbb R_+}1\{u\le \Upsilon_k^N(r^-)\}\,Q_k(dr,du).
$$
This is a genuinely stochastic homogeneous-mixing epidemic with reinfection and non-Markovian episode-specific trajectories [2210.04667].

The main asymptotic result is a functional law of large numbers for empirical average susceptibility and average infectivity. The deterministic limit $(\overline{\mathfrak S},\overline{\mathfrak F})$ solves a coupled nonlinear Volterra system:
$$
\overline{\mathfrak S}(t)=E\!\left[\gamma_0(t)e^{-\int_0^t\gamma_0(r)\overline{\mathfrak F}(r)\,dr}\right]
+\int_0^t E\!\left[\gamma(t-s)e^{-\int_s^t\gamma(r-s)\overline{\mathfrak F}(r)\,dr}\right]\overline{\mathfrak S}(s)\overline{\mathfrak F}(s)\,ds,
$$
$$
\overline{\mathfrak F}(t)=\overline I(0)\overline\lambda_0(t)+\int_0^t \overline\lambda(t-s)\overline{\mathfrak S}(s)\overline{\mathfrak F}(s)\,ds.
$$
The proof uses i.i.d. auxiliary processes and a propagation-of-chaos argument adapted to age-since-infection dynamics [2210.04667].

A notable asymmetry emerges in the limit: it depends on the full law of the susceptibility random functions but only on the mean infectivity functions. This is because infectivity contributes additively to the population force of infection, while susceptibility enters through exponential no-reinfection factors. The framework includes SIS, SIR, and SIRS as special cases, and under deterministic infectivity and susceptibility functions it recovers a Kermack–McKendrick SIRS PDE of the type studied by Inaba [2210.04667].

The same model yields a nonclassical endemic threshold. With
$$
R_0=\int_0^\infty \overline\lambda(t)\,dt
$$
and $\gamma_*=\lim_{t\to\infty}\gamma(t)$, extinction occurs when
$$
R_0<E\!\left[\frac1{\gamma_*}\right],
$$
while supercritical behavior is governed by the same harmonic-mean quantity rather than by $E[\gamma_*]$. This result makes heterogeneity in long-term post-infection susceptibility structurally important for persistence [2210.04667].

## 6. Adjacent models, deterministic analogues, and recurring misconceptions

Several adjacent models are relevant to the topic but are not, in the strict sense, full stochastic Kermack–McKendrick systems. One example is the birth–death process with immigration for the infected count $I(t)$, with rates $\lambda_n=n\lambda+\nu$ and $\mu_n=n\mu$. Its transient law is exactly negative binomial, with
$$
I(t)\sim NB(r,\beta(t)),\qquad r=\frac{\nu}{\lambda},
$$
and it explains overdispersion, a mode at zero when $r\le 1$, and large discrepancies between deterministic means and typical realizations. However, it is best understood as a stochastic linearized early-phase SIR model with importation, not as a full stochastic replacement of the Kermack–McKendrick infection-age system [2006.01586].

Another example is the kinetic particle formulation in which agents move in physical space, infection corresponds to the reactive collision $S+I\to I+I$, and recovery corresponds to $I\to R$. The microscopic model is stochastic, but the main equations are Boltzmann-type PDEs for phase-space densities $f_S$, $f_I$, and $f_R$. After spatial homogenization and, crucially, Maxwell-molecule closure, the macroscopic system reduces to
$$
\dot S=-\beta SI,\qquad \dot I=\beta SI-\gamma I,\qquad \dot R=\gamma I.
$$
This gives a stochastic microscopic underpinning of classical SIR while remaining distinct from standard stochastic compartment processes [2005.09967].

The gathering-based jump Markov SIR model occupies an intermediate position. Gatherings occur at rate $\mu N$, the gathering size $\Theta$ is random, and each susceptible attendee is infected with probability $1-(1-p)^{\tilde I^N}$ given $\tilde I^N$ infected attendees. The finite-$N$ process is a continuous-time jump Markov chain, and its mean-field limit yields
$$
s'=-\mu s B(i),\qquad i'=\mu s B(i)-\gamma i,\qquad r'=\gamma i,
$$
with
$$
B(i)=E\!\left[\Theta\left(1-(1-pi)^{\Theta-1}\right)\right].
$$
Its basic reproduction number is
$$
R_0=\frac{\mu p\,E[\Theta(\Theta-1)]}{\gamma}.
$$
This is a stochastic microscopic foundation for generalized incidence in a Markovian SIR model, not a non-Markovian infection-age Kermack–McKendrick theory [2203.08260].

The same boundary appears in deterministic age-of-infection inference work. The model
$$
\left\{
\begin{array}{l}
S'(t)=-\tau(t)S(t)\int_0^\infty \beta(a)i(t,a)\,da,\\
\partial_t i+\partial_a i=-\nu i(t,a),\\
i(t,0)=\tau(t)S(t)\int_0^\infty \beta(a)i(t,a)\,da
\end{array}
\right.
$$
is deterministic, and the corresponding Volterra equation for incidence is deterministic as well. The stochastic component is confined to an individual-based simulation model used to test robustness, not to the main analytical epidemic formulation [2205.15634]. Likewise, the deterministic review of Kermack–McKendrick-type models emphasizes that standard ODE formulations hide exponential waiting-time assumptions and that more realistic passage times are Gamma or Beta distributed, but it does not itself formulate a fully stochastic Kermack–McKendrick process [2104.07793]. The recent reinfection-and-infection-age paper is also deterministic: its central object is a nonlinear Volterra equation for incidence
$$
N(t)=\tau(t)S(t)\left[\Lambda(t)+\int_0^{t-t_0}\beta(a)e^{-(\nu+D)a}N(t-a)\,da\right],
$$
with the age-specific reproductive power
$$
\mathcal R(t,a)=\tau(t)S(t)\beta(a)e^{-(\nu+D)a},
$$
which is a renewal kernel rather than a stochastic intensity model [2512.05889].

A persistent misconception is therefore to treat every model with random simulation output, every stochastic surrogate for $I(t)$, or every deterministic renewal system with hidden probabilistic ingredients as a stochastic Kermack–McKendrick model in the same sense. The literature instead supports a sharper distinction. Some models are exact finite-population stochastic epidemics compared against Kermack–McKendrick bounds [1602.01730]; some are stochastic network or individual-based systems whose law-of-large-numbers limits are Kermack–McKendrick-type equations [1605.03555], [2210.04667]; some are stochastic derivations of non-Markovian infection-age or fractional models [1511.02545], [1505.02492]; and some are deterministic models or early-phase stochastic approximations that are relevant to Kermack–McKendrick theory without constituting full stochastic analogues [2205.15634], [2006.01586], [2104.07793], [2512.05889].

Source: https://www.emergentmind.com/topics/stochastic-kermack-mckendrick-model