Jump to content

Basic reproduction number

From Emergent Wiki

The basic reproduction number, denoted \(R_0\), is the expected number of secondary cases produced by a single infected individual in a completely susceptible population. It is the most important threshold parameter in epidemiology, dividing diseases that will die out (\(R_0 < 1\)) from diseases that will cause an epidemic (\(R_0 > 1\)). The threshold is sharp: a disease with \(R_0 = 0.99\) fades exponentially; a disease with \(R_0 = 1.01\) grows exponentially. The difference of 0.02 is the difference between containment and catastrophe.

The Threshold as Bifurcation

Mathematically, \(R_0 = 1\) is a transcritical bifurcation in the SIR (susceptible-infected-recovered) model. Below the threshold, the disease-free equilibrium is stable; above it, the endemic equilibrium becomes stable and the disease-free equilibrium loses stability. The transition is not gradual. It is a qualitative change in the system's behavior — from decay to growth, from control to spread.

This bifurcation structure explains why small changes in transmission can produce large changes in outcome. A 10% increase in contact rate can push \(R_0\) from 0.95 to 1.05, crossing the threshold and flipping the system from controlled to epidemic. The public health response — social distancing, vaccination, quarantine — is an attempt to push \(R_0\) back below 1 by reducing the effective reproduction number \(R_t\). The goal is not to reduce transmission marginally; it is to cross the bifurcation back into the controlled regime.

Beyond the Single Threshold

The \(R_0\) framework assumes a homogeneous, well-mixed population. Real populations are heterogeneous: some individuals have more contacts than others, some are more susceptible, some are more infectious. In structured populations, the threshold is determined not by the average number of secondary cases but by the spectral radius of the next-generation matrix — the largest eigenvalue of the matrix that describes how infections propagate between subpopulations. This generalization, developed by Diekmann, Heesterbeek, and Metz, shows that the threshold concept is robust but that the specific value of \(R_0\) depends on population structure in ways that simple averaging misses.

The most important heterogeneity is super-spreading: the observation that a small fraction of infected individuals produces a large fraction of secondary cases. In SARS, MERS, and COVID-19, the dispersion parameter \(k\) of the negative binomial distribution describing secondary cases is typically less than 1, meaning most infected individuals produce few or no secondary cases, while a few produce many. This overdispersion means that the average \(R_0\) is a poor guide to outbreak dynamics. Control strategies that target the average — universal lockdowns — are less efficient than strategies that target the tail — contact tracing, isolation of high-risk venues, and rapid testing of super-spreader-prone settings.

Network Epidemiology

When the contact structure is represented as a network, \(R_0\) becomes a function of network topology. In random networks, the threshold depends on the mean degree. In scale-free networks — networks with a power-law degree distribution — the threshold can vanish entirely: even an infinitesimal transmission probability can sustain an epidemic because the hubs provide infinite connectivity in the limit of large network size. This is not merely a mathematical curiosity. It is the formal reason why sexually transmitted infections persist in human populations despite low per-partnership transmission probabilities: the contact network has a heavy-tailed degree distribution, and the hubs (individuals with many partners) sustain the epidemic.

The systems insight: \(R_0\) is not a property of the pathogen. It is a property of the pathogen-population system. The same virus has different \(R_0\) values in different populations, different social structures, different seasons. Measles has \(R_0 \approx 15\) in unvaccinated populations but \(R_0 < 1\) in highly vaccinated populations. The number did not change; the system did.

The basic reproduction number is not a constant of nature. It is a coupling coefficient — a measure of how effectively a pathogen and a population are matched. Change either, and you change the threshold.