COVID-19 and Compartmental Models in Epidemiology
The severity of the current SARS-CoV-2 epidemic is undeniable: since the latest months of 2019, the COVID-19 outbreak is having a significant impact in the world at the macro level, starting its spread from China, then to the Asia-Pacific and then around the rest of the globe.
Photo by Edward Howell on Unsplash
More than 41 articles were published about COVID-19 just from the beginning of 2020, more than 10,000 GitHub forks of CSSEGISandData were made, various open data sets (including on GitHub and on Kaggle) shared the number of cases and the number of people affected as well as other metrics.
Many scientists around the globe participated in a global hackathon, trying to address current issues (lockdowns, lack of resources, etc.), trying to fit, predict and estimate the impact. Which models from these ones are close to the reality and which ones have a chance to be implemented as solutions for humanitarian response?
Basic Mathematical Models
There are various types of models in epidemiology, which are used for prediction of prevalence (total number of people infected), the duration of an epidemic spreading processes and other variables.
One of the most commonly used types are compartmental models, developed in the 1920s from the Kermack–McKendrick theory. In compartmental models, the population is divided into compartments of people with the same properties: people who are susceptible to the virus; those who have recovered; etc. There are three main types of compartmental models in epidemiology based on which many other models can be built upon:
- Susceptible-Infectious (SI), where each person can be in one of two states, either susceptible (S) or infectious (I);
- Susceptible-Infectious-Susceptible (SIS), where a person can be in one of two states, susceptible (S) or infectious (I), but can then become susceptible again;
- Susceptible-Infectious-Recovered (SIR), where a person can be in one of three states: susceptible (S), infectious (I) or recovered (R).
In the SIR model, one can study the number of people in each of three compartments: susceptible, infectious and recovered, denoted by the variables S, I and R correspondingly. The SIR system can be expressed by the set of ordinary differential equations proposed by O. Kermack and Anderson Gray McKendrick.
Another model, which has been proven to be much more realistic for some epidemics spreading models is the SEIR model. The main idea of this model is that for many important infections, there is a significant incubation period during which individuals have been infected but are not yet infectious themselves. During this period, the individual is in a special state E (exposed), an additional state to three states in the SIR model.
Some models described here are also explained on this page.
"Fitting the Curve" - Issues to Consider
During COVID-19 many open data sources were gathered, which explains the temptation of fitting the curve to empirical data or inferring parameters of spreading models. Here, we list some of the main issues with simply trying to do "curve fitting" approach:
- One of the first and simplest epidemiological model is the SI model. If you are in an infected state within this model, you stay infected forever. This one-way transition S → I is already not so realistic, at the same time the increasing number of infected cases can be obviously fit to the curves with the increasing number of cases, for example from this time series.
- Issues with more complex models, such as SIR or SEIR are more delicate. If we assume that in the beginning of an epidemic the fraction of susceptible individuals is still around 1 and that the number of infected people is exponentially growing as , we need to estimate the basic reproduction number of the virus as . The intuitive explanation of the basic reproduction number is the following: if each person is able to spread the disease to at least one other person before recovering, then the epidemic can continue, otherwise, if , the epidemic dies off. The problem here is that in reality, the population is very heterogeneous: there is no such thing as a single type of infectious individual since immune systems are varying from one person to another. Hence even if estimated from the empirical data, the true value will still need to be estimated with the confidence interval. For COVID-19 the estimates show that .
- In reality, one needs to take into account the information about so-called incubation period, the time elapsed between exposure to a pathogen, and when symptoms and signs start to appear. At the same time the effects from the mitigation policies, such as lockdown procedures, can be miscalculated using the simplest SIR models, since there total lockdown would correspond to , which gives unrealistic results.
For more publications on COVID-19 please see recent peer-reviewed articles, such as Fergusson et al.
Find more illustrations and code by B.Goncalves at Epidemiology101.
Predictions of next waves and role of mobility data
Prediction of next waves of COVID19 can be made using various models, including classical simple SIR model or its variations. One can already see the effect of curve flattening by changing the parameters, such as mobility index or average number of interactions between people. For spatial statistics one may need to use spatially embedded models, such as more generic and complex models Gleamviz model. For calibrating these models, one need to use expert-curated healthcare datasets available to COVID-19 researchers and data scientists. Mobility trends are important components for prediction of epidemics spreading, therefore we are working on making finding open mobility datasets. Below we collected some references of datasets and peer-reviewed publications:
- Mobility data on the effect of confinement on mobility: mobility data on mobility of people in residential, non-residential areas;
- Mobility data and predictions from the perspective of so-called "mobility index" can be found on the website of Cuebiq and illustrated in Gleamproject;
- Recent work and overview of predictive modeling of spreading COVID19 made and explained simply by B.Gonsalves, network scientist and physicist here;
- References resource of more than 128,000 scholarly articles about the novel coronavirus for use by the global research community: coronawhy.