Mathematics in Applied Sciences and Engineering https://ojs.lib.uwo.ca/mase Volume 1, Number 4, December 2020, pp.425-440 https://doi.org/10.5206/mase/10876 DETERMINING THE EFFECTIVENESS OF PRACTICING NON-PHARMACEUTICAL INTERVENTIONS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING AGENT-BASED MODELLING STEVEN KYLE D.C. VILLANUEVA AND CHRISTIAN ALVIN H. BUHAT Abstract. To determine the effectiveness of non-pharmaceutical interventions on an epidemic, we develop an agent-based model that simulates the spread of an infectious disease in a small community and its emerging phenomena. We vary parameters such as initial population, initial infected, infection rate, recovery rate, death rate, and asymptomatic rates, as inputs. Our simulations show that (i) random mass testing and quarantines decreases the number of deaths, infections, and time duration; (ii) social distancing lengthen outbreak period to an extent and helps flatten the epidemic curve; (iii) the most effective combination of NPIs to minimize death, infection, and duration is no mass testing, no social distancing, and a total lockdown but is not ethical; and (iv) the most feasible intervention is to implement an enhanced community quarantine while doing random testing on a maximum of 100 individuals. Results of this study can aid policymakers in determining interventions for their communities during a pandemic. 1. Introduction In late 2019, an unidentified coronavirus had emerged from Wuhan, China, and had initially spread to Southeast Asian countries. In February 2020, the disease was officially named Coronavirus Disease 2019 (COVID-19) and had been acknowledged by the World Health Organization as a global pandemic due to its rapid human to human transfer. Though it has an estimated 2 - 5% mortality rate, the death toll has reached half a million people worldwide since June 2020 [23, 24, 2]. In the Philippines, the first suspected COVID-19 case was confirmed last January 2020 [7, 10]. In early March 2020, the first local transmission was confirmed by the countrys fifth case [11]. Due to the quick surge in the number of suspected and confirmed cases, different interventions such as lockdown procedures were initially implemented in Metro Manila, which eventually expanded to the whole island of Luzon and ultimately the whole nation [21, 5]. The house bill Bayanihan to Heal as One Act was written and signed to provide the President further power to combat the pandemic [22]. However, the number of infected still rose due to factors such as difficulty in testing, thus the quarantine protocol continued [4]. Then, it was alleviated in May to minimal risk areas [9]. Since then, multiple interventions/ variations have been attempted by each LGUs (through the support of research) to try and fight the pandemic. One way experts study patterns of certain phenomena and gain insights in these pandemic situations is through modeling. Some of the advantages in modeling include lesser risk to subjects being infected compared to experimentation and fieldwork, easier to solve and analyze through computer algorithms, and can be done repeatedly with different conditions [15]. In this study, we use agent-based modeling (ABM) to simulate and analyze the underlying patterns and trends of such virus spread. Received by the editors 6 July 2020; accepted 23 December 2020; published online 29 December 2020. Key words and phrases. agent-based model, COVID-19, epidemics, non-pharmaceutical interventions . This study is supported by the University of the Philippines Resilience Institute, University of the Philippines. 425 https://ojs.lib.uwo.ca/mase https://dx.doi.org/https://doi.org/10.5206/mase/10876 426 SKDC. VILLANUEVA AND CAH. BUHAT Here, we first determine the effect of a certain epidemic by varying city information and virus dynam- ics. Then, we check if the different non-pharmaceutical interventions (NPIs) that have been implemented are necessary for controlling the spread of the said virus. After that, we want to determine if simulta- neously practice of such interventions is effective in combating the pandemic. We then determine what levels of preventive measures would minimize the number of deaths and the duration of the epidemic. Lastly, we investigate other parameters that may have an effect on the outbreak. 2. Model Framework We use an agent-based model (ABM) to simulate the transmission of the virus among individuals and to determine the effects of factors and NPIs in the process of infection. 2.1. Agent-Based Model. ABM is a type of microscale model used to simulate interactions among individuals agents and their complex behavior [3]. It aims to find and study emergent trends and phe- nomena that may arise on a given set of autonomous agents that follows specific predefined procedures [1, 16]. This has been utilized by different fields ranging from tourism, marketing, and epidemiology due to its flexibility, simplicity and affordability [16, 14, 20, 12]. Here, we consider a model that simulates the infection among individuals given various factors in an environment. We then determine the effects of various NPI to the number of infections, deaths, etc. The model determines significant developments whenever there are altercations in the population of virus dynamics. It also identifies whenever there is an effect in the the number of cases, deaths if agents practice the NPIs. The agents in this simulation have the following boolean attributes: infected?, con- firmed?, hospitalized?, dead?, and permissible? which is simulated per tick (14-day infectious period). Furthermore, parameters are used in the simulation, and it can be easily revised through the use of sliders, switches, or choosers of the program. They are categorized into two groups namely: city information and virus dynamics. 1 tick represents 1 infectious period for COVID-19 which is 14 days [6], thus, in the simulation environment, all rates are multiplied by 14 signifying the 14 days per tick effect in the simulation. We also incorporate scientific control in the model for the comparisons between different runs. We assume the control to take on the following default values from Table 1. Moreover, global variables are also declared to update the monitors and plots for every tick of the simulation. They are as follows: population, healthy, infected, confirmed, unconfirmed, hospitalized, dead, permissible, and max infected. Note that the NetLogo software has two built-in elements namely patches and ticks, which represent the space and time variables respectively. Since we are only concerned with the underlying behavior and trends of the agents including the repercussions after their interactions, we may not specify the variables in actual space and time measurements. However, if researchers and policymakers deemed it necessary for their analysis, they may assign an appropriate conversion table between the simulation and real-time values. Lastly, these are the procedures that the simulation does to affect all agents behaviors, attributes, and interactions. With all the variables and procedures declared, we use NetLogo programming environment for our ABM. The simulation environment of the program and its development is shown in Figs 1,2,3. The setup button generates the agents into random locations inside the environment. The number of agents placed depends on the initial population. Consequently, a 4 by 4 hospital, which has been colored blue, is placed in the bottom-left corner of the environment. All agents parameters are then varied depending on the initial infected and asymptomatic rate. All those who are infected and confirmed EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 427 Table 1. Parameters and their Descriptions Parameter Description Default value Reference initial population The initial number of agents in the simulation en- vironment 1000 Assumed initial infected The initial number of agents that are infected 10 Assumed testing capacity An integer variable which will be used in the test procedure - Varied quarantine The type of quarantine the city is implement- ing, which have 4 options: none, general, en- hanced. and total community quarantine which allows 100%, 50%, 25% and 0% of the present agents to move respectively. Note that the quar- antine cannot be changed in the middle of a run. - Varied social distancing A boolean variable which will be used in the check sd procedure - Varied infection rate The chance an agent gets infected if another agent is near them during an infectious period 70/14 [13] recovery rate The chance (in %) an agent don’t get infected pro- vided that they are hospitalized during an infec- tious period 20/14 [8] death rate The chance (in %) an agent die provided that they are infected during an infectious period 20/14 [8] asymptomatic rate The chance (in %) an infected agent will not be confirmed 18/14 [19] Figure 1. NetLogo user interface of the agent-based modelling program before a sim- ulation run becomes red in color. Agents who are infected but not yet confirmed turns to orange. The remaining agents stay green. Monitors will also be updated depending on the simulation. 428 SKDC. VILLANUEVA AND CAH. BUHAT Table 2. Procedure Names and their Descriptions Procedure Name Description clear Clears the world of agents, plots, monitors, and notes setup Runs the setup-hospital, setup-turtles, update-monitors, update-display pro- cedures, and reset the number of ticks setup-hospital Sets a small 4 by 4 unit of patches on the lower-left of the world to be colored blue, which will be designated as the hospital setup-turtles Creates a number of turtles in the world based on the initial population in random locations, sets the agents attributes based on the initial infected, and to run the check quarantine procedure check-quarantine Sets all present agents permission to move based on the chosen quarantine type update-monitors Updates all monitors on the program based on the global variables update-display Updates all agents to be hidden if dead sets all agents colors based on their attributes (red if infected and confirmed, orange if infected but not confirmed, and green if not infected) and checks whether they are inside the hospital or not start Runs the check sd, test , move, infect, recover, death, update-monitors, and update-display procedures, does a tick on the tick counter and runs the update-output once the simulation ends check sd Checks the social distancing variable and if true, all agents will move back whenever another agent is near them test Randomly picks agents based on the testing capacity variable and confirms the tested agents if they are infected move Sets all hospitalized agents movement based on their attributes (all infected towards or within the hospital and all not infected towards outside the hospital), and runs the outside movements procedure outside movements Sets all permitted agents movement outside the hospital to walk randomly infect Checks all infected agents and other non-infected agents near them and runs the get sick procedure based on the infection rate variable get sick Sets agent to be infected and also sets agent to be confirmed or not based on the asymptomatic rate variable recover Checks all hospitalized agents and runs the get healthy procedure based on the recovery rate variable death Checks all infected agents and sets them dead based on the death rate vari- able update-output Prints the results of the simulation in the output Once the setup is done, the start button proceeds to run the simulation. The procedures, which were defined in Table 2, for agents movement, testing measures, infections, deaths, and recoveries are done until it hits the stopping criterion. The run stops if (1) an extinction occurs where all the agents are dead or (2) the pandemic had subsided meaning there are no more infected agents. EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 429 Figure 2. NetLogo user interface of the agent-based modelling program during a running simulation Figure 3. NetLogo user interface of the agent-based modelling program after a sim- ulation run 3. Results and Discussion 3.1. Simulation outcome. Given various sets of variables, we determine if a certain epidemic situation either subside and ease up on the agents infections and deaths, or lead to extinction based on the stopping criterion through our model. Here are two runs that show the simulation results based on the same city information but with different virus dynamics, immediately followed by a comparison table of the simulation results. The simulation has three important results after every run: deceased, max infected cases, and ticks taken. Deceased pertains to those agents that have died while the ticks taken refers to the duration of the epidemic regardless if it ended in the survival of the agent set or pure extinction. Moreover, the max infected cases refer to the peak number of infected cases throughout the process. 430 SKDC. VILLANUEVA AND CAH. BUHAT Figure 4. Simulation Run A of the Program (default) Figure 5. Simulation Run B of the Program 3.2. Behavior of Simulation Results according to Number of Runs using Default Values. Consequently, we also observe the minimum number of runs required to observe at which point the simulation results approaches in terms of deceased, infected, and ticks taken. Also, we use the same number of runs to forecast all other remaining experiments. After proceeding on the experiments with a EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 431 Table 3. A Comparative Table of Runs A and B Variables Run A Run B initial population 1000 1000 initial infected 10 10 testing capacity 0 0 quarantine None None social distancing No No infection rate 70% 95% recovery rate 20% 1% death rate 20% 10% asymptomatic rate 18% 50% deceased 155 983 max infected cases 42 275 ticks taken 40 83 varying number of runs, we need to have at least 1000 runs for every experiment made to achieve a less than 1% error on average from 100000 runs. To ensure the stability of the results from the experiments, we consider 5000 runs for all succeeding simulations. Table 4. The Effects of Increasing Testing Capacity on Different Simulation Results testing cap deceased % change % sum max infected % change % sum ticks taken % change % sum 0 262.75 - - 60.32 - - 44.58 - - 10 256.67 -2.32 -2.32 59.90 -0.69 -0.69 43.60 -2.20 -2.20 20 250.37 -2.45 -4.71 59.29 -1.01 -1.70 42.46 -2.61 -4.75 30 246.76 -1.44 -6.09 59.12 -0.30 -1.99 41.83 -1.48 -6.17 40 242.95 -1.54 -7.54 59.35 0.39 -1.61 40.91 -2.21 -8.24 50 239.20 -1.54 -8.96 59.17 -0.29 -1.90 40.25 -1.61 -9.71 60 235.25 -1.65 -10.47 58.75 -0.71 -2.60 39.81 -1.11 -10.71 70 232.54 -1.15 -11.50 58.59 -0.27 -2.86 39.16 -1.61 -12.16 80 229.86 -1.15 -12.52 58.27 -0.55 -3.39 38.77 -1.02 -13.05 90 226.63 -1.40 -13.75 58.17 -0.17 -3.56 38.44 -0.85 -13.79 100 225.86 -0.34 -14.04 58.46 0.50 -3.08 37.93 -1.33 -14.93 3.3. Effects of Increments on Testing Capacity. Note that all deceased, max infected, and ticks taken results have been decreasing whenever we increase the testing capacity of the hospital. In a small community with a population of 1000 where 10 are infected by a virus, by randomly testing a maximum of 100 individuals per tick, we are able to lower the death toll by up to 14%, and the maximum number of infected cases by up to 3%. Also, the duration of the epidemic can be shortened by up to 14% as compared to runs without testing. However, as we increase the testing capacity, its effectiveness significantly reduces as well, as we can see in the percent change columns of Table 4. 3.4. Effects of Quarantine Types. From Table 1, the quarantine types have an irregular increment, which is as follows: 100% of the present agents are permissible to move if there is no quarantine, 50% 432 SKDC. VILLANUEVA AND CAH. BUHAT for general, 25% for enhanced, and 0% for a total quarantine. Thus, the percent differences displayed in Table 5 are all referenced from the control run, which is no community quarantine is in place. Table 5. The Effects of Quarantine Types on Different Simulation Results quarantine deceased % change max infected % change ticks taken % change none (100%) 263.2618 - 60.562 - 44.5142 - general (50%) 179.8694 -31.677 46.589 -23.072 37.172 -16.494 enhanced (25%) 135.3152 -48.601 38.3046 -36.751 32.5722 -26.827 total (0%) 88.0672 -66.548 29.359 -51.522 26.096 -41.376 The same trend can be seen here as well, as the testing capacity. As we implement stricter quarantine policies, we can reduce the death toll, maximum infected cases, and duration. By implementing a total community quarantine, we may be able to reduce virus-related deaths by up to 66%, infected cases by up to 51%, and the situation period by up to 41%. Table 6. The Effects of Social Distancing on Different Simulation Results social distancing deceased max infected ticks taken false 261.0946 60.0242 44.4006 true 379.6116 72.6072 55.1476 3.5. Effects of Social Distancing. Interestingly, the advantages of social distancing have increased the number of deaths, the total maximum of infected cases, and the duration of the epidemic. There is a 45% growth in death numbers, 21% in max infected cases, and 24% in ticks taken. This phenomenon might be a result of the stochasticity of the agent-based model and computational programming errors, such as the conditions set by the check sd procedure and zigzagging occurrence. For the procedure, all agents that are near other agents are checked. Then, an agent X will face the nearby agent Y and take a unit step backward. However, this might result in agent X being near to another agent, say Z. This can lead to a recursive problem due to the spatial capacity of the environment. A zigzagging occurrence also happens when a confirmed infected agent does winding movements on its way to the hospital since it also avoids nearby agents, in accordance with the social distancing procedure. This leads to delays in providing medical assistance to infected agents and further infections due to their unnecessary movements. A study on social distancing has also been proven that it can lengthen the duration of virus outbreaks. This intervention was also said to flatten infection curves [18]. However, by flattening curves, it does not directly imply that the peak, or the maximum infected cases, will be lowered as well. Furthermore, [17] have mentioned that even if social distancing is observed, the outbreak rebounds once it is removed. The researchers suggested that as much as social distancing can delay the spread of the virus, consequent testing of suspected individuals must also take place, as well as developments for improving healthcare systems. 3.6. Effects of Simultaneous NPIs. We now check the different combinations of the non-pharmaceutical interventions and identify which set of the said policies provides us with the minimal average number of fatalities, maximum number of infected cases, and the duration of the epidemic. Using various ar- rangements of such policies might vary its effects as well and might help us optimize on what policies EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 433 Figure 6. The Effects of Simultaneous Non-pharmaceutical Interventions on Deceased we can practice. With all previous simulations that have been done, we simulate 5000 runs on each of the following experiments. For a community, populated with 1000 individuals where 10 of those are infected, no mass testing, no social distancing, and a total lockdown must take place in order to minimize deaths, infections, and overall duration. Note that this combination offers the least mobility of agents. Thus, the outbreak slows down since there are no crowding situation in the hospital and no unnecessary movements. Only those who are infected and confirmed are allowed to be admitted to the hospital. However, this leads to ethical and humanitarian issues since asymptomatic individuals are left off to recover for themselves. Moreover, this is unrealistic since the entire population cannot be quarantined. Looking for an alternative, the next best option, which is the most ethical, is to implement an enhanced community quarantine with no social distancing while doing random testing on a maximum of 100 individuals. By looking at the figures and on sections 3.3 to 3.5, we can observe that the doing stricter quarantine procedures (excluding total) and having greater testing capacity, the deaths, infections, and duration can be minimized. Again, the social distancing protocol may be subject to change provided that the environment space is improved so that there are no zigzagging occurrences. Another way we can combat the virus is by delaying its spread or simply flattening the infection curve. By practicing no quarantine, no mass testing, and an active social distancing policy, we are able to allocate the maximum time for our scientists and researchers to study and formulate either a vaccine or cure. If successful, this completely eradicates the virus over a guaranteed period of time without having any rebound or subsequent waves. 434 SKDC. VILLANUEVA AND CAH. BUHAT Figure 7. The Effects of Simultaneous Non-pharmaceutical Interventions on Max Infected 3.7. Initial Population and Initial Infected Cases. Let us now determine if there will be significant insights that can be found by gradually increasing the initial number of infected cases on the simulation results. As the initial number of infected cases increases, the number of deaths and maximum infected cases increases as well, while the duration of the outbreak slows down. However, the rate of change in all aspects also decreases as the number of initial infected cases increases. Moreover, looking at Fig 9, if the number of initial infected cases is at least 70% of the general population, the population retains the same number of infected cases. This would imply that throughout the pandemic, the population has at least 70% infected individuals. Thus, quick and effective response against the pandemic must be taken into consideration. 3.8. Initial Population and Spatial Capacity of Environment. As for this experiment, we im- prove the program by considering a wider environment and check if there will be changes in the results of the simulation. Note that in both cases, we use the default (control) values. The expanded environ- ment is scaled 2 times larger than the original environment and same as before, we simulate this for 5000 runs each. Given the same population of 1000 people where 10 are infected, Table 7 have shown the drastic effect of the spatial capacity of the environment. Since there is a wider place the agents can move around, the chance of interaction, either inside or outside of the hospital, is lesser compared to the previous environment. This may play a very important role in dense cities and capitals. Governing bodies may opt to space out individuals by slowly decongesting highly populated areas to wider locations EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 435 Figure 8. The Effects of Simultaneous Non-pharmaceutical Interventions on Ticks Taken Figure 9. The Effects of Gradual Increments on Initial Infected Cases on Deceased and Max Infected with less population. Unfortunately, this may have a detrimental effect on the economic, social, and organizational aspects of society. Furthermore, this solution may also contaminate virus-free places. 436 SKDC. VILLANUEVA AND CAH. BUHAT Figure 10. The Effects of Gradual Increments on Initial Infected Cases on Ticks Taken Table 7. A Comparative Table of the Original and Modified Environments Variables and Results Original Program Modified Program % change Environment Size 30 x 30 60 x 60 - Hospital Size 4 x 4 8 x 8 - Average Deceased 259.6444 25.6146 -90.1347 Average Max Infected 60.1062 11.1038 -81.5264 Average Ticks Taken 44.1326 22.144 -49.8239 3.9. Virus Dynamics. We now vary these parameters such as infection rate, death rate, recovery rate, and asymptomatic rates to further analyze the effects of each of them in the transmission of the disease. On Figs 11, 12, and 13, we can see that the simulation results almost follow similar trends with each other. As the infection rate increases, the number of deaths, maximum infection and ticks taken increases as well. This trend also occurs in increments on asymptomatic rate. However, the growth in deceased and max infected results on infection rate seems to follow an exponential curve, while the rest looks linear. The results in varying asymptomatic rates may be due to the fact that it is harder to trace such patients thus lengthening and worsening the outbreak. On the other hand, there is an inverse relationship on recovery and death rate and in the simulation results. While the recovery rate build up, the deaths, infections and duration lessens in a linear manner. As the virus has an increasing death rate, the number of deaths, infections and duration decreases exponentially. The outcome from varying death rate may be caused by the chance of death overriding the chance of infection, thus reducing transmissions between agents. Also, we can consider an epidemic with a 0% death rate as a special case as we can see in Fig 14. After 5 runs of the specific experiment, the program does not reach any of the stopping criteria since this situation will not lead to extinction nor the end of the epidemic. It simply reaches an equilibrium EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 437 Figure 11. The Effects of Variations in Virus Dynamics on Deceased Figure 12. The Effects of Variations in Virus Dynamics on Max Infected with about 85% of the initial population being infected as seen on the first monitor. However, looking closer to the infected cases, the number of asymptomatic patients becomes greater than the number of confirmed ones as time passes. Unfortunately, this case will not be of further significance to the study of the virus since there are already confirmed death cases due to it. Thus, the virus has a nonzero death rate. 438 SKDC. VILLANUEVA AND CAH. BUHAT Figure 13. The Effects of Variations in Virus Dynamics on Ticks Taken Figure 14. A Simulation Run of the Program with Default Values but with 0% Death Rate Note that the governing bodies and institutions can make no particular policies that can help control the following virus rates. Nonetheless, these experiments are very beneficial and valuable to research and development teams as these show how crucial further studies on the virus affects the future of the community. 4. Conclusion The use of agent-based modeling has never been more important in our current pandemic situation. Without the use of heavy experimentation and costly operations, our model was able to determine how an epidemic/pandemic will behave based on its demographics and the dynamics of the infection/virus. Given an assumed control for the simulations, we have considered three vital outcomes of an epidemic, EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 439 namely the death toll, peak infected cases, and outbreak duration. We computed with at least 1000 runs for every experiment made to accurately estimate results from 100000 runs. We then tested non-pharmaceutical interventions individually among 1000 individuals, where 10 of those are infected. Mass testing can lower the death toll by up to 14%, the maximum number of infected cases by up to 3%, and the time extent by up to 14%. By practicing community quarantines, we may be able to reduce virus-related deaths by up to 66%, infected cases by up to 51%, and the outbreak period by up to 41%. Meanwhile, social distancing has lengthened the duration of the outbreak, by up to 24%, to provide ample time for the epidemic curve to flatten and for healthcare facilities to cope up with the number of cases. By combining these NPIs, the most optimal intervention was to implement no mass testing, no social distancing, and a total lockdown to be implemented in a small virus-infested community. This arrangement minimized the mobility of agents to its fullest. This resulted in the lowest deaths, infections, and duration in comparison to other combinations. Although this is unethical and detrimental to other fields such as, but not limited to, the economy. Thus, a more feasible, recommendable, and realistic option will be to implement an enhanced community quarantine with no social distancing while doing random testing on a maximum of 100 individuals. We may also consider another choice that is to have no mass testing, no quarantine but with active social distancing if the goal was to flatten the curve. Note that this may lead to the longest time taken by the outbreak, which can stretch the duration long enough for scientists to formulate either a vaccine and further control the infections. Additionally, the infected individuals must not reach 70% of the general population as this situation reaches equilibrium. This made the virus very hard to exterminate and we might need to resort to drastic measures. Also, a wider spatial capacity of the environment has made a significant development on the simulation results since the chances of interaction are minimized. This lowered deaths by up to 90%, infections by up to 81%, and time took by up to 49%. Lastly, increments on virus dynamics have been analyzed as well. An increase in infection and asymptomatic rates also led to an increase in deaths, infections, and duration. Conversely, as the recovery and death rate increase, the said results decrease. These results from our simulation model play a vital role in how governing bodies can optimally make guidelines such that we can properly utilize our resources, at the same time mitigating the adverse effects of the virus. This can also be a foundation for policymakers and for all other similar studies that dwell on epidemiology and important insights in mitigating virus spread. 5. Limitations The non-pharmaceutical interventions are assumed to be followed and practiced by all of the agents in the simulation environment. Furthermore, there is no delay time considered in such NPIs and no temporal analysis has been made, which may also play an important part in mitigating the outbreak. Moreover, the NPIs that were analyzed are all collective actions, hence, individuals practicing other solitary preventive measures were not considered. Finally, all agents are assumed to have the same protection or immunity rate against the virus. 6. Recommendations Note that assumptions such as movement in the study are random, which should be taken into con- sideration in considering the study. This study can be further extended by determining the minimum testing capacity needed to achieve minimal infections and deaths, in accordance with various parameter testing. Moreover, taking other precautionary measures and physiological factors can be made as proce- dures in the program as well. Among these are wearing of a mask, frequent handwashing, reproduction 440 SKDC. VILLANUEVA AND CAH. BUHAT rates, or limited foreign travels. Lastly, having actual and accurate data from a reliable and trustworthy source to assign on variables will make all the simulation results as close to the real-world situation. 7. Supplementary Files The agent-based model can be accessed by visiting the website: https://github.com/styledcv/EpidemicNPIsSimulations. The Program 1.nlogo is the original program that was used in the simulation testings done on the study. Meanwhile, Program 2.nlogo is the expanded version of Program 1 that was used on section 3.8. References [1] A. Bazghandi, Techniques, advantages and problems of agent based modeling for traffic simulation, International Journal of Computer Science Issues (IJCSI) 9(2012), 115-119. [2] BBC, Coronavirus: Covid-19 death toll hits 500,000 worldwide, https://www.bbc.com/news/world-53212430, ac- cessed 6 July 2020. [3] C. A. H. Buhat, J. F. Rabajante and V. G. V. Paller, , Spatiotemporal modeling of parasite aggregation among fish hosts in a lentic ecosystem, Modeling Earth Systems and Environment, in press (2020), DOI: https://doi.org/10.1007/s40808-020-00983-8. [4] C. A. H. Buhat, J. C. C. Duero, E. F. O. Felix, J. F. Rabajante, and J. B. Mamplata, Optimal Allocation of COVID- 19 Test Kits Among Accredited Testing Centers in the Philippines, Journal of Healthcare Informatics Research, in press (2020), DOI: https://doi.org/10.1007/s41666-020-00081-5. [5] Manacaan Palace, Proclamation No. 922 Declaring a State of Public Health Emergency Throughout The Philippines, https://www.officialgazette.gov.ph/downloads/2020/02feb/20200308-PROC-922-RRD-1.pdf, accessed 6 July 2020. [6] L. Dy and J. F. Rabajante, A COVID-19 infection risk model for frontline health care workers, Network Modeling Analysis in Health Informatics and Bioinformatics 9 (2020): 57. [7] E. M. Edrada, E.B. Lopez, J. B. Villarama, EP. Salva Villarama, B.F. Dagoc, C. Smith, AR. Sayo, JA. Verona, J. Trifalgar-Arches, J. Lazaro, and others, First COVID-19 infections in the Philippines: a case report, Tropical medicine and health 48(2020), 1-7. [8] S. Eikenberry, M. Mancuso, E. Iboi, T. Phan, K. Eikenberry, Y. Kuang, E. Kostelich, and A. Gumel, To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic, Infectious Disease Modelling 5(2020), 293-308. [9] Executive Order No. 112: Imposing an Enhanced Community QUarantine In High-Risk Geographic Areas of The Philippines and a General Community Quarantine in the Rest of the Country from 01 to 15 May, 2020, Adopting the Omnibus Guidelines on the Implementation Thereof, and for Other Purposes, http://www.covid19.gov.ph/wp- content/uploads/2020/05/2020030-EO-112-RRD.pdf, accessed 6 July 2020. [10] Department of Health, DOH Confirms First 2019-Ncov Case in the Country; Assures Public of Intensified Contain- ment Measures, https://www.doh.gov.ph/doh-press-release/doh-confirms-first-2019-nCoV-case-in-the-country, ac- cessed 6 July 2020. [11] Department of Health, DOH Confirms Local Transmission of Covid-19 In PH; Reports 6th Case, https://www.doh.gov.ph/doh-press-release/doh-confirms-local-transmission-of-covid-19-in-ph, accessed 6 July 2020. [12] E. Hunter, B. Mac Namee,and J. Kelleher, An open-data-driven agent-based model to simulate infectious disease outbreaks, Public Library of Science San Francisco, CA USA 13(2018): e0208775. [13] J. Ioannidis, The infection fatality rate of COVID-19 inferred from seroprevalence data, Bulletin of the World Health Organization, accepted (2020), Article ID: BLT.20.265892. [14] E. Keles, Emrah and FHE. Keles, Agent Based Marketing: An Inspiring Review, LCBR European Marketing Con- ference, Munich, Germany, August 2014. DOI: https://doi.org/10.13140/2.1.3429.5365. [15] E. Klipp, R. Herwig, A. Kowald, C. Wierling and H. Lehrach, Systems biology in practice: concepts, implementation and application, John Wiley & Sons (2005). [16] J. Loomis, C. Bond, D. Harpman, The potential of Agent-Based modelling for performing economic analysis of adaptive natural resource management, Journal of Natural Resources Policy Research 1(2008), 35-48. [17] L. Matrajt and T. Leung, Evaluating the Effectiveness of Social Distancing Interventions to Delay or Flatten the Epidemic Curve of Coronavirus Disease, Emerging Infectious Diseases 26(2020), 1740-1748. [18] G. J. Milne, and S. Xie, The effectiveness of social distancing in mitigating COVID-19 spread: a modelling analysis, medRxiv, DOI: https://doi.org/10.1101/2020.03.20.20040055. https://github.com/styledcv/EpidemicNPIsSimulations https://www.bbc.com/news/world-53212430 https://doi.org/10.1007/s40808-020-00983-8 https://doi.org/10.1007/s41666-020-00081-5 https://www.officialgazette.gov.ph/downloads/2020/02feb/20200308-PROC-922-RRD-1.pdf http://www.covid19.gov.ph/wp-content/uploads/2020/05/2020030-EO-112-RRD.pdf http://www.covid19.gov.ph/wp-content/uploads/2020/05/2020030-EO-112-RRD.pdf https://www.doh.gov.ph/doh-press-release/doh-confirms-first-2019-nCoV-case-in-the-country https://www.doh.gov.ph/doh-press-release/doh-confirms-local-transmission-of-covid-19-in-ph https://doi.org/10.13140/2.1.3429.5365 https://doi.org/10.1101/2020.03.20.20040055 EFFECTIVENESS OF PRACTICING NPIS IN IMPROVING VIRUS CONTROL IN A PANDEMIC USING ABM 441 [19] K. Mizumoto, K. Kagaya, A. Zarebski, and G. Chowell, Estimating the asymptomatic proportion of coronavirus dis- ease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020, Eurosurveillance 25(2020): 2000180. [20] S. Nicholls, B. Amelung, and J. Student, Agent-based modeling: A powerful tool for tourism researchers, Journal of Travel Research 56(2017), 3-15. [21] Office of the Presidential Spokeperson, On Code Red Sublevel 2, https://pcoo.gov.ph/OPS-content/on-code-red- sublevel-2/, accessed 6 July 2020 [22] Republic Act No. 11469: Bayanihan to Heal as One Act, https://www.senate.gov.ph/Bayanihan-to-Heal-as-One- Act-RA-11469.pdf, accessed 6 July 2020. [23] M. A. Shereen, S. Khan, A. Kazmi, N. Bashir, Nadia and R. Siddique, COVID-19 infection: Origin, transmission, and characteristics of human coronaviruses, Journal of Advanced Research 24(2020), 91-98. [24] Y. C. Wu, C. S. Chen, and Y. J. Chan, The outbreak of COVID-19: An overview, Journal of the Chinese Medical Association 83(2020), 217-220. Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Laguna 4031, Philippines E-mail address: sdvillanueva3@up.edu.ph Corresponding author, Institute of Mathematical Sciences and Physics, University of the Philippines Los Baños, Laguna 4031, Philippines University of the Philippines Resilience Insitute, University of the Philippines, Philippines E-mail address: chbuhat@up.edu.ph https://pcoo.gov.ph/OPS-content/on-code-red-sublevel-2/ https://pcoo.gov.ph/OPS-content/on-code-red-sublevel-2/ https://www.senate.gov.ph/Bayanihan-to-Heal-as-One-Act-RA-11469.pdf https://www.senate.gov.ph/Bayanihan-to-Heal-as-One-Act-RA-11469.pdf 1. Introduction 2. Model Framework 2.1. Agent-Based Model 3. Results and Discussion 3.1. Simulation outcome 3.2. Behavior of Simulation Results according to Number of Runs using Default Values 3.3. Effects of Increments on Testing Capacity 3.4. Effects of Quarantine Types 3.5. Effects of Social Distancing 3.6. Effects of Simultaneous NPIs 3.7. Initial Population and Initial Infected Cases 3.8. Initial Population and Spatial Capacity of Environment 3.9. Virus Dynamics 4. Conclusion 5. Limitations 6. Recommendations 7. Supplementary Files References