 Research article
 Open Open Peer Review
 Published:
Does the effect of gender modify the relationship between deprivation and mortality?
BMC Public Healthvolume 12, Article number: 574 (2012)
Abstract
Background
In this study we propose improvements to the method of elaborating deprivation indexes. First, in the selection of the variables, we incorporated a wider range of both objective and subjective measures. Second, in the statistical methodology, we used a distance indicator instead of the standard aggregating method principal component analysis. Third, we propose another methodological improvement, which consists in the use of a more robust statistical method to assess the relationship between deprivation and health responses in ecological regressions.
Methods
We conducted an ecological smallarea analysis based on the residents of the Metropolitan region of Barcelona in the period 1994–2007. Standardized mortality rates, stratified by sex, were studied for four mortality causes: tumor of the bronquial, lung and trachea, diabetes mellitus type II, breast cancer, and prostate cancer. Socioeconomic conditions were summarized using a deprivation index. Sixteen sociodemographic variables available in the Spanish Census of Population and Housing were included. The deprivation index was constructed by aggregating the abovementioned variables using the distance indicator, DP_{2}. For the estimation of the ecological regression we used hierarchical Bayesian models with some improvements.
Results
At greater deprivation, there is an increased risk of dying from diabetes for both sexes and of dying from lung cancer for men. On the other hand, at greater deprivation, there is a decreased risk of dying from breast cancer and lung cancer for women. We did not find a clear relationship in the case of prostate cancer (presenting an increased risk but only in the second quintile of deprivation).
Conclusions
We believe our results were obtained using a more robust methodology. First off, we have built a better index that allows us to directly collect the variability of contextual variables without having to use arbitrary weights. Secondly, we have solved two major problems that are present in spatial ecological regressions, i.e. those that use spatial data and, consequently, perform a spatial adjustment in order to obtain consistent estimators.
Background
Deprivation has been defined by Townsend as ‘a state of observable and demonstrable disadvantage relative to the local community or the wider society or nation to which an individual, family or group belong to’ [1]. It is a concept that has two distinguishable domains, material and social circumstances. The former relates to diet, health, clothing, housing, household facilities, environment and work. The latter is more difficult to measure as it relates to different forms of relationships, such as community integration, leisure, and formal participation in social institutions [1].
A deprivation index (DI) is recognized as a composite measure, where no single variable can be said to measure it but rather a number of variables contribute in some way [2]. Arealevel indicators are commonly used to evaluate the geographical distribution of socioeconomic inequalities in health [3–6]. Arealevel indicators have been used as proxies for individuallevel data when individual measures are not available (i.e. data from census is only provided at aggregated levels for reasons of confidentiality) [7].
Approaches to ecological smallarea analysis propose using the smallest practicable spatial scale for any study on inequalities in health [8, 9]. This diminishes the ecological bias, as the analysis is closer to the individual level [8, 10, 11]. Most studies that have focused on arealevel social factors use geographical boundaries, developed for censuses, as proxies for actual communities or neighborhoods [6]. Administrative boundaries allow for routinely collected data to be available at the smallest scale. In Spain, census tracts, the smallest administrative unit, tend to have a mean population of 1,000 subjects [12].
There are three main methods for constructing Deprivation Indices (DIs) [13, 14]. First, standardized zscores or log transformations of a priori selected variables. This was a very popular method up until the late 1980´s [15, 16]. The second method is the use of Principal Component Analysis (PCA) which has been the principal method used over the last twenty years. PCA transforms a number of possibly correlated variables into uncorrelated variables and permits the extraction of a smaller number of uncorrelated variables called principal components that collect a large percentage of the total variance contained in the original data which are then used to measure socioeconomic status [17–20]. The third method (the least common) uses feedback from health experts to assign weights to the selected variables [21–24].
The increased use of DIs has allowed for the identification of deprived areas and their association to morbidity and mortality, living in deprived areas is an effect modifier. In fact, some studies confirm that people living in deprived areas have higher rates of ill health and mortality [4, 5, 8, 25, 26]. However, material deprivation is manifested in higher rates of mortality differently by gender [27, 28]. This might be due in part to the fact that men and women perceive their environment differently [29–31], they may be exposed to different local environments [32], and their vulnerability to aspects in the local environment may vary [33].
In this study we propose improvements to the method of elaborating DIs. First, in the selection of the variables, we incorporated a wider range of both objective and subjective measures which makes for a more complete DI. Several studies that analyzed deprivationassociated excess mortality found that men had higher rates of excess mortality than women [27, 28]. A large proportion of excess deaths were found in dietrelated causes of death in females and smoking and alcohol consumption in males [27]. However, inferences about gender differences in relation to deprivation and mortality are sensitive to the choice of inequality measures used [34–36]. Making it important to incorporate indicators that are gender sensitive (i.e. residential/ environmental indicators as they are strongly related to women’s health and individual economic activity indicators as they are associated to men’s health) [32]. The strength of association between socioeconomic and health measures tend to be greater for men than for women [35]. Thus, contextual indicators are crucial if we want to have a comparative measure between men and women’s health.
Studies have shown that neighborhood context influences selfrated health, especially in women, beyond the effects of individual factors [33, 37]. Neighborhoods with worse socioeconomic conditions have a negative effect on health [38]. Additionally residential environment may be more important for women’s health while individual economic activity is more important for men’s health [32]. Other studies have found that occupational factors are more important for men’s health whereas the home environment is more important for women’s health [39].
Second, in the statistical methodology, we used a distance indicator (Pena Distance or DP_{2}) instead of the leading aggregating method PCA. DP_{2} overcomes several limitations of PCA, for instance, aggregating variables expressed in different units of measurement, arbitrary weights, and duplicate information [40–45]. The greatest advantage of DP_{2} with respect to PCA consists on how each handles redundant information. Thus, while in PCA the first component, i.e. the DI, is composed of those variables that capture greater variability it leaves some variables out of the index. In DP_{2} the index is composed of all the variables, leaving out of the index only the redundant information but including each and every one of the variables. Thus, in principle, DP_{2} would collect more variability [45].
In addition, we propose another methodological improvement, which consists in the use of a more robust statistical method to assess the relationship between deprivation and health responses in ecological regressions.
In this sense, we focused on mortality from trachea, bronchial and lung cancer and diabetes, because there is certain evidence that gender could act as an effect modifier in the relationship with deprivation, and breast and prostate cancer, given their known association with deprivation. Several studies have shown a positive association between lung cancer mortality and deprivation for men and a negative association for women [18, 28, 46]. This is explained because in men, lower socioeconomic status implies a higher prevalence of smoking, which leads to an increased risk of lung cancer mortality in individuals with lower socioeconomic status. On the other hand, in women, those with higher socioeconomic status are those with a higher prevalence of smoking, leading to an increased risk of dying from that cause in individuals with higher socioeconomic status. Regarding diabetes, those individuals with high educational levels have a significantly lower risk of dying from this cause (odds ratio, OR = 0.66) compared to those with low educational levels [47]. In the Whitehall study [48], we observed a similar gradient for mortality from diabetes mellitus than with the other pathologies studied; higher mortality from diabetes mellitus in individuals who were in the lowest levels of the hierarchy. This excess mortality appeared to be due to the increased frequency of cardiovascular risk in this type of study group. In addition, variations in mortality have also been observed for diabetes mellitus type II (DM2) at the residence area level. Thus, individuals living in more deprived areas tend to have worse lifestyle (higher values of body mass index –BMI, greater prevalence of smoking) which could lead to worse glycemic control and greater number of complications associated with diabetes mellitus. As in the case of lung cancer mortality, this association could be modified by sex.
However, the results from the Spanish literature, at least in those that use a spatial adjustment [18, 28, 46] do not systematically coincide with each other; particularly for prostate cancer and breast cancer. Diabetes mortality has shown a positive association with deprivation for both men [28, 49] and women [49]. Borrell et al.[18] and PuigpinósRiera et al.[46] found a positive association for breast cancer in Alicante (relative risk, RR = 1.55 95% confidence interval, 95% CI:1.042.19) yet the majority of other cities showed a negative association, not statistically significant, except for Vigo (RR = 0.54 95% CI 0.330.84). Benach et al.[27] found a negative association for breast cancer and deprivation with Index 1, which measured deprivation through unemployment, illiteracy and low social class; however, found no association with Index 2, measured deprivation through overcrowding and a small component of unemployment and illiteracy. For prostate cancer, some studies found no clear association between mortality risk and deprivation [46]. Benach et al.[27] found a negative association for prostate cancer and deprivation with Index 1; however, found no association with Index 2.
Our hypothesis is that the differences found in the literature are caused mainly by the lack of robustness of the statistical methodology applied and the lack of robustness of the spatial adjustment applied. Therefore, our objective in this paper is to evaluate a more robust methodology, in the way of developing DIs, the selection of component variables and how to build them, and, above all, in the statistical method applied to analyze how the effect of gender modifies the relationship between deprivation and mortality.
Methods
Design and study population
Most countries in the world periodically carry out population censuses that gather information on sociodemographic characteristics of the population resident in a country. In Spain, the Spanish Census of Population and Housing (CPH) conducted in 2001 reported nearly 41 million inhabitants. The province capital of Barcelona had over 1.5 million inhabitants and an average number of 1,201.33 inhabitants per census tract (standard deviation, SD, 526.66; Median 1,085.50). For this study, we conducted an ecological smallarea analysis based on the residents of the Metropolitan region of Barcelona. A total of 2,978 census tracts were examined.
Variables
Standardized mortality rates (SMR), stratified by sex, were studied for four mortality causes: tumor of the bronquial, lung and trachea (ICD10: C33C34), diabetes mellitus (ICD10:E10E14), breast cancer (ICD10:C50), and prostate cancer (ICD10:C61). The mortality data was provided from the Catalan Mortality Registry. We only used death certificates of residents of the Metropolitan Region of Barcelona who died between 1994 and 2007.
Deprivation index
Socioeconomic conditions were summarized using a DI in a census tract level. Sixteen sociodemographic variables available in the Spanish CPH were included. Four of these: unemployment rate, percent with low educational level, manual workers and temporary workers, have been used in several DIs for smallarea studies in Spain [46, 50]. The other variables (university education, monoparental homes, activity rate, immigrants, homes without heating, no toilet/bathroom in the home, bad communications, vandalism/crime, few green areas, exterior noise, contamination/bad smells, and dirty street) were chosen based on literature that focused on the association between contextual indicators and their effect on health, particularly those that analyzed gender differences [38]. We chose contextual indicators, most of them subjective, which were representative of the neighborhood characteristics (Table 1).
The DI was constructed by aggregating the abovementioned variables using another indicator, DP_{2}, instead of the standard PCA method. DP_{2} is an iterative procedure that weights partial indicators depending on their correlation with the global index [41, 45, 51]. It is able to capture all the nonredundant variance of the indicators (i.e. avoiding multicollinearity); thus, allowing for the inclusion of a greater number of variables [41, 42]. Given that DP_{2} uses all the valuable information contained in the partial indicators the more complete the final DI will be, since each variable contains unique information not present in others [44]. Magnifying the goodness of DP_{2} versus the extraction of an index through principal component analysis, especially when introducing a greater number of variables would not suffice for identifying a single component for the index. More than one component should be collected when applying PCA thus resulting in an arbitrary choice of weights for aggregating the components to obtain the index. Pena [41] and Zarzosa [42] showed that DP_{2} fulfills all the properties of a good composite indicator; monotony, unicity, invariance, homogeneity, transitivity, exhaustivity, existence and determination, and additivity.
The process for aggregating the variables using DP_{2} is a fourstep process previously described in detail [51].
Data analysis
When spatial data is available, the variability in the observed response is usually greater than the expected, resulting in over dispersion. In fact, it is important to distinguish between two sources of extravariability [52]. The first most important source is the socalled 'spatial dependence’; it is a consequence of the correlation between the spatial unit and the neighboring spatial units, generally contiguous geographic areas. Nearby areas are more similar than those far apart. Part of this dependence is not really a structural dependence, but mainly due to variables not included in the analysis. The second source is the independent extra variability, spatially uncorrelated, called heterogeneity (not spatial); it is the result of unobserved variables, without spatial structure, which could influence the response [52, 53].
To account for this extra variability, it is necessary to introduce some structure in the model. Otherwise, unless the model is linear, the estimates of the parameters of interest and the standard errors of the estimators will be biased; therefore any inference based on them, will be wrong [54].
In this paper we have chosen to use hierarchical Bayesian models, more specifically a model based on Besag, York and Mollié [55, 56] (BYM) to analyze the relationship between mortality and smallarea deprivation.
As we know, the idea is to introduce two random effects in a generalized linear model with Poisson response, in order to capture the extravariability [53].
Where O_{i} denotes the observed cases of the response variable in the census tract i, μ _{ i } is the relative risk in section i, Pop_{i} is the population (men or women) of the census tract, υ _{ i } is the random effect which reflects the heterogeneity; S_{i} is the random effect that reflects the spatial dependence, α is the intercept, interpreted as the logarithm of the baseline risk; IndexQ_{ki} is the quintile of the deprivation index (standardized) in the census tract i (the first quintile was taken as a reference); the β are the parameters of the model, which can be interpreted as the logarithms of the relative risks associated with the explanatory variables; Pop4564_{i} is the percentage of men (or women) from 45 to 64 years in the census tract i, and Pob65m_{i} is the percentage of men (or women) aged 65 years or older in the census tract i.
The DI was categorized into quintiles because the influence of deprivation on the (spatial) variation of mortality could be nonlinear.
The nonspatial random effect, also called heterogeneity, assumed to be normally distributed, with a mean of zero and constant variance. For the spatial random effect a conditional autoregressive model (CAR) is used [57, 58]. This approach, which is the most utilized and has the lowest computational cost, approximates the spatial dependence as an average spatial effect of neighboring areas [53]. The areas considered are the census tracts of each municipality in the Metropolitan Region of Barcelona and surrounding areas are defined as adjacent census tracts, in other words they share a border.
Note that the specified model does not use, as an offset the number of cases expected in the census tract, but the population (men or women) of the same. This is because, unlike the standard BYM model [53], here we use the crude death rate (from the census tract) and not the standardized mortality ratio as an indicator of mortality. The reason is to avoid the problem called 'mutual standardization' [59, 60]. Rosenbaun and Rubin [59] show how the use of standardized rates as the response variable in ecological regression models leads to biased results if only the answer, not predictors, are adjusted for the same confounder, usually the age distribution. When the predictor is not adjusted, it is implicitly assumed that the effect of predictor is constant for all strata of the confounding variable. This may be true for the effect of an air pollutant, in principle it is the same for all ages, but not for the deprivation index. Grisotto et al.[60], in line with Rosenbaun and Rubin [59] show that unbiased estimators can be obtained by adjusting the response and the predictors (the index of deprivation in this case) with the same variable (age distribution) or, even easier to use are crude rates as the response variable and entering age (as an average or structured) as an explanatory variable of the model. This is why we have introduced the age structure of the census tract (proportion of men and women aged 45 to 64 years and 65 years or more). The introduction of age also lets you control the age effect in the model. This modification allows to overcome methodological limitations of the standard implementation of the BYM model [18]. Furthermoe, and also getting no significance for the coefficients of the variables that represent the population from 45 to 64 years and older population or more at 65, none of the cases studied, we can say that the risks found not differ by age group. Moreover, because the estimates associated with variables representing the population aged 45 to 64 years and older population less than 65, were not significant in any of the cases studied, we can say that risks not found differ by age group.
A second difference from the standard BYM model, not so obvious, is the use of standardized explanatory variables. The reason is that the spatial random effect, approximated by CAR, can be correlated with some (or all) explanatory variables that have a similar spatial dependence (known as concurvity). If this were the case, there would be an overadjustment, unlike the phenomenon of multicollinearity in linear models, which could bias the estimates [61]. Simulations performed by us suggest that the problem could be solved by completely standardizing the potentially problematic explanatory variable. On the other hand, the introduction of the deprivation index (standardized) into quintiles, in addition to collecting a possible nonlinear effect, could mitigate much of the problem.
Spatial models were built as Bayesian hierarchical models with two stages [62] and estimated using the integrated nested Laplace approximation [62–64] (see Appendix: Annex).
Models were compared using the DIC (Deviance Information Criterion) [65] and the conditional predictive ordinate (CPO) for each observation (in fact –mean(log(cpo)) [66, 67]. CPO is a crossvalidated predictive approach i.e., predictive distributions conditioned on the observed data with a single data point deleted. Asymptotically the CPO statistic has a similar dimensional penalty as DIC. In this perspective, the CPO statistic may be similar to DIC. In both cases, the lower the DIC or the CPO, the better the model.
There was not experimental research in this work. All computations were carried out using the interface INLA [68], running directly in R (version R 2.11.0) [69].
Results and discussion
According to the Spanish CPH of 2001 the average number of inhabitants per census tract in the metropolitan region of Barcelona was 1201.33 (standard deviation 526.55; median 1085.50).
Table 1 describes de distribution of the socioeconomic indicators, used to construct the DI, by gender in the Metropolitan Region of Barcelona in 2001. We observed a high percent of temporary workers (20.9% women, 18.0% men) and manual workers (51.6% women, 56.0% men) for both men and women. Women had a higher rate of insufficient education (17.1% vs. 12%) and a lower activity rate (57.1% vs. 67.3%) than men. Over a third of the individuals claimed their neighborhoods had dirty streets, exterior noise, few green areas and vandalism/crime.
Tables 2 and 3 show the associations between causespecific mortality and deprivation by quintiles of the index, for men and women controlled for age. In the case of women (Table 2), we observed a positive association for diabetes mortality. However, both lung cancer mortality (RR in Q5 = 0.78; 95%CI: 0.650.93) and breast cancer mortality (RR in Q5 = 0.81; 95%CI: 0.710.92) had an inverse relationship with socioeconomic deprivation; women with greater deprivation had less mortality risk (statistical significance was present only in the fifth quintile).
For men (Table 3), we observed a positive association for lung cancer mortality; where, at greater socioeconomic deprivation there was greater mortality risk (RR in Q5 = 1.31; 95%CI: 1.191.43). For diabetes mortality we also observed a positive association with mortality risk and deprivation; however, statistical significance was found only in second and third quintiles. On the other hand, for prostate cancer mortality there is no systematic relationship between deprivation and mortality risk. Only the second quintile is significant with a RR = 1.13 (95% CI: 1.001.26).
Overall, our results are consistent with those provided by the literature. At greater deprivation, there is an increased risk of dying from diabetes for both sexes (although in the case of men the relative risks associated with the third quintile of deprivation onwards were not statistically significant) and of dying from lung cancer for men. On the other hand, at greater deprivation, there is a decreased risk of dying from breast cancer and lung cancer for women (although in both cases only the relative risk associated with the top quintile of deprivation was statistically significant). We did not find a clear relationship in the case of prostate cancer (presenting an increased risk but only in the second quintile of deprivation).
This study has been able to improve the statistical methodology applied in building deprivation indices as well as the robustness of spatial adjustments in ecological studies. Our results were consistent with the existing literature that analyzes the association between deprivation and mortality (lung cancer, diabetes mellitus, breast cancer and prostate cancer); however, our focus was in urban smallareas. As some literature suggests, mortality differs in rural versus urban settings, particularly in women [70]; thus, it would be important to explore geographic differences in more detail in future studies. In addition, future studies should also focus on researching these same causes of death in other European cities and additional causes of death in the Barcelona Metropolitan Region.
Conclusions
We believe our results were obtained using a more robust methodology. First off, we have built a better index that allows us to directly collect the variability of contextual variables without having to use arbitrary weights, as is assumed by the aggregation of components through PCA extraction. Building an index using PCA is adequate, the results do not differ from those obtained by DP_{2}, when you add only a few variables because the first component allows you to collect enough variability. When using a large number of variables, like shown in this article, it is clear that a single component is not enough to collect the variability contained in the variables and the index should thus be constructed by aggregating two or more components. The arbitrariness in the choice of aggregation weights is what we avoided by using the DP_{2} methodology as the algorithm used for constructing the index naturally allows you to collect all the variability contained in the variables excluding only the shared part and not the variable itself.
Secondly, we have solved two major problems that are present in spatial ecological regressions, i.e. those that use spatial data and, consequently, perform a spatial adjustment in order to obtain consistent estimators. In particular, the problem of mutual standardization and the problem of concurvity. As stated previously, the problem of concurvity is present when the spatial random effect is correlated with some (or all) of the explanatory variables that have similar spatial dependence. Using crude rates as indicators of mortality for the response variable and the introduction of age as an explanatory variable in the model allows us to ensure the collection of unbiased estimators. Furthermore, the introduction of age in the model allows us to monitor its effect on the results. On the other hand, the introduction of the index (standardized) into quintiles, in addition to collecting a possible nonlinear effect, could mitigate much of the problem of concurvity.
One limitation of this study is the limited variable selection due to lack of statistical information in censuses. Additionally, the deprivation index was computed using information from the 2001 Census of Population and Housing (the only information available) while mortality data was recorded using information from the period 1994–2007. Some characteristics recorded in 2001 may not capture the real environment of the areas in the pass, when the deaths occurred. However, the data collection period for the construction of the index corresponds to the half point of the period in which the mortality data was collected and the relative position of the census track with respect to the deprivation index (i.e. the quintile of the index where it was located) remains very stable in time. It is also important to point out that given that the data is over a decade old the present patterns of associations may have changed, primarily due to the increased amount of immigration since 2001 and the activity rate downfall due to the economic crisis.
Appendix
Annex
Spatial models were built as Bayesian hierarchical models with two stages [62, 63]. The first stage was the observational model $p\left(yx\right)$, where y denoted the vector of observations and x are the unknown parameters following a Gaussian Markov random field (GMRF) denoted as $p\left(x\theta \right)$. The second stage was given by the hyperparameters θ and their respective prior distribution $p\left(\theta \right)$. The desired posterior marginals
of the GMRF were approximated using the finite sum
where $\tilde{p}\left({x}_{i}\theta ,y\right)$and $\tilde{p}\left(\theta y\right)$ denoted approximations of $p\left({x}_{i}\theta ,y\right)$ and $p\left(\theta y\right)$, respectively. The finite sum (A1) was evaluated at support points ${\theta}_{k}$ using appropriate weights ${\Delta}_{k}$.
The posterior marginal $p\left(\theta y\right)$ of the hyperparameters is approximated using a Laplace approximation [71].
where the denominator ${\tilde{p}}_{G}\left(x\theta ,y\right)$ denoted the Gaussian approximation of $p\left(x\theta ,y\right)$ and ${x}^{*}\left(\theta \right)$ was the mode of the full conditional $p\left(x\theta ,y\right)$[72].
According to Rue et al.[63], it is sufficient to “numerically explore" this approximate posterior density using suitable support points ${\theta}_{k}$ in (A1). In this paper, these points were defined in the hdimensional space, using the strategy called central composite design. Here, centre points were augmented with a group of star points, which allowed for estimating the curvature of $\tilde{p}\left(\theta y\right)$.
Here, to approximate the first component of (A1) a simplified Laplace approximation (less expensive from a computational point of view with only a slight loss of accuracy) was used [62–64].
Abbreviations
 DI:

Deprivation index
 PCA:

Principal component analysis
 DP2:

Pena distance
 OR:

Odds ratio
 DM2:

Diabetis mellitus type 2
 BMI:

Body mass index
 RR:

Relative risk
 95% CI:

95% confidence interval
 CPH:

Census population and housing
 SD:

Standard deviation
 SMR:

Standardized mortality rates
 ICD10:

International classification of disease v. 10
 BYM:

Besag, York and Mollié model
 CAR:

Conditional autoregressive model
 DIC:

Deviance information criterion
 CPO:

Conditional predictive ordinate.
References
 1.
Townsend P: Deprivation. J Soc Policy. 1987, 16: 125146. 10.1017/S0047279400020341.
 2.
Folwell K: Single measure of deprivation. J Epidemiol Community Health. 1995, 49 (suppl 2): S51S56. 10.1136/jech.49.Suppl_2.S51.
 3.
Tunstall HVZ, Shaw M, Dorling D: Places and health. J Epidemiol Community Health. 2004, 58: 610. 10.1136/jech.58.1.6.
 4.
Townsend P, Phillimore P, Beattie A: Health and deprivation: inequality and the North. 1988, Routledge, London
 5.
Carstairs V, Morris R: Deprivation and health in Scotland. 1991, Aberdeen University Press, Aberdeen
 6.
Picket KE, Pearl MJ: Multilevel analysis of neighbourhood socioeconomic context and health outcomes: a critical review. J Epidemiol Community Health. 2001, 55 (2): 111122. 10.1136/jech.55.2.111.
 7.
Galobardes B, Lynch J, Davey Smith G: Measuring socioeconomic position in health research. Br Med Bull. 2007, 81 and 82: 2137.
 8.
OcañaRiola L, Saurina C, FernándezAjuria A, Lertxundi A, SánchezCantalejo C, Saez M, et al: Area deprivation adn mortality in the provincial capital cities of Andalusia and Catalonia (Spain). J Epidemiol Community Health. 2008, 62: 147152. 10.1136/jech.2006.053280.
 9.
Noble M, Wirght G, Dibben C, Smith GAN, McLennan D, Anttila C, et al: The English Indices of Deprivation 2004 (revised). 2004, Neighbourhood Renewal Unit, London, Office of the Deputy Prime Minister
 10.
Elliot P, Wartenberg D: Spatial Epidemiology: current approaches and future challenges. Environ Health Persp. 2004, 112: 9981006. 10.1289/ehp.6735.
 11.
Lawson AB: Statistical methods in spatial epidemiology. 2006, John Wiley and Sons, New Jersey, 2
 12.
Barcelò MA, Saez M, Saurina C: Spatial variability in mortality inequalities, socioeconomic deprivation, and air pollution in small areas of the Barcelona Metropolitan Region, Spain. Sci Total Environ. 2009, 407: 55015523. 10.1016/j.scitotenv.2009.07.028.
 13.
Carstairs V: Socioeconomic factors at area level and their relationship with Health. Sparial epidemiology. Methods and applications. Edited by: Elliott P, Wakefield JC, Best NG, Briggs DJ. 2000, Oxford University Press, New York, 5167.
 14.
Bell N, Schuurman N, Hayes MV: Using GISbased methods of multicriteria analysis to construct socioeconomic deprivation indices. Int J Health Geogr. 2007, 6: 1710.1186/1476072X617.
 15.
Townsend P: Life and Labour in London. 1987, CPAG, London
 16.
Carstairs V: Deprivation: explaining differences in mortality between Scotland and England and Wales. Br Med J. 1989, 299: 886889. 10.1136/bmj.299.6704.886.
 17.
Frohlich KL, Mustard C: A regional comparison of socioeconomic and health indices in a Canadian province. Soc Sci Med. 1996, 42 (9): 12731281. 10.1016/02779536(95)002200.
 18.
Borrell C, MaríDell’olmo M, Serral G, MartínezBeneito M, Gotsens M, MEDEA members: Inequalities in mortality in small areas of eleven Spanish cities (the multicenter MEDEA project). Health Place. 2010, 16 (4): 703711. 10.1016/j.healthplace.2010.03.002.
 19.
Pampalon R, Raymond G: A Deprivation Index for Health and Welfare Planning in Quebec. Chronic Dis Can. 2000, 21 (3): 104113.
 20.
DomínguezBerjón F, Borrell C, CanoSerral G, Esnaola S, Nolasco A, Pasarín MI, Ramis R, Saurina C, EscolarPujolar A: Constructing a deprivation index based on census data in large Spanish cities (the MEDEA Project) [In Spanish]. Gac Sanit. 2008, 22 (3): 179187. 10.1157/13123961.
 21.
Jarman B: Identification of underprivileged areas. Br Med J. 1983, 286: 17051708. 10.1136/bmj.286.6379.1705.
 22.
Talbot RJ: Underprivileged areas and health care planning: implications of use of Jarman indicators of urban development. Br Med J. 1991, 302: 283286. 10.1136/bmj.302.6771.283.
 23.
CarrHill R, Sheldon R: Designing a deprivation payment for general practitioners: The UPA (8) Wonderland. Br Med J. 1991, 302: 393396. 10.1136/bmj.302.6773.393.
 24.
Davey Smith G: Second thoughts on the Jarman Index – Claims being made for the underprivileged area score are unproved. Br Med J. 1991, 302: 359360. 10.1136/bmj.302.6773.359.
 25.
Regidor E, Calle ME, Navarro P, Domínguez V: Trends in the association between average income, poverty and income inequality and life expectancy in Spain. Soc Sci Med. 2003, 56: 961971. 10.1016/S02779536(02)001077.
 26.
Mackenbach JP: Health inequalities: Europe in profile. 2006, University Medical Center Rotterdam, Rotterdam
 27.
Benach J, Yasui Y, Borrell C, Pasarín MI, Martínez JM, Daponte A: The public health burden of material deprivation: excess mortality in leading causes of death in Spain. Prev Med. 2003, 36: 300308. 10.1016/S00917435(02)000142.
 28.
Benach J, Yasui Y, Borrell C, Saez M, Pasarin MI: Material deprivation and leading causes of death by gender: evidence from a nationwide smallarea study. J Epidemiol Commun Health. 2001, 55: 239245. 10.1136/jech.55.4.239.
 29.
Molinari C, Ahern M, Hendryx M: The relationship of community quality to the health of women and men. Soc Sci Med. 1998, 47 (8): 11131120. 10.1016/S02779536(98)001142.
 30.
Mohai P: Gender differences in the perception of most important environmental problems. Race Gender Class. 1997, 5: 153169.
 31.
Stafford M, Bartley M, Mitchell R, Marmot M: Characteristics of individuals and characteristics of areas: investigating their influence on health in the Whitehall II study. Health Place. 2001, 7 (2): 117129. 10.1016/S13538292(01)000041.
 32.
Stafford M, Cummins S, MacIntyre S, Ellaway A, Marmot M: Gender differences in the associations between health and neighbourhood environment. Soc Sci Med. 2005, 60: 16811692. 10.1016/j.socscimed.2004.08.028.
 33.
Ellaway A, Macintyre S: Women in their place. Gender and perceptions of neighbourhoods and health in the West of Scotland. Geographies of women’s health. Edited by: Dyck I, Davis Lewis N, McLafferty S. 2001, Routledge, London, 265281.
 34.
Mustard CA, Etches J: Gender differences in socioeconomic inequality in mortality. J Epidemiol Comm Health. 2003, 57: 974980. 10.1136/jech.57.12.974.
 35.
Macintyre S, McKay L, Der G, Hiscock R: Socioeconomic position and health: what you observe depends on how you measure it. J Public Health Med. 2003, 25 (4): 288294. 10.1093/pubmed/fdg089.
 36.
Kawachi I, Kennedy BP, Gupta V, ProthrowStith D: Women’s status and the health of women and men: a view from the States. Soc Sci Med. 1999, 48: 2132. 10.1016/S02779536(98)00286X.
 37.
Van Praag L, Bracke P, Christians W, Lavecque K, Pattyn E: Mental health in a gendered context. Gendered community effects on depression and problem drinking. Health Place. 2009, 15 (4): 990998. 10.1016/j.healthplace.2009.04.003.
 38.
Santos SM, Chor D, Werneck GL, FreireCoutinho ES: Association between contextual factors and selfrated health: a systematic review of multilevel studies [In Portuguese]. Cad Saude Publica. 2007, 23 (11): 25332554. 10.1590/S0102311X2007001100002.
 39.
Griffin J, Further R, Stansfeld S, Marmot M: The importance of low control at work and home on depression and anxiety: do these effects vary by gender and social class?. Soc Sci Med. 2002, 54 (5): 738798.
 40.
Windenberger F, Rican S, Jougla E, Rey G: Spatiotemporal association between deprivation and mortality: trends in France during the nineties. Eur J Public Health. 2011, 22 (3): 347353.
 41.
Pena JB: Problems in measuring welfare and related concepts (an application to the Spanish case) [in Spanish]. 1977, INE, Madrid
 42.
Zarzosa P: Poverty estimate in the Spanish autonomous regions, using the P2Distance indicator [in Spanish]. Estudios de Economia Aplicada. 2009, 27 (2): 397416.
 43.
Mishra SK: Construction of maximin and nonelitist composite indices – alternatives to elitist indices obtained by the principal component analysis. MPRA Paper, 3338. 2007, University Library of Munich, Germany
 44.
Montero JM, Chasco C, Larraz B: Building an environmental quality index for a big city: a spatial interpolation approach combined with a distance indicator. J Geogr Syst. 2010, 12: 43559. 10.1007/s1010901001086.
 45.
Salcedo N, Saurina C, Saez M: Deprivation Indices. Gender differences and geographical variations of mortality inequalities at the small area level. 2010, GRECS, University of Girona, Girona, mimeo
 46.
PuigpinósRiera R, MaríDell'Olmo M, Gotsens M, Borrell C, Serral G, Ascaso C, Calvo M, Daponte A, DomínguezBerjón FM, Esnaola S, Gandarillas A, LópezAbente G, Martos CM, MartínezBeneito MA, MontesMartinez A, Montoya I, Nolasco A, Pasarín MI, RodríguezSanz M, Saez M, SánchezVillegas P: Cancer mortality inequalities in urban areas: a Bayesian small area analysis in Spanish cities. Int J Health Geogr. 2011, 10: 610.1186/1476072X106.
 47.
Dray R, Gary T, Brancati F: Educational disparities in mortality among adults with Diabetes in the US. Diabetes Care. 2010, 33 (6): 12001205. 10.2337/dc092094.
 48.
Chaturvedi N, Jarrett J, Shipley MJ, Fuller JH: Socioeconomic gradient in morbidity and mortality in people with diabetes: cohort study findings from the Whitehall Study and the WHO Multinational Study of Vascular Disease in Diabetes. BMJ. 1998, 316 (7125): 100105. 10.1136/bmj.316.7125.100.
 49.
Espelt A, Borrell C, Roskam AJ, RodriguezSanz M, Stirbu I, DalmauBueno A, Regidor E, Bopp M, Martikainen P, Leinsalu M, Artnik B, Rychtarikova J, Kalediene R, Dzurova D, Mackenbach J, Kunst AE: Socioeconomic inequalities in diabetes mellitus across Europe at the beginning of the21st century. Diabetologia. 2008, 51 (11): 197179. 10.1007/s0012500811461.
 50.
Feja C, Alcalá JT, Rabanaque MJ, Saez M, MarcosGragera R, Martos CM: Geographical inequalities in mortality and incidence in larynx cancer in men: socioeconomic and environment factors [In Spanish]. Rev Esp Salud Publica. 2010, 84 (6): 74556.
 51.
Ivanovic B: How to prepare a list of indicators of development. Revue de Statistique Apliquée. 1974, 22 (2): 3750.
 52.
Lawson AB, Browne WJ, VidalRodeiro CL: Disease mapping with WinBUGS and MLwiN. 2003, John Wiley & Sons, Chichester
 53.
Barceló A, Saez M, CanoSerral G, MartínezBeneito MA, Martínez JM, Borrell C, OcañaRiola R, Montoya I, Calvo M, LópezAbente G, RodríguezSanz M, Toro S, Alcalá JR, Saurina C, SánchezVillegas P, Figueiras A: Methods to smooth mortality indicators: application to the analysis of inequalities in mortality in Spanish cities (the MEDEA Project) [In Spanish]. Gac Sanit. 2008, 22 (6): 596608. 10.1016/S02139111(08)753627.
 54.
Greene WH: Econometric Analysis. 2003, Pearson Education, Upper Saddle River, New Jersey, 191201. 5
 55.
Besag J, York J, Mollié A: Bayesian image restoration, with two applications in spatial statistics. Ann I Stat Math. 1991, 43: 159. 10.1007/BF00116466.
 56.
Mollié A: Bayesian mapping of disease. Markov Chain Monte Carlo in Practice. Edited by: Gilks WR, Richardson S, Spiegelhalter DJ. 1996, Chapman & Hall, New York, 359379.
 57.
Besag J: Spatial interaction and the statistical analysis of lattice systems (with discussion). J R Stat Soc Ser B. 1974, 36: 192236.
 58.
Clayton DG, Bernadinelli L, Montomoli C: Spatial correlation in ecological analysis. Int J Epidemiol. 1993, 22: 11931202. 10.1093/ije/22.6.1193.
 59.
Rosenbaum P, Rubin D: Difficulties with regression analyses of ageadjusted rates. Biometrics. 1984, 40: 437443. 10.2307/2531396.
 60.
Grisotto L, Catelan D, Accetta G, Biggeri A: Material deprivation as marker of health needs. Statistica. 2010, 70 (3): 343352.
 61.
Catelan D, Biggeri A, Lagazio C: On the clustering term in ecological analysis: how do different prior specifications affect results?. Stat Methods Appl. 2009, 18: 4961. 10.1007/s102600070089x.
 62.
Schrödle B, Held L: A primer on disease mapping and ecological regression using INLA. Comput Stat. 2011, 26 (2): 241258. 10.1007/s0018001002082. [http://www.rinla.org/papers]
 63.
Rue H, Martino S, Chopin N: Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations (with discussion). J R Stat Soc Ser B. 2009, 71: 319392. 10.1111/j.14679868.2008.00700.x. [http://www.rinla.org/papers]
 64.
Martino S, Rue H: Case studies in Bayesian computation using INLA. Complex Data Modeling and Computationally Intensive Statistical Methods. Edited by: Mantovan P, Secchi P. 2010, Springer, Milano, 99112. [http://www.rinla.org/papers]
 65.
Markov Chain Monte Carlo in Practice. Edited by: Gilks WR, Richardson S, Spiegelhalter DJ. 1996, Chapman & Hall, New York
 66.
Gelfand AE, Dey DK, Chang H: Model determination using predictive distributions with implementation via samplingbased methods (with discussion). Bayesian Statistics 4. Edited by: Bernardo JM, Berger JO, Dawid AP, Smith AFM. 1992, Oxford University Press, Oxford
 67.
Pettit LI: The conditional predictive ordinate for the normal distribution. J R Stat Soc Ser B. 1990, 52 (1): 175184.
 68.
,: RINLA project. 2010, http://www.rinla.org/home,
 69.
R Development Core Team. R: A language and environment for statistical computing. 2010, R Foundation for Statistical Computing, Vienna, Austria, http://www.Rproject.org,
 70.
Dobson A, McLaughlin D, Vagenas D, Wong KY: Why are death rates higher in rural areas? Evidence from the Australian Longitudinal Study on Women´s Health. Aust N Z J Public Health. 2010, 34 (6): 624628. 10.1111/j.17536405.2010.00623.x.
 71.
Tierney L, Kadane JB: Accurate approximations for posterior moments and marginal densities. J Am Stat Assoc. 1986, 81 (393): 8286. 10.1080/01621459.1986.10478240.
 72.
Rue H, Held L: Gaussian Markov Random Fields. 2005, Chapman & Hall/CRC, Boca RatonLondonNew YorkSingapore
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712458/12/574/prepub
Acknowledgements
We would like to the Catalan Mortality Register (Registre de Mortalitat de Catalunya, Servei d’Informació i Estudis, Direcció General de Recursos Sanitaris), Health Department, Catalan Government (Generalitat de Catalunya), for providing mortality data. This work was partially supported by the European Union, DGSANCO, Second Programme of Community action in the field of Health (2008–2013), project A/101156 and by the FIS (Health Research Fund), Spanish Ministry of Science and Innovation, project FIS08/0142.
Author information
Additional information
Competing interests
There are no conflicts of interest for any of the authors. All authors disclose any financial and personal relationships with other people or organizations that could inappropriately influence and/or bias their work.
Authors’ contributions
All authors contributed equally to the development of the manuscript. MS had the idea of writing this manuscript, based on the Masters thesis of NS. NS and CS were responsible for choosing the variables and collect data. CS and MS specified the models. MS and BB were responsible for the statistical analyzes. Interpretation of the results were done mainly by CS, NS and MS. Drafting of the manuscript were done by CS, MS and NS. All authors revised and approved the final version of the manuscript.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 Deprivation indexes
 Distance indicator DP_{2}
 Ecological regression
 Standardized mortality ratio
 Hierarchical Bayesian models
 Robust