🦠 Susceptible–Infectious–Recovered–Deceased (SIRD) Model

The Susceptible–Infectious–Recovered–Deceased (SIRD) model is an extension of the classic SIR compartmental framework, with an explicit compartment for disease-induced mortality. This refinement addresses a key limitation of the basic SIR model, in which the single “Removed” (R) category lumped together individuals who recovered and those who died. By separating deaths into their own compartment (D), the SIRD model allows epidemiologists to quantify fatalities more clearly and estimate measures like the case fatality ratio during an epidemic. This proved especially important in outbreaks such as COVID-19, where tracking the deceased separately from recoveries helped in assessing the disease’s severity.

📂 Compartmental Structure

Like SIR, the SIRD model partitions a fixed population N into distinct compartments, but now with four mutually exclusive groups:

• Susceptible (S): Individuals not yet infected but potentially at risk upon contact with an infectious case.
• Infectious (I): Individuals who are currently infected and capable of transmitting the disease.
• Recovered (R): Individuals who have survived the infection and acquired immunity (assumed permanent for the outbreak’s duration).
• Deceased (D): Individuals who have died from the disease.

All individuals must be in one of these four states at any given time, and the compartments sum to the total population N (which may decline if deaths occur and no new births are considered).

⚙️ Mechanistic Dynamics and Equations

The transitions between compartments are governed by a system of ordinary differential equations (ODEs) that extend the standard SIR formulation by adding a mortality term (μ) for the infectious compartment. The flow chart of transitions is:

S → I → R, I → D

In words, susceptible individuals become infected at a rate proportional to contacts with infectious individuals (rate β), infected individuals recover at rate γ, and infected individuals die from the disease at rate μ. The ODE system for the SIRD model can be written as:

• dS/dt = – (β × S × I) / N
• dI/dt = (β × S × I) / N – γ × I – μ × I
• dR/dt = γ × I
• dD/dt = μ × I

These equations conserve total population in the sense that d(S+I+R+D)/dt = 0 if we ignore natural births and non-disease deaths. The term (γ × I + μ × I) is the total removal rate from infectious status, split into recoveries (γ × I) and deaths (μ × I). The basic reproduction number R₀ for the SIRD model is given by:

R₀ = β / (γ + μ)

This reflects that an infected person causes new cases at rate β and remains infectious on average for 1 / (γ + μ) time units. The case fatality ratio (CFR) in this model is:

CFR = μ / (γ + μ)

📊 Parameter Estimates for Influenza

Below are typical parameter ranges for seasonal influenza:

Parameter	Symbol	Range (per day)	Description
Transmission	β	0.2 – 0.6	Rate of infection transmission
Recovery	γ	0.14 – 0.33	Rate of recovery (3–7 days infectious period)
Mortality	μ	0.0001 – 0.001	Disease-induced death rate (very low for influenza)
Basic Reproduction Number	R₀	1.2 – 1.5	Depends on β, γ, and μ
Case Fatality Ratio (CFR)	μ / (γ + μ)	<0.1%	Fraction of infected individuals who die

These values are illustrative and depend on context such as age, location, strain virulence, and intervention strategies.

💡 Significance and Extensions

The SIRD framework provides a clearer picture of disease outcomes by distinguishing between recovery and death. This is crucial for public health planning, as it enables direct computation of the case fatality rate during the epidemic. Early in the COVID-19 pandemic, for example, SIRD-type models were used to jointly estimate the per-day infection, recovery, and death rates from data, thereby informing projections of fatalities. In more complex models, additional compartments can be added to capture other stages or factors. Giordano et al. introduced the eight-compartment SIDARTHE model to analyze COVID-19 in Italy, explicitly splitting the “removed” class into Healed (H) and Extinct (E) compartments (for recovered and deceased, respectively). Such elaborations underscore the value of separating the fate of infected individuals: modeling recoveries and deaths separately yields better insight into the severity and likely impact of an outbreak.

In summary, the SIRD model is a fundamental epidemiological tool that enhances the classic SIR approach by accounting for disease mortality. It balances simplicity and realism, making it a baseline for studying diseases where deaths, though not negligible, occur alongside recoveries. By providing separate tallies for survivors and fatalities, SIRD helps public health officials and researchers evaluate both the spread and the lethality of an infectious disease outbreak.

📚 References

1. Anastassopoulou, C., Russo, L., Tsakris, A., & Siettos, C. (2020). Data-based analysis, modelling and forecasting of the COVID-19 outbreak. PLoS ONE, 15(3): e0230405.
2. Bailey, N. T. J. (1975). The Mathematical Theory of Infectious Diseases and its Applications (2nd ed.). Hafner Press.
3. Giordano, G., Blanchini, F., et al. (2020). Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nature Medicine, 26(6): 855–860.
4. Tang, L., Zhou, Y., Wang, L., et al. (2020). A review of multi-compartment infectious disease models. International Statistical Review, 88(2): 462–513.