Next Article in Journal
Topological Data Analysis with Spherical Fuzzy Soft AHP-TOPSIS for Environmental Mitigation System
Next Article in Special Issue
A Masked Self-Supervised Pretraining Method for Face Parsing
Previous Article in Journal
A Multi-View Ensemble Width-Depth Neural Network for Short-Term Wind Power Forecasting
Previous Article in Special Issue
Vasicek Quantile and Mean Regression Models for Bounded Data: New Formulation, Mathematical Derivations, and Numerical Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Equity-Based Optimization Model to Solve the Location Problem for Healthcare Centers Applied to Hospital Beds and COVID-19 Vaccination

1
Faculty of Natural Sciences and Mathematics, Escuela Superior Politécnica del Litoral ESPOL, Guayaquil 090902, Ecuador
2
Pacific International Center for Disaster Risk Reduction, Escuela Superior Politécnica del Litoral ESPOL, Guayaquil 090902, Ecuador
3
Faculty of Engineering, Universidad Espíritu Santo, Samborondón 0901952, Ecuador
4
School of Industrial Engineering, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile
5
Eating Behavior Research Center, School of Nutrition and Dietetics, Faculty of Pharmacy, Universidad de Valparaíso, Valparaíso 2360102, Chile
6
Center of Micro-Bioinnovation, Faculty of Pharmacy, Universidad de Valparaíso, Valparaíso 2360102, Chile
*
Author to whom correspondence should be addressed.
Submission received: 31 March 2022 / Revised: 16 May 2022 / Accepted: 19 May 2022 / Published: 26 May 2022

Abstract

:
Governments must consider different issues when deciding on the location of healthcare centers. In addition to the costs of opening such centers, three further elements should be addressed: accessibility, demand, and equity. Such locations must be chosen to meet the corresponding demand, so that they guarantee a socially equitable distribution, and to ensure that they are accessible to a sufficient degree. The location of the centers must be chosen from a set of possible facilities to guarantee certain minimum standards for the operational viability of the centers. Since the set of potential locations does not necessarily cover the demand of all geographical zones, the efficiency criterion must be maximized. However, the efficient distribution of resources does not necessarily meet the equity criterion. Thus, decision-makers must consider the trade-off between these two criteria: efficiency and equity. The described problem corresponds to the challenge that governments face in seeking to minimize the impact of the pandemic on citizens, where healthcare centers may be either public hospitals that care for COVID-19 patients or vaccination points. In this paper, we focus on the problem of a zone-divided region requiring the localization of healthcare centers. We propose a non-linear programming model to solve this problem based on a coverage formula using the Gini index to measure equity and accessibility. Then, we consider an approach using epsilon constraints that makes this problem solvable with mixed integer linear computations at each iteration. A simulation algorithm is also considered to generate problem instances, while computational experiments are carried out to show the potential use of the proposed mathematical programming model. The results show that the spatial distribution influences the coverage level of the healthcare system. Nevertheless, this distribution does not reduce inequity at accessible healthcare centers, as the distribution of the supply of health centers must be incorporated into the decision-making process.

1. Introduction and State-of-the-Art

In this section, we discuss introductory aspects of our investigation, detail the bibliographical review, give an illustrative example to facilitate understanding of the problem under analysis, and describe the contribution of this study, as well as detailing the contents of the sections of the paper.

1.1. Introductory Aspects

In public policy, social investment is one major concern for governments, especially when improving public healthcare systems. Access to these systems embraces a number of dimensions, including acceptability, accommodation, affordability, availability, and geographical location [1,2,3]. Some authors have addressed other barriers to be considered when measuring accessibility, such as dimensions related to age, behavior, ethnicity, financial status, information access, sex, and social class [4,5]. These dimensions allow the clustering of different accessibility measures according to whether spatial barriers are included in their quantification or not.
Certain decisions related to adequate health access focus on locating healthcare centers/facilities. This location is chosen, for example, to situate a set of additional centers or to provide a good supply of services, such as hospital beds, which improves the relationship between supply and demand. Unlike other facility location problems that focus on minimizing costs, decision-makers face complex challenges when deciding on locating the best sites for healthcare centers due to their possible inadequate spatial distribution or an insufficient number of them. These possibly counterproductive decisions can have a direct impact on the healthcare workers, increase the mortality and morbidity of the community [6], as well as limiting access to medical treatment [5,7,8,9,10].
Although the location of healthcare centers may meet an efficiency criterion, it must also guarantee an equitable distribution of resources [8], which can increase the inequity of access to healthcare centers [11,12]. If the focus is on a zone-divided region requiring the localization of healthcare centers, in each zone, it is necessary to select a geographical site (discrete center) where the demand is concentrated. Then, the decision-makers must determine where to situate the healthcare centers based on a set of possible locations that guarantee minimum and maximum thresholds for the operational viability of these centers. Since the set of potential locations does not necessarily cover the demand of all zones, decision-makers must maximize an efficiency criterion. However, efficient distribution of resources does not necessarily meet the equity criterion. Therefore, decision-makers must also consider the trade-off between these two criteria: efficiency and equity.
There are two types of accessibility: (i) revealed, which focuses on measuring the employment of the healthcare system; and (ii) potential, which measures the probable use of healthcare centers and not their utilization [5]. Some spatial measures have been studied on healthcare accessibility [13], many of which are based on the gravity model given by A i = j a j / f ( c i j ) , where i and j denote the healthcare demand i and supply j points, respectively; a j is the healthcare supply of the point j; c i j is the total cost between points i and j; and f is an impedance function. This function is decreasing and weights the supply from point j captured by the demand point i, that is, as the cost between points i and j increases, the supply of point j captured by point i decreases.
Some variants of the gravity model have been analyzed [1]. One of them is the two-step floating catchment area (2SFCA) method [5]. This method involves two stages: (i) the first calculates the ratio between the supply of healthcare centers of location j and the population that is distanced (or has a travel time) to a maximum threshold D of point j; and (ii) in the second, the accessibility of each demand point i is obtained as the sum of the supply-to-population ratio (calculated at the first stage) of all providers that are separated (or have a travel time) at most D from point i. Several studies adopted the 2SFCA method, or a variant of it, for analyzing accessibility to healthcare centers [14,15,16]. This accessibility is a non-medical factor that influences inequities in healthcare outcomes [11]; that is, there is a relation between inequity access to healthcare and inequitable distribution of illness [11,12,17,18]. For this reason, the World Healthcare Organization has issued recommendations for government actions. These actions are related to strengthening equity by improving investment in healthcare infrastructure, among other initiatives, to close the gap in healthcare outcomes [12].

1.2. Bibliographical Review

Several optimization models have been formulated for locating a set of healthcare facilities. Some of these models optimize an efficiency criterion, for example, to minimize the total travel distance or time [19,20,21], or to minimize the total operating and installation costs [22,23]. Other studies have focused on characterizing the location of healthcare facilities as a maximal covering location problem, which aims to situate a fixed number of healthcare facilities to maximize the demand coverage by these facilities [24]. The notion of covering a demand point refers to whether there is at least one facility at a maximum travel distance or threshold (or travel time) D from the demand point.
Various applications have adopted different coverage thresholds for demand points by healthcare facilities. For example, a travel distance threshold of fifteen miles was stated in [25] when studying spatial accessibility to healthcare centers in Illinois, United States (US), using a geographic information system. Similarly, a travel time threshold of three hours was assumed in [26] when analyzing spatial accessibility to major cancer care facilities in the US. In addition, a threshold defined as the third quartile of the data set was employed in a mobility survey in Florida, US [15]. A variant of the maximal covering location problem that included capacity and budget constraints was addressed in [27,28]. These authors assumed the number of facilities as a decision variable.
Other studies have adopted a multi-objective approach when locating a set of facilities. These objectives include costs, environment, equity concerns, and service level. For example, a two-criteria mathematical model was proposed in [29] for solving a multi-objective facility location problem. This model deals with efficiency (to minimize total costs) and equity concerns (to minimize the maximum distances from the facilities to the demand points). For an extensive literature review of the facility location problem under a multi-criteria approach, see [30], and for other p-center variants in a 2D Pareto front, see [31].
In the context of healthcare facility locations, under a multi-objective approach, a mathematical model was proposed in [23] for locating and sizing a set of medical departments based on finding a trade-off among costs, total travel distance, uncovered demand, and the number of changes to be implemented in the system. A similar approach was adopted in [21] for an application in a hierarchical healthcare system. Other investigations on this topic can be found in [32,33,34]. Similarly, a three-objective mathematical model was presented in [35] to locate preventive healthcare facilities, which included financial, coverage, and service level (waiting time) concerns.
Various metrics have been adopted to measure resource distribution and inequity. Among the measures of equity used in the context of facility location analysis, we have the standard deviation, variance, mean absolute deviation, and sum of absolute deviation [36,37]. These statistical metrics cannot be employed in any comparative study since they are not normalized metrics. Among normalized equity metrics, we have the Gini index, the Schutz index, and the coefficient of variation [38,39].
In the context of equity, the Gini index has been used in the selection of socioeconomic variables, in the establishment of relations of econometric models [40], and in the determination of conditioning factors for environmental equity [41], as a helpful method to establish equity in network programming problems [42].

1.3. Illustrative Example

To illustrate the problem that we address, let us suppose that a study region has been divided into six zones, each characterized by a discrete center in which a person is concentrated; see Figure 1a,b. Two healthcare centers, namely A and B, are located, with a healthcare supply of 1 and 2 units, respectively. The figures also show the coverage zones of each healthcare center (to which this center potentially offers its service).
Analyzing the scenario represented in Figure 1a, note that healthcare center A potentially offers one center to the four patients in its coverage region (zones 1, 2, 4, and 5). Hence, each zone receives 1 / 4 health units of center A, while healthcare center B potentially offers two centers to the four patients in its coverage region (zones 2 , 3 , 5 and 6).
Thus, each patient receives from center B 2 / 4 = 1 / 2 health units. Given that zones 2 and 5 can potentially access health centers offered by both centers A and B, the potential spatial accessibility to such centers of these zones is the sum of what they potentially receive from the centers A and B. Then, the potential spatial accessibility from zones 1, 2, 3, 4, 5 and 6 to healthcare centers is 1 / 4 , 3 / 4 , 1 / 4 , 1 / 4 , 3 / 4 , and  1 / 4 , respectively, yielding a Gini index of 0.267 , as computed in [38,39]. Moreover, if the locations of the health centers are those shown in Figure 1b, center A offers to each patient in zones 1, 2, 4 and 5 a value of 1 / 4 units of healthcare centers, while center B offers to each patient in zones 3 and 6 a value of 2 / 2 = 1 unit of healthcare center. Therefore, the potential spatial accessibility of zones 1, 2, 3, 4, 5 and 6 to healthcare centers is 1 / 4 , 1 / 4 , 1, 1 / 4 , 1 / 4 and 1, respectively, which produces a Gini index of 0.333 . Although all the zones of the hypothetical study region are covered by the health system in both scenarios, it is evident that an inadequate location of the centers could increase the inequity in the distribution of resources, evidenced by the difference in the value of the Gini index.

1.4. Contributions and Description of Sections

Since an efficient resource distribution does not necessarily guarantee equity and vice-versa [8], we propose a multi-objective approach. Our main contribution here is to solve the problem of location for healthcare centers considering, in addition to the costs of opening such centers, the trade-off among three elements: accessibility, demand, and equity, based on a mathematical programming model. To the best of our knowledge, these three elements have not been considered together until now in a multi-objective nonlinear programming model to solve this problem based on a coverage formulation using the Gini index to measure equity and accessibility. Our proposal seeks to set an efficiency criterion when covering the locations of healthcare centers subject to the criteria of equity and accessibility. The problem described fits perfectly with the challenge that governments face to minimize the impact of the pandemic on citizens, where healthcare centers can be either public hospitals that care for COVID-19 patients or vaccination points [43,44,45,46]. In this study, we focus on a type of accessibility named potential, due to confidentiality, and on the problem of locating healthcare centers in a zone-divided region.
From an algorithmic point of view, the problem to be solved here is relevant, since this is a discrete p-center problem, which has been proven to be NP-hard [47,48]. Although some variants of this problem can be solved in polynomial time [49,50], unfortunately, in the variant that we want to solve, this does not happen. Then, approximation algorithms are necessary to solve instances in a reasonable time.
Our study is structured in the following way. First, in Section 2, we formulate the proposed mathematical model. Then, in Section 3, an algorithm is presented for solving the optimization model. Next, in Section 4, numerical experiments are carried out based on simulations and examples. Finally, we end the article in Section 5 with discussion of some advantages and limitations of our study, as well as providing recommendations for further research.

2. Mathematical Model

In this section, we describe some notations and propose a three-objective mathematical programming model for the problem considered in this study that embeds efficiency and equity in the decision-making process of locating a set of healthcare centers.

2.1. Notations and Symbols

Before formulating the model, Table 1 presents a summary of some notations to be used throughout the document.

2.2. Objective Functions

We define a three-objective mathematical model according to:
(i)
Coverage. Maximize the coverage of the demand points, that is, maximize an efficiency criterion given by
max f 1 = i = 1 n w i x i ,
with n, w i and x i being defined in Table 1.
(ii)
Equity. Minimize inequity in access to healthcare centers. For this purpose, we use the Gini index as an inequity measure [51], which is defined as one-half of the relative mean absolute difference between all pairs of potential spatial accessibility (a measure of the potential use of the healthcare center) a i , with  i N , that is, one-half of the mean absolute difference between all pairs of potential spatial accessibility a i divided by the average of the values of a. Specifically, to minimize inequity in access to healthcare centers, we formulate the objective function given by
min f 2 = i = 1 n l = 1 n | a i a l | / 2 n i = 1 n a i ,
with n and a i being defined in Table 1.
(iii)
Accessibility. The third objective incorporated into our model is focused on maximizing the minimum potential spatial accessibility difference to zero and is stated as
max f 3 = min i N : a i 0 { a i } ,
with N and a i being defined in Table 1.

2.3. Constraints

To meet the three objectives described previously in (i)–(iii), a set of constraints must be satisfied, which defines its domain set detailed below.
The following inequalities allow us to identify, through variable x i , whether a demand zone i is covered by the healthcare center; that is, if there is at least one operational healthcare center at most distance D of zone i. In fact, if there is not a healthcare center located at most D of zone i, then j = 1 : d i j D m y j = 0 , which leads to x i = 0 . Reciprocally, if  x i = 0 , then it leads to y j = 0 . Therefore, we have the constraints established as
j = 1 : d i j D m y j x i , i N , x i y j , i N , j M , d i j D ,
with d i j , D, m, M, and  y j being defined in Table 1. The constraints also ensure the opening of a fixed number of facilities specified in advance, which are given by j = 1 m y j = k . To compute the potential spatial accessibility of a zone i according to the 2SFCA method, the constraints are stated as
r j = s j i = 1 : d i j D n w i , j M ,
a i = j = 1 : d i j D m r j , i N ,
with r j being defined in Table 1. Specifically, the constraint given in (5) calculates the supply-to-population ratio (within the coverage zone of a healthcare center j) of the healthcare center j, while the constraint defined in (6) computes the potential spatial accessibility of zone i to the healthcare centers.
To guarantee the minimum and maximum supply of a health center, based on the constraint given in (7), and the minimum and maximum global index of supply at the analyzed region, based on the constraint expressed in (8), the corresponding constraints are specified as
s min y j s j s max y j , j M ,
HBI min j = 1 m s j HBI max ,
with s j being defined in Table 1 and HBI denoting the hospital bed index also stated in this table. The basic domain constraints of our proposed model are presented as
y j , x i { 0 , 1 } , i N , j M , s j Z + { 0 } , j M , a i , r j 0 , i N , j M ,

3. Solution Algorithm for the Three-Objective Optimization Model

In this section, we describe the ε -constraint approach for obtaining the Pareto optimal set of the multi-objective mathematical model and the corresponding algorithm to optimize the model.

3.1. ε -Constraint Approach

To meet the different objectives described in the expressions given in (1)–(3), we adopt an ε -constraint method [52] for obtaining the Pareto optimal set and establishing the equilibrium point where we find the highest achievable well-being for society. For this purpose, we define as our primary objective to maximize the coverage of the demand points stated in (1). In contrast, the objectives of inequity defined in (2) and accessibility established in (3) are set as additional constraints. Thus, we define ε acc and ε GI as lower and upper bounds of inequity in minimum potential spatial accessibility and access to healthcare centers, respectively. Based on the above, the constraints derived when defining upper and lower bounds to the objectives defined in (2) and (3) are stated as
f 2 ε GI , f 3 ε acc ,
respectively. Note that the constraints defined in (10) are non-linear, so we proceed to linearize them by introducing auxiliary variables and constraints.
Given that there is a minimum supply of healthcare centers, i = 1 n a i > 0 namely, the expression given in (10) can be formulated as
i = 1 n l = 1 n | a i a l | ε GI 2 n i = 1 n a i .
Let h i l be a positive variable for each i , l N . Then, the constraints stated in (11) can be linearized by the expressions defined as
i = 1 n l = 1 n h i l = ε GI 2 n i = 1 n a i ,
where h i l a i a l h i l , i , l . Note that we must choose the minimum potential spatial accessibility among all zones. Hence, we define, for each zone i, a positive variable A i that takes the value a i providing that zone i is covered by a healthcare center, which is characterized by x i = 1 or x i = 0 . Let B be a large number and z be a positive variable such that z = min i N : a i 0 { a i } . Hence, the inequality stated in (10) can be replaced by z ε acc as long as for each zone i, A i = a i + B ( 1 x i ) , and  z A i , that is,
ε acc z a i + B ( 1 x i ) , i N .
Based on the above, the Pareto optimal set can be obtained by solving it iteratively. Then, the single-objective model is given by
max { f 1 } , subject to j = 1 : d i j D m y j x i , i N , x i y j , i N , j M , d i j D , j = 1 m y j = k , r j = s j / i = 1 : d i j D n w i , j M , a i = j = 1 : d i j D m r j , i N , s min y j s j s max y j , j M , HBI min j = 1 m s j HBI max , i = 1 n l = 1 n h i l = ε GI 2 n i = 1 n a i , ε acc z a i + B ( 1 x i ) , i N , y j , x i { 0 , 1 } , i N , j M , s j Z + { 0 } , j M , a i , r j , h i l , z 0 , i , l N , j M ,
for different values of ε GI and ε acc , where f 1 is defined in (1).
One of the disadvantages of the ε -constraint approach for solving a multi-objective mathematical model is the computational time spent when exploring ε GI and ε acc values that generate infeasible solutions. Hence, we must define a range for both values. By definition, the Gini index has a maximum value of one, whereas the minimum index Gini can be obtained by solving the mathematical model stated as
ε GI lower = max { f 2 } , subject to j = 1 : d i j D m y j x i , i N , x i y j , i N , j M , d i j D , j = 1 m y j = k , r j = s j / i = 1 : d i j D n w i , j M , a i = j = 1 : d i j D m r j , i N , s min y j s j s max y j , j M , HBI min j = 1 m s j HBI max , y j , x i { 0 , 1 } , i N , j M , s j Z + { 0 } , j M , a i , r j 0 , i N , j M ,
where f 2 is defined in (2).
Note that the lower minimum potential spatial accessibility is zero, whereas its upper minimum value can be obtained by solving the mathematical model expressed as
ε acc upper = max { f 3 } , subject to j = 1 : d i j D m y j x i , i N , x i y j , i N , j M , d i j D , j = 1 m y j = k , r j = s j / i = 1 : d i j D n w i , j M , a i = j = 1 : d i j D m r j , i N , s min y j s j s max y j , j M , HBI min j = 1 m s j HBI max , y j , x i { 0 , 1 } , i N , j M , s j Z + { 0 } , j M , a i , r j 0 , i N , j M ,
so that ε GI [ ε GI lower , 1 ] and ε acc [ 0 , ε acc upper ] , where f 3 is defined in (3).

3.2. Algorithm

Since ε GI and ε acc belong to a continuous set, we discretize their respective domains in intervals δ 1 and δ 2 defined in advance and we choose ε GI and ε acc as the limits for each interval. Thus, the procedure for obtaining the Pareto optimal set of the problem, Ω namely, is defined as { max { f 1 } , min { f 2 } , max { f 3 } } , subject to the constraints stated in (4) to (9), and formulated in Algorithm 1.
Algorithm 1 ε -constraint approach for obtaining the Pareto optimal set of the multi-objective mathematical model
Mathematics 10 01825 i009

4. Numerical Results

In this section, we propose an algorithm to generate instances of the location problem of healthcare centers and another algorithm to solve them. In addition, application examples of the proposed mathematical model are included and discussed. Furthermore, the computational burden to evaluate the efficiency of the proposed procedure is reported.

4.1. Simulation Algorithm

We propose Algorithm 2 to generate an instance of the problem defined in Section 2. To build the algorithm, we need to establish:
  • h : Upper limit on the horizontal axis and the vertical axis.
  • p : Number of inhabitants who may require healthcare.
  • l 1 : Real number between 0 and 1 that allows us to determine the radius of the interval to generate the demand of the zones.
  • l 2 : Real number between 0 and 1 that permits us to state the radius of the interval for the total zone demand.
  • K : Set of possible values for k.
  • s u min : Upper limit for s min .
  • s l max : Lower limit for s max .
  • U : A random variable with a continuous uniform distribution in ( α , β ) , denoted as U ( α , β ) .
Regarding other input and output parameters of Algorithm 2, see Table 1. It is important to mention that, on the one hand, Step 3 of Algorithm 2 allows us to assign a uniformly distributed demand to each zone in an interval centered on the ratio p / n with a radius of l 1 p / n . On the other hand, its Step 4 permits us to ensure that the value of the total demand of the zones is in an interval centered on p and with a radius of l 2 p . Similarly, s min must satisfy the condition s min s u min , while s max must satisfy the condition s max s l max , with  s min < s max .
Algorithm 2 Simulation algorithm for generating instances of the location problem for healthcare centers
Mathematics 10 01825 i010
To solve instances of the location problem for healthcare centers, we design Algorithm 3 that uses a sensor in a database. This algorithm achieves its objective by relying on a simulation algorithm, a database, a sensor located in it, and a solver responsible for obtaining the solution of the corresponding mathematical model.
Algorithm 3 Approach for solving instances of the location problem for healthcare centers using a database-sensor.
1:
Generate an instance of the problem based on Algorithm 2 and store it in a database.
2:
Detect the generation of the instance by means of a sensor in the database and build the corresponding mathematical model with the data from the instance at run-time.
3:
Send the built mathematical model to a solver by using a sensor.
4:
Receive the mathematical model from the solver and compile it.
5:
Store the results in the database and use them for later analysis and discussion.

4.2. Application Example

