Modeling an Apocalypse

How we use mathematics to describe a pandemic

Modeling an Apocalypse
Flattening the curve: S is the total susceptible population; I is the proportion infected; R refers to those previously infected and no longer susceptible. Based on the SIR epidemiological model.

This article provides general information about disease modeling. It does not include details specific to COVID-19.

“For God and Country!” has become passée. The new battle cry is “Flatten the curve!” We math nerds are clutching our pocket protectors in anticipation. These are words we can get behind.

So, what do they mean?

The “curve” we are called to flatten is part of an epidemiological model. Such models take our best guesses about a disease and swirl these guesses together. The guesses interact with one another and evolve. We then tweak the guesses and see how our tweaks affect the big picture. We’ll look at one of these models.

The new battle cry.

The SIR Model

This model is as simple as it gets. Each member of the population falls into one of three categories (or compartments):

  • Susceptible (S) — Most of the population begins in this category. Each member has potential either to join the Infected category or to ride out the situation unaffected. The size of this category can only shrink as members become infected.
  • Infected (I) — A few members of the population begin Infected. The only direction to go is to the Removed category. The size of this category may grow or shrink, depending on rates of infection and removal.
  • Removal (R) — Once an individual is no longer infected, that person joins the final category. This model doesn’t distinguish between a healthy recovery and a death. The size of this category can only grow as the disease progresses through the population.

The SIR model assumes the disease flows in one direction: S → I → R. Two educated guesses are necessary to calculate the changes in these variables over time:

  • Infection Rate (I-rate) — The infection rate tells us, on average, with how many individuals an infected person shares the disease. This variable is the one we are targeting when we “flatten the curve.” A rate of 1 means the disease is stable. Each person passes the disease to only one other person. Once that person recovers (or dies), they leave behind a replacement. This is like a population where the birth rate is two children per family. The population neither increases nor decreases in size. A rate greater than 1 means the disease is spreading; less than 1 means the disease has lost momentum.
  • Removal Rate (R-rate) — The removal rate is a number between 0 and 1, which represents how quickly the infected population recovers.
In the SIR model, the disease flows one way. Once an individual is no longer infected, that individual is removed from the model. A removed individual cannot infect others nor become reinfected.

The model is run as a dynamic computer simulation. The size of any one category (S, I, R) continuously updates, depending upon the other two. Three overlapping curves on a graph illustrate this interplay. The so-called curve flattening refers to the effect of reducing the number of individuals each person infects.

Compare the I curves between the two graphs. On the left, the steep incline shows a higher infection rate. On the right, disease progress has slowed; the curve has been “flattened.” The disease takes longer to run its course, but it infects fewer individuals. The dashed line helps compare the heights of the two curves.

Carriers

A modification to the SIR model adds the role of carriers. A carrier recovers from the disease in the sense that she is no longer symptomatic. She has the potential, however, to re-enter the infectious state. Such was the case with the infamous Typhoid Mary.

A carrier doesn’t recover from the disease, but remains asymptomatic with the potential to infect others.

The SIR model is the most basic of compartmental disease models. Other, more complex models build on it. The next time you hear the battle cry to “Flatten the curve,” you’ll have an idea what that looks like. At least from a mathematical perspective.