A Robust Numerical Framework for Multi-Strain Epidemics
Abstract
Epidemiological models often take the form of systems of differential equations of varying levels of complexity. In multi-strain models, complexity arises from the interaction of diverse pathogen types that may exhibit different patterns of cross-immunity due to the host's immune system. Where previous work has tended to focus on theoretical analysis of this complexity, or model simplification, less effort has been given to efficient numerical exploration of full multi-strain disease dynamics. Here, we propose a computational framework for representing multi-strain disease transmission models, through a mathematically transparent hard-coded mapping between sets and integers. By mapping infection histories into Boolean matrices, the framework eliminates explicit iteration over strain combinations while preserving key biological processes, including sequential infection and partial cross-immunity. Benchmark comparisons of numerical solutions of our matrix formulation and a set-indexed implementation demonstrate substantial performance gains. Across models with increasing numbers of strains, the matrix formulation consistently reduces runtime, with a speedup exceeding two orders of magnitude for systems with 11 strains. We present an application of this approach to the case of COVID-19 in Indonesia, focusing on the replacement of the Delta variant by Omicron. We hope that our approach represents a robust computational approach to automatically generate diverse multi-strain epidemic models that can be solved in a computationally efficient manner for a variety of scientific and public-health questions.
Related
- RoDPAA: A Robust Dissimilarity-Based Principal Angle Analysis Framework for Multi-View Microbiome Data Integration
- Gut microbiota and multi-omics analysis identify candidate gut-muscle associations underlying breast muscle quality differences between Mimas pigeon and pigeon Tianxiang 1 strain
- SeverYX: An ensemble framework for robust variant severity classification
- Hybrid local–global vision transformer framework for robust plant disease detection in precision agriculture
- How far can microbial monocultures predict growth in multi-strain communities?
BibTeX
@misc{ilahi2026robust,
title = {{A Robust Numerical Framework for Multi-Strain Epidemics}},
author = {Fadilah Ilahi and Ian Hall and Thomas House},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.15796},
eprint = {2610.15796},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.15796},
note = {CC0},
}
RIS
TY - UNPB
TI - A Robust Numerical Framework for Multi-Strain Epidemics
AU - Fadilah Ilahi
AU - Ian Hall
AU - Thomas House
PY - 2026
DA - 2026/10/11
PB - arXiv
DO - 10.48550/arxiv.2610.15796
UR - https://arxiv.org/abs/2610.15796
AB - Epidemiological models often take the form of systems of differential equations of varying levels of complexity. In multi-strain models, complexity arises from the interaction of diverse pathogen types that may exhibit different patterns of cross-immunity due to the host's immune system. Where previous work has tended to focus on theoretical analysis of this complexity, or model simplification, less effort has been given to efficient numerical exploration of full multi-strain disease dynamics. Here, we propose a computational framework for representing multi-strain disease transmission models, through a mathematically transparent hard-coded mapping between sets and integers. By mapping infection histories into Boolean matrices, the framework eliminates explicit iteration over strain combinations while preserving key biological processes, including sequential infection and partial cross-immunity. Benchmark comparisons of numerical solutions of our matrix formulation and a set-indexed implementation demonstrate substantial performance gains. Across models with increasing numbers of strains, the matrix formulation consistently reduces runtime, with a speedup exceeding two orders of magnitude for systems with 11 strains. We present an application of this approach to the case of COVID-19 in Indonesia, focusing on the replacement of the Delta variant by Omicron. We hope that our approach represents a robust computational approach to automatically generate diverse multi-strain epidemic models that can be solved in a computationally efficient manner for a variety of scientific and public-health questions.
ER -