Functional Central limit theorems for epidemic models with varying infectivity and waning immunity
Arsene-Brice Zotsa-Ngoufack
TL;DR
This work establishes a functional central limit theorem for a non-Markovian epidemic model with random, time-varying infectivity and gradual loss of immunity, building on a prior FLLN result. The fluctuations around the deterministic limit are governed by a two-dimensional Gaussian process and a coupled system of stochastic nonlinear integral equations, capturing the joint dynamics of the average force of infection and average susceptibility. Central to the analysis are new moment bounds for stochastic integrals with Poisson random measures, a quarantine-based approximation to bound descendant counts, and careful tightness arguments that exploit Regime decompositions and the proposed integral mappings. The results extend non-Markovian and interacting epidemic limit theories and provide a rigorous description of the stochastic fluctuations around endemic/exit dynamics with waning immunity, with potential implications for uncertainty quantification in similar infectious-disease models.
Abstract
We study an individual-based stochastic epidemic model in which infected individuals become susceptible again following each infection (generalized SIS model). Specifically, after each infection, the infectivity is a random function of the time elapsed since the infection, and each recovered individual loses immunity gradually (equivalently, becomes gradually susceptible) after some time according to a random susceptibility function. The epidemic dynamics is described by the average infectivity and susceptibility processes in the population together with the numbers of infected and susceptible/uninfected individuals. In \cite{forien-Zotsa2022stochastic}, a functional law of large numbers (FLLN) is proved as the population size goes to infinity, and asymptotic endemic behaviors are also studied. In this paper, we prove a functional central limit theorem (FCLT) for the stochastic fluctuations of the epidemic dynamics around the FLLN limit. The FCLT limit for the aggregate infectivity and susceptibility processes is given by a system of stochastic non-linear integral equation driven by a two-dimensional Gaussian process.
