Multi-group (or multi-patch) compartmental models are indispensable for accurately simulating infectious disease dynamics when the population structure is highly heterogeneous. By segmenting the total population into distinct interacting subgroups—such as age classes, regions, or behavioral cohorts—these models move beyond the homogeneous mixing assumption of classical SIR models to capture differential risks of infection and transmission across complex social structures. This framework is essential for designing targeted interventions reflecting group-specific contact patterns or susceptibility.
📐 Compartmental Structure and Flow Explanation
We consider a general framework for two interacting groups (Group 1 and Group 2), where movement between epidemiological states (S, I, T, R) occurs within each group, while transmission involves mixing between groups.
| Compartment | Notation | Description and Transition |
|---|---|---|
| Susceptible | Sᵢ | Uninfected individuals in Group i who may be infected by any interacting group. |
| Infectious | Iᵢ | Individuals in Group i capable of transmitting the disease. |
| Treated/Isolated | Tᵢ | Individuals in Group i under treatment or isolation, usually with reduced infectivity (σ). |
| Removed | Rᵢ | Individuals recovered or otherwise removed from the infectious process. |
The key mechanism is the Force of Infection (Λᵢ) applied to the susceptible population in each group, governed by weighted contacts and infectious pressure across all groups. This structure accommodates realistic assumptions such as proportionate mixing, where group-specific contacts depend on activity levels of all interacting groups.
🧮 Mathematical Formulation (Simplified Two-Group Model)
For two groups with constant population sizes Nᵢ and mass-action incidence modified by mixing and isolation effects, the Susceptible–Infectious–Treated–Removed (SITR) dynamics are:
dSᵢ/dt = − Sᵢ [ (aᵢ p₁ (I₁ + σ T₁) / N₁) + (aᵢ p₂ (I₂ + σ T₂) / N₂ ) ]
dIᵢ/dt = Sᵢ [ (aᵢ p₁ (I₁ + σ T₁) / N₁) + (aᵢ p₂ (I₂ + σ T₂) / N₂ ) ] − (α + φᵢ) Iᵢ
dTᵢ/dt = φᵢ Iᵢ − α_T Tᵢ
dRᵢ/dt = α Iᵢ + α_T Tᵢ
where:
• pⱼ is the proportion of contacts directed toward Group j.
• aᵢ is the average activity level (contact intensity) of Group i.
• σ (0 < σ < 1) is the reduced infectivity of treated individuals.
• φᵢ is the treatment/isolation rate in Group i.
• α and α_T are recovery/removal rates for Iᵢ and Tᵢ.
⚗ Parameter Definitions and Typical Ranges
| Parameter | Definition | Typical Range | Interpretation |
|---|---|---|---|
| α, α_T | Recovery/removal rates | 0.14 – 0.33 per day | Correspond to mean infectious periods of 3–7 days. |
| aᵢ | Group-specific contact/activity level | 0.2 – 1.0 | Higher for highly active groups (e.g., school-age children). |
| φᵢ | Isolation/treatment rate in Group i | 0.1 – 2.0 per day | High φᵢ means rapid case isolation. |
| σ | Reduction factor for infectivity under treatment | 0.1 – 0.7 | Reflects isolation effectiveness. |
| pⱼ | Proportionate mixing weight toward Group j | 0.0 – 1.0 | Sum over all groups equals 1. |
| R_c | Control reproduction number | 0.5 – 4.0 | Determined by activity levels, infectivity, and treatment rates. |
🎯 Applicability and Limitations
Applications
- Targeted Intervention Strategy
Multi-group models are essential for determining whether interventions—such as vaccination, treatment (φᵢ), or antivirals—should be concentrated on highly active groups (e.g., schoolchildren) or vulnerable populations (e.g., elderly). These models often reveal that focusing on the most active group yields near-optimal control benefits. - Heterogeneous Risk Quantification
They allow calculation of attack rates and final epidemic sizes by group, demonstrating how disease burden distributes unevenly across populations and identifying which groups act as primary transmission drivers. - Insight into Mixing Patterns
With structured contact matrices, these models quantify the effects of interventions like school closures or workplace distancing, revealing differential impacts across age or behavioral groups.
Key Assumptions and Weaknesses
- Mixing Homogeneity Assumption
Proportionate mixing assumes that contacts depend solely on group activity levels, ignoring preferential mixing (e.g., assortative mixing in sexual networks). - High Data Requirements
Requires empirical estimates of group-specific activity levels (aᵢ) and mixing probabilities (pᵢⱼ), often sourced from detailed contact surveys that may be incomplete or noisy. - Complexity–Realism Trade-off
Including additional epidemiological states or more groups increases realism but dramatically expands dimensionality, making model calibration and validation more difficult.
📚 Selected References
- Anderson, R. M., & May, R. M. Infectious Diseases of Humans: Dynamics and Control. Oxford University Press.
- Ferguson, N. M., et al. Strategies for mitigating an influenza pandemic. Nature (2006).
- Hethcote, H. W. The mathematics of infectious diseases. SIAM Review (2000).
- Mossong, J., et al. Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases. PLoS Medicine (2008).
- Wang, W., & Ruan, S. Bifurcations in an epidemic model with constant removal rate of the infectives. Journal of Mathematical Analysis and Applications (2004).