Abstract
Objectives To evaluate the transmission dynamics and health systems burden of COVID-19 with and without interventions, using an Agent Based Modeling (ABM) approach on a localized synthetic population.
Study design A synthetic population of Rangareddy district, Telangana state, India, with 5,48,323 agents and simulated using an ABM approach for three different scenarios.
Methods The patterns and trends of the COVID-19 in terms of infected, admitted, critical cases requiring intensive care and/ or ventilator support, mortality and recovery were examined. The model was simulated over a period of 365 days for a no lockdown scenario and two Non-Pharmaceutical Intervention (NPI) scenarios i.e., 50% lockdown and 75% lockdown scenarios.
Results Results revealed that the peak values and slope of the curve declined as NPI became more stringent. The peak values could facilitate policymakers to plan the required capacity to accommodate influx of hospitalizations.
Conclusions ABM provides better insight into projections compared to compartmental models. The results could provide a platform for researchers and modelers to explore using ABM approach for COVID-19 projections with inclusion of interventions and health system preparedness.
1. Introduction
The first reported case of the novel coronavirus (COVID-19 or SARS-CoV-2) in India dates back to January 30, 2020 when it was also announced as pandemic by WHO1,2. Since then, epidemic has spread across India infecting 44,94,389 people with 9,27,545 active cases, 34,90,908 recovered cases and 75,328 deaths as on Sep 10, 2020 3. Globally, COVID-19 has spread across 213 nations, infecting 28,035,700 people worldwide and claiming 908,991 lives as on Sep 10, 2020, posing a global health emergency 4–6.
In a country like India having a denser population, the situation poses a serious challenge 1. There are several underlying factors such as age, comorbidities, exposure to air pollution, amount of exposure to virus, etc., that determine transmission dynamics 8. Travel restrictions and 9,10other Non-Pharmaceutical Interventions (NPIs) such as restrictions on public gatherings, intra-city movements, etc., may flatten the curve. It is important to understand the time required between exposure and complete recovery, using these interventions, to take timely responsive actions against COVID-19 11.
Modelling is an effective tool to provide real-time projections and obviates the burden of investing time, cost and risks 12. Several models on COVID-19 projections exisit and most of them are compartmental models. Since compartmental models assume homogeneity in the compartments they do not account for individual level variations and interactions within the system, Agent Based Modeling (ABM) is a tool which considers these factors and may provide better insights 13. ABM considers the interaction between agents whilst also distinguishing them based on their individual parameters 14,15. Advancements in computational capabilities have increased interest among researchers across various verticals, including public health in recent times on ABM 16–18.
ABM incorporates actions of agents within the system, helping the model comprehend infection spread dynamics better 19,20. ABM follows a bottom-up approach wherein result of behaviour of individuals, defined as agents within the system 21,22. ABM allows to define unique characteristics to the agents to make each of them behave distinctly 22. Mixing patterns among the agents within a system play a vital role in dynamic transmission models for close contact infections23.
In the past, ABMs have been employed to address various infectious diseases such as, a bioterrorist introduction of smallpox 24, design of vaccination strategies for influenza 25, curtail transmission of measles through contact tracing and quarantine 26, control of tuberculosis 27, implementation of distancing measures and antiviral prophylaxis to control H5N1 influenza A (bird flu)28 and devise evacuation strategies in the event of airborne contamination 29.
An evaluation of use of ABM on COVID-19 globally suggests its use to measure the effects of lockdown on transmission dynamics 30–33, post-lockdown measures 34, use of control measures (face mask, social distancing)34,35, isolation of vulnerable proportion of population 34,36, contact tracing, intelligence of agents (based on awareness level 31 or protection level), contact tracing measures 31,35, good practices such as sneezing into one’s hands 32, both direct (upon contact) and indirect transmission (through suspended particles) 32, scheduled-based contacts 32–34,37 with close circle and in work, transport and public places 32–36,38; viral-load based transmissibility35, examination of genomic sequencing to determine the spread 39, etc.
In India, several COVID-19 models have been conducted based on Susceptible (S), Exposed (E), Infective (I) and Recovered (R) (SEIR) 40–43, Susceptible (S), Exposed (E), Symptomatic (I), Purely Asymptomatic (P), Hospitalized or Quarantined (H), Recovered (R) and Deceased (D) (SIPHERD)44, mathematical models 45,46, etc., to compare the spread during lockdown and no lockdown scenarios. Nation-wide models restrict the policymakers locally to devise strategies based on the results as they might not fit properly to the locality 47. In countries like India where people are diverse in all respects like population dynamics, contact network, migrating population, nature of work, etc., local models might prove effective and would assist policymakers to take local decisions for disease mitigation. The present study aims at an ABM to examine the patterns and trends of spread and the effect of NPIs in a synthetic population of a region in Telangana state of India.
2. Methods
A synthetic population of Telangana state, India, has been developed, the details of it are presented elsewhere 48. We used 5,48,323 agents from this population from the Rangareddy district.. The modeling follows an ABM approach using AnyLogic 8.5.2 University edition to model the interaction environment 49. The entire simulation and reporting follows the International Society for Pharmacoeconomics and Outcomes Research (ISPOR-SMDM) Modeling Good Research Practices and ethical good practice in modelling 50–52. These guidelines were used so that the assumptions, scope and shortcomings of the model are transparent to the readers and policymakers.
2.1. Synthetic population
Synthetic population is one of the commonly used approaches to represent a group of people, preserving the confidentiality of individuals. Synthetic population has statistical equivalence with the original population being represented and is indistinguishable from the census data 53. For this study, we used a synthetic population developed for Rangareddy, Telangana state consisting of 5,48,323 people (representing 10.35 % of Rangareddy’s population (n-52,96,741) as per Census of India 2011), to demonstrate the ABM 48,54. The population was categorized based on age as less than 5, 5 to 59, and above 60 with 47,039, 4,59,372 and 41,912 agents respectively8,23.
2.2. Transmission rates of COVID-19
Based on the WHO report on COVID-19 (16 to 24 Feb 2020), the transmission rates were set to vary from 1 to 10 percent55.
2.3. Contact network
The contact network plays a vital role in transmission dynamics. For contact rate estimation, we used a study from Ballabgarh, India, which determined the contact rates for close contact infections 56. The dataset representing the number of people met by each individual was input into the ‘Input Analyzer’ tool of Arena (Version 16.00.00002). Input Analyzer provides fitting distributions with associated errors of fitting. The contact rates of each group were found to follow normal distributions using "Input Analyzer" (Table 1).
Contact rate distribution of Ballabgarh 56
The population densities of various townships of India were calculated based on Indian population proportion (table 2)57,58. The population density of Ballabgarh was 551 people per square kilometer 56. This was used to proportionately determine the contact rates of towns based on their respective population densities, assuming Density-Dependent (DD) contact rate59,60
For estimation of overall average contact rate, weights equal to the corresponding proportions of people living in different townships were multiplied to their corresponding multiplication factors (slope in equation (1)) (table 2). Contact rate was thus derived as the product of this slope and the ratio of the population density of a particular township to the population density of Ballabgarh.
Classification of Indian towns with population proportions and densities 57
The number of contacts made by people under each town category were multiplied by their respective multiplication factors as reported in Table 2. These values were integrated to Input Analyzer and their respective contact rate distributions were then determined. Contact rate for the present study followed Lognormal distributions (Lognormal(µ, σ, Min)) (table 3).
2.4. State chart
A state chart represents the various states in which an agent would exist (figure 1) 61. The initial state of all agents healthy, during the start of simulation. The agents would interact with other agents in the population and transmit the infection. Infected agents undergo an incubation period and turn out to be either symptomatic or asymptomatic. They continue to contact other agents and transmit till they get admitted or recover. Once admitted, agents undergo treatment and either decease or recover whilst in any of the three levels of infections represented by ‘admitted’, ‘ICU’ and ‘ventilator’ states 21.
State chart for agent(s) (people)
2.5. Model parameters
Several sources including Models of Infectious Disease Agent Study (MIDAS) has been used as a source of acquiring parameters through the pre-prints and manuscripts available 62. Table 3 details the various parameters used for the model.
Parameters for the model
2.6. Model scenarios
In order to study the effect of minimization of number of contacts among the agents, three different scenarios were simulated for 365 days, results of which are discussed subsequently.
First scenario, the "no lockdown" scenario includes no NPI put in place. In order to study the effect of lockdown, the contact rate of people needs to be reduced. To achieve this, the results of a study that presented the proportion of contacts made by individuals of different age groups were utilized. The number of contacts made at different locations namely home, school, work and others, as designated by the authors of the study was used to enact the lockdown 68. To simulate the NPI scenarios, the number of people met in work and other places were reduced by 75 percent and 50 percent for the two scenarios whereas the contacts in school was completely discarded owing to the closure of schools. These resulted in reduced diffusion of the infection across the population, results of which are discussed subsequently.
3. Results
Simulations were run for different age groups as per the categorization for all the three scenarios. Detailed day-wise data of the number of people in each health state are provided in the supplementary files. The results of the same and their interpretations are discussed below.
It was observed that (Figure 2), number of uninfected people declines as the stringency of the imposed lockdown increases. After a duration of one year, proportion of people who remain uninfected was 28.53, 76.33 and 93.8 percent in No lockdown, 50% lockdown and 75% lockdown scenarios respectively. It was observed that (Figure 2 b) the spread of the infection observed as peak infections in scenarios 2 & 3 were 129779 and 33973 which would be reached in a period of just 33 and 25 days in a no lockdown condition. The peak values 191907, 37790 and 7986 in figure 2 c corresponds to 35%, 6.89% and 1.46% respectively of the initial uninfected population.
a) Number of uninfected people - all age groups b) Number of infected people - all age groups (Cumulative) c) Number of infected people (for a given instant)
The peak values in figure 3 a) corresponds to 31.71%, 6.28% and 1.33% of the initial uninfected population respectively for the three scenarios. There is an occurrence of a maximum equal to 173892 on 44th day, 34414 on 84th day and 7281 on 90th day for the no lockdown, 50% and 75% lockdown scenarios respectively. The peak admissions correspond to 6.42%, 1.33% and 0.28% of the initial healthy population respectively for the three scenarios. The peak number of patients in ICU, 2390, 496 and 94 respectively for the three scenarios, which indicate the minimum number of intensive care setups required for the three scenarios. Similarly, figure 3 d) indicates the minimum number of ventilator setups, 1929, 405 and 78 for the three scenarios respectively.
a) Number of asymptomatic people – all age groups b) Number of admitted people– all age groups c) Number of people in ICU – all age groups d) Number of people using ventilators – all age groups
Number of people deceased - all age groups Maximum deaths as seen in figure 4 correspond to 0.42%, 0.14% and 0.04% of the initial uninfected population respectively for three of the scenarios.
Table 4 indicates the peak values for various states possessed by the agents. The peak values also decrease for all these states as the stringency of lockdown is increased, indicating the effectiveness of lockdown measures. The values (in %) indicate the percentage with respect to the initial healthy population. A significant drop in peak number of ICUs required from 2390 for a no lockdown condition to 94 for a 75% lockdown is evident. Concurrently, the peak number of ventilators decline from 1929 for a no lockdown condition to 78 for a 75% lockdown condition.
Peak values for various health conditions
4. Discussion
This approach based on the synthetic population of 5,43,823 agents in Rangareddy District, for three different NPI scenarios (no lockdown, 50% lockdown and 75% lockdown) projects that the transmission rate of COVID-19 could be effectively brought down by stringency of lockdown measures. The study was performed at the district-level is a major strength of the study as it facilitates decision-making easier to policymakers at specific regions 47. The simulation results are presented using ISPOR-SMDM Modeling Good Research Practices and ethical good practice in modelling.
Synthetic populations are most often generated from open source data such as Open data from London Imperial College, etc.,34, US Census Bureau data,32,36, Australian Census data39, etc. The geographical scope of the study governs the number of agents. For example, 10 million stochastic agents for the State of Delaware, US32, a scaled-down simulation of New York with 10000 agents31, synthetic population of NYC with 500,00036, 5000 agents in the premises of a University in Italy37, 24 million agent representing Australia 39, 750,805 agents representing Urmia, Iran 33, etc. In this study, we have scaled-down the synthetic population of Rangareddy district to 10.35% and used for the ABM.
Contact networks augment variations in behavioral aspects of agents in each model. Classification of a group of people aged greater than 65 and/ or with underlying illness as obesity, chronic cardiac or respiratory illness, and diabetes 34, awareness level (that enhances protection), use of contact tracing mechanisms 35, schedule-based contacts with house members 31–32,36; close contacts 36, closed spaces such as office spaces, university 31,37, indirect contact with suspended viral particles, public gatherings at café, gym, hospitals, transport, 32, touching contaminated surfaces, washing hands 31–35,69; etc., have all been modeled. Present study modeled the contacts of agents based on the contact rates derived based on population densities.
There are also a spectrum of scenarios analyzed by modelers with an aim to determine the ones that outperform others, such as the no lockdown scenario that is included in almost all studies to be used as a base for comparison, control measures such as face masks, physical distancing, shielding of vulnerable population 34, lockdowns for varying durations 34, reducing contacts in external settings whilst maintaining the close contacts constant32, lifting lockdowns based on age-groups 36, and effect of contact tracing of symptomatic individuals 31,36. NPIs such as 50 and 75 percent lockdowns have been used as scenarios in the present study.
In the present study, a decrease in number of contacts at various locations such as in schools, works, etc., was incorporated to enact lockdown scenarios whilst maintaining the contacts made at house 41. It is evident from the results of present study that as the percentage of lockdown imposed was increased, the magnitudes of peak infections reduced with a delay in their corresponding occurrences, which provides more time for the policymakers to increase their capacities to meet the influx of cases. A team of researchers from The Center For Disease Dynamics, Economics & Policy (CDDEP) and Princeton University using ABM estimated the state-wise capacity requirements to accommodate the influx of hospitalizations to help the policymakers to increase their capacities to match the influx based on estimates in India 70,71.
Considering some other parameters such clustering in contact networks, especially in the context of spread of infections would provide more accurate results 72–75. GIS information and migration routes could be included to improve the projections in specific areas 20,21. Wearable devices could be integrated with mobiles to provide real-time monitoring of COVID-19 patients 76. Exploring the contact network and dynamics of different regions would help us to represent the region-specific disease spread better 47
There are certain limitations to the study as parameters such as underlying health conditions, migration routes, adoption of control measures (face mask, social distancing, etc.), longitudinally varying lockdown phases, etc., have not been considered. The parameters which were used in the model were form different countries and may not represent the India or Rangareddy district scenario. The results of simulation model clearly indicate that the peak values could significantly be reduced by increasing the lockdown imposed. Thus, the importance of reducing the number of contacts, i.e., social distancing, is apparent through the results of this study and flattening the disease curve.
5. Conclusions
Majority of the ABM studies focus on specific regions that is a major strength of ABM as it allows defining characteristics at individual level 36. We present this ABM using AnyLogic on a synthetic population in Rangareddy district, Telangana state, India. Further, data specific to India to parametrize such ABM will be critical. Having a synthetic population of a country can provide several options to create ABMs for several disease conditions apart from COVID-19 and may prove efficient for decision-making.
Data Availability
Detailed day-wise data of the number of people in each of the state (age-group wise), the AnyLogic model file, synthetic population (text file) and input template (Spreadsheet) are available in the links provided.
https://cloud.anylogic.com/model/7cd10c0c-f1c1-4b8f-9aac-0bf37a45379a?mode=SETTINGS
Data Availability
Detailed output (age-wise), AnyLogic model file, synthetic population and input spreadsheet are available at: https://cloud.anylogic.com/model/7cdl0c0c-flcl-4b8f-9aac-0bf37a45379a?mode=SETTINGS and https://osf.io/utmhg/files/
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this article.
Funding information
No funds were received for conducting this modeling study.
References
