Household epidemic models revisited
Frank Ball, Tom Britton, Peter Neal
TL;DR
The paper advances a generalized stochastic household epidemic framework defined by $(X_G,X_L)$, capturing global and local infection contacts, and introduces a swapped-contact variant with probability $p$ for local-to-global transitions. It derives rigorous large-population limits, including a central limit theorem for the final size of a major epidemic and fixed-point characterizations of the limiting outbreak proportion, while giving explicit variance formulas and practical computation guidance. The authors prove monotonicity results: the probability of a major outbreak increases with both household size $h$ and swap probability $p$ for any $(X_G,X_L)$, and the limiting final size $z^{(h,p)}$ increases in $h$ and $p$ when the pgf of $X_L$ is log-convex (as in traditional mixed-Poisson household models); they also present counterexamples when this condition fails. The methodology hinges on enhanced embedding via Sellke constructions, branch- ing-process approximations, and careful variance analysis, complemented by numerical illustrations that validate asymptotic results and reveal nuanced behaviors. The work also discusses extensions to unequal households and multi-type structures, showing the broader relevance of the embedding approach to complex two-level mixing epidemics.
Abstract
We analyse a generalized stochastic household epidemic model defined by a bivariate random variable $(X_G, X_L)$, representing the number of global and local infectious contacts that an infectious individual makes during their infectious period. Each global contact is selected uniformly among all individuals and each local contact is selected uniformly among all other household members. The main focus is when all households have the same size $h \geq 2$, and the number of households is large. Large population properties of the model are derived including a central limit theorem for the final size of a major epidemic, the proof of which utilises an enhanced embedding argument. A modification of the epidemic model is considered where local contacts are replaced by global contacts independently with probability $p$. We then prove monotonicity results for the probability of the major outbreak and the limiting final fraction infected $z$ (conditioned on a major outbreak). a) The probability of a major outbreak is shown to be increasing in both $h$ and $p$ for any distribution of $X_L$. b) The final size $z$ increases monotonically with both $h$ and $p$ if the probability generating function (pgf) of $X_L$ is log-convex, which is satisfied by traditional household epidemic models where $X_L$ has a mixed-Poisson distribution. Additionally, we provide counter examples to b) when the pgf of $X_L$ is not log-convex.
