Are we there yet? Decision-making under extreme uncertainty in the COVID-19 pandemic
Peter Solymos, Subhash Lele, Khalid Lemzouji, analythium.io, May 15, 2020
As we’ve all come to know, flattening the curve is the process of reducing the number of infections and easing the burden on the health system during the COVID-19 pandemic. Protective measures, such as travel restrictions, social distancing, staying home and basic hygiene, have been implemented with success to flatten the curve. Comparisons across countries with different approaches, timelines and levels of government intervention indicate that slowing the spread of the disease is possible via these intervention strategies. In this post, we explore how we can provide the best scientifically-based evidence for such decisions.
There is no shortage of data for COVID-19. Indeed, global counts of confirmed cases and deaths are seemingly everywhere. We hear daily updates on the news and we see digital dashboards showcasing global maps and curves of cases, which have proliferated following the data dashboard popularized by the Johns Hopkins University team (find a non-exhaustive list of apps here, see a list of global and regional data interfaces here and here). How does one use this data to better understand underlying mechanisms of the disease?
Number of total cases (top) and percent daily change (bottom) in the five health zones of Alberta. Data source: Provincial Surveillance Information, image source: https://hub.analythium.io/covidapp/.
To do this we must move beyond graphing the curve and implement an epidemiological model. Epidemiological experts base their analysis on SEIR epidemic models. These models are deterministic differential equation models, which classify the population in susceptible, exposed, infected, and recovered sub-groups. A nice discussion of those models is given by Dan Coombs in an earlier blog post, and Alison Hill has provided a detailed implementation of the SEIR model with extensive documentation and references here.
This SEIR model is composed of compartments based on the progression of the disease. The parameters of the SEIR model describe the rate at which individuals move into different compartments, i.e. as susceptible individuals get exposed, or infected individuals recover or die. Clinical observations are used to determine some of these parameters, for example time to recovery, or duration of ICU admission. The recently released scenarios by the Government of Alberta are timely examples of such models in action for emergency preparedness and capacity building.
Unfortunately, the parameters concerning the infected individuals with less severe symptoms or without any symptoms are uncertain, and estimating the prevalence and incubation periods for these groups is more difficult. There is also variation among individuals depending upon pre-existing conditions, gender and age, which can result in biased estimates, particularly when testing selectively focuses on high risk groups.
Another consequence of individual variation is that the incubation periods vary drastically between infected individuals, some individuals being still infectious after the mandatory 14-day quarantine period. Focusing on a strict deterministic rule and ignoring the individual variations may lead to an incomplete understanding of the risk associated with such rare events.
As an alternative to previous approaches, we formulate a hierarchical statistical model to link observation to epidemiological mechanisms. Hierarchical models are statistical abstractions for situations when, for example, the true numbers of new infections are not directly observable. The deterministic SEIR model captures how the infections are changing over time. However, this model doesn’t describe how the true infections are related to the observed number of cases; this is where a hierarchical model can help. A hierarchical model provides the means to connect observed and latent processes, and to estimate model parameters and their associated uncertainties. This is important, as failure to account for observation error can lead to biased and misleading estimates of epidemic parameters.
Our hierarchical model requires detailed data on the daily number of new cases within small geographic areas. We chose Alberta as a test case because it is close to home and such data has been published by Alberta Health for 132 local geographical areas. Yet the same model can easily be applied to other jurisdictions as well. Now let Yi,t be the observed number of daily cases in area i on day t, and Ni,t be the true number of daily infections in area i on day t. Let Ni,t be a random variable that follows a Poisson distribution with rate λi,t, or in shorthand: (Ni,t | λi,t) ~ Poisson(λi,t).
The variable λi,t is also the expected number of true cases per unit time (days in our example) from a Poisson process. The rate takes positive values and can depend on covariates: log(λi,t) = Xi,tTβ, where X is a design matrix with the covariates as columns and β is the coefficient vector describing the relationship between the covariates and the rate parameter. The rate parameter can be a function of area specific demographic and census data, e.g. proportion of age classes, occupations, the size density of the population.
The rate can also depend on time varying effects, such as interventions, weather, or the true number of cases on previous days. These temporal dynamics can capture community transmission within areas with already positive cases. The proximity to such hotspots can introduce spatial dependence and travel-related infections into the model. For example, the number of infections in the neighbourhood two weeks ago might influence the true cases today. We would also expect interaction between temporal and spatial effects, e.g. neighbourhood effects would have been larger before travel restrictions.
How does this relate to the observed cases? We can describe the number of observed cases (Yi,t) as another Poisson distributed random variable: (Yi,t | Ni,t) ~ Poisson(αi,t Ni,t), where αi,t is the proportion of infected cases reported. The αi,t variable takes value between 0 and 1 and can be modelled on the logit scale as a function of covariates, such as changes in testing and reporting over time, the number of tests performed, demographic composition of the sample tested, lag time involved in testing. We can also use analogues to inform this part of the model, e.g. past influenza infection rates. If the spreading mechanism for influenza and COVID-19 is similar (infection rates might be different), then past flu spread should help us predict COVID spread.
Statistical models allow to estimate the model parameters (point estimates) and their associated uncertainty. This uncertainty in our cases is not normally distributed because of the exponential transformation used to link the covariates to the support for the random variable. This formulation inherently allows for longer tailed distributions and therefore provides a more realistic depiction of risk due to extremes from longer tailed distributions. Such characteristics of statistical modelling are most useful when making predictions about unobserved cases for during forecasting.
The ability to forecast the true number of cases for the different geographic areas with associated prediction error is the true advantage of the hierarchical state-space model outlined above. Understanding the spatio-temporal dynamics of the disease spread would help to inform future policy. For example, we can use the model to evaluate different scenarios based on the estimated number of actual cases rather than the observed number of cases.
Imagine, for example, that in a rural area with 0 current confirmed cases, there might already be an unknown number of infections, which number can be a function of nearby areas. This estimate can give early warnings for nearby hospitals to expect more patients and be an important tool in emergency preparedness. This prediction can be made more spatially accurate and precise by incorporating demographic data, like age, occupation, existing health conditions.
The Model in Action
To build this model, we started gathering the data: we have daily case numbers from the 132 local areas in Alberta, the 2017 census data for population sizes, and the intervention dates by the Government of Alberta. We also include data from all countries around the world. We built a website (https://hub.analythium.io/covidapp/) where we collected relevant information about the COVID-19 pandemic.
Users can look at different countries and see how the cumulative cases compare to each other (see linked video at the end of the post). Furthermore, they can also drill down to the Canadian data to explore cases between provinces and territories to compare daily cases, rates, and doubling times with interactive visualization. We also added the Alberta data by health regions and we are now working on even finer resolution data. We relied heavily upon the data by the Johns Hopkins team and the data provided by the Government of Alberta. Our scripts automatically update the data every day. The app in action can be seen at the end of this post.
Number of total cases in 132 local geographic areas in Alberta (left) and over time (right). Data source: Provincial Surveillance Information, image source: https://hub.analythium.io/covidapp/.
The model can be refined by including key data, like demography, occupations, age distributions, etc. Other proxy information would also be very valuable, e.g. cell phone based average driving distances by regions. Another piece of information that would help to connect the true and observed cases is the reporting error for diseases with similar spreading mechanisms, i.e. influenza.
If you are interested in working with us to improve the website, help out with data collection, or modelling, please reach us at firstname.lastname@example.org. We would love to take the hierarchical model to the next level to inform decision making with the best evidence during extreme uncertainty in the COVID-19 pandemic.
YouTube video: https://youtu.be/AjzygcnuY0o