A Primer on COVID‐19 Mathematical Models

Abstract

The emergence of severe acute respiratory syndrome coronavirus 2 (or coronavirus disease [COVID]-19) has led to a widespread global pandemic ((1)). COVID-19 symptoms and mortality are disproportionately more severe in people with obesity and obesity-related comorbidities ((2, 3)). This is of concern for the United States, where ~42% have obesity, and of these, 85% have type 2 diabetes. Without effective vaccines or therapies, there is an urgent need for prevention and treatment strategies for COVID-19 infections. Developing such strategies requires a thorough understanding of how COVID-19 infections spread in a population. Mathematical modeling predicts the spread of COVID-19 and permits testing different strategies to control and reduce population disease impacts. The New York Times published COVID-19 predictions from five mathematical models with widely varying results ((4)). A recent review ((5)) also highlighted model differences and noted that a fluid situation dependent on many external variables should not be viewed as a “crystal ball.” Why do COVID-19 models yield such variable predictions? How are the different COVID-19 models developed? How should models be applied? Here, we provide a brief tutorial addressing these questions. The majority of models fall into the following two categories, which are applied for different purposes: projections versus statistical forecasts. Projections are deterministic and they explain what could happen under a set of underlying hypotheses, while statistical forecasts use observed data to predict what will happen ((6)). In 1927, Kermack and McKendrick ((7)) developed the first continuous variable projection model of epidemic population dynamics, often referred to as a “susceptible, infected, and removed” (SIR) model. They compartmentalized a constant population into three states. The first state represents individuals susceptible to the disease with the number of susceptible individuals on day of the epidemic, denoted by The second state is the number of infected individuals on day of the epidemic, denoted by Finally, individuals removed from the infectious disease dynamics through death or recovery with immunity on day of the epidemic are denoted by The key property of projection models like the Kermack-McKendrick system is that they are based on logical assumptions of the underlying mechanics of a process and they can be developed without data to immediately address “what if” questions. Here, r declines as the number of infected individuals gets higher, reflecting increased social distancing during peak infectivity. On the other hand, increases when the number of infected individuals decrease. The values , K, and n are parameters that could be fit to data once available. The changes in assumptions alter the projections as depicted in Figure 2. Once data is fit to model parameters, SIR models can predict when infections peak or how high the peak may be, but as pointed out in Jewell et al. ((5)) and observed in Figure 2, these are dependent on the underlying model assumptions. Forecasting models are more useful after data has been gathered. For instance, most SIR COVID-19 models use fixed parameters, which lead to a constant referred to as or “ naught.” is defined as the number of secondary infections that will result from a single infected individual being introduced into a completely susceptible population. A calculated reproduction number during an ongoing epidemic is the “effective” reproduction number, , which is the observed average number of secondary cases from a single primary case per unit time. accounts for changes to through control measures such as social distancing. Within the SIR framework, , where is the time to recover and is the population size. Absent of interventions, the estimated of COVID-19 is between 1.5 and 6.7 ((9)). However, this value is not constant but changing daily. Renewal equations allow us to estimate values ((10)) of that can be fit to a statistical model. Similar to Massad et al. ((6)), we fit an exponential decay model to New York’s data yielding , though other statistical models such as one with an asymptote could be considered. Figure 3 depicts the fitted effective reproduction number versus the calculated values for New York. This forecast model will inform improved SIR model projections and generate prediction intervals. Projection models for the SARS-1 epidemic overestimate compared with statistical forecasting estimates of reproduction numbers ((6)). It is important that researchers use projection models for hypothetical scenario development and statistical forecasting models for forecasting, which are two different goals. Including the differential effects of COVID-19 on individuals with obesity into both projection and forecasting models by separating each compartment into normal weight and obesity populations will be key to understanding the population-wide impacts of COVID-19 on obesity. Obesity researchers can apply both types of models to evaluate prevention and treatment strategies for COVID-19 in persons with obesity.