A multi-criteria spatial deprivation index to support health inequality analyses
- Pablo Cabrera-Barona^{1}Email author,
- Thomas Murphy^{1},
- Stefan Kienberger^{1} and
- Thomas Blaschke^{1}
https://doi.org/10.1186/s12942-015-0004-x
© Cabrera Barona et al.; licensee BioMed Central. 2015
Received: 12 December 2014
Accepted: 19 February 2015
Published: 20 March 2015
Abstract
Background
Deprivation indices are useful measures to analyze health inequalities. There are several methods to construct these indices, however, few studies have used Geographic Information Systems (GIS) and Multi-Criteria methods to construct a deprivation index. Therefore, this study applies Multi-Criteria Evaluation to calculate weights for the indicators that make up the deprivation index and a GIS-based fuzzy approach to create different scenarios of this index is also implemented.
Methods
The Analytical Hierarchy Process (AHP) is used to obtain the weights for the indicators of the index. The Ordered Weighted Averaging (OWA) method using linguistic quantifiers is applied in order to create different deprivation scenarios. Geographically Weighted Regression (GWR) and a Moran’s I analysis are employed to explore spatial relationships between the different deprivation measures and two health factors: the distance to health services and the percentage of people that have never had a live birth. This last indicator was considered as the dependent variable in the GWR. The case study is Quito City, in Ecuador.
Results
The AHP-based deprivation index show medium and high levels of deprivation (0,511 to 1,000) in specific zones of the study area, even though most of the study area has low values of deprivation. OWA results show deprivation scenarios that can be evaluated considering the different attitudes of decision makers. GWR results indicate that the deprivation index and its OWA scenarios can be considered as local estimators for health related phenomena. Moran’s I calculations demonstrate that several deprivation scenarios, in combination with the ‘distance to health services’ factor, could be explanatory variables to predict the percentage of people that have never had a live birth.
Conclusions
The AHP-based deprivation index and the OWA deprivation scenarios developed in this study are Multi-Criteria instruments that can support the identification of highly deprived zones and can support health inequalities analysis in combination with different health factors. The methodology described in this study can be applied in other regions of the world to develop spatial deprivation indices based on Multi-Criteria analysis.
Keywords
Deprivation Analytical Hierarchy Process (AHP) Ordered Weighted Averaging (OWA) Geographically Weighted Regression (GWR) HealthResumen
Antecedentes
Índices de privación son medidas útiles para analizar inequidades en salud. Existen varios métodos para construir estos índices, sin embargo pocos estudios han usado Sistemas de Información Geográfica (SIG) y métodos Multi-Criterio para esta construcción. Este estudio aplica Evaluación Multi-Criterio para calcular los pesos de los indicadores del índice de privación, y también un enfoque SIG de lógica difusa para crear distintos escenarios de este índice.
Métodos
El Proceso Analítico Jerárquico (AHP) es usado para obtener los pesos de los indicadores del índice. La Sumatoria Lineal Ordenada Ponderada (OWA) que usa cuantificadores lingüísticos es aplicada para crear diferentes escenarios de privación. La Regresión Ponderada Geográficamente (GWR) y el índice Moran’s I son empleados para explorar relaciones espaciales del índice de privación y sus escenarios, con dos factores relacionados a salud: distancia a servicios de salud y porcentaje de personas que nunca han tenido un nacido vivo. Este último indicador fue considerado como la variable dependiente de la GWR. El caso de estudio es la Ciudad de Quito, en Ecuador.
Resultados
El índice basado en el método AHP muestra media y alta privación (0,511 a 1,000) en zonas específicas del área de estudio, no obstante, la mayoría del área de estudio tiene bajos niveles de privación. Los resultados de OWA muestran escenarios de privación que pueden ser evaluados considerando diferentes actitudes de los tomadores de decisión. Los resultados de GWR indican que el índice de privación y sus escenarios OWA pueden ser considerados como estimadores locales de fenómenos relacionados a la salud. Los cálculos de Moran’s I demuestran que varios escenarios de privación, en combinación con el factor de ‘distancia a servicios de salud’, podrían ser variables explicativas del porcentaje de personas que nunca han tenido un nacido vivo.
Conclusiones
El índice basado en el método AHP y los escenarios OWA de privación son instrumentos de análisis Multi-Criterio que pueden apoyar a la identificación de zonas con pobreza, y en combinación con otros factores de salud, pueden apoyar al análisis de inequidades en salud. La metodología descrita puede ser aplicada en otras regiones del mundo para desarrollar índices de privación basados en análisis Multi-Criterio.
Palabras clave
Privación Proceso Analítico Jerárquico Sumatoria Lineal Ordenada Ponderada Regresión Ponderada Geográficamente SaludBackground
Approaches to developing deprivation indices are diverse [1-4], and area-based deprivation indices have been proven to be useful in identifying patterns of inequalities in health outcomes [1-11]. Deprivation can be defined as any disadvantage of an individual or human group, related to the community or society to which the individual or human group belongs, and these disadvantages can be of social or material nature [4,5]. Social deprivation can be linked to concepts of social fragmentation [11], and material deprivation can be related to the concept of poverty in terms of the lack of basic goods. These two kinds of deprivation are closely linked to public health and wellbeing [12]. Measuring deprivation requires the identification of two main issues: which indicators to be used to construct a deprivation index, and how to combine these indicators. The criteria for choosing the different indicators that compose deprivation indices can vary. In general, they depend on the availability of information in census and the objective of the study [2-4,8,9]. There are referential studies on constructing multiple deprivation indices, such as the Townsend Deprivation Index, which uses four indicators of material and social deprivation [4]; the Under Privileged Area score, also known as the Jarman Deprivation score, which considers eight deprivation indicators, and has been used to determine remuneration for physicians in United Kingdom [13,14]. Another known measure is the Carstairs deprivation index [15], which is very similar to the Townsend index but is a Scottish reality-based index. Common indicators for these three indices are overcrowding and unemployment. Townsend and Carstairs indices also include a very specific variable available in the British Census, namely the indicator of “Non car ownership”. More recent efforts have used other kinds of indicators from different domains, including health, housing and vulnerability of the population, for the construction of deprivation indices [1-3,6-9]. However, the most common deprivation domains that can support studies of health are related to occupation, education and household conditions, including overcrowding [3]. Once the indicators for a deprivation index are chosen, the next important step is to define how they are going to be combined. Deprivation indicators can be combined using (i) simple additive techniques, using (ii) weights for each indicator, or using (iii) multivariate techniques [16]. The first technique just adds the deprivation indicators [4,16], the second technique can include expert-based weights [17], and the third technique commonly uses indicators weights created using statistical analysis such as the Principal Component Analysis [2,18].
Deprivation indices are constructed by integrating indicators generally extracted from census areas data [18,19]. In many parts of the world where census data are available, such indices can be geo-referenced using GIS. Subsequently, such geo-referenced data allow further spatial analyses, such as investigating spatial correlations [20], performing accessibility analysis [21], analyzing geographical patterns [22] or studying multiple scale evaluations [10] of deprivation measures.
However, regarding the capacity of deprivation indices to be represented spatially explicit, there has been surprisingly little discussion so far about the spatial perspectives of these indices [10,17,22]. There is also not much documented experience - at least not through systematic comparisons of different scenarios - on how to construct these indices spatially.
Based on this background, this paper shows the development of a deprivation index using techniques from Multi-Criteria decision making [8,17,23] and GIS-based fuzzy methods [17,24]. This methodology will show how an Analytical Hierarchy Process (AHP) is applied to obtain the weights for the different indicators that make up the deprivation index. AHP is a Multi-Criteria evaluation method that takes information from experts’ judgments [23]. We then apply Ordered Weighted Averaging (OWA) in order to create different deprivation scenarios [17]. The indicators used to construct our spatial deprivation index follow a rights-based approach [25,26], and are extracted from the 2010 Ecuadorian Population and Housing Census. This rights-based approach prioritizes latent problems in Latin America, where basic needs problems (for example, not having sewerage systems) are more common than, for example, in European countries. The indicators used represent education, health, employment and housing conditions in census blocks of our study area, Quito City, Ecuador. This area has a total of 4034 census blocks, and the census block is considered the smallest area from which census information could be extracted.
A spatial explorative analysis using Geographically Weighted Regression (GWR) and Moran’s I is applied to the deprivation index and its scenarios to evaluate how they are spatially related to the following health factors: distance to health services, and the percentage of people that have never had a live birth. The distance to health services is considered a variable of health accessibility that could be considered in relation to deprivation measures in order to identify its effects on health [21]. The health factor of the percentage of people that have never had a live birth is related to a Population Census variable called “number of people that have never had a live birth”. This indicator can represent health inequalities: when a woman’s child is not born alive, this could be considered to be a consequence of a health condition, such as reproductive or maternal health problem [27]. This indicator can be calculated using information available in the 2010 Ecuadorian Population and Housing Census, and therefore could be considered a useful health-related indicator that can be analyzed together with deprivation indices to be obtained from future Census data. This variable is obtained from women’s answers about how many live births they have had. At the time of a child’s birth, he or she is considered to be a “live birth” if he or she shows vital life signals such as breathing and movement.
Methods
Study area and materials
The information to construct the deprivation index was derived from the 2010 Ecuadorian Population and Housing Census [28]. The advantages of using Population and Housing Census information to construct indices are that census data are commonly open access, and follow a standardization that allows a comparison of information between different places and time. A geocoded shape file of census blocks was also used in order to link the 2010 Census data to the 4034 census areas that make up the study area.
For the calculation of the distance to health services, a data set of the geo-referenced health services in Quito City was used. This data set was provided by Ecuador’s Ministry of Health.
Multi-Criteria Evaluation
A Multi-Criteria Evaluation (MCE) includes knowledge derived from different resources that can be integrated with GIS methods in order to support different kinds of analyses [23]. MCE combines information obtained from different criteria to produce an evaluation index [31] and a weight is allocated to each criterion, to represent the importance of the criterion. In this study, the Analytical Hierarchy Process (AHP) was applied. AHP is a MCE method developed by Saaty [32] that offers practical support for decision making and a straightforward way to obtain weights from criteria [23]. MCE methods, including AHP, also have the capacity to be integrated into GIS-based environments [23,24,33-36] and these GIS-based MCE approaches have been widely and successfully applied in environmental analysis [23,24,31,33,35,37].
The first step in an AHP is to identify a set of criteria. For the deprivation index developed in this study, the criteria are factors or variables that are considered to determine deprivation. A rights-based approach was used to choose the factors that make up the deprivation index [25,26], taking into account the framework of Buen Vivir (Good Living), that is based on human rights and nature rights. The Buen Vivir concept considers that in order to achieve a better quality of life, including time for leisure and harmony with nature, basic needs should first be satisfied [26]. Buen Vivir cannot be achieved if people do not have access to services that ensure their wellbeing and allow them to develop capabilities that create equal opportunities for everyone [26]. To have a good education, health, and to live in conditions of dignity, encourages actions that allow people to construct cohesive societies of Good Living. Human rights are universal. Therefore, the Buen Vivir concept can be applied in other countries, and it is not a concept which focuses only on Ecuador.
Criteria to construct the deprivation index
Dimension | Variable | Indicator/Criteria |
---|---|---|
Education | Level of instruction | A: % of people without any level of instruction |
Health | Health insurance | B: % of people with no health insurance |
Disabled people | C: % of people that area disable for more than a year | |
Employment | Workers with no payment | D: % of people that work with no payment |
Housing | House structure | E: % of houses with damaged roofs |
Overcrowding | F: % of houses with more than four people per dormitory | |
House services | G: % of houses without sewerage system | |
H: % of houses without connection to the electric grid | ||
I: % of houses without connection to the drinking water supply system | ||
J: % of houses without cobbled or paved access roads | ||
K: % of houses without garbage collection service | ||
House general condition | L: % of houses that are shacks |
The indicators used to construct the index represent socio-economic problems: people with no education and people that work for no payment. People who are physically disabled for over a year will be limited in their normal work and daily activities, and those without insurance will be extremely disadvantaged when it comes to health care services. The housing indicators used represent limitations of access to services and a lack of quality of life in the households.
Variances Inflation Factors (VIF) were calculated for all indicators used in order to identify multi-colinearities. VIF shows how much the variance of an estimated regression coefficient is increased as a result of the colinearities between two variables. All VIF obtained were less than 5, which means all selected indicators can be used for the construction of the index.
The key step in any AHP is the creation of a pairwise comparison matrix to compute weights for each criterion while reducing the complexity of the phenomenon in question, because only two criteria are compared at one time [38]. For the comparison of the resulting pairwise matrix, a unified scale is used. The grade of importance of each indicator is evaluated in relation to all other indicators. The importance scale ranges from 1 to 9, whereby 1 means equal importance, 3 means moderate importance, 5 means strong or essential importance, 7 represents very strong importance and 9 indicates extreme importance. Values of 2, 4, 6 and 8 can also be used and are considered as intermediate values.
Results of the AHP method
Indicator (Criteria) | A | B | C | D | E | F | G | H | I | J | K | L | Weights w _{ j } |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
A | 1 | 0,0757 | |||||||||||
B | 1 | 1 | 0,0757 | ||||||||||
C | 1/2 | 1/2 | 1 | 0,0408 | |||||||||
D | 2 | 2 | 3 | 1 | 0,1410 | ||||||||
E | 1 | 1 | 2 | 1/2 | 1 | 0,0757 | |||||||
F | 1 | 1 | 2 | 1/2 | 1 | 1 | 0,0757 | ||||||
G | 2 | 2 | 3 | 1 | 2 | 2 | 1 | 0,1410 | |||||
H | 1 | 1 | 2 | 1/2 | 1 | 1 | 1/2 | 1 | 0,0757 | ||||
I | 2 | 2 | 3 | 1 | 2 | 2 | 1 | 2 | 1 | 0,1410 | |||
J | 1/2 | 1/2 | 1 | 1/3 | 1/2 | 1/2 | 1/3 | 1/2 | 1/3 | 1 | 0,0408 | ||
K | 1 | 1 | 2 | 1/2 | 1 | 1 | 1/2 | 1 | 1/2 | 2 | 1 | 0,0757 | |
L | 1/2 | 1/2 | 1 | 1/3 | 1/2 | 1/2 | 1/3 | 1/2 | 1/3 | 1 | 1/2 | 1 | 0,0408 |
The normalized value for each cell of N is obtained by calculating the ratio of each importance value α _{ ij } of the pairwise comparison matrix and the values sum of each column of this matrix.
Afterwards, all the row values of the normalized matrix are added, and then the sum is divided by the number of the indicators used to construct the deprivation index. The result of this operation is a vector that contains the weights for each indicator (criterion), the eigenvector.
Where n represents the number of criteria and λ _{ max } is obtained as follows: a second vector is obtained by multiplying the eigenvector and the pairwise comparison matrix. Then a third vector is obtained by dividing the values of the second vector by the values of the eigenvector. λ _{ max } is the average of all the components of this final vector [39].
Random indices
n | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
RI | 0 | 0 | 0,58 | 0,90 | 1,12 | 1,24 | 1,32 | 1,41 | 1,45 | 1,49 | 1,51 | 1,48 | 1,56 | 1,57 | 1,59 |
The CR obtained was 0,0019, a value lees than 0,10. This value means that the pairwise comparison matrix is satisfactory [39], which is to say that there is a reasonable level of consistency in the experts’ judgments [38,40]. The weights obtained for each indicator and the CR are also showed in Table 2.
A first representation of the deprivation index was calculated based on the AHP weights by adding the deprivation weighted indicators. Linear min-max normalization was applied to this deprivation index. Values closer to 1 represent higher deprivation. We call the result of this calculation the AHP-based deprivation index.
Ordered Weighted Averaging (OWA)
The Ordered Weighted Averaging (OWA) provides an extension of the Boolean and weighted aggregation operations [39,41]. It ranks the criteria in a MCE and addresses the uncertainty from criteria interaction [24]. OWA works not only with criteria weights (w _{ j }. j = 1,2,3,…, n) but principally with order weights (v _{ j }. j = 1,2,3,…, n). Criteria weights are assigned to each criterion and indicate the level of importance of each criterion [42]. We applied AHP to calculate criteria weights. On the other hand, order weights depend on the ranking of each criterion rather than on its attributes. Order weights are assigned differentially in each location, depending on the respective criterion rank order [43]:
For example, if v _{1}, v _{2} and v _{3} are order weights that have to be applied to the AHP-based weighted criteria X, Y and Z, for instance, if at one location the rank order is YXZ and in another location it is ZYX, the order weights are assigned as v _{1} * Y + v _{2} * X + v _{3} * Z and v _{1} * Z + v _{2} * Y + v _{3} * X, respectively.
Where u _{ j } is the criterion weight reordered according to each criterion attribute value, v _{ j } is the order weight and Z _{ ij } is the sequence obtained by reordering the attribute values. When using different order weights, different results can be produced. From a GIS-based perspective, therefore, using different Boolean operations such as union (OR) and intersection (AND), or weighted linear combination [43-46] will result in different spatial patterns.
Where the parameter ∝ changes depending on which linguistic quantifier it belongs to, and can vary from, “at least one” to “all” quantifiers [38,42,44].
Properties of Regular Increasing Monotone (RIM) quantifiers
∝ | Linguistic Quantifier (Q) | OWA weights ( v _{ j } ) | Decision strategies |
---|---|---|---|
→0 | At least one | v _{1}=1; v _{ j }=0, for all other order weights | Extremely optimistic |
0,1 | Few | * | Very optimistic |
0,5 | Some | * | Optimistic |
1 | Half | v _{ j }=1/n for all j | Neutral |
2 | Many | * | Pessimistic |
10 | Most | * | Very pessimistic |
→∞ | All | v _{ n }=1; v _{ j }=0, for all other order weights | Extremely pessimistic |
In the OWA procedure, it is very important to evaluate the decision strategies. These strategies range between extremely optimistic and extremely pessimistic. These strategies are to be interpreted according to the following logic: in the extremely optimistic strategy, the decision maker’s attitude leads to weighting the highest possible outcome value (for this study the outcome value is the value of deprivation). From a probabilistic perspective, an extremely optimistic strategy is a situation in which a probability of 1, the highest probability, is assigned to the highest value at each location [45]. In other words, the highest ordered weight is assigned to the highest value at each location. The linguistic quantifier for the extremely optimistic strategy is “At least one”, and this linguistic quantifier is equivalent to the logic OR (union) [44], meaning that something is true if at least one logic operand is true.
The other extreme is the extremely pessimistic strategy, where the decision maker’s attitude leads to weighting the lowest possible outcome value. From a probabilistic perspective, in this strategy, the probability of 1 is assigned to the lowest value at each location [45]. The linguistic quantifier for the extremely pessimistic strategy is “All”, and this linguistic quantifier is equivalent to the logic AND (intersection) [44], meaning that something is true if all logic operands are true.
The neutral decision strategy represents a full-tradeoff between criteria, where equal order weights are applied to all possible values at each location. When increasing the degree of optimism from the neutral strategy, greater order weights are assigned to the higher criterion values and smaller weights to the lower criterion values.
Illustration of OWA calculation for four criteria values, for the linguistic quantifier 2
j | Criterion values | Criterion weights w _{ j } | Ordered criterion values z _{ ij } | Ordered criterion weights u _{ ij } | \( {\left({\displaystyle \sum_{\boldsymbol{k}=1}^{\boldsymbol{j}}}{\boldsymbol{u}}_{\boldsymbol{k}}\right)}^{\propto } \) | v _{ j } | v _{ j } * z _{ ij } |
---|---|---|---|---|---|---|---|
1 | 0,20 | 0,30 | 0,80 | 0,35 | (0,35)^{2} = 0,1225 | (0,1225 - 0) = 0,1225 | 0,098 |
2 | 0,80 | 0,35 | 0,50 | 0,10 | (0,45)^{2} = 0,2025 | (0,2025 -0,1225) = 0,08 | 0,04 |
3 | 0,50 | 0,10 | 0,30 | 0,25 | (0,70)^{2} = 0,49 | (0,49-0,2025) = 0,2875 | 0,08625 |
4 | 0,30 | 0,25 | 0,20 | 0,30 | (1)^{2} = 1 | (1–0,49) = 0,51 | 0,102 |
∑ | 1 | 1 | 1 | OWA _{ i }=0,33 |
The process described in our illustration was applied to all the 4034 census blocks of our study area, for all the 12 chosen indicators, for each one of the seven quantifiers: At least one (Extremely optimistic), Few (Very optimistic), Some (Optimistic), Half (Neutral), Many (Pessimistic), Most (Very pessimistic) and All (Extremely pessimistic).
We applied the tool to compute the OWA for the 4034 census blocks. While the tool is running, the linguistic quantifier selected is translated to a numeric parameter ∝ where the values assigned are: 0,001 for At least one, 0,1 for Few, 0,2 for Some, 1 for Half, 2 for Many, 10 for Most and 1000 for All [44]. Using the tool, the OWA was computed for all seven decision strategies, yielding seven scenarios for the deprivation index. These seven scenarios were normalized on a scale from 0 to 1 using linear min-max normalization.
Spatial relationships between the different deprivation measures and health factors
Two health indicators were chosen to evaluate the relation of the OWA deprivation scenarios and the AHP-based deprivation index with health: the distance of each census block to the nearest health service and the percentage of people in each census block that have never had a live birth. These two health indicators represent two different direct measures of the health dimension: a spatial measure of distances and a social measure of a health outcome. For the indicator of distance to health services, first, the centroids for each of the 4034 census blocks were calculated. Sizes of the census blocks differed all over the study area (sizes from around 3200 square meters to sizes of more than 400 000 square meters), therefore, centroids are a good representation of each census block location. Then, 128 health services were identified in the study area and the distances from each census block centroid to the nearest health service were calculated.
The indicator of percentage of people that have never had a live birth was calculated for each census block, using information available in the 2010 Ecuadorian Population and Housing Census: the “number of people that have never had a live birth” and the population of each census block.
This equations means that at every location s, all coefficients β _{ k } need to be estimated, and ε(s _{ i }) is a random error with a mean of zero and a constant variance [50].
Where w _{ ij } is the weight of location s _{ j } that is used to estimate a parameter β _{ k } at the location s _{ i }, and \( {h}_{ij}^{\alpha } \) is the distance between observations s _{ j } and s _{ i } [50].
The aim of applying GWR in this study is to explore how the AHP-based deprivation index and its OWA scenarios relate to health factors by determining the spatial correlations of these relationships. The GWR technique is complemented with the application of the Global Moran’s I.
Where n is the number of spatial units to be taken into account, x is a value of a unit, \( \overline{x} \) is the mean of all values across all n units, and w _{ij} is the spatial weight matrix that is a function of the distance that describes the neighborhood of spatial units. A positive Moran’s I indicates the existence of clusters of similar values, while a negative Moran’s I indicates clusters of dissimilar values. Moran’s I closer to 0 indicates weak autocorrelation [52].
Results and Discussion
Deprivation index and its OWA scenarios
The deprivation scenario with the linguistic quantifier “All” (logic AND) is considered the “worst-case scenario” [44] and in the case of our study, the worst-case scenario means that no action needs to be taken regarding the socio-economic deprivation in almost the entire Quito City territory. Nevertheless, this scenario could be useful to detect the most deprived areas, and can discern areas where taking immediate action is required to reduce socio-economic deprivation. On the other hand, the deprivation scenario with the linguistic quantifier “At least one” (logic OR), representing extreme optimism, shows larger deprivation areas. With this strategy, a larger number of deprivation areas should be considered for socio-economic recuperation, but this may not be feasible for decision makers due to time- and financial constraints. The “Half” scenario means that if the decision makers’ risk-taking is neutral, only the AHP deprivation index constructed based on the experts’ judgments can be considered. The “Few” and “Some” scenarios could support decision making that identifies areas where an extensive social-improvement program could work for most of the city, while the scenarios “Many” and “Most” could support decision making that focuses on taking action in highly socio-economically deprived areas without excessive financial/time investment.
Spatial relationships between the different deprivation measures and health factors
GWR Statistics for all regressions performed
AHP-based | At least one | Few | Some | Half | Many | Most | All | |
---|---|---|---|---|---|---|---|---|
AIC | 18470,02 | 18743,45 | 18634,98 | 18483,08 | 18470,02 | 18505,63 | 18684,81 | 19653,30 |
R ^{ 2 } | 0,59 | 0,50 | 0,53 | 0,58 | 0,59 | 0,59 | 0,57 | 0,36 |
Moran’s I statistics for the residuals of all regression performed
AHP-based | At least one | Few | Some | Half | Many | Most | All | |
---|---|---|---|---|---|---|---|---|
Moran’s Index | −0,005 | 0,051 | 0,035 | −0,000 | −0,006 | −0,008 | 0,006 | 0,215 |
z-score | −1,083 (Random) | 10,198 (Clustered) | 7,130 (Clustered) | −0,049 (Random) | −1,076 (Random) | −1,478 (Random) | 1,175 (Random) | 42,908 (Clustered) |
p-value | 0,279 | 0,000 | 0,000 | 0,961 | 0,282 | 0,139 | 0,239 | 0,0000 |
Conclusion
Our AHP-based deprivation index is a multidimensional index that considers a rights-based conceptual approach useful to representing deprivation in dimensions of education, health, employment and housing. We conclude that our deprivation index has the potential to explain the socio-economic deprivation in the study area accurately because i) the important rights-based indicators used, ii) the consistency of experts’ opinions in the AHP method, and because iii) the several alternative deprivation scenarios allow decision makers to identify urgent zones that can be addressed efficiently and also to the identification of a broader spectrum of zones that can be addressed using more resources. These OWA deprivation scenarios can be considered useful tools for decision makers and health planners. The different decision strategies offer different options when dealing with socio-economic deprivation in the study area. If decision makers decide not to use the AHP-based deprivation index, they can opt for a variety of tradeoff deprivation scenarios (“Few”, “Some”, “Many” and “Most”) that can guide them to where their work will yield better results by saving time and financial resources. The “All” scenario is also interesting when it comes to identifying very deprived zones. These zones represent bigger gaps in the quality of life, and people living there should be considered a priority by health planners and city authorities. The GWR models show that the deprivation index and its scenarios can be related to health factors, and that several deprivation scenarios in combination with the ‘distance to health services’ factor, could be considered explanatories variables to predict the percentage of people that have never had a live birth.
One limitation of this study is that no analysis of uncertainty was elaborated for the OWA scenarios. Even though this is not an objective of this article, we consider that a future study can incorporate uncertainty analysis for different OWA deprivation scenarios. Another limitation is that this study does not develop a complete statistical deprivation-health factor model. We reiterate that the GWR and Moran’s I analyses should only be seen as an exploratory analysis, and more research regarding this issue is needed. Future research could include the incorporation of more explanatory health variables that could interact with the AHP-based deprivation index and the OWA deprivation scenarios. The identification of additional health problems that can be explained to some degree by the methods implemented in this study is also important. Further work can include variations of the Multi-Criteria evaluation used, for example, the use of different techniques to obtain criteria weights and order weights for the deprivation indicators.
This study has several strengths, and can be considered as one of the first instances where Multi-Criteria evaluation methods such as AHP and OWA are utilized to create a deprivation index and deprivation scenarios. A strength of this study is that the AHP-OWA approach captures quantitative and qualitative information to produce different scenarios that are useful for decision makers when faced with different decision strategies due to constraints in time and financial resources. A further value of our study is that the OWA method is spatial in the sense that it aggregates the criteria for each census block depending on their values, and this aggregation is done for all linguistic quantifiers. Another strength of this work is the fact that the OWA procedure was automated with the development of the Python toolbox, which allows more efficient calculation of OWA deprivation scenarios for future studies.
The methodology described in this study can be applied in other regions of the world to develop spatial deprivation indices based on Multi-Criteria analysis.
An important contribution of this study is that the mixed method of applying AHP to calculate deprivation criteria weights and OWA to create different deprivation scenarios is a methodology that can be carried out in other studies beyond Latin America. The indicators considered in this study are common Population and Housing Census variables. However, as AHP and OWA methods are techniques that can be adapted to specific problems and phenomena, future studies can use the methodology presented here considering different deprivation indicators. Furthermore, the methods and results showed in this study can be considered as important tools to support health planners and decision makers.
Declarations
Acknowledgements
The presented work has been funded by the Government of Ecuador through the Ecuadorian Secretary of Higher Education, Science, Technology and Innovation (SENESCYT) and the Ecuadorian Institute of Educational and Scholarship Credits (IECE) (Scholarship contract No. 375–2012). It has also partially been funded by the Austrian Science Fund (FWF) through the Doctoral College GIScience (DK W 1237 N23) at the University of Salzburg.
Authors’ Affiliations
References
- Niggebrugge A, Haynes R, Jones A, Lovett A, Harvey I. The index of multiple deprivation 2000 access domain: a useful indicator for public health? Soc Sci Med. 2005;60:2743–53.View ArticlePubMedGoogle Scholar
- Pampalon P, Pamel D, Gamache P, Raymond G. A deprivation index for health planning in Canada. Chronic Dis Canada. 2009;29(4):178–91.Google Scholar
- Pasetto R, Sampaolo L, Pirastu R. Measures of material and social circumstances to adjust for deprivation in small-area studies of environment and health: review and perspectives. Ann Ist Super Sanita. 2010;46(2):185–97.PubMedGoogle Scholar
- Townsend P. Deprivation. J Soc Policy. 1987;16:125–46.View ArticleGoogle Scholar
- Testi A, Ivaldi E. Material versus social deprivation and health: a case study of an urban area. Eur J Health Econ. 2009;10:323–8.View ArticlePubMedGoogle Scholar
- Adams J, Ryan V, White M. How accurate are Townsend deprivation scores as predictors of self-reported health? A comparison with individual level data. J Public Health. 2004;27(1):101–6.View ArticleGoogle Scholar
- Boyle P, Gatrell A, Duke-Williams O. Do area-level population change, deprivation and variations in deprivation affect individual-level self-reported limiting long-term illness? Soc Sci Med. 2001;53:795–9.View ArticlePubMedGoogle Scholar
- Cabrera Barona P. A multiple deprivation index and its relation to health services accessibility in a rural area of Ecuador. In: Vogler R, Car A, Strobl J, Griesebner G, editors. GI_Forum 2014. Geospatial Innovation for Society: 1–4 July 2014; Salzburg. Wichmann/OAW; 2014. p. 188-91.Google Scholar
- Havard S, Deguen S, Bodin J, Louis K, Laurent O, Bard D. A small-area index of socioeconomic deprivation to capture health inequalities in France. Soc Sci Med. 2008;67:2007–16.View ArticlePubMedGoogle Scholar
- Schuurman N, Bell N, Dunn JR, Oliver L. Deprivation indices, population health and geography: an evaluation of the spatial effectiveness of indices at multiple scales. J Urban Healt. 2007;84(4):591–603.View ArticleGoogle Scholar
- Stjärne MK, Leon A, Hallqvist J. Contextual effects of social fragmentation and material deprivation on risk of myocardial infarction—results from the Stockholm Heart Epidemiology Program (SHEEP). Int J Epidemiol. 2004;33:732–41.View ArticlePubMedGoogle Scholar
- Pampalon R, Raymond G. A deprivation index for health and welfare planning in Quebec. Chronic Dis Canada. 2000;21(3):104–13.Google Scholar
- Jarman B. Identification of underprivileged areas. Br Med J. 1983;286:1705–9.View ArticleGoogle Scholar
- Jarman B. Underprivileged areas: validation and distribution of scores. Br Med J. 1984;289:1587–92.View ArticleGoogle Scholar
- Carstairs V, Morris R. Deprivation: explaining differences in mortality between Scotland and England and Wales. Br Med J. 1989;299:886–9.View ArticleGoogle Scholar
- Folwell K. Single measures of deprivation. J Epidemiol Community Health. 1995;49(2):S51–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Bell N, Schuurman N, Hayes M. Using GIS-based methods of multicriteria analysis to construct socio-economic deprivation indices. Int J Health Geogr. 2007;6(17):1–19.Google Scholar
- Lalloué B, Monnez JM, Padilla C, Kihal W, Le Meur N, Zmirou-Navier D, et al. A statistical procedure to create a neighborhood socioeconomic index for health inequalities analysis. Int J Equity Health. 2013;12(21):1–11.Google Scholar
- Messer LC, Laraia BA, Kaufman JS, Eyster J, Holzman C, Culhane J, et al. The development of a standardized neighborhood deprivation index. J Urban Healt. 2006;83(6):1041–62.View ArticleGoogle Scholar
- Hogan JW, Tchernis R. Bayesian factor analysis for spatially correlated data, with application to summarizing area-level material deprivation from census data. J Am Stat Assoc. 2004;99(466):314–24.View ArticleGoogle Scholar
- Jordan H, Roderick P, Martin D. The index of multiple deprivation 2000 and accessibility effects on health. J Epidemiol Community Health. 2004;58:250–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Benach J, Yasui Y. Geographical patterns of excess mortality in Spain explained by two indices of deprivation. J Epidemiol Community Health. 1999;53:423–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Feizizadeh B, Blaschke T. Landslide risk assessment based on GIS multi-criteria evaluation: a case study in Bostan-Abad County, Iran. J Earth Sci Eng. 2011;1:66–71.Google Scholar
- Feizizadeh B, Blaschke T, Nazmfar H. GIS-based ordered weighted averaging and Dempster–Shafer methods for landslide susceptibility mapping in the Urmia Lake Basin. International Journal of Digital Earth: Iran; 2012.View ArticleGoogle Scholar
- Mideros A. Ecuador: defining and measuring multidimensional poverty, 2006–2010. Cepal Rev. 2012;108:49–67.Google Scholar
- Ramírez R. La vida (buena) como riqueza de los pueblos. Hacia una socio ecología política del tiempo. Economía e Investigación IAEN; 2012.Google Scholar
- Pan American Health Organization. Salud en las Américas Publicaciones Científicas y Técnicas 636. Capítulo de Ecuador; 2012.Google Scholar
- Instituto Nacional de Estadísticas y Censos. Censo de Población y Vivienda 2010. 2014. http://www.ecuadorencifras.gob.ec/banco-de-informacion/. Accessed 10 Octo. 2014.Google Scholar
- Lozano Castro A. Quito, Ciudad Milenaria, Forma y símbolo. 1st ed. Abya Yala; 1991.Google Scholar
- Carrión F, Erazo Espinosa J. La forma urbana de Quito: una historia de centros y periferias. Bulletin de l‘Institut Français d‘Études Andines. 2012;41(3):503–22.View ArticleGoogle Scholar
- Yu J, Chen Y, Wu J. Cellular automata based spatial multi-criteria land suitability simulation for irrigated agriculture. Int J Geogr Inf Sci. 2011;25(1):131–48.View ArticleGoogle Scholar
- Saaty TL. A scaling method for priorities in hierarchical structure. J Math Psychol. 1977;15(3):34–9.View ArticleGoogle Scholar
- Feizizadeh B, Blaschke T. GIS-multicriteria decision analysis for landslide susceptibility mapping: comparing three methods for the Urmia lake basin, Iran. Nat Hazards. 2013;65:2105–28.View ArticleGoogle Scholar
- Marinoni O. Implementation of the analytical hierarchy process with VBA in ArcGIS. Comput Geosci. 2004;30:637–46.View ArticleGoogle Scholar
- Joerin F, Theriault M, Musy A. Using GIS and outranking multicriteria analysis for land-use suitability assessment. Int J Geogr Inf Sci. 2001;15(2):153–74.View ArticleGoogle Scholar
- Carver S. Integrating multi-criteria evaluation with geographical information systems. Int J Geogr Inf Sci. 1991;5(3):321–39.View ArticleGoogle Scholar
- Ramanathan R. A note on the use of the analytic hierarchy process for environmental impact assessment. J Environ Manag. 2001;63:27–35.View ArticleGoogle Scholar
- Boroushaki S, Malczewski J. Implementing an extension of the analytical hierarchy process using ordered weighted averaging operators with fuzzy quantifiers in ArcGIS. Comput Geosci. 2008;34:399–410.View ArticleGoogle Scholar
- Gómez Delgado M, Barredo Cano JI. Sistemas de Información geográfica y evaluación multicriterio en la ordenación del territorio. RA-MA Editorial. 2005.Google Scholar
- Saaty TL. The analytic hierarchy process: planning. Resource Allocation. McGraw-Hill: Priority Setting; 1980.Google Scholar
- Malczewski J. GIS‐based multicriteria decision analysis: a survey of the literature. Int J Geogr Inf Sci. 2006;20(7):703–26.View ArticleGoogle Scholar
- Meng Y, Malczewski J, Boroushaki S. A GIS-based multicriteria decision analysis approach for mapping accessibility patterns of housing development sites: a case study in Canmore, Alberta. J Geogr Inf Syst. 2011;3:50–61.Google Scholar
- Jian H, Eastman R. Application of fuzzy measures in multi-criteria evaluation in GIS. Int J Geogr Inf Sci. 2000;14(2):173–84.View ArticleGoogle Scholar
- Malczewski J. Ordered weighted averaging with fuzzy quantifiers: GIS-based multicriteria evaluation for land-use suitability analysis. Int J Appl Earth Observation Geoinf. 2006;8:270–7.View ArticleGoogle Scholar
- Malczewski J, Chapman T, Flegel C, Walters D, Shrubsole D, Healy MA. GIS-multicriteria evaluation with Ordered Weighted Averaging (OWA): case study of developing watershed management strategies. Environ Plan A. 2003;35(10):1769–84.View ArticleGoogle Scholar
- Yager RR. Quantifier guided aggregation using OWA operators. Int J Intell Syst. 1996;11:49–73.View ArticleGoogle Scholar
- Liu X, Shilian H. Orness and parameterized RIM quantifier aggregation with OWA operators: a summary. Int J Approx Reason. 2007;48:77–97.View ArticleGoogle Scholar
- Clement F, Orange D, Williams M, Mulley C, Epprecht M. Drivers of afforestation in Northern Vietnam: assessing local variations using geographically weighted regression. Appl Geogr. 2009;29:561–76.View ArticleGoogle Scholar
- Gao J, Li S. Detecting spatially non-stationary and scale-dependent relationships between urban landscape fragmentation and related factors using Geographically Weighted Regression. Appl Geogr. 2011;31:292–302.View ArticleGoogle Scholar
- Krivoruchko K. Spatial regression models: concepts and comparison. In: In spatial statistical data analysis for GIS users. Firstth ed. Redlands, California: Esri Press; 2011. p. 483–537.Google Scholar
- Braz Junior G, Cardoso de Paiva A, Corrêa S. Classification of breast tissues using Moran’s index and Geary’s coefficient as texture signatures and SVM. Comput Biol Med. 2009;39:1063–72.View ArticlePubMedGoogle Scholar
- Cai X, Wang D. Spatial autocorrelation of topographic index in catchments. J Hidrol. 2006;328:581–91.View ArticleGoogle Scholar
- Malczewski J. Integrating multicriteria analysis and geographic information systems: the ordered weighted averaging (OWA) approach. Int J Environ Technol Manage. 2006;6(1/2):7–19.View ArticleGoogle Scholar
- Amiri MJ, Mahiny AS, Hosseini SM, Jalali SG, Ezadkhasty Z, Karami SH. OWA Analysis for ecological capability assessment in watersheds. Int J Environ Res. 2013;7(1):241–54.Google Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.