To analyze the impact of this novel approach when locating healthcare centers, an instance was obtained when applying Algorithm 2 with the values n = 10 , h = 1.7 , p = 1 million of inhabitants, l 1 = 0.4 , l 2 = 0.1 , and k { 3 , 4 , 5 , 6 } . The values of these parameters have been set for experimental purposes. However, they depend on the criteria of decision-makers when looking for a balance between the service level and available budget.
To test with an adequate minimum and maximum global index of healthcare center supply ( HBI min and HBI max ), we analyzed data from the World Bank Group [53] related to hospital beds in several countries worldwide. For example, in 2018, Colombia reported that the HBI was 1.7 hospital beds per 1000 inhabitants, while it was 1.4 in Ecuador in 2014. Based on the above, we adopted 1 and 2 hospital beds per 1000 inhabitants as the minimum and maximum global index of healthcare centers, respectively, a more realistic condition in South America. The execution of Algorithm 2 does not only provide us with the spatial location of the study zones and their respective demands, but it also offers upper and lower bounds for s min and s max , respectively. Based on the bounds calculated by the algorithm, the values s min = 250 and s max = 400 are obtained.
Figure 2 shows the spatial distribution of ten zones, denoted from 1 to 10, each one represented by a discrete center whose diameter is proportional to the demand for healthcare centers in each zone. In other words, zone 8 is the area with the least demand, while zone 5 has the highest demand. Each line in Figure 2 indicates a connection between a pair of zones whose length is given by the Euclidean metric. Otherwise, the distance between each pair of nodes is provided by the shortest path between the nodes.
As described in Algorithm 1, the Pareto optimal set of the instance is obtained from the application of the ε -constraint method. We vary both ε acc and ε GI in a fine grid. In this application example, we assume δ 1 = δ 2 = 500 as the number of intervals in which the domain of ε GI and ε acc is divided. We implement the models in GAMS 24.0.2 and use ILOG CPLEX 12.5 on a computer equipped with an Intel Core i7-5500U 2.66 Ghz processor and 8 GB of RAM.
The application results of the proposed instance for k { 3 , 4 , 5 , 6 } are shown in Table 2. For each value of k, the Pareto-frontier set is shown, in which each point in the set is identified as a scenario. In addition, Table 3 presents the values of objective functions for each scenario described in Table 2.
The results of the proposed instance for k = 6 show that the Pareto optimal set corresponds to the mesh sets (Figure 3a), where the Pareto-frontier set is made up of twelve non-dominated solutions (Figure 3b). Figure 4 presents the location of healthcare centers when k = 6 for different coverage values, where a black square represents the healthcare center.
An aspect worth noting is that the increase in the number of healthcare centers to be installed does not necessarily increases the level of system coverage and additionally it could increase inequity. Indeed, taking the results of Table 2, note that the maximum coverage level obtained when k = 3 (scenario 3–7) is less compared to that obtained when k = 4 (scenario 4–7). From there, when k = 5 and k = 6 , the maximum coverage level remains constant because it is the maximum value that can be achieved in instance 1. Nevertheless, when increasing from k = 5 to k = 6 healthcare facilities, the Gini index might improve, although the minimum potential spatial accessibility decreases (see Table 3 scenarios 5–7 and 6–11). However, the opposite effect happens if we compare scenarios 5–7 and 6–12 in Table 3.
The location of the facilities and healthcare center supply play an essential role in reducing inequity. From Figure 5c,d, note that, in both cases, the same level of coverage exists. Nonetheless, the inequity level in scenario 6–12 is greater compared to scenario 6–11 (Table 3). Thus, we can see that, although a zone has the largest population, a health center should not necessarily be installed (Figure 5b), even if its sanitary supply is the largest allowed when it is installed (see scenario 6–12 on Table 2 and Figure 5d).
A similar pattern can be seen in another random instance (Figure 6a), which was generated according to Algorithm 2 with its parameters such as instance 1. To solve that, we applied Algorithm 1. By analyzing a non-dominated Pareto solution set (Figure 6b), observe that, when k = 7 , the coverage level increases from 722,229 (square marks) to 751,945 (plus marks), producing a decrease in both inequity and accessibility of healthcare systems (increasing in the Gini index) and also a decrease in the minimum potential accessibility in some cases (see Table 4 and Table 5).
An exhaustive analysis of the Pareto-frontier sets shows that, in some scenarios, an increase in the total healthcare center supply of the system does not necessarily cause a decrease in inequity. For example, the total healthcare center supply in scenario 5–1 is 1723 hospital beds (Table 4) and its Gini index is 0.439 (Table 5). However, in scenario 5–7, its Gini index decreases to 0.289 although its total healthcare center supply goes down to 1549 hospital beds. A similar situation can be seen in scenarios 6–1 and 6–8 reported in Table 4 and Table 5. In addition, we can observe that, although the total healthcare center supply of the system is the same, its efficiency (measured by coverage) and inequality (measured by the Gini index) are different, which clearly highlights the need to adequately manage the available resources (see all scenarios where k = 7 in Table 4 and Table 5).

4.3. Computational Burden

To evaluate the computational performance of the proposed procedure, we apply Algorithm 2 to generate instances of different sizes. According to the inputs given in this algorithm, in all the generated instances, we considered HBI min = 1 , HBI max = 2 , h = 2 , p = 0.2 , l 1 = 0.2 , l 2 = 0.2 , and k { 5 , 6 , 7 , 8 , 9 } .
Table 6 reports the results of the application of our model in instances with 10, 15, 20 and 50 zones. Considering in each case different values of k, healthcare facilities are installed. Each of these instances was run within a maximum runtime of 2 days (48 h).
As can be seen in Table 6, the implementation of the model is adequate in small instances (less than 20 zones). However, it is evident that, in instances that exceed 50 zones, more than two days will be required to obtain the Pareto frontier, so the use of meta-heuristic-based approaches could be considered.

5. Conclusions, Limitations, and Future Research

In this section, we provide the conclusions of our study, as well as some of its limitations, and ideas for further investigation.

5.1. Concluding Remarks

In this paper, we addressed the location problem for healthcare centers under a multi-criteria approach that embeds efficiency and equity concerns. For this purpose, we have adopted the coverage of healthcare systems as an efficiency criterion, and the Gini index as a proxy variable that measures the inequity level on accessibility to healthcare centers.
The results obtained in this study show that the application of an efficiency criterion in the location of facilities does not guarantee equity in resource distribution and vice versa. Therefore, the adopted approach in this study was appropriate. Note that the spatial distribution influenced the coverage level of the healthcare system. Nevertheless, this is not decisive for reducing inequity in accessible healthcare centers since the distribution of healthcare center supply must be incorporated into the decision-making process. In fact, there could be situations with a low total healthcare center supply (see Table 2, scenario 6–11) that are more equitable than a scenario with a high healthcare center supply (see Table 2, scenario 6–12) for the same coverage level. However, the minimum potential spatial accessibility would be sacrificed. In contrast, a decrease in the total healthcare center supply could increase inequality and decrease minimum potential spatial accessibility, although it might increase the coverage of the healthcare system.
For example, results from the numerical application showed 10 non-dominated Pareto solutions with a coverage level of 915,367 inhabitants (circle marks in Figure 7), which performed poorly both in equity and minimum potential spatial accessibility compared to non-dominated Pareto solutions with a coverage level of 949,469 inhabitants (plus marks). Therefore, the performance of the suggested actions depends to a great extent on the number of facilities to be installed and on the total healthcare center supply as well. This gives us an opportunity to expand the study by incorporating budget constraints. In this context, from the results reported in Table 2 and Table 3, we have noticed that an increase in total healthcare center supply could increase the coverage of the system. Nonetheless, this does not necessarily reduce inequity. For example, scenario 4–4 had a total healthcare center supply of 1175 hospital beds distributed in zones 3, 7, 9 and 10, which covered 915,367 inhabitants and produced a Gini index of 0.111. In addition, increasing the total healthcare center supply up to 1302 hospital beds (scenario 4–5) distributed into zones 1, 5, 7 and 10, the coverage increases up to 933,294 inhabitants, but the Gini index was 0.151. However, this fact cannot be considered a rule since there are cases in which an increase in the total healthcare center supply produced an increase in the coverage and a decrease in the Gini index, that is, the inequality declined. For example, this happened when comparing scenarios 3–1 and 3–2.

5.2. Limitations and Future Research

A limitation of the present investigation is that the simulation model could not be used outside the limits of the parameters studied in which it was constructed. This may cause a false appreciation of the multi-objective linear programming problem when considering other scenarios. Undoubtedly, future research should analyze other efficiency criteria that contribute to reducing the inequality gap. In addition, some variants to this problem can be addressed in future research, such as the inclusion of decay functions in the expression measuring potential spatial accessibility. This could reflect a gradual supply capture according to the distance that separates a demand zone from each healthcare center. It is also important to consider restrictions on the budgetary availability of resources to install healthcare centers, which may vary according to geographical zones.
Since the method we use to solve our problem, that is, the ε -constraint method, deals with mixed integer linear problems at each iteration, it is important to mention that a study of techniques that employ, for example, reformulations or relaxations [2,54,55] might be applied to improve the performance of the method and its execution times.
Finally, due to what was stated in Section 4.3 about runtimes, the utilization of heuristic techniques is something that should be explored, even when there is no guarantee of obtaining an optimal solution when this approach is applied.

Author Contributions

Conceptualization, E.J.D., X.C., C.M.-B.; data curation, E.J.D., X.C., C.M.-B.; formal analysis, E.J.D., X.C., C.M.-B., V.L., F.R.; investigation, E.J.D., X.C., C.M.-B., F.R., V.L.; methodology, E.J.D., X.C., C.M.-B., F.R., V.L.; writing—original draft, E.J.D., X.C., C.M.-B., F.R.; writing—review and editing, V.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by Decanato de Investigación de la Escuela Superior Politécnica del Litoral ESPOL, Ecuador (E.J.D., X.C., C.M.-B); FONDECYT grant number 1200525 (V.L.) and FONDECYT grant number 11190004 (F.R.), both from the National Agency for Research and Development (ANID) of the Chilean government under the Ministry of Science, Technology, Knowledge, and Innovation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Codes and data are available from the authors upon request.

