key: cord-0119066-y11z3zi2 authors: Mourad, Ayman; Mroue, Fatima title: Modeling and Simulation of the spread of coronavirus disease (COVID-19) in Lebanon date: 2020-08-07 journal: nan DOI: nan sha: 712099e693ce3d93aca29862e0cc4b58fc840c1b doc_id: 119066 cord_uid: y11z3zi2 In this paper, we develop a probabilistic mathematical model for the spread of coronavirus disease (COVID-19). It takes into account the known special characteristics of this disease such as the existence of infectious undetected cases and the different social and infectiousness conditions of infected people. In particular, it considers the social structure and governmental measures in a country, the fraction of detected cases over the real total infected cases, and the influx of undetected infected people from outside the borders. Although the model is simple and allows a reasonable identification of its parameters, using the data provided by local authorities on this pandemic, it is also complex enough to capture the most important effects. We study the particular case of Lebanon and use its reported data to estimate the model parameters, which can be of interest for estimating the spread of COVID-19 in other countries. We show a good agreement between the reported data and the estimations given by our model. We also simulate several scenarios that help policy makers in deciding how to loosen different measures without risking a severe wave of COVID-19. We are also able to identify the main factors that lead to specific scenarios which helps in a better understanding of the spread of the virus. After the identification of a novel strain of coronavirus in China in December 2019, later labeled COVID-19, the virus spread in China and most of the countries of the world. The World Health Organization (WHO) declared the outbreak as a public health emergency of international concern [31] by January 2020 and the situation was declared as a pandemic by March 2020 [30] . The global spread of the virus has lead to an unprecedented infectious disease crisis worldwide with 15,666,671 confirmed cases and 636,787 deaths by July 30, 2019. With the uncertainties about the COVID-19 virus, governments must put in place measures to minimize mortality rates as well as the economic impact of the viral spread [3] . In this context, predictive mathematical models for epidemics are fundamental, for policy makers to take the right decisions in order to minimize mortality and morbidity and flatten the epidemic curve for the sake of reducing the strains on the healthcare system, in a timely way taking into account the socioeconomic consequences [3] . Moreover, the effects of mitigation strategies significantly vary depending on population characteristics such as age structure, the stage of the epidemic such as the number of infected and immunized people and the level of population compliance with the measures. Therefore, model-based predictions of the eventual effect of different mitigation measures on the epidemic, taking into consideration such local parameters, are crucial to support evidence-based decisions. Various models have been proposed in the literature for disease spreads. They can be categorized into agent-based models (ABM) [4, 1] and compartmental models [14, 8, 5] . Compartmental models are built on differential equations and assume that the population is perfectly mixed with people moving between compartments such as susceptible (S), infected (I) and recovered (R) [21, 28, 18, 19] . These models revealed the threshold nature of epidemics and triumphed in explaining "herd immunity". However, they fail to capture complex social networks and the behavior of individuals who may adapt depending on disease prevalence. On the other hand, agent-based models can capture irrational behavior and complex social networks and are used to simulate the interactions of autonomous agents that can be either individuals or collective entities [9] . Furthermore, ABM approaches in epidemiology can simulate such complex dynamic systems with less simplification of the rich variation among individuals. Also, the statistical variance is more evident than in compartmental models, whose smooth curves often misleadingly express more certainty than justified, due to the randomization at each run [17, 16, 32, 20] . Several models have been developed for the COVID-19 pandemic. Some are compartmental models such as [25, 2, 6, 12] while others are stochastic transmission models [13, 23] and agent-based models [7, 20] . In the present work, we propose a stochastic individual agentbased model for the spread of COVID-19. The model takes into account the characteristics of the virus and attributes to them probability distributions. We mention, for instance, the incubation period, infectiousness, testing sensitivity, ... It also incorporates probability distributions for social and individual conditions such as: family size within and outside the same household, number of people contacted per day and their subdivision into known versus unknown contacts... These distributions are fitted for the case of Lebanon but they can be modified depending on the country. In Section 2, we describe the parameters used in the model and the rational of using them. Section 3 contains the calibration of the model for the case of Lebanon, the validation of the results obtained by the simulations by comparing to the real data and the predictions of the strength of the epidemic by the end of August and September depending on the mitigation measures imposed by policy makers such as the closure or opening of the airport and the strictness of measures in the country. Finally, we conclude with some remarks in Section 4. In this model, we consider each individual as an agent with variety of variables that are considered relevant to model COVID-19 spreading. In particular, to each infected individual are attributed the following parameters: -the times of infection and detection: it is important to reduce the delay between the infection time and the detection time because during this period, infected individuals remain in the community while possibly infecting others [22] . -the probability of infection (transmitting the virus to others depending whether the person is wearing a mask, obeying hygiene measures, keeping a safe distance, ...): this probability differs also depending on the nature of contact between the infected person and his contacts where we differentiate family members in the same household from contacts outside the house such as at work or in public transportation. The probability of infecting family members sharing the same house is assumed to be higher during the period extending from the instant of infection to the time of detection. -the incubation period (time from exposure to illness onset): the general incubation period for COVID-19 varies for different individuals and ranges from 2 to 14 days. It is 5-6 days on average with the 95 th percentile of the distribution at 12.5 days and the log-normal distribution provides the best fit to the data for incubation period estimates [24, 26] (see Figure 1 ). -the connections: those are either family members or contacts of the individual such as coworkers, neighbors,... To account for family members we consider a log-normal probability distribution function for the average family size or household size (see -contact-tracing: the identification of the source of infection, tracing and isolating its contacts are crucial for breaking the chain of transmission and for the control of the epidemic [13] . Two parameters influence the effectiveness of contact-tracing and isolation: the transmissibility of the pathogen (measured by the basic reproductive number) and the proportion of presymptomatic transmission [11] . Very high levels of contact tracing are required in the case of presymptomatic infectiousness [13] . In the present model, we account for contact tracing efficiency by including a probability distribution for the proportion of newly infected people whom source of infection is unknown. We also consider the probability of a detected infected person not to restrict to isolation measure and keep roaming around irresponsibly. -the localization: in order to account for the effect of airport and border opening, we split individuals into travelers and locals. For travelers, we distinguish between those having a positive screening result and those with a negative result. Those with a positive test result are assumed to be quarantined and are no longer infectious. For the negative result group, we consider the probability of having a false negative test result in addition to the percentage of those who do not go into a period of isolation. On the other hand, the local infected residents are subdivided into two groups: symptomatic with positive test result versus those screened due to their contact with an infected person before onset of symptoms. The probability distribution function for travelers is based on data collected from local health authorities. We note that the ABM approach needs realistic data, typically obtained from a census, to create agents that effectively capture human behavior [10] . . Travelers with negative PCR results upon arrival were isolated for a maximum of 3 days. However, people were not abiding by the preventive measures, and travelers were not respecting the isolation period [15] . The virus rebounded to reach a total number of 4205 by July 29 and the authorities reinstated lock-down from July 30 to August 3 and from August 6 to August 10. In this context, it is obvious that there is an urgent need for evidence-based decisions particularly that the country is facing an unprecedented economic crisis. Calibration. In order to obtain a realistic model for the case of Lebanon, first we used part of the real data to estimate the parameters that are involved in the model. Then we validated the model for the second part of the real data. Indeed, using numerical optimization techniques, we calibrated the model and computed the parameters in order to reproduce results that are close to real data for the period before June 30, published on the website of the ministry of public health of Lebanon (https://coronanews-lb.com/). During this period, strict measures regarding the travelers were imposed in Lebanon, and the number of social events (such as weddings) was still small. On the other hand, as of July 1 st , measures were relaxed and people were less respectful to the measures so the average number of new people met per day was increased, the quarantine period for travelers was restricted to 3 days and the percentage of travelers respecting the quarantine duration decreased. For the incubation period, we used a discrete distribution obtained from [24, 26] as in Figure 1 . On the other hand, due to the lack of census in Lebanon, we used the reports [27, 33] to deduce a distribution for family size in the same household (see Figure 2 ). In the absence of official reports, the distribution of new people met per day was also deduced from our familiarity with the Lebanese context where public transport is not prevalent, the number of family members met per day is high, ... Figure 2 . Distribution for the family size in the same household. Validation. The simulation results that we obtained for the suggested model are very close to the real data for the period up to June 30 due to the parameters estimation. Moreover, these results are close to the real data for the period from July 1 st to July 19 which validates the model. Figure 3 shows the daily new confirmed cases as reported on the website of the ministry of public health in Lebanon. Figure 4 shows the daily new cases from February 21 till July 19 obtained by the simulation with the aforementioned parameters and Figure 5 shows the cumulative number of cases per day for both the real data and the simulations. In order to obtain a representative probabilistic result, the simulations were obtained from 60 runs of the probabilistic model, so we represent the result as the average number of daily new cases within the range of the corresponding standard deviation. Prediction. The prediction of the time course of the epidemic being crucial for policymakers. We simulate in this paragraph four scenarios that differ according to the measures (strict or loose) and to the airport (open or closed) for the periods of the first two weeks and the last two weeks of August. First, we consider that during August, strict measures are adopted by citizens but the airport kept open. This scenario (Scenario I) corresponds to the government policy in response to the bound in the number of cases during the last week of July. Figure 7 shows the curve of the total number of cases until the end of August. This number reaches by the end of August 7,893 cases on average ± 2,193. Figure 6 shows the effectiveness of strict measures in decreasing the number of people infected by the same person as shown in the higher probability of a person not to infect any other person when strict measures are imposed during August. Second, we consider that during August, less stringent measures are adopted by citizens and they are mostly behaving as in the period following the lift of the lock-down. Moreover, we assume that the airport is kept open. Figure 8 shows the result of the simulation of this scenario (Scenario II) where the total number of cases reaches by the end of August 8,793 cases on average ±2,549. We also consider two more scenarios: Scenarios III and IV where strict measures and no strict measures are imposed during all of August respectively, but for both scenarios the airport closes for the last two weeks of August. In Table 1 , we summarize the total number of cases for the four scenarios by the end of August. We proposed, in the present work, a stochastic individual-based model for the spread of COVID-19 in Lebanon that takes into account some characteristics of the virus and the Lebanese population. On one hand, we considered probability distributions for the incubation period of the virus and for the available testing sensitivity. On the other hand, we took into Figure 9 ). Consequently, a plateau is observed in the total number of cases by September (if the airport is closed) compared to a continuous increase in the number of cases during September if the airport is kept open (see Figure 10 ). Despite the lack of census in Lebanon and the absence of public official reports containing real contact-tracing data the suggested model captured effectively the pattern of propagation of the virus in the country and is shedding light on the efficiency of mitigation measures imposed by policy makers. Comparing large-scale computational approaches to epidemic modeling: agent-based versus structured metapopulation models Athanasios Tsakris, and Constantinos Siettos. Data-based analysis, modelling and forecasting of the covid-19 outbreak How will countrybased mitigation measures influence the course of the covid-19 epidemic? The Lancet Agent-based modeling: Methods and techniques for simulating human systems. Proceedings of the national academy of sciences Mathematical models in population biology and epidemiology Can the covid-19 epidemic be managed on the basis of daily data? Modelling transmission and control of the covid-19 pandemic in australia Mathematical tools for understanding infectious disease dynamics Modelling to contain pandemics Controlling pandemic flu: the value of international air travel restrictions Factors that make an infectious disease outbreak controllable Modelling the covid-19 epidemic and implementation of population-wide interventions in italy Feasibility of controlling covid-19 outbreaks by isolation of cases and contacts. The Lancet Global Health The mathematics of infectious diseases Lebanon reinstates lockdown measures after virus rebound An open-data-driven agent-based model to simulate infectious disease outbreaks A taxonomy for agent-based models in human infectious disease epidemiology The basic si model An sirs model with a nonlinear incidence rate Universal masking is urgent in the covid-19 pandemic: Seir and agent based models, empirical validation Containing papers of a mathematical and physical character The effectiveness of contact tracing in emerging epidemics Early dynamics of transmission and control of covid-19: a mathematical modelling study. The lancet infectious diseases Early transmission dynamics in wuhan, china, of novel coronavirusinfected pneumonia A conceptual model for the outbreak of coronavirus disease 2019 (covid-19) in wuhan, china with individual reaction and governmental action Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data The development and state of the art of adult learning and education (ale): National report of lebanon. National Committee for Illiteracy and Adult Education The quasi-stationary distribution of the closed endemic sis model Director-generals opening remarks at the media briefing on covid-19 -11 World Health Organization et al. Statement on the second meeting of the international health regulations detail/30-01-2020-statement-onthe-second-meeting-of-theinternational-health-regulations Agent-based modeling in public health: current applications and future directions. Annual review of public health Population & housing in lebanon