Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Development of an interactive, agent-based local stochastic model of COVID-19 transmission and evaluation of mitigation strategies illustrated for the state of Massachusetts, USA

  • Alexander Kirpich ,

    Contributed equally to this work with: Alexander Kirpich, Vladimir Koniukhovskii

    Roles Conceptualization, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliations Department of Population Health Sciences, Georgia State University, Atlanta, Georgia, United States of America, School of Public Health, Georgia State University, Atlanta, Georgia, United States of America

  • Vladimir Koniukhovskii ,

    Contributed equally to this work with: Alexander Kirpich, Vladimir Koniukhovskii

    Roles Conceptualization, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft

    Affiliation EPAM Systems, Saint Petersburg, Russia

  • Vladimir Shvartc,

    Roles Conceptualization, Data curation, Methodology, Writing – original draft

    Affiliation EPAM Systems, Saint Petersburg, Russia

  • Pavel Skums,

    Roles Conceptualization, Methodology, Visualization, Writing – original draft

    Affiliation Department of Computer Science, Georgia State University, Atlanta, Georgia, United States of America

  • Thomas A. Weppelmann,

    Roles Conceptualization, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Department of Internal Medicine, University of South Florida, Tampa, Florida, United States of America

  • Evgeny Imyanitov,

    Roles Conceptualization, Methodology, Writing – original draft

    Affiliation N.N. Petrov Research Institute of Oncology, Saint Petersburg, Russia

  • Semyon Semyonov,

    Roles Data curation, Project administration, Resources, Software, Visualization

    Affiliation Quantori, Cambridge, Massachusetts, United States of America

  • Konstantin Barsukov,

    Roles Data curation

    Affiliation EPAM Systems, Saint Petersburg, Russia

  • Yuriy Gankin

    Roles Conceptualization, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    yuriy.gankin@quantori.com

    Affiliation Quantori, Cambridge, Massachusetts, United States of America

Abstract

Since its discovery in the Hubei province of China, the global spread of the novel coronavirus SARS-CoV-2 has resulted in millions of COVID-19 cases and hundreds of thousands of deaths. The spread throughout Asia, Europe, and the Americas has presented one of the greatest infectious disease threats in recent history and has tested the capacity of global health infrastructures. Since no effective vaccine is available, isolation techniques to prevent infection such as home quarantine and social distancing while in public have remained the cornerstone of public health interventions. While government and health officials were charged with implementing stay-at-home strategies, many of which had little guidance as to the consequences of how quickly to begin them. Moreover, as the local epidemic curves have been flattened, the same officials must wrestle with when to ease or cease such restrictions as to not impose economic turmoil. To evaluate the effects of quarantine strategies during the initial epidemic, an agent based modeling framework was created to take into account local spread based on geographic and population data with a corresponding interactive desktop and web-based application. Using the state of Massachusetts in the United States of America, we have illustrated the consequences of implementing quarantines at different time points after the initial seeding of the state with COVID-19 cases. Furthermore, we suggest that this application can be adapted to other states, small countries, or regions within a country to provide decision makers with critical information necessary to best protect human health.

Introduction

The epidemic of a novel coronavirus was first detected in the city of Wuhan in the Chinese province Hubei on December of 2019 [14]. Despite the unprecedented efforts from Chinese authorities including the complete lockdown of the entire city of Wuhan on January 22, 2020 the virus has rapidly spread to all continents except Antarctica. The World Health Organization (WHO) officially declared the coronavirus a global pandemic on March 11, 2020 [5], only three months after the first case was detected. The novel coronavirus is now officially named SARS-CoV-2 and the disease caused by it has been called COVID-19 [6] to distinguish from SARS-CoV and the corresponding severe acute respiratory syndrome (SARS) pandemic from 2003 [7, 8]. Despite the much lower case-fatality rate, SARS-CoV-2 has caused morbidity and mortality orders of magnitude higher than severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) combined [9]. As of May 13, 2020, more than 4.1 million infections have been reported worldwide, with more than 287, 000 deaths due to complications of COVID-19 [10]. As of May, 2020 there is neither an effective virus-specific treatment, nor Food and Drug Administration (FDA) approved vaccine available for SARS-CoV-2 [1118]. As such, social distancing and quarantine are the only available measures to reduce the transmission and prevent overwhelming the capacity of existing healthcare systems. Social distancing in this context means practice of physical distancing between individuals to prevent transmissions. Quarantine involves more targeted actions and restrictions for those individual who are or are suspected to be infected with the goal of keeping them away from the others. Starting at the epicenter of Hubei province [19] in January, 2020, governments around the world have implemented society lockdown measures of varying degrees [2024]. Since such measures remain the only available tools to control the spread of the pandemic, it is critical to understand the transmission dynamics of SARS-CoV-2 in the population. This would allow for the prediction of COVID-19 cases and deaths over time under different mitigation strategies, which could be implemented to reduce morbidity and mortality; as well as the allocation of limited resources to medical providers.

To achieve this goal, multiple approaches can be implemented that are typically driven by the quality and precision of the available data. The most commonly reported data for epidemics are the incidence of new cases and deaths represented as a time series over the fixed intervals (e.g. days or weeks) aggregated across multiple regions and reporting sources [2527]. This aggregated data can be used for incidence curve reconstruction, modeling, and prediction when more detailed information about each infected individual is not available [2832]. Such models are called compartmental models [33] where the study population is divided into groups while individuals within each group are assumed to have the same characteristics of interest (e.g. susceptible, infected, vaccinated, or immune). The compartmental models are formulated via a defined system of differential equations that allows for both deterministic and stochastic formulations to quantify the uncertainty of the model fit. Those models provide insight into the underlying epidemic dynamics and allow for the prediction of future trends in incidence under different transmission scenarios, including interventions such as social distancing, quarantines, and vaccinations. Another approach are the agent-based models [3439], which utilize a synthetic population, that attempts to realistically represent the social interactions in time and space between individuals with different characteristics and infection status [36]. Compared to compartment models, agent-based models rely on detailed data about the study population, making them more computationally intensive while allowing for more realistic simulations of human transmission pathways.

An example of an agent-based model that has been developed to study influenza but has been successfully applied to study other respiratory diseases, including the efficacy of quarantine measures for containment of SARS-CoV-2 is the FluTE model [35, 40] which was specifically adopted with the name Corvid [38]. The FluTE model is based on the assumption that the synthetic study population is parsed in social subgroups and the interactions between them are modeled at different community levels (household, neighborhood, work), age groups, time of the day, and other characteristics. In this work, we propose an agent based model that can be applied by regional governments to local epidemiological settings, using incidence data from Massachusetts, USA. The proposed model incorporates the infected individuals that are reported in the beginning of an epidemic for a given area and their personal characteristics such as the date the infection was confirmation, geographic location, and demographic characteristics. Thus, by using a limited amount of information, such as the populations of each zip code and number of reported cases, a robust simulation is generated with time series data of predicted disease morbidity and mortality. Furthermore, by incorporating different quarantine strategies, the reduction in the number of new cases and deaths can be estimated for each locale depending on their unique characteristics. Since this model uses limited input data that are publicly available and is implemented in the form of an interactive web application, we believe that this tool could represent a widely adaptable format for state and local governments and health officials to make informed decisions as they consider easing or ceasing mandatory quarantines. Once effective treatments and vaccines become available, this framework could also be used to allocate treatment resources and plan vaccination campaigns tailored to fit different geographic regions. The reported demographics data can also be depersonalized in accordance to HIPAA regulations [41] to make the use of the model versatile and not to violate the privacy of individuals.

Materials and methods

The model structure

The model stochastic simulations are generating the infected individuals at different times and stages. Those individuals are indexed by k and have individual characteristics presented below: (1)

The characteristics of each individual Qk are determined and updated within the simulation process and have the following details:

  • (xk, yk)—the Cartesian coordinates (in pixels) of the individual Qk that do not change within the simulation process after they are introduced;
  • tinfk—the time of infection onset for the individual Qk which is measured in relation to the simulation baseline time denoted as 0;
  • detk—the detection time variable measured in days that corresponds to a period from an infection acquisition until the proper diagnosis and reporting of individual Qk;
  • stgk = (stg1(k), stg2(k), stg3(k))—the vector of durations of three disease infection stages measured in days that characterize the infectivity of the given individual during those stages. It is assumed that stg1(k) + stg2(k) + stg3(k) = detk;
  • agek—the age of the individual Qk at the time of the infection onset;
  • radk—the distance (in meters) up to which the individual Qk is able to infect the nearby individuals;
  • pcont(k)—the probability that during each day the individual Qk has any contacts which lead to new infections;
  • —the individual-specific parameters that define the distribution of the number of successful infection transmissions to other individuals within a given day. This number is generated randomly for each day t, provided that the individual has any transmissions on the given day (according to the contact probability pcont(k)).
  • R0(k)—the individual’s reproduction number. This variable stores the number of individuals that is infected by Qk during the infection period detk. The average of those reproduction numbers across individuals and simulations are used to estimate the basic reproduction number [42] which is a characteristic of the entire epidemic;
  • severk—the disease severity variable for the individual Qk that takes three values, where 1 corresponds to lethal, 2 corresponds to severe, and 3 corresponds to mild; the disease severity does not change for a given individual after it is determined randomly from a trinomial distribution;
  • durk—the disease duration from the infection onset to cure (or death) in days, which is generated randomly based on the severk parameter;
  • stk(t)—the status of the individual Qk at a given day t. The status of the individuals within the simulation is expected to change over time and is expected to take the following values:
    • stk(t) = 0—the individual is detected based on the external information i.e. from the reported data that are used as the model input;
    • stk(t) = 1—the individual is infected but has not been identified as such yet;
    • stk(t) = 2—the individual has been infected and detected as such, which has also implied the individual’s isolation (quarantine).
    • stk(t) = 3—the individual has recovered and is immune;
    • stk(t) = 4—the individual has deceased.

The asymptomatic infections are incorporated inside the model via three severities that are generated for each infected individual. This incorporation of asymptomatic infections into those with a “mild” severity status is made to improve the model computational tractability. Since the model parameters are calibrated based on the detected cases and the output of the model presents the detected cases as well the prediction abilities for detected and reported cases within the model are preserved. The severities are discussed in details in the manuscript supplement S1 Appendix. In the beginning of the simulations the model utilizes multiple local epidemic epicenters . Those epicenters serve as the model initial conditions and represent the introductory geographic points for the index cases that are introduced into the susceptible population. The epicenters can either correspond to the actual address coordinates for those places where the initial outbreaks were detected or to the centers of the corresponding aggregated geographic units. The latter may be the case, if either the exact infection acquisition locations are not known, or the privacy concerns prevent the inclusion of such data into the model. In the latter case the centers of the aggregated geographic units are taken as epicenters Ei for each .

The local epicenters in the model are defined by a pair of geographic coordinates (Lat, Long) and by an epicenter-specific region radius Ri which is defined in meters. Therefore, for the epicenter region is defined by a triplet: (2)

The epicenter regions are defined from the surveillance epidemiological data [43]. As the initial conditions in addition to the local epicenters the model incorporates the areas of high density for , where each Pj represents a large city or a densely populated area and which is also defined by a triplet: (3)

In the model the reporting times (days) for the initial index cases for each epicenter i precede the modeled epidemic starting time which corresponds to the baseline time slot t = 0. Therefore, the reported time slot indexes across the epicenters Ei are denoted as with the corresponding times . The earliest reported cases and their dates are used for the model input with indexes such that and the corresponding times . The corresponding number of confirmed and reported infections for each local epicenter Ei up to and including the time for is denoted as . The corresponding set of infected and reported (i.e. with the status stk(t) = 1) individuals across all epicenter is denoted as: (4) where is the global index for initial cases across all times and is the total number of the initial index cases that is simulated within the model based on the input data. The tilde notation for -s in D emphasizes the link to the model input data. The newly infected individuals are generated spatially in relation to the population area centers in the beginning and later in relation to the previously infected individuals. The individual’s coordinates are generated randomly using the distribution mechanism described in details in the manuscript supplement S1 Appendix.

The time index that corresponds to individual day within the model is denoted as t and is equal to 0 at the model baseline. The simulation baseline time t = 0 corresponds to the latest reporting time of the earliest reported cases that are used for the model input. The actual infection times for those index cases precede the selected baseline simulation time t = 0 due to the infectivity periods generated for those index cases prior to their reporting. The actual simulation starting time that accounts for the infectivity periods is denoted as t = Tmin and is smaller than the baseline time t = 0. This simulation starting time t = Tmin is generated within the model, while the baseline time t = 0 is defined by the data and is defined by the largest index within the set of calibration indexes . The largest simulation time t = Tmax is determined by the model user based on the desired length of prediction. The initial set of index cases D from (4) defines the model initial cases that are allocated across the local epicenters (2) at times up to the baseline time t = 0.

Based on the model geographic characteristics (2) and (3) and the initial set of reported individuals D from (4) the new lists of of modeled individuals are simulated for time slots t ∈ [Tmin;Tmax] where TmaxTmin + 1 is the total number of the simulated time slots. The simulated lists have the following format: (5) where the value of is defined by the simulation at every simulation time step t ∈ [Tmin;Tmax]. During this procedure the input set individuals D defined in (4) is allocated between the different epicenters and time slots within the lists defined by (5). The allocations of the set D is performed during the time slots where .

The overall model flow

The entire modeling process can be summarized via the following steps:

  • The model input time interval is determined by fixing the first reporting indexes out of the total where . Those indexes correspond to the reporting time slots . This completely defines the list of reported index cases D from (4) that are used as the model initial conditions. The baseline time of the model t = 0 is assumed to correspond to .
  • The individuals from the reported set D that are defined in (4) are assigned to the local epicenters of the future epidemic Ei for based on the available (from the input data) geographic distribution.
  • The geographic data about the areas of high density Pj for are incorporated into the model.
  • The model is initialized with the index cases from D. Based on those index cases that are defined in (4) the initial infection time Tmin is determined. This step is necessary to incorporate the infection times that have been present before the first reporting time into the model.
  • The final time point of the stochastic simulations Tmax is define by the user based on the desired study and prediction goals.
  • The initial list of infected individuals is initialized at time Tmin only with the earliest model input cases from the list D(Tmin).
  • The infected list of individuals for the time slot t + 1 is generated sequentially for all t ∈ [Tmin;Tmax − 1] based on the list of individuals from previous time slot t and the individual’s characteristics within the list are updated at this time step t + 1. The details of the new infection generations are provided in S1 Appendix.

Based on the lists at every time slot t ∈ [Tmin;Tmax] the infected modeled population summaries can be computed and summarized. In particular, the total number of currently infected but not identified individuals (i.e. those with the status st(t) = 1) is saved into Inf(t) variable for every t. The total number of treated or quarantined individuals (i.e. with the status st(t) = 2) is saved into Treat(t) variable for every t. The total number of recovered individuals (i.e. with the status st(t) = 2) is saved into Recov(t) variable for every t. The total number of deceased individuals (i.e. with the status st(t) = 4) is saved into Dead(t) variable for every t. Those numbers are used in the model calibration procedures, epidemiological summaries and in the model predictions. The model input utilizes only the first reported indexes with the corresponding reported times for with the total number of reported indexes equal to and . The remaining reported indexes for some integer τ are divided into the two groups: (6)

The first group of the reported indexes from (6) is used for the model calibration and estimation of the unknown parameters. The second group of the reported indexes from (6) is used to evaluate the quality of the model predictions. The Massachusetts surveillance data that are used for the model calibration, validation and predictions are freely available at the Massachusetts Department of Public Health web site [43]. The first reported date which corresponds to the time index in the model is March 13, 2020. The latest reported date that is used for the model input is March 26, 2020 which corresponds to the time index in the model. The time indexes that correspond to and are April 14, 2020 and April 22, 2020 respectively. The parameter optimization is performed by minimizing the sum of squared differences between the model-produced outputs and the calibration data by using the Nelder–Mead numerical minimization method [44]. The additional details about the model formulation, parameterization, and calibration are summarized in the S1 Appendix.

After the model calibration is performed various quarantine and transmission intervention strategies are investigated within the calibrated model. Those intervention strategies are based on the assumption that the probability of contacts between the individuals decreases after the quarantine measures are enforced. Within the model this is implemented by an immediate change in the contact probability parameter starting from a certain calendar date. The earlier implementation of the quarantine in the model is represented by an earlier calendar date when the change in parameter value occurs. Those earlier quarantine dates result in the smaller transmission probability and fewer infected individuals in comparison to the later quarantine dates. The quarantine date in the end of the simulation period corresponds to no quarantine scenario. As an example, three scenarios with different quarantine start dates are presented and the corresponding reduction in the number of cases is discussed. There is also an option to consider other quarantine enforcement dates interactively within the current model implementation.

The model availability

The model has been implemented in multiple environments which include the application tool for Microsoft Windows [45] and the web prediction tool [46, 47] (summaries only). The model application tool for Microsoft Windows is freely available under the terms of the MIT license [48]. The tool source code, the application, and the relevant documentation are available on GitHub [45]. The current model tool has been calibrated based on the state of Massachusetts (United States) incidence data [43]. In addition to that the user has an option to adjust interactively the tool parameters which include, in particular, the transmission parameters and the quarantine implementation dates(s). Overall, the proposed framework and the code are fairly general and can be adopted for other areas and territories where the demographics of the incidence cases and population characteristics are known with at least some geographic precision, and where the rapid evaluations of social distancing measures have to be quantified.

Results

Within the model multiple epidemic progression scenarios can be considered. In particular, three different quarantine strategies are presented in this manuscript as an illustration of the model. The alternative quarantine scenarios can be produced and customized interactively within the model application tool if necessary. The difference between the presented scenarios is in the quarantine date at which the transmission probability parameter changes to smaller one. The smaller transmission parameter values result in fewer infections therefore the earlier quarantine dates result in less infections overall in comparison to the later dates. The predicted numbers of infections for each scenario are defined by the quarantine implementation date and those number can be compared. The first scenario corresponds to the quarantine date on March 29, 2020 i.e. the early reduction in contact probabilities and social distancing between individuals. The second scenario assumes the implementation of the quarantine measures on April 6, 2020, and the third scenario assumes the implementation of the quarantine measures on April 13, 2020. The point estimates and the prediction bands have been produced by replicating each of the three scenarios and taking the median values across 500 model runs for the point estimates and 5-th and 95-th percentiles for the 90% prediction intervals. The results are summarized in Tables 1 and 2 for the model-predicted cases and deaths, respectively. For example, the summaries from Table 1 can be compared after one month of the baseline date i.e. on April 26, 2020. For the first scenario the model predicts 24, 039 cumulative cases (with the 90% PI (20, 665;27, 296)), for the second scenario the model predicts 56, 587 cumulative cases (with the 90% PI (46, 944;66, 401)), and for the thirds scenario the model predicts 123, 351 cumulative cases (with the 90% PI (100, 113;144, 018)). Compared to the quarantine start date in the second scenario, the first scenario results in 58% reduction in cumulative cases on April 26, 2020, in 63% reduction in cumulative cases on May 26, 2020, and in 65% reduction in cumulative cases on June 26, 2020. Compared to the quarantine start date in the third scenario, the first scenario results in 81% reduction in cumulative cases on April 26, 2020, in 87% reduction in cumulative cases on May 26, 2020, and in 88% reduction in cumulative cases on June 26, 2020. Based on the model outputs the earliest quarantine measures and the reduction in contact probabilities can be extremely beneficial in mitigation of the outbreak consequences. The analogues summaries for the model-predicted death across the three scenarios are summarized in Table 2. Compared to the quarantine start date in the second scenario, the first scenario results in 52% reduction in cumulative deaths on April 26, 2020, in 63% reduction in cumulative deaths on May 26, 2020, and in 65% reduction in cumulative deaths on June 26, 2020. Compared to the quarantine start date in the third scenario, the first scenario results in 70% reduction in cumulative deaths on April 26, 2020, in 86% reduction in cumulative deaths on May 26, 2020, and in 88% reduction in cumulative deaths on June 26, 2020.

thumbnail
Table 1. The predicted number of cumulative cases produced by the model over time for three different quarantine scenarios and three time periods together with the corresponding 90% prediction intervals.

https://doi.org/10.1371/journal.pone.0247182.t001

thumbnail
Table 2. The predicted number of cumulative death produced by the model over time for three different quarantine scenarios and three time periods together with the corresponding 90% prediction intervals.

https://doi.org/10.1371/journal.pone.0247182.t002

The model is presented via the graphic user interface (GUI) application for Microsoft Windows as well as the as web prediction tool [46, 47] (summaries only) that can be used for the geographic visualization of various epidemiological curves and geographic visualization of cases in the state of Massachusetts. Users can work with the tool and utilize the available model customizations. The appearance of the GUI for MS Windows application is presented in Fig 1. The state of Massachusetts MassGIS data were used to produce the tool map. The data are available for download and are public records from the Bureau of Geographic Information (MassGIS), Commonwealth of Massachusetts, Executive Office of Technology and Security Services [49].

thumbnail
Fig 1. The tool graphic user interface.

The visualization tool GUI for Microsoft Windows. The tool provides the geographic visualizations of epidemic on the state of Massachusetts map and constructs the epidemiological summaries and curves. The tool allows interactive model calibrations and step-by-step simulations.

https://doi.org/10.1371/journal.pone.0247182.g001

The presented application does not require installation and can be launched directly by running the executable file. The tool allows the parameters optimizations and the visual animated simulations of the model outputs and their comparison with the reported data. The user can interactively customize the most important simulation parameters, change the duration of prediction, and adjust manually the locations of high population density. This allows the user to consider multiple scenarios of the epidemic spread. In addition, the user has flexibility to re-run the model multiple times either step by step or entirely for all time slots. The user can scroll via each day of the recently completed simulation to see the visualized results of that specific day and choose which epidemiological curves to include in the summary graphs. The modeled cases and other summaries are saved into the comma-separated values (csv) files after the end of each simulation. The tool also provides the estimate of the population basic reproduction number for each simulation run together with the corresponding 90% confidence intervals. The estimates for are provided based on the quantiles of the individual’s reproduction numbers R0(k) across multiple individuals k both before and after the quarantine date that is defined by the user. The resulting distribution of the individual’s reproduction numbers R0(k) from the tool before and after the quarantine are summarized in Fig 2.

thumbnail
Fig 2. The distribution of the individual’s reproduction numbers R0(k).

The example of the tool output for the distribution of the individual’s reproduction numbers R0(k). The output graph contains the estimated probability density function of the individual’s reproduction numbers R0(k) together with the cumulative distribution functions both before and after the quarantine implementation date.

https://doi.org/10.1371/journal.pone.0247182.g002

The example of the summary graphs for the model-produced outputs for the second scenario from Tables 1 and 2 are presented in Fig 3, which contains the four combined graphs available in the “Statistics” tab in the top right corner of the tool. Those graphs within the tool can be produced by setting the “Max Simulation Time” and “Forecast Day” fields to July 15, 2020 and by running the model 500 times by using the “Daily Forecast Evaluation” button. The 500 runs are necessary to produce the median predictions and the corresponding 90% uncertainty prediction bands across those runs by taking the 5-th and the 95-th percentiles across those situations for each time slot. Those graphs include the cumulative numbers of reported cases and deaths, together with the currently hospitalized patients and unreported cases. The graphs also include the reported data in blue. The calibration time period is highlighted in blue and is bounded by vertical bars.

thumbnail
Fig 3. The model-produced predictions.

The median of the model-produced 500 runs together with the corresponding 90% uncertainty prediction bands for different model outputs. The top left graph includes the cumulative numbers of reported cases. The top-right graph summarizes the cumulative subset of the reported cases that have deceased. The bottom graphs summarize the number of hospitalized and unreported patients in the given moment of time. The reported data are displayed in blue for visual comparison.

https://doi.org/10.1371/journal.pone.0247182.g003

Discussion

In this work the local agent-based modeling framework for respiratory diseases has been presented. This framework incorporates the reported geographic incidence data that are typically available from surveillance, which include individual’s age, infection status, and the severity of the disease. The model accounts for the latent period of the individual’s infection before detection and proper reporting as well as for different disease severity levels. The framework also allows to incorporate the exact geographic addresses of individuals (if available) or the random geographic distribution of individuals within those aggregated districts where they are reported in case of privacy concerns. The model allows to perform predictions with different levels of social isolation between individuals and quarantine measures. Those measures are implemented at different times to compare different quarantine scenarios. As expected, there was a decrease in the cumulative incidence and deaths inversely proportional to the date quarantine was implemented; which resulted in approximately 50-80% reduction in cases and deaths depending on the scenario. This agrees with already published results that strict social distancing combined with proper testing will keep the disease at level that does not overwhelm the capacity of heath care system [18].

Compared to complex agent-based models, the compartmental models are based on the assumptions of homogeneous mixing and can be parameterized by a relatively small set of rates and initial conditions. The main challenges for the compartmental models [30, 31, 50] are the determination of the compartment types that are used in the model, the assignment of individuals between compartment i.e. the specifications of the set of rules that assign each particular individual to each type of compartments, and the determination of the parameters of interest which can either be postulated from external sources or estimated from data.

The agent based models, due to their inherited complexity, incorporate separate individuals with multiple different characteristics and parameters per individual. This adds another layer of parameterization flexibility, but also introduces another layer of modeling challenges, since the number of individual’s characteristics within the model is determined by the modeler [18, 3436]. This typically implies some additional assumptions about the model behavior and the overall model parametrization. Those assumptions are expected to be region-specific, since the transmission times, patterns, and other characteristics typically wary from region to region. The respiratory infections, which include SARS-CoV-2, add another level of complexity to any model due to quality of reported data. The asymptomatic cases are typically neither identified nor properly recorded and surveillance systems for symptomatic cases are never perfect either. As a result a substantial (but unknown) amount of cases is not reported and the modeler has to account for those asymptomatic and unreported symptomatic cases who still participate in the disease transmission process.

Ideally, the model has to be: 1) flexible enough to incorporate the possible social and geographic characteristics of individuals and to provide the way to realistically represent the social interactions and the disease transmission mechanisms; 2) simple enough to avoid the problems with parameterization, but able to capture the actual transmission patterns with the goal of predictions and intervention studies; 3) utilize the available surveillance and public health data in the best possible way. The best possible way in this context means, that all the information from the data that can be used to answer the questions of interest are utilized, while the number of assumptions within the model beyond the information available from the data is the smallest possible that is necessary to implement the model.

In the case of COVID-19, an epidemic which has quickly evolved into a pandemic, the local epidemic developments in every region are expected to have different dynamics influenced by multiple region-specific factors. Thus, an agent based model which utilizes local settings is likely superior to a global agent-based model in this setting and can be implemented with minimal inputs as long as local data are available. In this example, we chose regional data for the state of Massachusetts, however we believe this framework and interactive tool could be adopted and useful for small or middle size countries or other administrative districts within a larger country, that have comparable reporting and data quality across different administrative regions.

Conclusion

In this paper, we have presented a novel, localized agent-based model that can be used within minimal input data, which is publicly available and tailored to the population distributions of Massachusetts, USA. After calibration the model provided a good estimation of the actual incidence, hospitalizations, and death rates, with the added benefit of estimating the number of undetected infections in the population. Given the necessity for making decisions of easing or ceasing quarantines that are specific to a state or county based on their reported case counts, adaptation of this framework could prove to be very useful with efforts to reopen the economy, while quantifying the disease burden posed by such decisions. In addition, this model could be used for future outbreaks of other novel respiratory diseases to protect public health and possibly designed tailored interventions of treatment and vaccination campaigns.

Supporting information

S1 Appendix. Model details.

The details about the model formulation, parameterization, and calibration.

https://doi.org/10.1371/journal.pone.0247182.s001

(PDF)

References

  1. 1. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The Lancet. 2020;395(10223):497–506.
  2. 2. Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. The Lancet. 2020;395(10223):507–513. pmid:32007143
  3. 3. Munster VJ, Koopmans M, van Doremalen N, van Riel D, de Wit E. A novel coronavirus emerging in China—key questions for impact assessment. New England Journal of Medicine. 2020;382(8):692–694.
  4. 4. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A novel coronavirus from patients with pneumonia in China, 2019. New England Journal of Medicine. 2020;.
  5. 5. WHO;. https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-COVID-19---11-march-2020.
  6. 6. WHO;. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/naming-the-coronavirus-disease-(covid-2019)-and-the-virus-that-causes-it.
  7. 7. Peiris JS, Yuen KY, Osterhaus AD, Stöhr K. The severe acute respiratory syndrome. New England Journal of Medicine. 2003;349(25):2431–2441.
  8. 8. WHO;. https://www.who.int/ith/diseases/sars/en/.
  9. 9. Mahase E. Coronavirus: COVID-19 has killed more people than SARS and MERS combined, despite lower case fatality rate. BMJ. 2020;. pmid:32071063
  10. 10. WHO;. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports.
  11. 11. Jin YH, Cai L, Cheng ZS, Cheng H, Deng T, Fan YP, et al. A rapid advice guideline for the diagnosis and treatment of 2019 novel coronavirus (2019-nCoV) infected pneumonia (standard version). Military Medical Research. 2020;7(1):4. pmid:32029004
  12. 12. Yu F, Du L, Ojcius DM, Pan C, Jiang S. Measures for diagnosing and treating infections by a novel coronavirus responsible for a pneumonia outbreak originating in Wuhan, China. Microbes and infection. 2020;. pmid:32017984
  13. 13. Chen WH, Strych U, Hotez PJ, Bottazzi ME. The SARS-CoV-2 vaccine pipeline: an overview. Current Tropical Medicine Reports. 2020; p. 1–4.
  14. 14. Amanat F, Krammer F. SARS-CoV-2 vaccines: status report. Immunity. 2020;. pmid:32259480
  15. 15. Bicher MR, Rippinger C, Urach C, Brunmeir D, Siebert U, Popper N. Agent-Based Simulation for Evaluation of Contact-Tracing Policies Against the Spread of SARS-CoV-2. medRxiv. 2020.
  16. 16. Silva PC, Batista PV, Lima HS, Alves MA, Guimarães FG, Silva RC. COVID-ABS: An Agent-Based Model of COVID-19 Epidemic to Simulate Health and Economic Effects of Social Distancing Interventions. Chaos, Solitons & Fractals. 2020; p. 110088.
  17. 17. MS N, Jammy GR, Pant R, Choudhury L, R A, Yeldandi V, et al. An Agent Based Model methodology for assessing spread and health systems burden for Covid-19 using a synthetic population from India. medRxiv. 2020.
  18. 18. Aleta A, Martin-Corral D, y Piontti AP, Ajelli M, Litvinova M, Chinazzi M, et al. Modeling the impact of social distancing, testing, contact tracing and household quarantine on second-wave scenarios of the COVID-19 epidemic. medRxiv. 2020;. pmid:32511536
  19. 19. Reuters;. https://www.reuters.com/article/us-china-health-who-idUSKBN1ZM1G9.
  20. 20. The Wall Street Journal;. https://www.wsj.com/articles/italy-bolsters-quarantine-checks-after-initial-lockdown-confusion-11583756737.
  21. 21. The New York Times;. https://www.nytimes.com/2020/03/22/world/europe/germany-coronavirus-budget.html.
  22. 22. The Guardian;. https://www.theguardian.com/world/2020/mar/28/covid-19-may-be-peaking-in-parts-of-spain-says-official.
  23. 23. The New York Times;. https://www.nytimes.com/2020/03/24/world/asia/india-coronavirus-lockdown.html.
  24. 24. The Atlantic;. https://www.theatlantic.com/ideas/archive/2020/03/why-theres-no-national-lockdown/609127/.
  25. 25. Google: Coronavirus (COVID-19);. https://news.google.com/covid19/map?hl=en-US&gl=US&ceid=US:en.
  26. 26. Johns Hopkins University: COVID-19 Case Tracker;. https://coronavirus.jhu.edu/map.html.
  27. 27. WHO: Daily Situation Reports;. https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports.
  28. 28. Flaxman S, Mishra S, Gandy A, Unwin H, Coupland H, Mellan T, et al. Report 13: Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries. 2020;.
  29. 29. Kucharski AJ, Russell TW, Diamond C, Liu Y, Edmunds J, Funk S, et al. Early dynamics of transmission and control of COVID-19: a mathematical modelling study. The lancet infectious diseases. 2020;. pmid:32171059
  30. 30. Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, Hyman J, et al. Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020. Infectious Disease Modelling. 2020;5:256–263. pmid:32110742
  31. 31. Roosa K, Lee Y, Luo R, Kirpich A, Rothenberg R, Hyman JM, et al. Short-term forecasts of the COVID-19 epidemic in Guangdong and Zhejiang, China: February 13–23, 2020. Journal of Clinical Medicine. 2020;9(2):596. pmid:32098289
  32. 32. Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. The Lancet Public Health. 2020;. pmid:32220655
  33. 33. Brauer F. Compartmental models in epidemiology. In: Mathematical epidemiology. Springer; 2008. p. 19–79.
  34. 34. Venkatramanan S, Lewis B, Chen J, Higdon D, Vullikanti A, Marathe M. Using data-driven agent-based models for forecasting emerging infectious diseases. Epidemics. 2018;22:43–49.
  35. 35. Chao DL, Halloran ME, Obenchain VJ, Longini IM Jr. FluTE, a publicly available stochastic influenza epidemic simulation model. PLoS computational biology. 2010;6(1). pmid:20126529
  36. 36. Eubank S, Guclu H, Kumar VA, Marathe MV, Srinivasan A, Toroczkai Z, et al. Modelling disease outbreaks in realistic urban social networks. Nature. 2004;429(6988):180–184. pmid:15141212
  37. 37. Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science. 2020;368(6489):395–400. pmid:32144116
  38. 38. Chao DL, Oron AP, Srikrishna D, Famulare M. Modeling layered non-pharmaceutical interventions against SARS-CoV-2 in the United States with Corvid. medRxiv. 2020.
  39. 39. Hoertel N, Blachier M, Blanco C, Olfson M, Massetti M, Rico MS, et al. A stochastic agent-based model of the SARS-CoV-2 epidemic in France. Nature Medicine. 2020; p. 1–5.
  40. 40. Koo JR, Cook AR, Park M, Sun Y, Sun H, Lim JT, et al. Interventions to mitigate early spread of SARS-CoV-2 in Singapore: a modelling study. The Lancet Infectious Diseases. 2020;. pmid:32213332
  41. 41. Edemekong PF, Haydel MJ. Health Insurance Portability and Accountability Act (HIPAA). In: StatPearls [Internet]. StatPearls Publishing; 2019.
  42. 42. Milligan GN, Barrett AD. An Essential Guide. Wiley Online Library; 2015.
  43. 43. www.mass.gov;. https://www.mass.gov/info-details/archive-of-covid-19-cases-in-massachusetts.
  44. 44. Nelder JA, Mead R. A simplex method for function minimization. The computer journal. 1965;7(4):308–313.
  45. 45. github.com;. https://github.com/quantori/COVID19-MA-Transmission.
  46. 46. tableau.com;. https://public.tableau.com/profile/quantori#!/vizhome/ModelingCovidMassachusets/COVID-19Massachusettsforecast.
  47. 47. tableau.com;. https://public.tableau.com/profile/quantori#!/vizhome/ModelingCovidMassachusets_Map/InteractiveMap.
  48. 48. MIT Licence;. https://opensource.org/licenses/MIT.
  49. 49. MassGIS;. https://www.mass.gov/get-massgis-data.
  50. 50. Chowell G, Tariq A, Hyman JM. A novel sub-epidemic modeling framework for short-term forecasting epidemic waves. BMC medicine. 2019;17(1):164.