Acknowledgments

The authors warmly thank the Editors and Reviewers for their helpful comments which have led to an improved version of our paper.

Conflicts of Interest

There are no conflicts of interest declared by the authors.

References

  1. Neutens, T. Accessibility, equity and health care: Review and research directions for transport geographers. J. Transp. Geogr. 2015, 43, 14–27. [Google Scholar] [CrossRef]
  2. Cabezas, X.; Martin-Barreiro, C.; Garcia, S.; Delgado, E.; Leiva, V. A two-stage location problem with order solved using a Lagrangian algorithm and stochastic programming for a potential use in COVID-19 vaccination based on sensor-related data. Sensors 2021, 21, 5352. [Google Scholar] [CrossRef] [PubMed]
  3. Rojas, F.; Leiva, V.; Wanke, P.; Lillo, C.; Pascual, J. Modeling lot-size with time-dependent demand based on stochastic programming and case study of drug supply in Chile. PLoS ONE 2019, 14, e0212768. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Aday, L.; Andersen, R. A framework for the study of access to medical care. Health Serv. Res. 1974, 9, 208–220. [Google Scholar]
  5. Luo, W.; Wang, F. Measures of spatial accessibility to health care in a GIS environment: Synthesis and a case study in the Chicago region. Environ. Plan. B 2003, 30, 865–884. [Google Scholar] [CrossRef] [Green Version]
  6. Daskin, M.; Dean, L. Location of health care facilities. Oper. Res. Health Care 2005, 70, 43–76. [Google Scholar]
  7. Carr, B.G.; Branas, C.C.; Metlay, J.P.; Sullivan, A.F.; Camargo, C.A.J. Access to emergency care in the United States. Ann. Emerg. Med. 2009, 54, 261–269. [Google Scholar] [CrossRef] [Green Version]
  8. Culyer, A.J. Efficiency, Equity and Equality in Health and Health Care; Working Papers 120cherp; Centre for Health Economics, University of York: York, UK, 2015. [Google Scholar]
  9. Dai, D. Black residential segregation, disparities in spatial access to health care facilities, and late-stage breast cancer diagnosis in metropolitan Detroit. Health Place 2010, 16, 1038–1052. [Google Scholar] [CrossRef]
  10. Hiscock, R.; Pearce, J.; Blakely, T.; Witten, K. Is neighborhood access to health care provision associated with individual-level utilization and satisfaction? Health Serv. Res. 2008, 43, 2183–2200. [Google Scholar] [CrossRef] [Green Version]
  11. Commission on Social Determinants of Health. Closing the Gap in a Generation: Health Equity through Action on the Social Determinants of Health; Technical Report; World Health Organisation: Geneva, Switzerland, 2008. [Google Scholar]
  12. Meyer, S.B.; Luong, T.C.; Mamerow, L.; Ward, P.R. Inequities in access to healthcare: Analysis of national survey data across six Asia-Pacific countries. BMC Health Serv. Res. 2013, 13, 238. [Google Scholar] [CrossRef] [Green Version]
  13. Higgs, G. A literature review of the use of GIS-based measures of access to health care services. Health Serv. Outcomes Res. Methodol. 2004, 5, 119–139. [Google Scholar] [CrossRef]
  14. Boisjoly, G.; Deboosere, R.; Wasfi, R.; Orpana, H.; Manaugh, K.; Buliung, R.; El-Geneidy, A. Measuring accessibility to hospitals by public transport: An assessment of eight Canadian metropolitan regions. J. Transp. Health 2020, 18, 100916. [Google Scholar] [CrossRef]
  15. Mao, L.; Nekorchuk, D. Measuring spatial accessibility to healthcare for populations with multiple transportation modes. Health Place 2013, 24, 115–122. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, X.; Yang, H.; Duan, Z.; Pan, J. Spatial accessibility of primary health care in China: A case study in Sichuan province. Soc. Sci. Med. 2018, 209, 14–24. [Google Scholar] [CrossRef]
  17. Nunes, B.P.; Thumé, E.; Tomasi, E.; Duro, S.M.S.; Facchini, L.A. Socioeconomic inequalities in the access to and quality of health care services. Rev. Saúde Pública 2014, 48, 968–976. [Google Scholar] [CrossRef]
  18. National Research Council (US) Committee on Future Directions for Behavioral and Social Sciences Research at the National Institutes of Health. The Influence of Inequality on Health Outcomes. In New Horizons in Health; Singer, B., Ryff, C., Eds.; National Academies Press: Washington, DC, USA, 2001; Chapter 7. [Google Scholar]
  19. Cocking, C.; Flessa, S.; Reinelt, G. Improving access to health facilities in Nouna district, Burkina Faso. Socio-Econ. Plan. Sci. 2012, 46, 164–172. [Google Scholar] [CrossRef]
  20. Galvao, R.D.; Acosta Espejo, L.G.; Boffey, B. A hierarchical model for the location of perinatal facilities in the municipality of Rio de Janeiro. Eur. J. Oper. Res. 2002, 138, 495–517. [Google Scholar] [CrossRef]
  21. Mestre, A.M.; Oliveira, M.D.; Barbosa-Póvoa, A. Organizing hospitals into networks: A hierarchical and multiservice model to define location, supply and referrals in planned hospital systems. Spectrum 2012, 34, 319–348. [Google Scholar] [CrossRef]
  22. Ghaderi, A.; Jabalameli, M.S. Modeling the budget-constrained dynamic uncapacitated facility location–network design problem and solving it via two efficient heuristics: A case study of health care. Math. Comput. Model. 2013, 57, 382–400. [Google Scholar] [CrossRef]
  23. Stummer, C.; Doerner, K.; Focke, A.; Heidenberger, K. Determining location and size of medical departments in a hospital network: A multiobjective decision support approach. Health Care Manag. Sci. 2004, 7, 63–71. [Google Scholar] [CrossRef]
  24. Gunes, E.; Yaman, H.; Cekyay, B.; Verter, V. Matching patient and physician preferences in designing a primary care facility network. J. Oper. Res. Soc. 2014, 65, 483–496. [Google Scholar] [CrossRef]
  25. Luo, W.; Wang, F. Spatial accessibility to primary care and physician shortage area designation: A case study in Illinois with GIS approaches. In Geographic Information Systems and Health Applications; Khan, O.A., Skinner, R., Eds.; IGI Global: Hershey, PA, USA, 2003; pp. 261–279. [Google Scholar]
  26. Shi, X.; Alford-Teaster, J.; Onega, T.; Wang, D. Spatial access and local demand for major cancer care facilities in the United States. Ann. Assoc. Am. Geogr. 2012, 102, 1125–1134. [Google Scholar] [CrossRef]
  27. Griffin, P.M.; Scherrer, C.R.; Swann, J.L. Optimization of community health center locations and service offerings with statistical need estimation. IIE Trans. 2008, 40, 880–892. [Google Scholar] [CrossRef]
  28. Smith, H.K.; Harper, P.R.; Potts, C.N. Bicriteria efficiency/equity hierarchical location models for public service application. J. Oper. Res. Soc. 2013, 64, 500–512. [Google Scholar] [CrossRef]
  29. Bhattacharya, U.; Rao, J.; Tiwari, R. Bi-criteria multi facility location problem in fuzzy environment. Fuzzy Sets Syst. 1993, 56, 145–153. [Google Scholar] [CrossRef]
  30. Farahani, R.Z.; SteadieSeifi, M.; Asgari, N. Multiple criteria facility location problems: A survey. Appl. Math. Model. 2010, 34, 1689–1709. [Google Scholar] [CrossRef]
  31. Dupin, N.; Nielsen, F.; Talbi, E.G. Unified polynomial dynamic programming algorithms for p-center variants in a 2D Pareto front. Mathematics 2021, 9, 453. [Google Scholar] [CrossRef]
  32. Khodaparasti, S.; Maleki, H.R.; Bruni, M.E.; Jahedi, S.; Beraldi, P.; Conforti, D. Balancing efficiency and equity in location-allocation models with an application to strategic EMS design. Optim. Lett. 2016, 10, 1053–1070. [Google Scholar] [CrossRef]
  33. Zhang, W.; Cao, K.; Liu, S.; Huang, B. A multi-objective optimization approach for health-care facility location-allocation problems in highly developed cities such as Hong Kong. Comput. Environ. Urban Syst. 2016, 59, 220–230. [Google Scholar] [CrossRef]
  34. Zhou, L.; Geng, N.; Jiang, Z.; Wang, X. Multi-objective capacity allocation of hospital wards combining revenue and equity. Omega 2018, 81, 220–233. [Google Scholar] [CrossRef]
  35. Dogan, K.; Karatas, M.; Yakici, E. A model for locating preventive health care facilities. Cent. Eur. J. Oper. Res. 2020, 28, 164–172. [Google Scholar] [CrossRef]
  36. Drezner, Z.; Thisse, J.F.; Wesolowsky, G.O. The minimax-min location problem. J. Reg. Sci. 1986, 26, 87–101. [Google Scholar] [CrossRef]
  37. Maimon, O. The variance equity measure in locational decision theory. Ann. Oper. Res. 1986, 6, 147–160. [Google Scholar] [CrossRef]
  38. Eiselt, H.; Laporte, G. Objectives in location problems. In Facility Location: A Survey of Application and Method; Drezner, Z., Ed.; Springer: New York, NY, USA, 1995. [Google Scholar]
  39. Ogryczak, W. Inequality measures and equitable approaches to location problems. Eur. J. Oper. Res. 2000, 122, 374–391. [Google Scholar] [CrossRef]
  40. Wang, F.; Zhou, K.; Wang, M.; Wang, Q. The impact analysis of land features to JL1-3B nighttime light data at parcel level: Illustrated by the case of Changchun, China. Sensors 2020, 20, 5447. [Google Scholar] [CrossRef]
  41. Lay, U.S.; Pradhan, B.; Yusoff, Z.B.M.; Abdallah, A.F.B.; Aryal, J.; Park, H.J. Data mining and statistical approaches in debris-flow susceptibility modelling using airborne LiDAR data. Sensors 2019, 19, 3451. [Google Scholar] [CrossRef] [Green Version]
  42. Luque-Nieto, M.Á.; Moreno-Roldán, J.M.; Otero, P.; Poncela, J. Optimal scheduling and fair service policy for STDMA in underwater networks with acoustic communications. Sensors 2018, 18, 612. [Google Scholar] [CrossRef] [Green Version]
  43. Rojas, F.; Leiva, V.; Huerta, M.; Martin-Barreiro, C. Lot-size models with uncertain demand considering its skewness/kurtosis and stochastic programming applied to hospital pharmacy with sensor-related COVID-19 data. Sensors 2021, 21, 5198. [Google Scholar] [CrossRef]
  44. Cortes-Carvajal, P.D.; Cubilla-Montilla, M.; Gonzalez-Cortes, D.R. Estimation of the instantaneous reproduction number and its confidence interval for modeling the COVID-19 pandemic. Mathematics 2021, 10, 287. [Google Scholar] [CrossRef]
  45. Alkadya, W.; ElBahnasy, K.; Leiva, V.; Gad, W. Classifying COVID-19 based on amino acids encoding with machine learning algorithms. Chemom. Intell. Lab. Syst. 2022, 224, 104535. [Google Scholar] [CrossRef]
  46. Xu, J.; Tang, Y. Bayesian framework for multi-wave COVID-19 epidemic analysis using empirical vaccination data. Mathematics 2021, 10, 21. [Google Scholar] [CrossRef]
  47. Hsu, W.L.; Nemhauser, G. Easy and hard bottleneck location problems. Discret. Appl. Math. 1979, 1, 209–215. [Google Scholar] [CrossRef] [Green Version]
  48. Hochbaum, D. When are NP-hard location problems easy? Oper. Res. 1984, 1, 201–214. [Google Scholar] [CrossRef]
  49. Megiddo, N.; Tamir, A.; Dupin, N.; Nielsen, F.; Talbi, E.G. New results on the complexity of p-centre problems. SIAM J. Comput. 1983, 12, 751–758. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, D.; Li, J.; Wang, H. Efficient algorithms for the one-dimensional k-center problem. Theor. Comput. Sci. 2015, 592, 135–142. [Google Scholar] [CrossRef]
  51. Gini, C. Variabilita e Mutuabilita. Contributo allo Studio delle Distribuzioni e delle Relazioni Statistiche; [Fasc. I.]; Tipogr. di P. Cuppini: Bologna, Italy, 1912. [Google Scholar]
  52. Marler, R.T.; Arora, J.S. Survey of multi-objective optimization methods for engineering. Struct. Multidiscip. Optim. 2004, 26, 369–395. [Google Scholar] [CrossRef]
  53. World Bank Group. Hospital Beds per 1000 Habitants. Available online: https://data.worldbank.org/indicator/SH.MED.BEDS.ZS?name_desc=true (accessed on 9 August 2021).
  54. Calik, H.; Tansel, B. Double bound method for solving the p-center location problem. Comput. Oper. Res. 2013, 40, 2991–2999. [Google Scholar] [CrossRef] [Green Version]
  55. Elloumi, S.; Labbé, M.; Pochet, Y. A new formulation and resolution method for the p-center problem. INFORMS J. Comput. 2004, 16, 84–94. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Hypothetical spatial distribution of health centers and their coverage zone in: (a) scenario 1; and (b) scenario 2; where “A” and “B” in both plots are the healthcare centers, whereas “1–6” are the geographical zones, with Mathematics 10 01825 i001 and Mathematics 10 01825 i002 sketching a zone discrete center and a center location, respectively.
