INFEKTA: A General Agent-based Model for Transmission of Infectious Diseases: Studying the COVID-19 Propagation in Bogotá - Colombia ===================================================================================================================================== * Jonatan Gomez * Jeisson Prieto * Elizabeth Leon * Arles Rodríguez ## Abstract The transmission dynamics of the coronavirus - COVID-19-have challenged humankind at almost every level. Currently, research groups around the globe are trying to figure out such transmission dynamics using different scientific and technological approaches. One of those is by using mathematical and computational models like the compartmental model or the agent-based models. In this paper, a general agent-based model, called INFEKTA, that combines the transmission dynamics of an infectious disease with agents (individuals) that can move on a complex network of accessible places defined over a Euclidean space representing a real town or city is proposed. The applicability of INFEKTA is shown by modeling the transmission dynamics of the COVID-19 in Bogotá city, the capital of Colombia. ## Introduction Infectious diseases have a substantial impact on public health, health care, macroeconomics, and society. The availability of options to control and prevent the emergence, expansion or resurgence of pathogens warrants continuous evaluation using different methods. Mathematical models allow characterizing both the behavior and the emergent properties of biological systems, such as the transmission of infectious disease. [1–3]. Many biological systems have been modeled in terms of complexity since their collective behavior cannot be simply inferred from the understanding of their components [4, 5]. Complex systems are computational approaches that make use of computer-based algorithms to model dynamic interactions between individuals agents (e.g. persons, cells) or groups and their properties, within, and across levels of influence [6, 7]. In general, agent-based modeling (ABM) can be used for testing theories about underlying interaction mechanics among the system’s components and their resulting dynamics. It can be done by relaxing assumptions and/or altering the interaction mechanisms at the individual agent level. ABMs can increase our understanding of the mechanisms of complex dynamic systems, and the results of the simulations may be used for predicting future scenarios [8]. In this paper, we introduce an ABM, called INFEKTA (Esperanto word for infectious), for modeling the transmission of infectious disease, applied to the coronavirus COVID-19. INFEKTA models the disease transition at the person level and takes into consideration individual infection disease incubation periods and evolution, medical preconditions, age, daily routines (movements from house to destination places and back, including transportation medium if required), and enforced social separation policies. ### Complex Systems Approaches for Epidemic Models The complex system model approach considers a system as a large number of entities (equally complex systems that have autonomous strategies and behaviors) that interact with each other in local and non-trivial ways [9–11]. This approach provides a conceptual structure (a multi-level complex network [12]) that allows characterizing the interrelation and interaction between elements of a system and between the system and its environment [13]. In this way, a system is composed of sub-systems of second order, which in turn may be composed of subsystems of the third-order [14]. Transmission dynamics of infectious diseases are not traditionally modeled at the individual level, but at the population-level with a compartmental model. However, some recent research use agent-based modeling for doing that [15]. ### Compartmental Model A compartmental model tracks changes in compartments without specifying which individuals are involved [16] and typically reflects health states relevant for transmission (e.g., susceptible, exposed, infectious and recovered). Basically, these kinds of models represent epidemics of the communicable disease using a population-based, non-spatial approach. The conceptual framework for this approach is rooted in the general population model which divides a population into different population compartments [17]. Compartmentalization typically reflects health states relevant for transmission (e.g., susceptible, exposed, infectious and recovered, in short **SEIR**), though more partitioning is possible according to age and/or other relevant host characteristics. Heterogeneous and temporal behavior is modeled through the incorporation of relevant time-dependent social mixing, community structures and seasonality, relevant for infectious disease dynamics [18, 19]. Process dynamics are captured in transition rates, representing the rate by which an average individual transitions between compartments. ### Agent-based models for infectious disease Agent-based models (ABMs) are a type of computer simulation for the creation, disappearance, and movement of a finite collection of interacting individuals or agents with unique attributes regarding spatial location, physiological traits and/or social behavior [20–22]. ABMs work bottom-up, with population-level behavior emerging from the interactions between autonomous individuals and their environment [22, 23]. They allow the history of every individual to be tracked and network structures to be explicitly represented. In general, ABMs allow: i) To introduce local interaction rules at the individual level, which closely coincide with physical and social interaction rules; (ii) To include behaviors that may be randomized at the observational level, but can be deterministic from a mathematical point of view; (iii) To incorporate a modular structure and to add information through new types of individuals or by modifying current rules; and (iv) To observe systems dynamic that could not be inferred from the examination of the rules of particular individuals [8]. Because of the complex nature of interactions among individuals, an ABM may not provide an understanding of the underlying mechanisms of the complex system. However, changes in the local rules repertoire of the individuals can be introduced into the ABM to observe changes in the complex system [24]. When ABM is used for epidemic modeling, infectious disease transmission dynamics is expected to emerge from the interaction between local interactions between the individuals. In this way, each virtual individual includes its health state relevant for the transmission of the infectious disease (severity and time in it), a kind of individual SEIR (or more complex) model. Health state transmission rates are usually approximated from the rates obtained by a compartmental model but are used at the individual level, i.e., when individuals interact with each other. ### INFEKTA Agent-Based Model Our agent-based model of infectious disease propagation, called INFEKTA, consists of five-layer components: **i) Space:** virtual space where individuals move and interact; **ii) Time:** virtual time used by individuals for moving and interacting; **iii) Individuals:** virtual persons being simulated, **iv) Infectious disease dynamics:** Specific disease parameters that modeled how an individual can get infected and how the infectious disease will evolve; and **v) Social separation policy:** Set of rules defined for restricting both access to places and mobility of individuals. #### Space The virtual space (for a city or town being studied) is a Euclidean complex network [12]: Nodes are places (located in some position of the 2D Euclidean space) where individuals can be at some simulation time and edges are routes (straight lines) connecting two neighbor places. #### Place (Node) A place may be of three kinds: home (where individuals live), public transportation station (PTS), and interest place (IP) i.e., school, workplace, market, and transportation terminal. IPs and PTSs are defined in terms of capacity (maximum number of individuals that can be at some simulation step time). IPs and PTSs may be restricted, during some period, to some or all individuals. Place restriction is established according to the social separation rule that is enforced during such a period. #### Neighbor (Edge) A PTS is a neighbor to another according to the public transportation system of the city or town being studied. Homes and IPs are considered neighbors to its closest1 PTS. No home is neighbor to any other home neither an IP is neighbor of any other IP. Finally, a home and an IP are considered neighbors if they are neighbors of the same PTS. ### Time Virtual time is defined in INFEKTA at two resolution levels: days for modeling the transmission dynamics of the infectious disease, and hours for modeling the moving and interaction of individuals. Therefore, if an individual gets infected more than once during the same day, INFEKTA considers all of them as a single infection event. Any individual movement is carried on the same one hour, it was started, regardless of the traveled Euclidean distance neither the length of the path (number of edges in the complex network). ### Individuals A virtual individual in INFEKTA is defined in terms of his/her demographic, mobility and infectious disease state information. #### Demographics The demographic information of a virtual individual consists of: **i) Age** of the individual; **ii) Gender** of the individual *female* o *male*; **iii) Location** of the individual at the current time step; **iv) Home** of the individual, **v) Impact level of medical preconditions** on the infectious disease state if the individual is infected, and **vi) IP interest** of going to certain type of IPs. #### Mobility The ability of an individual to move through space (we use the graph defining the space for determining the route as proposed in [25, 26]). Each individual has a **mobility plan for every day**, plan that is carried on according to the enforced social separation policy and her/his infectious disease state. The mobility plan is modeled in INFEKTA as a collection of simple movement plans to have **i) Policy**: social separation policy required for carrying on the mobility plan; **ii) Type**: may be *mandatory*, i.e., must go to the defined interest place) or *optional*, i.e., any place according to individual’s preferences; **iii) Day**: day of the week the plan is carried on, maybe *every, week, weekend, Monday*, …, *Sunday* ; **iv) Going Hour**; **v) Duration** in hours for coming back to home; and **vi) Place**: if plan type is mandatory, it is a specific place, otherwise it is an IP selected by the individual according to his/her IP preferences. ### Infectious Disease Dynamic Any individual can potentially be in one of seven different infectious disease states or health states: Immune (*M*), Susceptible (*S*), Exposed (*E*), Seriously-Infected (*I**S*), Critically-Infected (*I**C*), Recovered (*R*), and Dead (*D*). As can be noticed, we just adapt the terminology from the compartmental models in epidemiology – namely, from the SEIR (Susceptible-Exposed-Infectious-Recovered) model. In INFEKTA, the infectious state of the SEIR model is divided into exposed, seriously-infected and critically infected in order to capture how age, gender, IP preferences, medical preconditions (co-morbidity), and social separation policies can impact the evolution of the infectious disease in an individual. INFEKTA introduces both the *M* state since some individuals are naturally immune to or can become immune to (after recovering) to certain infectious diseases and the Dead (*D*) state to distinguish between recovered and dead individuals. Figure 1 shows the transition dynamics of the infectious disease at the individual level in INFEKTA. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F1) Fig 1. Transmission dynamics of the infectious disease at individual level in INFEKTA. The infectious disease states are, Susceptible (*S*), Exposed (*E*), Seriously-Infected (*I**S*), Critically-Infected (*I**C*), Recovered (*R*), Immune (*M*) and Dead (*D*). Probabilities are individual based and can be defined in terms of current location, age, gender, and so on. Symbol ![Graphic][1], on an state *X*, indicates that an individual must stay some period of time 𝒯*X* at such state *X* before being able to change to other state. Since rates are defined at the individual level, these rates can be defined by taking into consideration, for example, rates at the population level (obtained from a compartmental model), age, gender and co-morbidity presented in the individual. Remember that those rates are not defined at some time scale (as in compartmental models) but define the rule determining changes in health states of individuals being close enough for interacting at the infectious disease transmission level or after some period of time being in some state. In this way, an individual can change with probability *α* from *S* state to state *E* if close enough to an infected (*E, I**S*, and *I**C*) individual and will change from state *I**S* to state *I**C* if has been at state *I**S* with probability *θ* if has been on state *I**S* a period of time ![Graphic][2]. * *α*: is the transmission rate and incorporates the encounter rate between susceptible and infectious individuals together with the probability of transmission. * *β*: is the rate at which individuals move from the exposed (*E*) to the Seriously-Infected state (*I**S*). It’s reciprocal (1*/β*) is the average latent (exposed) period. * *θ*: is the rate at which individuals move from the Seriously-Infected (*E*) to the Critically-Infected state (*I**S*). Its reciprocal (1*/θ*) is the average critical period. * *ϕ*: is the death rate. * *ω*: is the immune rate that incorporates the probability of becoming immune. * 𝒯*E*: Time an individual will be at the Exposed state (*E*) before changing to the Recovered (*R*) or Seriously Infected (*I**S*) states. * ![Graphic][3] Time an individual will be at the Seriously Infected state (*I**S*) before changing to the Recovered (*R*) or Critically Infected (*I**C*) states. * ![Graphic][4] Time an individual will be at the Critically Infected (*I**C*) before changing to the Recovered (*R*) or Dead (*D*) states. * 𝒯*R*: Time an individual will be at the Recovered state (*R*) before changing to the Susceptible (*S*) or Immune (*M*) states. INFEKTA can consider that two individuals were close enough for interacting at the transmission of the infectious disease if they were at the same place (home, interest place or public transportation station) at the same time or if they were close enough (in the Euclidean space) while moving. In order to simplify this checking process, it is possible to consider that an individual just visited its home, final interest place, and both the initial and final PTSs when using the public transportation system. ### Social Separation Policy The social separation policy is described in INFEKTA as a finite sequence of rules, each rule having **i) Start Time**: an initial day for applying the social separation policy rule; **ii) End Time**: final day for ending the social separation policy rule; **iii) Level**: indicates the kind of restriction applied to the mobility of persons and accesses to places, and **iv) Enforce**: defines the specific mobility and access restrictions of the social separation policy. ### Modeling the Transmission Dynamics of the COVID-19 in Bogotá - Colombia INFEKTA is used for modeling the Transmission dynamic of the COVID-19 in Bogotá city, the largest and crowded city in Colombia. Bogotá is the capital city of Colombia, its urban perimeter population is 7.412.566, is composed by 20 *districts*, and its massive public transportation system is called *Transmilenio* (**TM**). TM is a bus-based system, which has 143 stations and moves near to 2.500.000 citizens every day. #### Virtual Space Setup Geographical information of Bogotá is used as the Euclidean space where the moving and interaction complex network is defined. Each one of the TM stations is located and added to the complex network according to the real TM system [27]. Also, the airport and the regional bus terminal are located and connected to the nearest TM station. Demographic information from nineteen districts of Bogotá (only urban districts) is used for generating in the Euclidean space interest places (Workplaces (W), markets (M), and schools (S)), homes (H), and people(P). Places are generated, in each one of the districts, following a 2D multivariate normal distribution *N* ∼ (*µ*, Σ) (*µ* is the geographic center of the district and Σ is the co-variance matrix defined by the points determining the perimeter of the district). The number of places in each district was generated based on the population density of each district according to the projections to 2015 [28]. Table 1 shows the amount of data generated for each type of place and for people, also the number of TM stations (Bus), and terminal transportation that we use in the simulation, and Table 2 shows a detailed information of the number of interest places generated by district. View this table: [Table 1.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/T1) Table 1. Data used in the simulation. (*) real places. View this table: [Table 2.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/T2) Table 2. Number of interest places generated by districts Figure 2 shows the virtual Euclidean map of Bogotá with the georeferenced real and simulated places; also, the figure shows the associated complex network of connected places (nodes are places and edges are routes between places), the graph was drawn with Gephi [29]. ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F2) Fig 2. Bogotá with georeferenced places (left) and Bogotá complex network (right). #### Individuals setup An heterogeneous (varying gender, age, district and home) group of one thousand one individuals (1001) is generated using a stratified sampling based on the demographic information of the city for each district according to the projections to 2015 [28]. An individual is classified, according to her/his age, as: *Child* = [0-12], *Adolescence* = [13-18], *Adult* = [19-49], *Senior* = [50-69], and *older* = 70+. Table 3 shows the total demographic information of virtual people. View this table: [Table 3.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/T3) Table 3. Demographic information of virtual people. District (D), Total(T), Male(M), Female(F). Also, a sequence of activities was assigned randomly to each individual to define a diary routine. This was done according to the person’s age and the hour of the day. For example, some agents *Adolescence* go to school, and some agents *adult* go to work. Some agents may move using the PTI system and some others while going directly to its destination place. The route an individual takes is defined according to the complex network. Figure 3 shows three examples of different routines (paths over the graph) for the individuals. ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F3) Fig 3. Example routines (paths) carried by individuals [Age, Gender, IP, District]. Individual 5 (left): [Child, F, School, 01]; Individual 55 (middle): [Senior, F, Workplace, 01]; Individual 999 (right): [Older, M, Market, 19]. The explicit impact level of age and medical preconditions on the state of the COVID-19 dynamic is not included in this preliminary modeling. We wrapped them in the transition rates and allow modelers to change and play with different rates. Therefore, we set the initial values of these rates as shown in Table 4. View this table: [Table 4.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/T4) Table 4. Parameters of INFEKTA and their estimations for COVID-19. #### Social Separation Rule Setting The level attribute of the social separation rule for the COVID-19 in the virtual Bogotá city is defined as follows: * None: No restrictions to the mobility neither to access to places. * Soft: Few places are restricted (depending on the type, capacity, etc). i.e., randomly close *U* ∼ (*size*(*places*) ∗ 0.3, *size*(*places*) ∗ 0.5). * Medium: Many places and few stations are restricted (depending on the type, capacity, etc). Some type of individuals is restricted to stay at home (except those with the required mobility level). i.e., randomly close *U* ∼ (*size*(*places*) ∗ 0.5, *size*(*places*) ∗ 0.7). * Extreme: Few places are accessible to persons while few stations are restricted. Almost every individual is restricted to stay at home (except those with the required mobility level). i.e., randomly close *U* ∼ (*size*(*places*) ∗ 0.7, *size*(*places*) ∗ 0.9). * Total: Scare places are accessible to persons, and almost all stations are not accessible. Individuals are restricted to stay at home (except those with the required mobility level). i.e., randomly close *U* ∼ (*size*(*places*) ∗ 0.9, *size*(*places*)). ## Results The INFEKTA implementation of the COVID-19 propagation model for Bogotá city is implemented using the simulation modeling software tool *AnyLogic*, see S1 Simulator. Also, the methodology (Data preprocessing; places, population, and routes assignation; network creation) is available in a Github repository, see S2 Repository. Figure 4 shows the initial configuration of virtual individuals in the simulation, represented as small points. All of them are initially assigned to homes, and almost all of them are in the Susceptible state (*S*) (green points). Just ten (10) individuals are considered at state Exposed (*E*). ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F4) Fig 4. Initial configuration of INFEKTA. The simulator presents the dynamic of the Individuals (agents) in Bogotá, the parameters of INFEKTA, the social separat**P**i**a**o**us**n**ed** rules that can be applied, and the Epidemic dynamics (agents state per hour). Figure 5 shows the state of the COVID-19 dynamics (propagation of the virus) at three different simulation time steps. Figure 5 (left) shows such dynamics after some simulated minutes of starting the simulation. Notice that, individuals are moving (mainly using Transmilenio) according to their assigned routines. Figure 5 (middle) shows the COVID-19 dynamics after 400 simulated hours. At this time, all individuals have been Exposed (yellow points), meaning that some individuals are seriously infected (red points), critically infected (dark red points), recovered (blue points), and few of them become immune (purple points). Figure 5 (right) shows the COVID-19 dynamics after 800 simulated hours. Notice that, the virus has already peaked once and the first wave is over. The propagation is now stable and 20% of the population is now immune to the virus. As expected, the number of exposed, seriously and critically infected cases grows exponentially when a social separation policy is not enforced. Clearly, the transmission dynamics of the infectious disease (in this case COVID19) emerge from the local interaction of the individuals. ![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F5) Fig 5. Evolution of the epidemic dynamics of INFEKTA. Initial configuration (left), 400 hours (middle), and 800 hours (right). We analyze the sensitivity of some parameters to check the robustness of the model. Sensitivity to infection rate is shown in Figure 6. Notice that by increasing or decreasing the infection disease rate (Figure 6 (left)), the peak of the transmission dynamics is reached sooner or later on time. When low infection diseases rates, the number of cases is also low, reducing the impact on the economy but a new peak is likely to appear after that. On the other hand, for high infection disease rates (Figure 6 (right)) the peak is reached in an early stage, and around half of the population is on one of the infected states (Exposed, Seriously-Infected, Critically-Infected). ![Fig 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F6.medium.gif) [Fig 6.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F6) Fig 6. Sensitive analysis for the infection probability *α. α* = 0.05 (left), *α* = 0.1 (middle), *α* = 0.2 (right). Sensitive to recovery time (𝒯*R*) is shown in Figure 7. If the recovery time is reduced to 2 days (Figure 7 (left)) there is not a significant difference in the epidemic behavior respect to the predefined 5 days. However, if recovery time is increased to 10 days (Figure 7 (right)) the transmission disease dynamics show a second peak after 500 hours (approx 21 days). ![Fig 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F7.medium.gif) [Fig 7.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F7) Fig 7. Sensitive analysis for the Recovery Time 𝒯*R*. 𝒯*R* = 2 days (left), 𝒯*R* = 5 days (middle), 𝒯*R* = 10 days (right). Also, we try different scenarios where each one of the social separation policies is enforced just after 800 simulation hours, see Figure 8. As can be noticed, its evident how social separation rules help to mitigate the exponential growth on the transmission disease dynamics (COVID19), reducing the number of infectious cases (Exposed, Seriously-Infected, and Critically-Infected). Interestingly, when the extreme social separation rule (access to approximately 75% of interest places is restricted) the transmission disease dynamic displays similar behavior to the one when a total social separation rule (access to approximately 95% of interest places is restricted) is enforced. ![Fig 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/22/2020.04.06.20056119/F8.medium.gif) [Fig 8.](http://medrxiv.org/content/early/2020/04/22/2020.04.06.20056119/F8) Fig 8. Social Separation after 800 hours and initialized in hour 0. None (top-left), Soft (top-middle), Medium (top-right), Extreme (bottom-left), and Total (bottom-right). ## Conclusions and Future Work Modeling the Transmission dynamic of an infectious disease such as the COVID-19 is not an easy task due to its highly complex nature. When using an agent-based model, several different characteristics can be modeled, for example, the demographic information of the population being studied, the set of places and the mobility of agents in the city or town under consideration, social separation rules that may be enforced, and the special characteristics of the infectious disease being modeled. INFEKTA is a general agent-based model that allows researches to combine and study all of those characteristics. Our preliminary results modeling the transmission dynamics of the coronavirus COVID-19 in Bogotá city, the largest and crowded city in Colombia, indicate that INFEKTA may be a valuable asset for researchers and public health decision-makers for projecting future scenarios when applying different social separation policy rules and controlling the expansion of an infectious disease. Although we are doing a rough and no so real approximation of the transmission dynamics of the COVID19, we are able to obtain similar behaviors, in our preliminary experiments, to those reported for the COVID19 in the real world. Therefore, INFEKTA may be able to provide more accurate results if its parameters are set to real ones: disease transmission rates, virus, incubation periods, comorbidity, houses, interest places, routines, population size (close to nine millions of virtual individuals for the Bogotá city). Our future work will concentrate on studying the transmission of COVID-19 in Bogotá city by considering different scenarios of social separation rules and by using more realistic information about: i) Relation between personal information and propagation rates of the COVID-19, ii) Places and routes, iii) Population size, and iv) Age and Medical preconditions. ## Data Availability No medical data from external sources is used in this paper. Demographic data of the Bogotá city is available on Colombian government sites. Transmilenio data is available at [http://www.transmilenio.gov.co](http://www.transmilenio.gov.co) [https://cloud.anylogic.com/model/9806d370-f0a8-48d1-b3e1-3537721b39ba](https://cloud.anylogic.com/model/9806d370-f0a8-48d1-b3e1-3537721b39ba) [https://github.com/japrietov/INFEKTA](https://github.com/japrietov/INFEKTA) ## Supporting information **S1 Simulator. INFEKTA Anylogic Simulator**. A hands-on version where anybody can try to change some of the parameters is available at INFEKTA anylogic. **S2 Repository. INFEKTA repository**. A repository containing the source code of the simulator and a technical report explaining the modeling methodology is available at INFEKTA github. ## Footnotes * 1 in the 2D Euclidean space. * Received April 6, 2020. * Revision received April 20, 2020. * Accepted April 22, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Lessler J, Edmunds W, Halloran M, Hollingsworth T, Lloyd A. Seven challenges for model-driven data collection in experimental and observational studies. Epidemics. 2015;10:78–82. 2. 2.1. Shapshak P, 2. Balaji S, 3. Kangueane P, 4. Chiappelli F, 5. Somboonwit C, 6. Menezes LJ, et al., editors Escobar-Ospina ME, Gómez J. In: Shapshak P, Balaji S, Kangueane P, Chiappelli F, Somboonwit C, Menezes LJ, et al., editors. Artificial Life and Therapeutic Vaccines Against Cancers that Originate in Viruses. Cham: Springer International Publishing; 2019. p. 149–305. 3. 3.An G, Mi Q, Dutta-Moscato J, Vodovotz Y. Agent-based models in translational systems biology. WIREs Systems Biology and Medicine. 2009;1(2):159–171. 4. 4.Meyers RA. Encyclopedia of Complexity and Systems Science. New York, NY: Springer; 2009. 5. 5.1. Pagel M Mitchell M, Newman M. In: Pagel M, editor. Complex Systems Theory and Evolution. Oxford University Press; 2005. 6. 6.Galea S, Riddle M, Kaplan GA. Causal thinking and complex system approaches in epidemiology. International journal of epidemiology. 2010;39(1):97–106. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyp296&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19820105&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274491000019&link_type=ISI) 7. 7.Ladyman J, Lambert J, Wiesner K. What is a complex system? European Journal for Philosophy of Science. 2013;3(1):33–67. 8. 8.Komosinski M, Adamatzky A. Artificial Life Models in Software. ex2nd ed. Springer Publishing Company, Incorporated; 2014. 9. 9.Luke DA, Stamatakis KA. Systems science methods in public health: dynamics, networks, and agents. Annual review of public health. 2012;33:357–376. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-publhealth-031210-101222&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22224885&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304202700022&link_type=ISI) 10. 10.Merelli E, Rucco M, Sloot P, Tesei L. Topological Characterization of Complex Systems: Using Persistent Entropy. Entropy. 2015;17(10):6872–6892. 11. 11.Sayama H. Introduction to the modeling and analysis of complex systems. Binghamton University, SUNY; 2015. 12. 12.Van Der Hofstad R. Random Graphs and Complex Networks Vol. I.. vol. I. Cambridge Series in Statistical and Probabilistic Mathematics; 2017. 13. 13.Balaji PG, Srinivasan D. An introduction to multi-agent systems. Studies in Computational Intelligence. 2010;310:1–27. 14. 14.Alcocer-Cuarón C, Rivera AL, Castan∼o VM. Hierarchical structure of biological systems. Bioengineered. 2014;5(2):73–79. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4161/bioe.26570&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24145961&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) 15. 15.Willem L, Verelst F, Bilcke J, Hens N, Beutels P. Lessons from a decade of individual-based models for infectious disease transmission: a systematic review (2006-2015). BMC Infectious Diseases. 2017;17:612–. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-017-2699-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) 16. 16.Garnett GP, Cousens S, Hallett TB, Steketee R, Walker N. Mathematical models in the evaluation of health programmes. The Lancet. 2011;378(9790):515–525. 17. 17.Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proceedings of the royal society of london Series A, Containing papers of a mathematical and physical character. 1927;115(772):700–721. 18. 18.Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS medicine. 2008;5(3). 19. 19.Salathé M, Jones JH. Dynamics and control of diseases in networks with community structure. PLoS computational biology. 2010;6(4). 20. 20.Jit M, Brisson M. Modelling the epidemiology of infectious diseases for decision analysis. Pharmacoeconomics. 2011;29(5):371–386. 21. 21.van Kleef E, Robotham JV, Jit M, Deeny SR, Edmunds WJ. Modelling the transmission of healthcare associated infections: a systematic review. BMC infectious diseases. 2013;13(1):294. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2334-13-294&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23809195&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) 22. 22.Willem L, Verelst F, Bilcke J, Hens N, Beutels P. Lessons from a decade of individual-based models for infectious disease transmission: a systematic review (2006-2015). BMC infectious diseases. 2017;17(1):612. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-017-2699-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F22%2F2020.04.06.20056119.atom) 23. 23.Russell S, Norvig P. Artificial intelligence a modern approach. 3rd ed. New Jersey: Prentice-Hall; 2010. 24. 24.Helbing De. Social Self-Organization: Agent-Based Simulations and Experiments to Study Emergent Social Behavior. Springer-Verlag Berlin Heidelberg; 2012. 25. 25.Rodriguez A, Gomez J, Diaconescu A. Towards Failure-Resistant Mobile Distributed Systems Inspired by Swarm Intelligence and Trophallaxis.In: Proceedings of the European Conference on Artificial Life 2015, At The University of York UK; 2015. p. 448–455. 26. 26.1. Pichardo-Lagunas O, 2. Miranda-Jiménez S Rodríguez A, Gómez J, Diaconescu A. Exploring Complex Networks with Failure-Prone Agents. In: Pichardo-Lagunas O, Miranda-Jiménez S, editors. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). vol. 10062 LNAI. Cham: Lecture Notes In Computer Science; 2017. p. 81–98. 27. 27.www transmilenio gov co Subgerencia Técnica y Servicios Bogotá D C. Trazados Troncales de TRANSMILENIO; 2019. Available from: [https://datosabiertos-transmilenio.hub.arcgis.com/datasets/trazados-troncales-de-transmilenio/data](https://datosabiertos-transmilenio.hub.arcgis.com/datasets/trazados-troncales-de-transmilenio/data). 28. 28.Secretaría Distrital de Planeación. Información, cartografía y estadística; 2016. Available from: [http://www.sdp.gov.co/gestion-estudios-estrategicos/](http://www.sdp.gov.co/gestion-estudios-estrategicos/) informacion-cartografia-y-estadistica/consulta. 29. 29.Bastian M, Heymann S, Jacomy M. Gephi: An Open Source Software for Exploring and Manipulating Networks; 2009. [1]: F1/embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif