Multi-objective Risk-based Resource Allocation for Urban Pandemic Preparedness: The COVID-19 Case in Bogotá, Colombia ====================================================================================================================== * Jeisson Prieto * Jonatan Gomez ## ABSTRACT Determining how best to allocate resources to be used during a pandemic is a strategic decision that directly affects the success of pandemic response operations. However, government agencies have finite resources, so they can’t monitor everything all of the time: they have to decide how best to allocate their scarce resources (i.e., budget for antivirals and preventive vaccinations, Intensive Care Unit (ICU), ventilators, non-intensive Care Unit (non-ICU), doctors) across a broad range of risk exposures (i.e., geographic spread, routes of transmission, overall poverty, medical preconditions). This paper establishes a comprehensive risk-based emergency management framework that could be used by decision-makers to determine how best to manage medical resources, as well as suggest patient allocation among hospitals and alternative healthcare facilities. A set of risk indexes are proposed by modeling the randomness and uncertainty of allocating resources in a pandemic. The city understudy is modeled as a Euclidean complex network, where depending on the neighborhood influence of allocating a resource in a demand point (i.e., informing citizens, limit social contact, allocate a new hospital) different network configurations are proposed. Finally, a multi-objective risk-based resource allocation (MoRRA) framework is proposed to optimize the allocation of resources in pandemics. The applicability of the framework is shown by the identification of high-risk areas where to prioritize the resource allocation during the current COVID-19 pandemic in Bogotá, Colombia. Keywords * Risk-based Resource Allocation * Risk assessment * Multi-objective Optimization * Pandemic control * COVID-19 * Urban Spatial analysis ## 1. Introduction Pandemics are large-scale outbreaks of infectious diseases that cause significant social, political, and economical disruption (PAHO, 2009; Koyuncu and Erol, 2010; Madhav et al., 2017). Policy attention has focused on the need to identify emerging outbreaks that might lead Pandemics, and to expand investment to build preparedness and health capacity (Lederberg et al., 2003). In the preparedness pandemic, effective allocation of limited health resources (i.e., budget for antivirals and preventive vaccinations, Intensive Care Unit (ICU), ventilators, non-intensive Care Unit (non-ICU), doctors) plays a critical role in order to reduce the number of cases, hospitalization, and deaths. Despite the great advances in prevention and treatment of infectious diseases, the world is unaware to respond to Pandemic or any similarly global public-health emergency (Who, 2011; Madhav et al., 2017). In most countries, health care systems operate at or above maximally designed capacity. Many hospitals just do not have sufficient pre-existing resources to respond to surge capacity in an outbreak (Biddison et al., 2019). Unlike with natural disasters, where the greatest need for resources often occurs early in the time course, pandemic resource requirements will build over months. Outbreaks that become pandemics generally do not take hold in multiple locations at exactly the same time, they are geographically and temporally patchy (Koyuncu and Erol, 2010). Many government agencies and health planners are responsible for overseeing and monitoring future outbreaks. However, they can’t monitor everything all of the time, they have to decide how best to allocate their scarce resources across a broad range of risk exposures. This is called “risk-based resource allocation.” (Farrell et al., 2013). Different types of government agencies face risk-based resource allocation decisions: agricultural land and water resources (Romero, 2000; Wolgin, 1975; Li et al., 2016), system design in a Distributed Environment (Yeddanapudi et al., 2008; Qiu et al., 2008), terrorism (Willis, 2007; Quadrifoglio, 2008; Ray et al., 2009), and Natural Hazards (i.e., tornados, hurricanes, earthquakes) (Vaziri et al., 2010; Zolfaghari and Peyghaleh, 2015; Murphy and Gardoni, 2007; Vaziri, 2008). In the Risk-based resource allocation for the pandemic response, a demand point has one (or more) associated risk (i.e., geographic spread, routes of transmission, risk factors for infection, overall poverty, medical preconditions) and the objective is to choose the amount to be invested in several interventions which minimize the overall risk exposed by the demand points according to budget constraints and health benefits. Due to the randomness and uncertainty of conditions, not only one but a set of risks may adversely affect the allocation of resources in the geographical space. Then, the objectives (one objective for each risk that a demand point may be exposed) must be optimized simultaneously (Yan and Haimes, 2011; Sun et al., 2014), but there exists a trade-off among objectives, i.e., an improvement gained for one objective is only achieved by making concessions to another objective. This paper aims to describe and illustrate a Multi-objective Risk-based Resource Allocation framework (MoRRA) that could be used by decision-makers to determine how best to manage medical resources, as well as suggest patient allocation among hospitals and alternative healthcare facilities. In MoRRA, different definitions of risk-based resource allocation are given depending on the geographical space and its neighborhood configuration. This study was carried out during the COVID-19 in Bogotá, Colombia to identify geographic areas with high-risk factors in where to prioritize surveillance to control the outbreak and to generate recommendations for future outbreaks. The remainder of this paper is organized as follows. In the next section, the background knowledge, the risk-based resource allocation (RRA) problem, and the formulation of the multi-objective RRA (MoRRA) are given. In Section 3, the experimental setup for MoRRA is described. Afterward, numerical results of a case study (the current COVID-19 Pandemic in Bogotá, Colombia) are presented in section 4 where it is demonstrated how MoRRA could help decisionmakers to determine the resource allocation and potential resource shortages in the healthcare system. Finally, conclusions and potential future developments are discussed in the last section. ## 2. Risk-based Resource Allocation in a Multiobjective Framework Following the risk-based resource allocation methodology proposed in (Farrell et al., 2013), the proposed framework involves three main stages. The first stage is the identification and definition of the risk. The second stage is the estimation of the level of risk posed in a demand point. Once the risk has been defined and measured, an optimal strategy is proposed to minimize the risk exposure. ### 2.1. Risk definition Although there are different definitions of risk, we use the one given by (Kelman, 2018). Risk is composed by two components, *hazard*, and *vulnerability*. Definition 1. **Hazard** is the probability that a disaster (i.e., COVID-19) occurs. Definition 2. **Vulnerability** is the possibility that damages (i.e., fatalities, injuries, property damage, or other consequences) occur at a demand point because a resource is not allocated. Risk is then defined as the expected damages due to a particular hazard for a given area and reference period. Based on mathematical calculations, the risk of the demand point *i* can be determined as a product of hazard (*H*) and vulnerability (*V*) (DHA, 1992). ![Formula][1] ### 2.2. Risk measuring #### 2.2.1. Hazard assessment The hazard assessment describes the identification of what hazards can be expected and how they might change in the short and medium-term as a result of environmental phenomena or processes (Kelman, 2018). First of all, all of the potential hazards are identified. Then the areas that could be affected by the hazard are marked, this is called Hazard Mapping. The magnitude, intensity, and frequency of the hazards are determined and the causes of the hazards are investigated. Hazards could include earthquakes, volcanic eruptions, floods, drought, cyclones, and epidemics. #### 2.2.2. Vulnerability assessment Vulnerability Assessment describes the degree to which socioeconomic systems and physical assets in geographic areas are either susceptible or resilient to the impact of a disaster (i.e., pandemic). Several models have been proposed to establish vulnerable urban areas over the infectious disease domain, i.e., vector-borne diseases (Hagenlocher et al., 2014), Dengue (de Mattos Almeida et al., 2007), malaria (Kienberger and Hagenlocher, 2014; Hagenlocher and Castro, 2015), Ebola (Moore et al., 2017), and COVID-19 (Mishra et al., 2020; Prieto et al., 2021). ### 2.3. Risk strategy #### 2.3.1. Resource Allocation Problem The resource allocation problem seeks to find an optimal allocation of a fixed amount of resources to activities to minimize the cost incurred by the allocation. Given a finite set of resources ℛ = {(*r**1*, *r**2*, … *r**a*) | *r*i ∈ ℝ+} whose total amount is equal to *T*, it is required to allocate it to *a* activities so that the objective value *f* (ℛ) is minimized, see equation 2. The objective value may be interpreted as the cost or loss, or the profit or reward, incurred by the resulting allocation (Katoh and Ibaraki, 1998). ![Formula][2] where *r**i* represents the amount of resource allocated to activity *i* and *f**i*(*r**i*) is the cost incurred by allocate the resource *r**i* at the *i*-th activity. If the resource is divisible, *r**i* is a continuous variable that can take any non-negative value. If it represents persons, processors, or trucks, on the other hand, variable *r**i* becomes a discrete variable that takes non-negative integer values (*discrete resource allocation problem*). #### 2.3.2. Urban space Let 𝒮 the geographical space under study (i.e. state, country, or city) defined in terms of a finite set of *P* smaller spatial units (i.e. countries, census tracts, or zip codes); that is 𝒮 = {1, 2, …, *P*}; is an Euclidean complex network *G* = (𝒮, *E*) (Van Der Hofstad, 2017), where the spatial units (or vertex) 𝒮 are located in some position of the 2D Euclidean space and edges *E* are connection between two nodes given the spatial relation *meet*. Definition 3. The spatial relation *meet*(*i, j*) occurs when *i* has at least one point in common with *j* in the boundary, but their interiors do not intersect (Egenhofer, 1990). Definition 4. An **adjacent vertex** *u* of a vertex *v* in a graph *G* is a vertex that is connected to *v* by an edge (i.e., *meet*(*u, v*)). Definition 5. The **neighborhood** of a vertex *v* (*N*(*v*)) in a graph *G* is the subgraph of *G* induced by all vertices adjacent to *v*. Definition 6. Regardless the metric space under consideration (points, spatial units, binary strings, DNA strands) (Prieto et al., 2019), we will call 𝒩 the **neighborhood class** of all neighborhoods in a graph *G*, i.e., 𝒩 = {*N*(*v*)|*v* ∈ 𝒮}. #### 2.3.3. Risk-based Resource Allocation Problem (RRA) Let Λ be a risk values associated for each spatial unit in 𝒮; that is Λ = {(*λ*1, *λ*2 … *λ*p)| *λ**i* ∈ ℝ+}, the **Risk-based Resource Allocation Problem**, looks for the optimal way to allocate the resources *R* to each demand point (spatial unit) *i* such that the overall risk over 𝒮 is minimized. Here, the cost incurred *f*Λ,*i*(*r**i*) by allocate the resource *r**i* at the *i*-th activity depends on the neighborhood influence of allocate a resource in *i* (i.e., informing citizens, limit social contact, allocate a new hospital), see equation 3. ![Formula][3] where *r**i* is the impact factor to allocate a resource to *i* demand point (0 ≤ *r**i* *≤*1), *N*(*i*) is the neighborhood of *i*-th demand point, and *a**i,k* is the influence factor in *k* when a resource is allocated in *i* (0 ≤ *a**i,k* *≤* 1). Then, the objective function *f*Λ is calculated among spatial units 𝒮. ![Formula][4] ![Formula][5] Here, depending on the network configuration, three configurations in RRA are proposed (Fig. 1). ![Fig. 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F1.medium.gif) [Fig. 1:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F1) Fig. 1: Different configurations of the Risk-based Resource Allocation problem (RRA). Here, the cost incurred *f*Λ,*i*(*r**i*) by allocate the resource *r**i* at the *i*-th activity depends on the neighborhood influence. Without neighbor influence (RRA-I) (left), with neighbor influence at same scale (RRA-II) (middle), and with neighbor influence at different scale (RRA-III) (right). Definition 7. The **RRA-I** configuration happens when there are not neighborhood influence (∀*i* ∈ 𝒮, *N*(*i*) = Ø). So, the cost incurred *f*Λ,*i*(*r**i*) is defined as. ![Formula][6] Definition 8. The **RRA-II** configuration happens when there are neighborhood influence at same scale (*α* *i,k* = *a*). So, the cost incurred *f*Λ,*i*(*r**i*) is defined as. ![Formula][7] Definition 9. The **RRA-III** configuration happens when there is neighborhood influence at a different scale. So, the cost incurred *f*Λ,*i*(*r**i*) is defined as. ![Formula][8] Using the adjacency matrix *A* of *G*, where the *α**i,k* are the the weight of the edge *w*(*i, k*), the objective function can be evaluated in terms of *A*, see Appendix A. ![Formula][9] ![Formula][10] where *χ* is the Indicator function that determines when a value of *A* is different to 0, *I* is the identity matrix, *diag* is the function that diagonalizes the vector of resources *R*, and Λ are the risks associated with the demand points. #### 2.3.4. Multi-objective RRA (MoRRA) The **multi-objective optimization problem (MoP)** can be mathematically defined as follows. ![Formula][11] where *x* = (*x*1, *x*2, …, *x**n*)*T* is the *n*-dimensional decision variable vector from the decision space Ω; *f* : Ω → Θ ℝ*m* consists a set of the *m* objective functions that map *x* from *n*-dimensional decision space Ω to *m*-dimensional objective space Θ. Definition 10. Given two decision vectors *x, y* ϵ Ω, *x* is said to Pareto **dominate** *y*, denoted by *x ≺ y*, iff *f**i*(*x*) ≤ *f**i*(*y*), for every *i* ∈ {1, 2, …, *m*}, and *f**j*(*x*) *< f**j*(*y*), for at least one index *j* ∈ {1, 2, …, *m*}. Definition 11. A decision vector *x** ∈ Ω is **Pareto optimal** iff there is no *x* ∈ Ω such that *x ≺ x**. Definition 12. The **Pareto set (PS)** is defined as ![Formula][12] Definition 13. The **Pareto Front (PF)** is defined as ![Formula][13] Since objectives in (11) conflicted with each other, no point in Ω simultaneously minimizes all the objectives. The best trade-offs among the objectives can be defined in terms of PF. Due to the randomness and uncertainty of conditions (environmental, operational), RRA also brings many risks that may adversely affect the allocation of resources in the geographical space. Therefore, it is necessary to introduce a comprehensive set of risk indexes by modeling the randomness and uncertainty of the RRA problems. Then, **Multi-objective Risk-based Resource Allocation** aims to optimal way to allocate *R* to each demand point (spatial unit) *i*, in a set of *M* risk indexes; that is ![Graphic][14]. ![Formula][15] ## 3. Resource allocation for COVID-19 in Bogotá, Colombia ### 3.1. Study area and data sources The applicability of the problem is shown for the curent COVID-19 Pandemic in Bogotá city, the largest and crowded city in Colombia. Bogotá is a metropolitan city with 7.412.566 inhabitants living in an area of 1775km (995km urban and 718km rural), at an altitude of 2640m, with an annual temperature ranging from 6 to 20°C, and annual precipitation of over 840mm. Bogotá is composed of 621 Urban Sectors (Urban Sector is a cartographic the division created by the National Administrative Department of Statistics (DANE) (DANE, 2018a).) Fig. 2 shows the distribution of the Urban sectors over Bogotá. ![Fig. 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F2.medium.gif) [Fig. 2:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F2) Fig. 2: Spatial distribution of Bogotá, Colombia using Urban sectors. Information was obtained from the National Department of Statistics (DANE), District Planning Secretary of Bogotá (SDP), and the District Mobility Secretary of Bogotá (SDM). Data comprised public information about demographic, transportation, socio-economic, and health conditions reported from 2011 to 2020. A summary of the datasets is presented as follows: * **MON_2017 (SDP, 2011, 2017):** Dataset provided by SDP containing a set monograph which provides a physical, demographic and socioeconomic description of Bogotá and its districts. * **SDM_2017 (SDM, 2018):** Dataset provided by SDM presenting detailed official information of mobility characterization in Bogotá. * **CNPV_2018 (DANE, 2018a)**: Dataset provided by DANE containing the national census made in 2018 which provides socio-demographic statistics of Colombia. * **DANE_2018 (DANE, 2018b):** Dataset provided by DANE containing the results of the Multidimensional Poverty Index which encompasses educational and health quality, work and housing conditions, and access to public services. * **DANE_2020 (DANE, 2020):** Dataset provided by DANE presenting a vulnerability index based on demographic and health conditions relevant for COVID-19 pandemic. The complex network for the different RRA configurations is built (Fig. 3). For RRA-II, the *α* value is fixed in 0.5, which means that influence in the neighborhood is half when a resource is allocated. For RRA-III, the *α**ik* values are fixed depending on the distance between the spatial units; that is less distance more influence. The *α**ik* are normalized over the range 0 (less influence) and 0.5 (more influence). ![Fig. 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F3.medium.gif) [Fig. 3:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F3) Fig. 3: Complex network representation of the three RRA configurations. RRA-I(left), RRA-II (middle), and RRA-III (right). ### 3.2. Risk definition in COVID-19 #### 3.2.1. Pandemic hazard Natural disasters (i.e., the COVID-19 pandemic) rarely xist, because disasters are social, arising from a combination of hazard and vulnerability, with vulnerability as the causative factor. The disaster occurs at multiple levels simultaneously, with responses to a hazard exposing as many vulnerability problems as the original hazard. The failure to heed to the pandemic plans alongside the lack of health-care accessible to everyone meant that the hazard could not be addressed effectively and vulnerability fundamentals were revealed (Kelman, 2020). Then, based on mathematical calculations, we assume hat the hazard (the new coronavirus) is constant for all spatial units; that is ∀*i* ∈ 𝒮, *H*(*i*) = 1. #### 3.2.2. Pandemic vulnerability Three domains are proposed on (Prieto et al., 2021) to deibe the vulnerability for the COVID-19 in Bogotá, Colombia. Those domains are: (i) Where and how he/she lives (life), (ii) Where and how he/she works (work), and (iii) Where and how he/she moves around (movement). Table 1 shows the domains proposed, the vulnerability factors associated with them, and the dataset used to calculate the values for each factor. View this table: [Table 1](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/T1) Table 1 Vulnerability domains for the COVID-19 case in Bogotá, Colombia (Prieto et al., 2021). * **Where and how he/she lives (life)**: Several demographic factors influence the degree of risk of the geographic areas in a pandemic. The literature emphasizes factors such as *urban density, age*, urban living (i.e., *socio-spatial segregation*), literacy and healthcare quality (i.e., *poverty index*). Further, most data on the COVID-19 pandemic suggest that people with underlying *comorbid* conditions are more vulnerable than people without them. * **Where and how he/she works (work)**: Urban sectors with high-density facilities (i.e, *educational buildings, cultural buildings, sport buildings, food markets*, all *formal labor, informal labor*) are more vulnerable to the spread of contagious diseases due to space limitations within and between households, growth and mobility, and limited water, sanitation, and hygiene (WASH) infrastructure. * **Where and how he/she moves around (movement)**: Understanding the *public transportation dependency, transmission routes*, and infection vulnerability factors (i.e., *geographic impact*) provides the baseline for epidemiological modeling that can inform the planning of response and containment efforts to reduce the likelihood of the disease spreading. Then, based on mathematical calculations, a spatial unit *i* has associated three vulnerable factors: life, work, and movement; that is ∀*i* ∈ 𝒮, *V* (*i*) = {*V**life*(*i*), *V**work*(*i*), *V**movement*(*i*)}. ### 3.3. Risk assessment in COVID-19 Taking the hazard as constant *H*(*i*) = 1, we have only to worry about the Vulnerability Assessment. Previous works (Prieto et al., 2021) proposed a framework for Urban Vulnerability Assessment (UVA) that condense a set of vulnerability factors (Table 1) into a vulnerability index that allowed us to establish and rank potentially vulnerable urban areas in Bogotá. To build the vulnerability index for the three vulnerable domains (life, work, and movement), the following steps are applied: 1. The raw data for each factor is normalized across all spatial units over the range 0 (best) to 1 (worst). The normalization is already calculated in (Prieto et al.,2021). 2. Synthesize the normalized information of all spatial units into *k* partitions which groups spatial units with similar vulnerability profiles. We use *k* = 10 to generate 10 ranges of vulnerabilities (from 0.5 to 0.95, with step of 0.1). 3. The clusters’ centroids of each group are used to sort the vulnerability factors in descending order. This sort is interpreted as vulnerability ranking which is used for the analysis. 4. Then, to aggregate the *L* ranks (one for each vulnerable factor, then for life *L* = 5, work *L* = 6, movement *L* = 3) in a unique vulnerability ranking the Borda’s method is used. 5. The unique vulnerability ranking is then transformed into a vulnerability index, where a higher rank indicates higher vulnerability. The output of this process is the three vulnerability index (one for each domain). Finally, to quantify the risk we follow the risk definition presented in 1. ![Formula][16] ![Formula][17] ![Formula][18] where *Life*(*i*), *W ork*(*i*) and *Movement*(*i*) are the life risk, work risk and movement risk, respectively; for the spatial unit *i*. Fig. 4 shows the final three risk indexes. ![Fig. 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F4.medium.gif) [Fig. 4:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F4) Fig. 4: Risk indexes generated for the current COVID-19 pandemic in Bogota, Colombia. ### 3.4. Risk strategy in COVID-19 #### 3.4.1. Experimental Setup To solve the formulated multi-objective risk-based resource allocation problem, a comparison with different multiobjective algorithms (MOEA/D (Zhang and Li, 2007), NSGA-III (Deb and Jain, 2013), RVEA (Cheng et al., 2016), and ARMOEA (Tian et al., 2017)) was made on different configurations of MoRRA (i.e., RRA-I, RRA-II, RRA-III) and different allocation percentages (i.e., *T* = 10%, *T* = 25%, *T* = 50%). Here, the allocation percentages (amount of resource) *T* is the total of spatial units where a resource should be allocated (i.e., *T* = 10% means that only 10% of the total space units will receive the resource). Also, the impact factor *r**i* is assumed to 0.5, which means that allocates a resource in the spatial unit *i* would reduce the risk in half. ##### Compared Algorithms The following four state-of-the-art algorithms for multi-objective functions are considered as peer algorithms. * MOEA/D (Zhang and Li, 2007): It is representative of the decomposition-based method, the basic idea of MOEA/D is to decompose a MOP into several single-objective optimization subproblems through aggregation functions and simultaneously optimizes them. * NSGA-III (Deb and Jain, 2013): It is based on decomposition with Pareto-adaptive weight vectors. This approach automatically adjusts the weight vectors by the geometrical characteristics of the Pareto front. * RVEA (Cheng et al., 2016): It is a scalarization approach, and termed angle penalized distance approach that dynamically adjusts the distribution of the reference vectors to balance convergence and diversity of the solutions in the PFs. * ARMOEA (Tian et al., 2017): It uses an enhanced inverted generational distance indicator, in which an adaptation method adjusts a set of reference points based on the indicator contributions of candidate solutions. ##### Performance Metrics To evaluate the performance of different MOEAs for the RRA problems, each algorithm was run for the same number of generations, and the resulting solutions (known as Pareto front approximations), are compared using functions that measure two qualities: (i) solution accuracy, i.e., to determine how similar the solution is to the true Pareto front and (ii) solution diversity, i.e., to evaluate how well distributed are the points in the solution. We selected two of the most used metrics called Δ*p* (Schutze et al., 2012) and Coverage over Pareto Front (*C P F*) (Tian et al., 2019) to compare the accuracy and the diversity of the solutions found by the different algorithms. ##### Parameter Settings In this subsection, we first present the general parameter settings for the experiments, and afterward, the specific parameter settings for each algorithm in comparison are given. 1. **MoRRA parameters:** We made the experiment with different configurations of RRA (RRA-I, RRA-II, RRA-III) and allocation percentages (*T* = 10%, *T* = 25%, *T* = 50%). The decision variables are equal to the number of urban sectors in Bogotá (*D* = 631) and the objective functions are fixed to 3 (*M* = 3) representing the different risks. 2. **Settings for the operators**: We select the simulated binary crossover (SBX) (Deb et al., 1995) and the polynomial mutation (PM) (Deb and Goyal, 1996) as genetic operators for our experiments. For the SBX, the distribution index is set to *η**c* = 30 and the crossover probability *p**c* = 1.0 is used in all algorithms; for PM the distribution index and the mutation probability are set to *η**m*= 20 and *p**m* = 1/n, respectively. 3. **Population size**: For all algorithms, the population size is determined by the simplex-lattice design factor *H* together with the objective number *M* (Das and Dennis, 1998). Then the population size using this approach is set to 105 individuals. 4. **Termination Condition**: Every algorithm stops when the number of function evaluations reaches the maximum number. For all configurations and allocation percentages the maximal number of generations is set to 10000. 5. **Specific Parameter Settings in Each Algorithm**: For MOEA/D, the weights vectors are calculate using the penalty-based boundary intersection (PBI), the neighborhood size *T* is set to 20, and the penalty parameter *θ* in PBI is set to 5, as recommended in (Zhang and Li, 2007). For RVEA, the parameter controlling the rate of change of penalty (*α*) and the frequency of employing reference vector adaptation (*fr*) are fixed in 2 and 0.1, respectively, as recommend in (Cheng et al., 2016). #### 3.4.2. Pareto front The statistical results of the Δ*p* and *C P F* metrics values obtained by the four algorithms for the different configurations and allocation percentages over 20 independent runs are summarized in Table 2, where the best results are highlighted. The Wilcoxon rank-sum test is adopted to compare the results obtained by the four compared algorithms at a significance level of 0.05 (here, the MOEA/D algorithm is taken as the reference’s algorithm). Symbol ‘+’ indicates that MOEA/D is significantly outperformed by the compared algorithm according to a Wilcoxon rank-sum test, while ‘−’ means that the compared algorithm is significantly outperformed by MOEA/D. Finally, ‘≈’ means that there is no statistically significant difference between the results obtained by MOEA/D and the compared algorithm. It can be seen that MOEA/D shows the best overall performance among the four compared algorithms over the Δ*p* metric in the experiments, while RVEA shows the best overall performance over the *C P F* metric in experiments. The results obtained by RVEA (good performance on Δ*p* and *C P F*) in the different configurations and allocation percentages will be used in the rest of the paper for the following results. View this table: [Table 2](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/T2) Table 2 Statistics Δ*p* and *C P F* metric values of the Pareto-optimal solutions founded by the four compared algorithms for the different RRA configurations and amount allocations. The numbers in parentheses are the standard deviations. The range of the *non-dominated* solutions found with VEA are shown in Fig. 5. The *Pareto front* behavior shows promising convergence performance as well as a good distribution on problems with different configurations and allocation percentages. ![Fig. 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F5.medium.gif) [Fig. 5:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F5) Fig. 5: The approximate Pareto optimal solutions obtained by RVEA on problems with different configurations (RRA-I, RRA-II, RRA-III) and allocation percentages (10%, 25%, 50%). #### 3.4.3. Decision making To visualize the solution in the Bogotá network map, a pseudo-weight vector approach proposed in (Deb, 2011) is used. This method calculates the normalized distance to the best solution regarding each objective. First, we select one solution with equal pseudo-weights *life* = 0.33, *w**work* = 0.33, *w**movement* = 0.33) for each different RRA configurations and an allocation percentages equal to %10 (Fig 6). The results shows an interesting scenario where the spatial correlation between urban sectors is not remarkable getting an unbiased risk-based resource allocation for COVID-19. ![Fig. 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F6.medium.gif) [Fig. 6:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F6) Fig. 6: Solution of MoRRA visualized in the Bogotá complex network with 10% of allocation percentage and using equal pseudo-weights (*w**life* = 0.33, *w**work* = 0.33, *w**movement* = 0.33). Further, in order to support the decision-maker to find the most preferred solution with a good balance between robustness and the nominal quality, different combinations of pseudo-weights are applied in solutions with RRA-III configuration and 10% of allocation percentage (Fig. 7). In the first combination, we give more weight to life risk (*w**life* = 0.8, *w**work* = 0.1, *w**movement* = 0.1). In the second combination, we give more weight to work risk (*w**life* = 0.1, *w**work* = 0.8, *w**movement* = 0.1). Finally, in the third combination, we give more weight to movement risk (*w**life* = 0.1, *w**work* = 0.1, *w**movement* = 0.8). The results indicate that the MoRRA framework proposed could be used to recommend actions for before, during, and after a pandemic that is, to planning and coordination efforts through leadership and coordination across sectors, to assess if the risk of a pandemic could increase in specific geographic areas. ![Fig. 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/03/01/2021.02.24.21252407/F7.medium.gif) [Fig. 7:](http://medrxiv.org/content/early/2021/03/01/2021.02.24.21252407/F7) Fig. 7: Solution of MoRRA visualized in the Bogotá complex network with different configurations of RRA and 10% of allocation percentage and using different combinations of pseudo-weights. ## 4. Conclusions and Future Work A Multi-objective Risk-based Resource Allocation (MoRRA) amework for Pandemic Preparedness is proposed. MoRRA could be used to build evidence for planning, modeling, and epidemiological studies to better inform the public, policymakers, and international organizations and funders as to where and how to improve surveillance, response efforts, and delivery of resources. The proposed MoRRA is tested in the current COVID-19 Pandemic in Bogotá city, the largest and crowded city in Colombia. MoRRA creates not only one, but a set of risk indices (i.e., life, work, movement) and uses them to apply the risk-based resource allocation. Although the risk factors involved in the framework are uctural, the proposed approach is flexible, does not require xpert support or knowledge, and allows policy-makers, and international organizations to prioritize resource allocation in short and long-term actions for affected populations in a city. For instance, using the solutions with more weight in the movement risk a short-actions (i.e, staying home, limit close contact, avoid crowds, limit non-essential travel) can be taken to reduce the risk in the city. Further, using the solutions with more weight in the life risk, it is possible to advance in long-term territorial reorganization plans (i.e., reduce socio-spatial segregation, decent housing, bio-secure protocols for high-density facilities) as our results indicate for the COVID-19 in the urban area of Bogotá. ### A. Matrix Notation of RRA To derive (9) for (4), the Euclidean Complex Network *G* is represented with its adjacency matrix *A*. Then, we want to demonstrate. ![Formula][19] when *f*Λ(*A, R*) is written as. ![Formula][20] Here, *χ* is the Indicator function that determines when a value of *A* is different to 0, *I* is the identity matrix, *diag* is the function that diagonalizes the vector of resources *R*, and A are the risks associated with the demand points. Expanding *f*Λ(*A, R*), we have. ![Formula][21] Solving. ![Formula][22] Grouping similar terms. ![Formula][23] Each term in the column vector can be written as (3). ![Formula][24] Applying 1-Norm. ![Formula][25] So, it is sufficient to prove that (4) could be written as (9) when the graph *G* is representing by its adjacency matrix *A*. ## Data Availability No medical data from external sources is used in this paper. Demographic data of Bogotá city is available on Colombian government sites. * Received February 24, 2021. * Revision received February 24, 2021. * Accepted March 1, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. Alfonso R Ó.A., 2016. Densidades divergentes y segregación socio-espacial en la bogotá de hoy: un análisis a partir de un índice de cali-dad de la densidad, in: VIII Seminario Internacional de Investigación en Urbanismo, Barcelona-Balneário Camboriú, Junio 2016, Departament d’Urbanisme i Ordenació del Territori. Universitat Politècnica. 2. Biddison, E.L.D., Faden, R., Gwon, H.S., Mareiniss, D.P., Regenberg, A.C., Schoch-Spana, M., Schwartz, J., Toner, E.S., 2019. Too many patients… a framework to guide statewide allocation of scarce mechan-ical ventilation during disasters. Chest 155, 848–854. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chest.2018.09.025&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F01%2F2021.02.24.21252407.atom) 3. Cheng, R., Jin, Y., Olhofer, M., Sendhoff, B., 2016. A reference vector D Estadística: guided evolutionary algorithm for many-objective optimization. IEEE Transactions on Evolutionary Computation 20, 773–791. 4. DANE, 2018a. Departamento Administrativo Nacional de COLOMBIA - Censo Nacional de Población y Vivienda - CNPV - 2018. URL: [http://microdatos.dane.gov.co/index.php/catalog/643](http://microdatos.dane.gov.co/index.php/catalog/643). 5. DANE, 2018b Departamento Administrativo Nacional de Es-Pobreza multidimensional en Colombia. URL: [https://www.dane.gov.co/files/investigaciones/condiciones\_vida/pobreza/2018/bt\_pobreza\_multidimensional\_18.pdf](https://www.dane.gov.co/files/investigaciones/condiciones\_vida/pobreza/2018/bt_pobreza_multidimensional_18.pdf). 6. DANE, 2020. Departamento Administrativo Nacional de Estadística: Índice de vulnerabilidad por manzana con el uso de variables de-mográficas y comorbilidades. URL: [https://www.dane.gov.co/files/comunicados/Nota\_metodologica\_indice\_de\_vulnerabilidad.pdf](https://www.dane.gov.co/files/comunicados/Nota\_metodologica_indice_de_vulnerabilidad.pdf). 7. Das, I., Dennis, J.E., 1998. Normal-boundary intersection: A new method for generating the pareto surface in nonlinear multicriteria optimization problems. SIAM journal on optimization 8, 631–657. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/S1052623496307510&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000074047900001&link_type=ISI) 8. Deb, K., 2011. Multi-objective optimisation using evolutionary algorithms: an introduction, in: Multi-objective evolutionary optimisation for prod-uct design and manufacturing. Springer, pp. 3–34. 9. Deb, K., Agrawal, R.B., et al., 1995. Simulated binary crossover for continuous search space. Complex systems 9, 115–148. 10. Deb, K., Goyal, M., 1996. A combined genetic adaptive search (geneas) for engineering design. Computer Science and informatics 26, 30–45. 11. Deb, K., Jain, H., 2013. An evolutionary many-objective optimization algo-rithm using reference-point-based nondominated sorting approach, part i: solving problems with box constraints. IEEE transactions on evolu-tionary computation 18, 577–601. 12. Dha, U., 1992. Internationally agreed glossary of basic terms related to disaster management. UN DHA (United Nations Department of Human-itarian Affairs), Geneva. 13. Egenhofer, M., 1990. A mathematical framework for the definition of topo-logical relations, in: Proc. the fourth international symposium on spatial data handing, pp. 803–813. 14. Farrell, D., Gebre, B., Hudspeth, C., Sellgren, A., 2013. Risk-based re-source allocation. McKinsey Center for Government. 15. Goldstein, E., Lipsitch, M., 2020. Temporal rise in the proportion of younger adults and older adolescents among coronavirus disease (covid-19) cases following the introduction of physical distancing measures, germany, march to april 2020. Eurosurveillance 25, 2000596. 16. Gomez, J., Prieto, J., Leon, E., Rodríguez, A., 2021. Infekta—an agent-based model for transmission of infectious diseases: The covid-19 case in bogotá, colombia. PloS one 16, e0245787. 17. Hagenlocher, M., Castro, M.C., 2015. Mapping malaria risk and vulnerabil-ity in the united republic of tanzania: a spatial explicit model. Population health metrics 13, 1–14. 18. Hagenlocher, M., Kienberger, S., Lang, S., Blaschke, T., 2014. Implications of spatial scales and reporting units for the spatial modelling of vulnerability to vector-borne diseases. GI_Forum 2014, 197. 19. Katoh, N., Ibaraki, T., 1998. Resource allocation problems, in: Handbook of combinatorial optimization. Springer, pp. 905–1006. 20. Kelman, I., 2018. Lost for words amongst disaster risk science vocabulary? International Journal of Disaster Risk Science 9, 281–291. 21. Kelman, I., 2020. Covid-19: what is the disaster? Social Anthropology. 22. Kienberger, S., Hagenlocher, M., 2014. Spatial-explicit modeling of social vulnerability to malaria in east africa. International journal of health geographics 13, 1–16. 23. Koyuncu, M., Erol, R., 2010. Optimal resource allocation model to mitigate the impact of pandemic influenza: a case study for turkey. Journal of medical systems 34, 61–70. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20192056&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F01%2F2021.02.24.21252407.atom) 24. Lederberg, J., Hamburg, M.A., Smolinski, M.S., et al., 2003. Microbial threats to health: emergence, detection, and response. National Academies Press. 25. Li, M., Guo, P., Singh, V.P., Yang, G., 2016. An uncertainty-based frame-work for agricultural water-land resources allocation and risk evaluation. Agricultural Water Management 177, 10–23. 26. Madhav, N., Oppenheim, B., Gallivan, M., Mulembakani, P., Rubin, E., Wolfe, N., 2017. Pandemics: risks, impacts, and mitigation, in: Disease Control Priorities: Improving Health and Reducing Poverty. 3rd edition. The International Bank for Reconstruction and Development/The World Bank. 27. de Mattos Almeida, M.C., Caiaffa, W.T., Assunçao, R.M., Proietti, F.A., 2007. Spatial vulnerability to dengue in a brazilian urban area during a 7-year surveillance. Journal of Urban Health 84, 334–345. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11524-006-9154-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17243024&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F01%2F2021.02.24.21252407.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246522000003&link_type=ISI) 28. Mishra, S.V., Gayen, A., Haque, S.M., 2020. Covid-19 and urban vulnerability in india. Habitat international 103, 102230. 29. Moore, M., Gelfeld, B., Adeyemi Okunogbe, C.P., 2017. Identifying future disease hot spots: infectious disease vulnerability index. Rand health quarterly 6. 30. Murphy, C., Gardoni, P., 2007. Determining public policy and resource allocation priorities for mitigating natural hazards: A capabilities-based approach. Science and Engineering Ethics 13, 489–504. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18066680&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F01%2F2021.02.24.21252407.atom) 31. PAHO, 2009. Leadership during a pandemic: What your municipality can do. 32. Prieto, J., Gomez, J., Leon, E., 2019. Multi-objective evolutionary algo-rithm for dna codeword design, in: Proceedings of the Genetic and Evo-lutionary Computation Conference, pp. 604–611. 33. Prieto, J., Malagón, R., Gomez, J., León, E., 2021. Urban vulnerability assessment for pandemic surveillance. medRxiv. 34. Qiu, Y., Ge, P., Yim, S.C., 2008. Risk-based resource allocation for collaborative system design in a distributed environment. Journal of Mechanical Design 130. 35. Quadrifoglio, L., 2008. A bottom-up risk-based resource allocation methodology to counter terrorism. International journal of society sys-tems science 1, 4–25. 36. Ray, J., Boggs, P., Gay, D., Lemaster, M., Ehlen, M., 2009. Risk-based decision making for staggered bioterrorist attacks: Resource allocation and risk reduction in “reload” scenarios. Sandia Report, SAND2009-6008. 37. Romero, C., 2000. Risk programming for agricultural resource allocation: A multidimensional risk approach. Annals of Operations Research 94, 57–68. 38. Schutze, O., Esquivel, X., Lara, A., Coello, C.A.C., 2012. Using the averaged hausdorff distance as a performance measure in evolutionary mul-tiobjective optimization. IEEE Transactions on Evolutionary Computa-tion 16, 504–522. 39. SDM, 2018. Secretaría Distrital de Movilidad: Observatorio de Movilidad Bogotá D.C. 2017. URL: [http://www.sdp.gov.co/sites/default/files/determinantes\_de\_la\_distribucion\_espacial\_de\_informalidad\_laboral\_en\_bogota.pdf](http://www.sdp.gov.co/sites/default/files/determinantes\_de\_la\_distribucion\_espacial\_de\_informalidad_laboral_en_bogota.pdf). 40. SDP, 2011. Secretaría Distrital de Planeación: Mono-grafías de las localidades Bogotá D.C. 2011. URL: [http://www.sdp.gov.co/gestion-estudios-estrategicos/informacion-cartografia-y-estadistica/repositorio-estadistico/monografia-de-bogota-2011%5D](http://www.sdp.gov.co/gestion-estudios-estrategicos/informacion-cartografia-y-estadistica/repositorio-estadistico/monografia-de-bogota-2011). 41. SDP, 2017. Secretaría Distrital de Planeación: Mono-grafías de las localidades Bogotá D.C. 2017. URL: [http://www.sdp.gov.co/gestion-estudios-estrategicos/informacion-cartografia-y-estadistica/repositorio-estadistico/monografia-de-bogota-2017%5D](http://www.sdp.gov.co/gestion-estudios-estrategicos/informacion-cartografia-y-estadistica/repositorio-estadistico/monografia-de-bogota-2017). 42. SDP, 2018. Secretaría Distrital de Planeación: Determinantes de la distribución espacial de la informalidad laboral en Bo-gotá. URL: [http://www.sdp.gov.co/sites/default/files/determinantes\_de\_la\_distribucion\_espacial\_de\_informalidad\_laboral\_en\_bogota.pdf](http://www.sdp.gov.co/sites/default/files/determinantes\_de\_la\_distribucion\_espacial\_de\_informalidad_laboral_en_bogota.pdf). 43. Sun, L., DePuy, G.W., Evans, G.W., 2014. Multi-objective optimization models for patient allocation during a pandemic influenza outbreak. Computers & Operations Research 51, 350–359. 44. Tian, Y., Cheng, R., Zhang, X., Cheng, F., Jin, Y., 2017. An indicator-based multiobjective evolutionary algorithm with reference point adaptation for better versatility. IEEE Transactions on Evolutionary Computation 22, 609–622. 45. Tian, Y., Cheng, R., Zhang, X., Li, M., Jin, Y., 2019. Diversity assessment of multi-objective evolutionary algorithms: Performance metric and benchmark problems [research frontier]. IEEE Computational In-telligence Magazine 14, 61–74. 46. Van Der Hofstad, R., 2017. Random Graphs and Complex Networks Vol. I.. volume I. Cambridge Series in Statistical and Probabilistic Mathematics. 47. Vaziri, P., 2008. Earthquake risk mitigation: hazard identification and resource allocation. 48. Vaziri, P., Davidson, R.A., Nozick, L.K., Hosseini, M., 2010. Resource allocation for regional earthquake risk mitigation: a case study of tehran, iran. Natural hazards 53, 527–546. [GeoRef](http://medrxiv.org/lookup/external-ref?access_num=2012093321&link_type=GEOREF) 49. Who, D.G., 2011. Report of the review committee on the functioning of the international health regulations (2005) in relation to pandemic (h1n1) 2009. Sixty-fourth World Health Assembly: World Health Organization, 49–50. 50. Willis, H.H., 2007. Guiding resource allocations based on terrorism risk. Risk Analysis: An International Journal 27, 597–606. 51. Wolgin, J.M., 1975. Resource allocation and risk: A case study of small-holder agriculture in kenya. American Journal of Agricultural Eco-nomics 57, 622–630. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1238880&link_type=DOI) 52. Yan, Z., Haimes, Y.Y., 2011. Risk-based multiobjective resource allocation in hierarchical systems with multiple decisionmakers. part i: Theory and methodology. Systems Engineering 14, 1–16. 53. Yeddanapudi, S.R.K., Li, Y., McCalley, J.D., Chowdhury, A.A., Jewell, W.T., 2008. Risk-based allocation of distribution system maintenance resources. IEEE Transactions on Power Systems 23, 287–295. 54. Zhang, Q., Li, H., 2007. Moea/d: A mult iobjective evolutionary algorithm based on decomposition. IEEE Transactions on evolutionary computa-tion 11, 712–731. 55. Zolfaghari, M.R., Peyghaleh, E., 2015. Implementation of equity in re-source allocation for regional earthquake risk mitigation using two-stage stochastic programming. Risk analysis 35, 434–458. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-11.gif [11]: /embed/graphic-12.gif [12]: /embed/graphic-13.gif [13]: /embed/graphic-14.gif [14]: /embed/inline-graphic-1.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-19.gif [17]: /embed/graphic-20.gif [18]: /embed/graphic-21.gif [19]: /embed/graphic-27.gif [20]: /embed/graphic-28.gif [21]: /embed/graphic-29.gif [22]: /embed/graphic-30.gif [23]: /embed/graphic-31.gif [24]: /embed/graphic-32.gif [25]: /embed/graphic-33.gif