Figure 1. Hypothetical spatial distribution of health centers and their coverage zone in: (a) scenario 1; and (b) scenario 2; where “A” and “B” in both plots are the healthcare centers, whereas “1–6” are the geographical zones, with Mathematics 10 01825 i001 and Mathematics 10 01825 i002 sketching a zone discrete center and a center location, respectively.
Mathematics 10 01825 g001
Figure 2. Spatial location of zones in a hypothetical instance 1, where 1 to 10 are the geographical zones, with Mathematics 10 01825 i003 sketching a discrete center whose diameter is proportional to the demand for healthcare centers in the indicated zone.
Figure 2. Spatial location of zones in a hypothetical instance 1, where 1 to 10 are the geographical zones, with Mathematics 10 01825 i003 sketching a discrete center whose diameter is proportional to the demand for healthcare centers in the indicated zone.
Mathematics 10 01825 g002
Figure 3. Solution for k = 6 with: (a) Pareto optimal set; and (b) Pareto-frontier set; where black planes to the left indicate the mesh sets and black dots to the right are the non-dominated solutions.
Figure 3. Solution for k = 6 with: (a) Pareto optimal set; and (b) Pareto-frontier set; where black planes to the left indicate the mesh sets and black dots to the right are the non-dominated solutions.
Mathematics 10 01825 g003
Figure 4. Spatial location of healthcare centers for different coverage levels with k = 6 and: (a) 858,449 inhabitants; (b) 915,367 inhabitants; (c) 949,469 inhabitants; and (d) 1,044,779 inhabitants; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i004 and Mathematics 10 01825 i005 sketching a zone discrete center and a healthcare center location, respectively.
Figure 4. Spatial location of healthcare centers for different coverage levels with k = 6 and: (a) 858,449 inhabitants; (b) 915,367 inhabitants; (c) 949,469 inhabitants; and (d) 1,044,779 inhabitants; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i004 and Mathematics 10 01825 i005 sketching a zone discrete center and a healthcare center location, respectively.
Mathematics 10 01825 g004
Figure 5. Spatial location of healthcare centers with a maximum coverage level (1,044,779) for different values of k in instance 1 with: (a) k = 3 ; (b) k = 4 ; (c) k = 5 ; and (d) k = 6 ; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i006 and Mathematics 10 01825 i007 sketching a zone discrete center and a healthcare center location, respectively.
Figure 5. Spatial location of healthcare centers with a maximum coverage level (1,044,779) for different values of k in instance 1 with: (a) k = 3 ; (b) k = 4 ; (c) k = 5 ; and (d) k = 6 ; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i006 and Mathematics 10 01825 i007 sketching a zone discrete center and a healthcare center location, respectively.
Mathematics 10 01825 g005
Figure 6. A randomly generated instance 2 for: (a) spatial zone distribution; and (b) non-dominated Pareto solutions when k = 7 ; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i008 sketching a discrete center whose diameter is proportional to the demand for healthcare centers in the indicated zone; and the indicated optimum coverage level f 1 * .
Figure 6. A randomly generated instance 2 for: (a) spatial zone distribution; and (b) non-dominated Pareto solutions when k = 7 ; where 1 to 10 are the geographical zones, with Mathematics 10 01825 i008 sketching a discrete center whose diameter is proportional to the demand for healthcare centers in the indicated zone; and the indicated optimum coverage level f 1 * .
Mathematics 10 01825 g006
Figure 7. Spatial location of health centers for the indicated optimum coverage level f 1 * with k = 6 .
Figure 7. Spatial location of health centers for the indicated optimum coverage level f 1 * with k = 6 .
Mathematics 10 01825 g007
Table 1. Acronyms, notations and symbols.
Table 1. Acronyms, notations and symbols.
Acronym
Notation/SymbolDefinition
Sets
N { 1 , , n } Set of demand zones.
M { 1 , , m } Set of potential locations for a healthcare center.
Parameters
kNumber of healthcare facilities to be installed.
DThreshold of travel distance (or travel time).
w i Demand of zone i.
d i j Distance between zone i and zone j.
s min ( s max )
Minimum (maximum) supply of a healthcare center measured in terms of
hospital bed units.
HBI min ( HBI max )
Minimum (maximum) hospital bed index measured in terms of
hospital bed units per 1000 inhabitants.
BA big number.
Decision variables
y j One (1) if a healthcare center is located at j, and zero (0) otherwise.
x i

One (1) if the demand zone i is covered by the healthcare center,
that is, there is at least one operational healthcare center located at
most D of zone i), and zero (0) otherwise.
s j 0
Supply of healthcare center located at zone j measured in terms of
hospital beds.
a i 0
Potential spatial accessibility of zone i measured in terms of
hospital beds.
r j 0
Supply-to-population ratio within the coverage zone of the
healthcare center j measured in terms of hospital beds per inhabitants.
A i 0 , h i l 0 , z 0 Auxiliary variables.
Table 2. Pareto frontier set for the listed k, scenario, zone, and other indicators in instance 1.
Table 2. Pareto frontier set for the listed k, scenario, zone, and other indicators in instance 1.
kScenario ε GI ε acc (Hospital Beds/ 1000 Inhabitants)Coverage (Inhabitants)Healthcare Supply (Hospital Beds by Zone)Total (Hospital Beds)
12345678910
33–10.7074002.78100315,74200004000395250001045
3–20.4013251.94052597,62500400000362039901161
3–30.4013251.79220633,72700347004000040001147
3–40.3301001.32870708,57240000003500297001047
3–50.2126751.25454803,88240000000250004001050
3–60.2146001.24218807,83240000400002500001050
3–70.1664750.97026915,36700400000251003941045
44–10.3012251.94052727,0370040002520362039901413
4–20.2011251.79220820,0570034400400337040001481
4–30.2377001.32870837,9844000002502982500001198
4–40.1106501.25454915,3670025000025002754001175
4–50.1510751.25454933,2944000002500252004001302
4–60.1530001.24218937,2444000040025002510001301
4–70.1145000.970261,044,7790000250025204004001302
55–10.3378001.80456729,037025034300263004004001656
5–20.3397251.92198729,037026239400297003994001752
5–30.3416501.94052729,037026940000285003994001753
5–40.3435752.14446729,037352003770260004003101699
5–50.3455002.21244729,037372003990276004003091756
5–60.1068001.79220949,469003430250400336039901728
5–70.0529001.254541,044,779003020250025002504001452
66–10.3397252.49054729,03737102673970274002504001959
6–20.2454002.17536858,44936204003882822680003092009
6–30.2473252.21244858,44937204004002872770003082044
6–40.1453002.10738915,36734100366025303934003092062
6–50.1472252.13210915,36734800374025804004003092089
6–60.1587751.34724949,46925003060278400040025001884
6–70.1607001.39050949,46925002500278400400032401902
6–80.1626251.43994949,46925002500277400400034301920
6–90.1645501.48938949,46925003630278400400025001941
6–100.1664751.79220949,46900400025240025025040001952
6–110.0471251.254541,044,77931000250260037404003991993
6–120.0875501.328701,044,77940027404002524003630002089
Where GI: Gini index.
Table 3. Objective values in the listed scenario, zone, and other indicators for instance 1.
Table 3. Objective values in the listed scenario, zone, and other indicators for instance 1.
ScenarioCoverageGI ( * * ) Potential Spatial Accessibility (Hospital Beds by Zone/1000 Inhabitants)MPSA ( * )
12345678910
3–1315,7420.7074000003.090903.46163.4616003.0909
3–2597,6250.400011.94261.94261.94260001.94281.94281.942601.9426
3–3633,7270.401321.81621.81621.8162001.7983001.81621.79831.7983
3–4708,5720.329761.334201.3342001.57351.59391.59391.33421.57351.3342
3–5803,8820.211661.334201.33421.258901.25891.34171.34171.33421.25891.2589
3–6807,8320.213941.33421.24341.33421.2434001.34171.34171.33421.24341.2434
3–7915,3670.166460.97250.97250.97251.240001.24001.34711.34710.97251.24000.9725
4–1727,0370.300211.94261.94261.942601.947301.94281.94281.942601.9426
4–2820,0570.200891.80891.80891.8089001.79831.80861.80861.80891.79831.7983
4–3837,9840.237561.334201.334201.93181.33971.34171.34171.33421.33971.3342
4–4915,3670.110631.27651.27651.27651.258901.25891.34171.34171.27651.25891.2589
4–5933,2940.150891.334201.33421.25891.93181.25891.35241.35241.33421.25891.2589
4–6937,2440.152951.33421.24341.33421.24341.931801.34711.34711.33421.24341.2434
4–71,044,7790.114210.97250.97250.97251.25891.93181.25891.35241.35240.97251.25890.9725
5–1729,0370.337751.80652.44222.44221.894602.4413002.44222.44131.8065
5–2729,0370.339711.92812.59432.59431.925102.5941002.59432.59411.9251
5–3729,0370.341621.94262.62662.62661.942902.5402002.62662.54021.9426
5–4729,0370.343502.14662.14452.14662.147602.1445002.14663.31652.1445
5–5729,0370.344512.21332.21292.21332.212802.2133002.21333.45362.2128
5–6949,4690.106771.80411.80411.804101.93181.79831.80331.80331.80411.79831.7983
5–71,044,7790.051281.34211.34211.34211.25891.93181.25891.34171.34171.34211.25891.2589
6–1729,0370.339722.49452.49112.49452.493002.4907002.49453.72482.4907
6–2858,4490.245362.18002.17872.18002.17862.17912.1773002.18003.38352.1773
6–3858,4490.246142.21332.21602.21332.21282.21772.2147002.21333.45812.2128
6–491,5360.145292.10992.11032.10992.110202.10992.10922.10922.10993.24762.1092
6–591,5360.146292.13332.13512.13332.135102.13242.14672.14672.13333.29502.1324
6–6949,4690.158772.18571.35182.185702.14821.79832.14672.14672.18571.79831.3518
6–7949,4690.160622.22951.39562.229502.14821.79832.14672.14672.22951.79831.3956
6–8949,4690.162612.27561.44182.275602.14051.79832.14672.14672.27561.79831.4418
6–9949,4690.164502.32431.49042.324302.14821.79832.14672.14672.32431.79831.4904
6–10949,4690.166341.94511.94511.945101.94731.79832.68342.68341.94511.79831.7983
6–111,044,7790.047042.00651.74972.00652.03292.00911.25572.00722.00722.00652.03291.2557
6–121,044,7790.087551.33421.94012.03091.94011.94731.79831.94821.94822.03093.04171.3342
Where ( * ) MPSA: Minimum potential spatial accessibility; GI ( * * ) : Gini index.
Table 4. Pareto frontier set for the listed k, scenario, zone, and other indicators in instance 2.
Table 4. Pareto frontier set for the listed k, scenario, zone, and other indicators in instance 2.
kScenario ε GI ε acc
(Hospital Beds/
1000 Inhabitants)
Coverage
(Inhabitants)
Healthcare Supply (Hospital Beds by Zone)Total
(Hospital Beds)
12345678910
33–10.7016503.09000305,75103150000040036501080
3–20.6419801.88490429,26700004000033229601028
3–30.5553621.76130522,14600400000395003291124
3–40.5572871.77366522,14600400000400003471147
3–50.4302471.52646637,08200308032440000001032
3–60.4321721.59444637,08200344033940000001083
3–70.4340971.67478637,08200384035540000001139
3–80.4360221.70568637,08200400036540000001165
3–90.4764441.21128657,00800003780025004001028
3–100.4167731.06914751,94540000000025003781028
44–10.6015583.09000390,8980274026400035531701210
4–20.4417961.73658637,0822970002504000003421289
4–30.4437211.86636637,0823450002504000003841379
4–40.4456461.98996637,0823910002504000004001441
4–50.3494041.70568751,9450040003624000250001412
4–60.4052241.21128759,4820000400002502504001300
4–70.3513291.06914854,4193990000002502504001299
55–10.4379472.50290637,082297040003624000002641723
5–20.4398722.65122637,082352040003624000003121826
5–30.4417962.77482637,082400039803614000003541913
5–40.3474791.97142751,945386000250400025404001690
5–50.3494042.00850751,945400000256398026204001716
5–60.2685611.70568854,419004000370400025025001670
5–70.2897341.06914942,833399250000002502504001549
66–10.4571952.71302637,08230802930400400250004002051
6–20.3455542.6265725,49638825033703364000003442055
6–30.3397802.28042751,94528203050323400026202501822
6–40.3417052.53998751,94531203980361400029202772040
6–50.3436302.5647751,94537903170328400029503362055
6–60.2550872.0085854,41940000025739802632504001968
6–70.1954171.70568942,83302504000370400025025001920
6–80.2185151.069141,027,98040025002500002552504001805
77–10.3628782.33604722,2292500250250353397250003052055
7–20.3648032.36694722,2292500250250360320250003752055
7–30.3628782.32986725,4962502502500354400250003012055
7–40.3648032.36076725,4962502502500359333250003632055
7–50.3648032.31132739,5562500250034836125002713252055
7–60.3667282.36076739,5562500250035930625002563842055
7–70.3686522.37312739,5562500250036129425002504002055
7–80.3667282.29278751,9452500250034531625029003542055
7–90.3686522.33604751,9452500250035427725027404002055
7–100.3705772.34222751,9452500250035528025027004002055
7–110.2454632.26188854,4193180250029939602602502822055
7–120.1780931.86636942,8333492500025040002502503062055
7–130.1395971.557361,027,9800250322250333400025025002055
Where GI: Gini index.
Table 5. Objective values in the listed scenario, zone, and other indicators for instance 2.
Table 5. Objective values in the listed scenario, zone, and other indicators for instance 2.
ScenarioCoverageGI ( * * ) Potential Spatial Accessibility
(Hospital Beds by Zone/1000 Inhabitants)
MPSA ( * )
12345678910
3–1305,7510.7015203.5628000003.48243.561903.4824
3–2429,2670.641971.88740001.8874002.89042.888501.8874
3–3522,1460.555351.763702.7600002.76001.7637001.76241.7624
3–4522,1460.557191.776302.8271002.82711.7763001.81691.7763
3–5637,0820.430192.118701.530701.52881.53071.5307001.53071.5288
3–6637,0820.432122.258401.599701.59961.59971.5997001.59971.5996
3–7637,0820.434082.410501.676301.67511.67631.6763001.67631.6751
3–8637,0820.435812.488301.706901.72231.70691.7069001.70691.7069
3–9657,0080.475871.783601.211301.78361.211302.1765101.21131.2113
3–10751,9450.415311.073602.218301.07361.14471.07362.1765101.14471.0736
4–1390,8980.6002903.099103.10050003.090643.0934703.0906
4–2637,0820.441711.976802.773701.97681.97651.7380001.97651.7380
4–3637,0820.443722.105603.029702.10562.10371.8668002.10371.8668
4–4637,0820.445612.229103.201602.22912.15221.9903002.15221.9903
4–5751,9450.349172.474201.706901.70811.70691.70692.176501.70691.7069
4–6759,4820.404641.887401.211301.88741.211302.17652.43961.21131.2113
4–7854,4190.351221.070902.282301.07091.21131.07092.17652.43961.21131.0709
5–1637,0820.437933.271403.303602.50532.50642.5041002.50642.5041
5–2637,0820.439873.419003.596502.65292.65172.6517002.65172.6517
5–3637,0820.441433.539303.848702.77702.77512.7767002.77512.7751
5–4751,9450.347472.215703.188202.21572.15221.97692.211302.15221.9769
5–5751,9450.349312.281603.221102.28162.14752.00982.281002.14752.0098
5–6854,4190.268442.511901.706901.74591.70691.70692.17652.43961.70691.7069
5–7942,8330.289051.07092.82762.282301.07091.21131.07092.17652.43961.21131.0709
6–1637,0820.456803.906604.171402.71413.34472.9601002.71332.7133
6–2725,4960.345533.27232.82763.669402.62682.62802.6277002.62802.6268
6–3751,9450.339772.865103.039002.28102.28212.28192.281002.28212.2810
6–4751,9450.341683.303103.379402.54082.54192.54052.542202.54192.5405
6–5751,9450.343613.172103.582702.56492.56552.56522.568302.56552.5649
6–6854,4190.255062.286303.221102.28632.14752.00982.28972.43962.14752.0098
6–7942,8330.195342.51192.82761.706901.74591.70691.70692.17652.43961.70691.7069
6–81,027,9800.218501.07362.82762.28502.93611.07361.21131.07362.22002.43961.21131.0736
7–1722,2290.362623.446803.63862.93612.33672.96762.7150002.33622.3362
7–2722,2290.364773.479903.66952.93612.36972.99852.5339002.36712.3671
7–3725,4960.362833.45162.82763.633602.34142.96262.7220002.33122.3312
7–4725,4960.364753.47522.82763.663702.36502.99272.5645002.36132.3613
7–5739,5560.364793.423303.614502.31312.94352.630302.64462.31212.3121
7–6739,5560.366703.475203.663802.36502.99282.500902.49822.36142.3614
7–7739,5560.367593.484603.684102.37443.01302.472702.43962.38162.3744
7–8751,9450.366713.409103.596502.29892.92552.52452.524702.29412.2941
7–9751,9450.368633.451603.644102.34142.97302.43272.385502.34172.3414
7–10751,9450.369073.456303.651102.34612.98012.43982.350602.34872.3461
7–11854,4190.245372.743203.117702.26442.26422.26382.26362.43962.26422.2636
7–12942,8330.177982.11642.82762.804302.11641.86751.87762.17652.43961.86751.8675
7–131,027,9800.138472.18802.82761.55752.93611.57131.55751.55752.17652.43961.55751.5575
Where ( * ) MPSA: Minimum potential spatial accessibility; GI ( * * ) : Gini index.
Table 6. Computational study results for the indicated n and k.
Table 6. Computational study results for the indicated n and k.
nkRuntime (in Hours)
105 9.33
6 7.25
7 8.30
8 8.01
9 9.20
155 11.15
6 13.14
7 12.42
8 11.23
9 12.31
205 19.76
6 18.25
7 21.24
8 22.28
9 23.77
305 32.21
6 33.02
7 36.12
8 36.53
9 37.10
505>48
6>48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Delgado, E.J.; Cabezas, X.; Martin-Barreiro, C.; Leiva, V.; Rojas, F. An Equity-Based Optimization Model to Solve the Location Problem for Healthcare Centers Applied to Hospital Beds and COVID-19 Vaccination. Mathematics 2022, 10, 1825. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111825

AMA Style

Delgado EJ, Cabezas X, Martin-Barreiro C, Leiva V, Rojas F. An Equity-Based Optimization Model to Solve the Location Problem for Healthcare Centers Applied to Hospital Beds and COVID-19 Vaccination. Mathematics. 2022; 10(11):1825. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111825

Chicago/Turabian Style

Delgado, Erwin J., Xavier Cabezas, Carlos Martin-Barreiro, Víctor Leiva, and Fernando Rojas. 2022. "An Equity-Based Optimization Model to Solve the Location Problem for Healthcare Centers Applied to Hospital Beds and COVID-19 Vaccination" Mathematics 10, no. 11: 1825. https://0-doi-org.brum.beds.ac.uk/10.3390/math10111825

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop