Next Article in Journal
Flood Disaster Risk Assessment of Rural Housings — A Case Study of Kouqian Town in China
Next Article in Special Issue
What Causes Environmental Inequalities and Related Health Effects? An Analysis of Evolving Concepts
Previous Article in Journal
Exposures of 129 Preschool Children to Organochlorines, Organophosphates, Pyrethroids, and Acid Herbicides at Their Homes and Daycares in North Carolina
Previous Article in Special Issue
The Socioeconomic Determinants of Health: Economic Growth and Health in the OECD Countries during the Last Three Decades
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Relationship Quantification between Environmental, Socioeconomic and Health Data at Different Geographic Levels

1
French National Institute for Industrial Environment and Risks, Parc Technologique Alata, BP 2, 60550 Verneuil-en-Halatte, France
2
University of Picardie Jules Verne, 33 rue St Leu, Amiens 80039, France
3
University Hospital of Amiens, Place Victor Pauchet Amiens 80054, France
4
Regional Observatory of Health and Social Issues in Picardie (OR2S), 3, rue des Louvels, Amiens 80036, France
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2014, 11(4), 3765-3786; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110403765
Submission received: 31 October 2013 / Revised: 18 March 2014 / Accepted: 19 March 2014 / Published: 3 April 2014
(This article belongs to the Special Issue Inequalities in Health)

Abstract

:
Spatial health inequalities have often been analyzed in terms of socioeconomic and environmental factors. The present study aimed to evaluate spatial relationships between spatial data collected at different spatial scales. The approach was illustrated using health outcomes (mortality attributable to cancer) initially aggregated to the county level, district socioeconomic covariates, and exposure data modeled on a regular grid. Geographically weighted regression (GWR) was used to quantify spatial relationships. The strongest associations were found when low deprivation was associated with lower lip, oral cavity and pharynx cancer mortality and when low environmental pollution was associated with low pleural cancer mortality. However, applying this approach to other areas or to other causes of death or with other indicators requires continuous exploratory analysis to assess the role of the modifiable areal unit problem (MAUP) and downscaling the health data on the study of the relationship, which will allow decision-makers to develop interventions where they are most needed.

1. Introduction

Analyzing the relationship between the environment and health has become a major issue for public health in France as forecasted by the national plans for health and environment (NPHE). Two priority areas were selected during the first NPHE: (1) preventing health risks related to the quality of resources and to chemicals and (2) developing environmental health through research, expertise, training and information. In 2009, the second NPHE was prepared from the perspective of the upcoming conference on health and the environment organized by the World Health Organization. Two main axes were prioritized: (1) identifying and managing geographic areas where hotspot exposures to substances present in air, soil, water, and foods resulting from anthropic activities suspected of generating potentially increasing risks to human health and (2) reducing environmental health inequalities. Thus, environmental health inequality has become a substantial topic that guides policy developments in France. To address this aim, there is an urgent need for tools that can quantify the spatial relationships between the environment, socioeconomics and health and that can highlight areas with strong inequalities.
Health inequalities are a quite recent study topic. Previous studies were essentially based, at an individual level, on specific surveys [1,2] and, at a spatially aggregated level (administrative unit), on specific regions [3,4]. At a regional scale, data are often available at a fine level or resolution. This allows for building environmental, socioeconomic and health indicators at different spatial scales; for example, Salmond et al. [5] built a new census-based index of deprivation based on the smallest possible geographical areas.
Regarding health data, there are strict privacy rules for individual-level health data that prohibit their public release. Aggregated data are only available at the geographic level, from which disclosure and reconstruction of patient identity are impossible. In France these census units could be regions or counties. This aggregation unfortunately results in incidence or mortality rates that can be unreliable over small and/or sparsely populated areas. This effect, known as the “small number problem” [6], should be corrected for an accurate evaluation of health-environment relationships.
Several authors have already addressed the spatial relationships between health data and environmental data. One of the issues faced by spatial epidemiologists and for exposure assessment is the combination of data measured for very different spatial scales and with different levels of reliability. In reality, the analysis of cancer mortality maps is often hindered by the presence of noise caused by unreliable extreme rates computed from sparsely populated geographic units. A number of approaches have been developed to improve the reliability of risk estimates [7,8]. The most commonly used are Bayesian methods [9], which are commonly referred to as the BYM model. Bayesian methods prohibit any change of scales, an operation that is easily conducted within the framework of kriging. Goovaerts and Gebreab [10] conducted a simulation-based evaluation of the performance of geostatistical and full Bayesian disease-mapping models, and they found that the geostatistical approach yielded smaller prediction errors and more precise and accurate probability intervals and that it allowed for better discrimination between counties with high and low mortality risks.
Poisson kriging, in this context, presents a spatial methodology that allows for filtering the noise caused by the small number problem and enables the estimation of mortality risk and the associated uncertainty at different spatial scales. This approach has been implemented to modeling cancer risk by a number of authors: Oliver et al. [11] studied cases of cancer in children under fifteen years of age, and Goovaerts and collaborators considered lung cancer [12,13], breast cancer [14,15], prostate cancer [16], cervical cancer [17], and pancreatic cancer [18], and all found it to be relevant for this particular problem.
Selection of scale is perhaps the most important factor in creating and analyzing a relationship between environmental exposure and health outcomes [19]. This issue is similar to the modifiable area unit problem (MAUP), a term introduced by Openshaw [20,21]. The MAUP can cause differences in the analytical results of the same input data compiled under different zoning systems [22,23].
The present study aims to evaluate spatial relationships at three levels of aggregation: the IRIS level, an intermediate scale (the grid level), and the county level between health outcomes (mortality attributable to cancer) initially aggregated to the county level, district socioeconomic covariates, and exposure data modeled on a regular grid. The approach is illustrated using age-adjusted lip, oral cavity and pharynx, and pleural cancer mortality rates over the period 2000–2009 for the Picardy region. The deprivation index and trace metal exposure indicators are used as putative risk factors.

2. Materials and Methods

2.1. Study Area

The region of Picardy covers an area of roughly 19,500 km2 and is located between the North Artois, the Ile-de-France in the south, the Bay of the Somme to the west and the East Champagne. It covers the departments of Somme, Oise and Aisne. The urbanization rate in this region is far below the national average (60.4% versus 74% for the entire country). The agricultural sector provides more than 4% of French agricultural production. It also has significant industrial activity through which fine chemicals and specialty chemicals account for nearly 15% of jobs and the vehicle industry comprises 40% of industrial employment (26.5% of assets employed in industry versus 19.5% at the national level). Three administrative or statistical spatial units, of different sizes, were considered: IRIS (the smallest administrative units in Picardy, 2,129 units) with irregular sizes and shapes, 64 km2 grid cells (308 units) that are all squares of same size, and counties (112 units) with irregular sizes and irregular shapes. Figure 1 shows the counties of the study area.

2.2. Data

2.2.1. Exposure Indicators

The environmental indicators (inhalation and ingestion) used were those described in Caudeville et al. for building GIS-based modeling platforms for quantifying human exposure to chemical substances [24]. The exposure indicators integrate soil, water, air, food, demographic and behavioral geo-referenced data to construct population exposure doses and associated risks at a fine resolution (1 km2 grid).
Figure 1. Map of the study area.
Figure 1. Map of the study area.
Ijerph 11 03765 g001

0201 Anizy-le-Château6001 Attichy8001 Abbeville Nord
0202 Aubenton6002 Auneuil8002 Abbeville Sud
0203 Bohain-en-Vermandois6004 Beauvais Sud-Ouest8003 Acheux-en-Amiénois
0204 Braine6005 Betz8004 Ailly-le-Haut-Clocher
0205 La Capelle6006 Breteuil8005 Ailly-ser-Noye
0206 Le Catelet6007 Chaumont-en-Vexin8006 Albert
0207 Charly6008 Clermont8007 Amiens Ouest
0208 Château Thierry6009 Compiègne Nord8008 Amiens Nord-Ouest
0209 Chauny6010 Le Coudray-Saint-Germer8009 Amiens-Nord-Est
0210 Condé-en-Brie6011 Creil-Nogent-sur-Oise8010 Amiens-Est
0211Coucy-le-Château Auffrique6012 Crépy-en-Valois8011 Ault
0212 Craonne6013 Crèvecoeur-le-Grand8012 Bernaville
0213 Crécy-sr-Serre6014 Estrées-Saint-Denis8013 Boves
0214 La Fère6015 Formerie8014 Bray-sur-Somme
0215 Fère-en-Tardenois6016 Froissy8015 Chaulnes
0216 Guise6017 Grandvilliers8016 Combles
0217 Hirson6018 Guiscard8017 Conty
0218 Laon Nord6019 Lassigny8018 Corbie
0219 Marle6020 Liancourt8019 Crécy-en-Ponthieu
0220 Moy-de-l’Aisne6021 Maignelay-Montigny8020 Domart-en-Ponthieu
0221 Neufchatel-sur-Aisne6022 Marseille-en-Beauvaisis8021 Doullens
0222 Neuilly-Saint-Front6023 Méru8022 Gamaches
0223 Le Nouvion-en-Thiérache6024 Mouy8023 Hallencourt
0224 Oulchy-le-Château6025 Nanteuil-le-Haudoin8024 Ham
0225 Ribemont2026 Neuilly-en-Thelle8025 Homoy-le-Bourg
0226 Rozoy-sur-Serre2027 Nivillers8026 Molliens-Dreuil
0227 Sains-Richaumont6028 Noailles8027 Montdidier
0229 Saint-Simon6029 Noyon8028 Moreuil
0230 Sissonne6030 Pont-Sainte-Maxence8029 Moyenneville
0231 Soissons-Nord6031 Ressons-sur-Matz8030 Nesle
0232 Vailly-sur-Aisne6032 Ribécourt-Dreslincourt8031 Nouvion
0233 Vermand6033 Saint-Just-en-Chaussée8032 Oisemont
0234 Vervins6034 Senlis8033 Péronne
0235 Vic-sur-Aisne6035 Songeons8034 Picquigny
0236 Villers-Cotterets6036 Chantilly8035 Poix-de-Picardie
0237 Wassigny6037 Compiègne-Sud-Est8036 Roisel
0238 Laon-Sud6039 Montataire8037 Rosières-en-Santerre
0239 Saint-Quentin-Nord6040 Beauvais-Nord-Ouest8038 Roye
0240 Saint-Quentin-Sud6041 Compiègne Sud-Oest8039 Rue
0241 Soissons-Sud6097 Compiègne8040 Saint-Valery-sur-Somme
0242 Tergnier6098 Creil8041 Villers-Bocage
0297 Laon6099 Beauvais8042 Amiens Sud-Est
0298 Saint-Quentin8043 Amiens Sud-Ouest
0299 Soisson 8044 Amiens Nord
8046 Friville-Escarbotin
8098 Abbeville
8099 Amiens
Laon(ville et cantons) comprend les cantons 0218,0238 et 0297Compiègne(ville et cantons) comprend les cantons 6009, 6037, 6041 et 6097Amiens(ville et cantons) comprend les cantons 8007,8008, 8009, 8010 8042,8044,8045 et 8099
Saint-Quentin(ville et cantons) comprend les cantons 0239,0240 et 0298Creil-Nogent-sur-Oise comprend les cantons 6011 et 6098 Abbeville(ville et cantons) comprend les cantons 8001,8002, et 8098
Soissons(ville et cantons) comprend les cantons 0231,0241 et 0299Beauvais(ville et cantons) comprend les cantons 6004, 6040 et 6099
Trace elements (nickel-Ni, cadmium-Cd, and lead-Pb) were modeled within the Picardy region [25] for various exposure pathways: atmospheric contaminant inhalation and ingestion of soil, vegetation, meat, eggs, milk, fish and drinking water.

2.2.2. Deprivation Index (SE)

The deprivation index used was developed by Rey [26] and was built at the French census block (IRIS) using the following socioeconomic data: the median household income, the percentage of high school graduates in the population aged 15 years and older, the percentage of blue-collar workers in the active population, and the unemployment rate. The deprivation index was also constructed for the county. For each county, the deprivation index was calculated as the population-weighted average score for all of the IRISes in the county.

2.2.3. Health Data

The health data came from the Regional Health Observatory of Picardy [27], where the age-adjusted mortality rates are calculated for each county from 2000 to 2009. Table 1 shows the cumulative, maximum and minimum number of mortality and age-adjusted rates per 100,000 person-years by county from 2000 to 2009.
Table 1. Cumulative, maximum and minimum number of mortality and age-adjusted rates per 100,000 person-years by county, 2000–2009.
Table 1. Cumulative, maximum and minimum number of mortality and age-adjusted rates per 100,000 person-years by county, 2000–2009.
Cancer MortalityNumbers of CasesAge-adjusted Rates Per 100,000 Person-years
Lip, oral cavity and pharynx cancer mortality
Cumulative1,32716.26
Minimum12.81
Maximum12837.4
Pleural cancer mortality
Cumulative2633.78
Minimum00
Maximum1811.94
Figure 2 shows the spatial distribution age-adjusted lip, oral cavity, pharynx, and pleura cancer mortality per 100,000 person-years. It should be noted that (1) the population is not evenly distributed throughout the study area and (2) the age-adjusted rate calculated from the less-populated counties tend to be less reliable. This implies that the interpretation of the map must be carried out with caution. The scatter plot at the bottom of Figure 1 illustrates this effect, well-known as the “small number problem”.
Table 2 presents the different spatial scales of measurement and the approaches used to homogenize spatial coverage.
Figure 2. (a) Map of log population density. Geographic distribution of age-adjusted mortality rates per 100,000 person-years recorded over the period 2000–2009 for: (b) lip, oral cavity and pharynx; (c) pleura cancer mortality. The bottom scatter plots illustrate: (d) the age-adjusted mortality rates for lip, oral cavity and pharynx cancers plotted against population density and (e) the age-adjusted mortality rates of pleura cancers plotted against population density.
Figure 2. (a) Map of log population density. Geographic distribution of age-adjusted mortality rates per 100,000 person-years recorded over the period 2000–2009 for: (b) lip, oral cavity and pharynx; (c) pleura cancer mortality. The bottom scatter plots illustrate: (d) the age-adjusted mortality rates for lip, oral cavity and pharynx cancers plotted against population density and (e) the age-adjusted mortality rates of pleura cancers plotted against population density.
Ijerph 11 03765 g002
Table 2. Spatially resolved data types and approaches used to homogenize spatial coverage.
Table 2. Spatially resolved data types and approaches used to homogenize spatial coverage.
IndicatorVariablesSourcesSpatial Scale or ResolutionSpatial Operation
SocioeconomicSE: Deprivation index French census Rey et al. [26]Vector data from the IRIS.Spatial population-weighted aggregation
ExposureF1: Exposure inhalation indicator
F2: Exposure ingestion indicator
Caudeville et al. [24,25]Raster data of 1 km2 gridSpatial aggregation
HealthLip, oral cavity and pharynx cancer mortalityRegional Health Observatory of Picardy [27]Vector data from the county databasePoisson kriging
Pleural cancer mortality

2.3. Methods

2.3.1. The Geostatistical Approach: Correcting Small Numbers and Estimating the Corresponding Risk at Different Spatial Scales

To correct for the instability attributable to the small number problem, a number of algorithms have been developed that aim at estimating risk. The geostatistical approach, in this context, presents an interesting alternative; it conducts the noise filtering and allows for risk estimation along with the associated uncertainty at different scales. This section provides a brief overview of the geostatistical methodology for estimating risk values. See Goovaerts [17] for more details about this approach.
The cancer mortality count d(να) within a county να is interpreted as the realization of a random variable D(να) that is Poisson distributed with a parameter (expected number of counts) that is the product of the population size n(να) by the local risk R(να). R(να) might be thought of as a noise-filtered mortality rate for area να, which we also refer to as the mortality risk. It is estimated by using a variant of kriging with nonsystematic errors known as Poisson kriging [28].
The mortality risk and the associated kriging variance for an area να are estimated as:
Ijerph 11 03765 i001
Kriging variance is computed as follows:
Ijerph 11 03765 i002
where x represents either an area να (ATA kriging). The kriging weights (λi) and the Lagrange parameter μ (να) are computed by solving the Poisson kriging system of equations:
Ijerph 11 03765 i003
Ijerph 11 03765 i004
where δij = 1 if i = j and 0 otherwise. The “error variance” term, m*/n(νi) leads to smaller weights for rates measured over smaller population sizes. The covariance CR (νi, νj) is approximated as the population-weighted average of the point-support covariance CR (h) computed between any two discrete locations between the areas νi and νj.

2.3.2. Spatial Autocorrelation

Global Moran’s I was calculated for all of the explanatory variables as well as for the dependent variables within three spatial structures to determine the role of spatial representation using global spatial autocorrelation. The Global Moran’s I spatial autocorrelation statistic measures the self-similarity of a spatial variable’s value as a function of adjacency [29], using a first-order Queen’s case spatial weight matrix and 999 permutations.

2.3.3. Exploring the Relationships between Health, Environment and Socioeconomic Factors

Analyses of correlations between health data and putative factors are traditionally performed using a global or “aspatial” regression model, under the implicit assumption that the impact of variables is constant over the entire study area. This assumption is likely unrealistic for large areas, which can display large geographic variations. Fotheringham and colleagues developed Geographically Weighted Regression (GWR) to explore spatial non-stationarity and map statistics to visualize the spatial patterns of the relationships between dependent and independent variables [30,31,32].

Aspatial Regression

The explanatory power of SE and exposure indicators was first investigated using the following multiple linear regression model:
Ijerph 11 03765 i005
where γ is the kriging risks estimate for observation i, β0 is the intercept, βk is the regression coefficient (slope) of each factor xk and εi is the error term. To account for the reliability of the kriged risks in the regression, each observation receives a weight that is the reciprocal of the kriging variance [33].

Geographically Weighted Regression.

In geographically weighted regression, the regression is conducted within local windows centered around each observation. The regression coefficients are thus location-dependent:
Ijerph 11 03765 i006
Within each window, observations are weighted according to their proximities to the center of the window. A variety of distance decay functions are available. In this paper, we used the XX function, which is characterized by a bandwidth that corresponds to the distance beyond which the weight rapidly approaches zero.
The bandwidth is estimated by minimizing the AICc value:
Ijerph 11 03765 i007
where n is the number of observations in the dataset, Ijerph 11 03765 i008 is the estimate of the standard deviation of the residuals, and tr(S) is the trace of the hat matrix. For more information on the theory and practical application of GWR, the reader is referred to Fotheringham et al. [34].

3. Results

3.1. Poisson Kriging for Health Indicator

Figure 3 and Figure 4 show the risk values with the corresponding prediction variance estimated by Poisson kriging at: (a) the county level; (b) the grid level and (c) the IRIS level. All maps are substantially smoother than the original rate map because the noise caused by small population sizes has been filtered. These maps allow a better visualization of areas of higher risks: the lip, oral cavity and pharynx cancer mortality rates vary between 2.81 and 37.40 per 100,000 person-years. After the application of Poisson kriging, the minimum rate increased from 2.81 to 8.87 deaths/100,000 person-years, and the maximum rate of 37.40 decreased to 25.14 deaths per 100,000 person-years. We can note, for instance, that the high rates recorded in sparsely populated counties such as Sains-Richaumont (37.40 deaths/100,000 person-years), north of the Aisne department, are strongly smoothed (24.46 deaths/100,000 person-years). The highest rate recorded in a densely populated county (i.e., Abbeville North county-26.60 deaths/100,000 person-years) remained nearly the same after smoothing (24.90 deaths/100,000 person-years). Zero pleural cancer mortality rates recorded in sparsely populated counties were also smoothed, leading to minimum values of 1.00 deaths/100,000 person-years.
The maps of the kriging variance indicate the higher reliability of risks estimated in densely populated areas such as Amiens, Beauvais, Saint Quentin, and Abbeville. The variance of the risk estimates decreased as the geographic unit area increased: from the IRIS level to the grid level and then to the county level (Table 3).
The risk estimates are characterized by positive spatial autocorrelation within the three spatial scales (p ≤ 0.05) but display low levels of statistically significant spatial autocorrelation at the IRIS level in comparison with the grid and county levels (Table 3). In this case, the counties are internally homogeneous in terms of mortality according to the risks estimated by kriging.
Figure 3. Maps of the lip, oral cavity and pharynx cancer mortality risk estimates and the corresponding prediction variance computed by Poisson kriging at three spatial scales: (a) county level; (b) grid level and (c) IRIS level.
Figure 3. Maps of the lip, oral cavity and pharynx cancer mortality risk estimates and the corresponding prediction variance computed by Poisson kriging at three spatial scales: (a) county level; (b) grid level and (c) IRIS level.
Ijerph 11 03765 g003
Figure 4. Maps of the pleural cancer mortality risk estimates and the corresponding prediction variances computed by Poisson kriging at three spatial scales: (a) county level; (b) grid level; and (c) IRIS level.
Figure 4. Maps of the pleural cancer mortality risk estimates and the corresponding prediction variances computed by Poisson kriging at three spatial scales: (a) county level; (b) grid level; and (c) IRIS level.
Ijerph 11 03765 g004
Table 3. Summary statistics for health indicators after applying Poisson kriging.
Table 3. Summary statistics for health indicators after applying Poisson kriging.
Lip. Oral Cavity and Pharynx Cancer Mortality
Estimation Type MeanMinMaxMorans’I
County Level
Kriging risk15.598.8825.140.65 (0.001)
Kriging variance8.361.8713.42
Grid Level
Kriging risk15.328.3125.920.78 (0.001)
Kriging variance16.062.8130.09
IRIS Level
Kriging risk15.357.3826.560.96 (0.001)
Kriging variance22.524.133.24
Pleural Cancer Mortality
VariablesMeanMinMaxMoran’s I
County Level
Kriging risk3.1618.720.52 (0.001)
Kriging variance1.920.433.07
Grid Level
Kriging risk2.990.878.480.62 (0.001)
Kriging variance3.650.636.67
IRIS Level
Kriging risk3.210.879.040.93 (0.001)
Kriging variance4.990.897.32

3.2. Spatial Aggregation for the Explanatory Variables

The mean values of the explanatory variables under each of the three spatial structures are similar (Table 4). The variables F1 (exposure inhalation indicator) and F2 (exposure ingestion indicator) display greatly reduced variability under different spatial structures. The variable SE has the highest levels of variability. In contrast, the SE index is affected by the use of different spatial structures. The lowest average was −4.5 at the county level in comparison with −7.32 at the IRIS level, and the variance decreased with increasing aggregation, unlike with F1 and F2, for which the lowest and strongest averages were somewhat similar for the three spatial scales.
All of the variables were characterized by positive spatial autocorrelation within the three spatial scales at levels of p ≤ 0.05. For F1 and F2, the Moran’s I values of the three spatial structures were similar (Table 4); however, the F1 and F2 were not affected by the use of different spatial structures. This is explained by the fact that the exposure indicator presented a homogeneous distribution within each county (the original resolution of exposure indicator data was a grid (15 × 10 km); see Caudville et al. [24,25]. Figure 5 and Figure 6 show maps of the SE index and the F1 aggregated at (a) county level, (b) grid level and (c) IRIS level.
Table 4. Summary statistics for the explanatory variables.
Table 4. Summary statistics for the explanatory variables.
VariablesMeanMinMaxVarianceMoran's I
County Level
SE: Deprivation index0.61−4.53.482.840.63(0.001)
F1: Exposure inhalation indicator0.080.060.130.00020.81(0.001)
F2: Exposure ingestion indicator0.270.270.390.0020.61(0.001)
Grid Level
SE: Deprivation index0.58−5.14.13.020.70(0.001)
F1: Exposure inhalation indicator0.080.060.130.00020.88(0.001)
F2: Exposure ingestion indicator0.270.280.480.0020.61(0.001)
IRIS Level
SE: Deprivation index0.48−7.3−84.620.55(0.001)
F1: Exposure inhalation indicator0.080.060.152,00E−040.91(0.001)
F2: Exposure ingestion indicator0.260.310.680.0030.65(0.001)
Figure 5. Maps of the deprivation index computed at three spatial scales: (a) county level; (b) grid level; (c) IRIS level.
Figure 5. Maps of the deprivation index computed at three spatial scales: (a) county level; (b) grid level; (c) IRIS level.
Ijerph 11 03765 g005
Figure 6. Maps of the exposure inhalation indicator (F1) aggregated at: (a) county level; (b) grid level; (c) IRIS level.
Figure 6. Maps of the exposure inhalation indicator (F1) aggregated at: (a) county level; (b) grid level; (c) IRIS level.
Ijerph 11 03765 g006

3.3. Aspatial Regression

Table 5 shows the correlation coefficient between each covariate and mortality rates before and after noise-filtering using Poisson kriging. Filtering the noise because of the small number problem clearly enhanced the explanatory power of the covariates: the proportion of variance explained (adjusted R2) increased by nearly one order of magnitude: lip, oral cavity and pharynx cancer mortality, 0.11 to 0.26, and pleural cancer mortality, 0.11 to 0.25. The uncertainty attached to the risk estimates can be accounted for by weighting each estimate according to the inverse of its kriging variance, leading to correlation coefficients and adjusted R2 values that were slightly larger for pleural cancer mortality and a slightly lower mortality rate for lip, oral cavity and pharynx cancers.
Application of the linear model for lip, oral cavity and pharynx cancer mortality data explained a moderate proportion of the total variance (adjusted R2 = 0.22 and 0.19) at the county level and grid level, respectively. This proportion was lower when the analysis was conducted at the IRIS level (adjusted R2 = 0.11). It is noteworthy that the correlation coefficients for the SE factor were always significant for the different aggregation levels and were higher at the county level than at the IRIS level, which was the expected result because aggregation is known to increase the strength of correlation [20,35]. Linear association between SE index and cancer mortality has been demonstrated in other works [36,37,38,39].
Table 5. Results of the correlation analysis.
Table 5. Results of the correlation analysis.
Lip. Oral Cavity and Pharynx Cancer Mortality
Regression ModelSEF1F2Adjusted R2
County-level correlation
Age-adjusted rate0.32 *−0.110.040.11
Kriging risk0.53 *−0.270.060.26
Kriging risk (weighted)0.49 *−0.260.030.22
Grid-level correlation (64 km)
Kriging risk0.49 *−0.280.030.24
Kriging risk (weighted)0.43 *−0.260.010.19
IRIS-level correlation
Kriging risk0.37 *−0.210.010.15
Kriging risk (weighted)0.32 *−0.13 *−0.030.11
Pleural Cancer Mortality
Regression ModelSEF1F2Adjusted R2
County-level correlation
Age-adjusted rate−0.130.35 *0.030.11
Kriging risk−0.180.51 *0.020.25
Kriging risk (weighted)−0.160.52 *−0.010.28
Grid-level correlation (64 km)
Kriging risk−0.180.47 *0.040.20
Kriging risk (weighted)−0.170.49 *0.030.24
IRIS-level correlation
Kriging risk−0.010.46 *0.060.22
Kriging risk (weighted)0.040.50 *0.050.28
Notes: * Significant at α = 0.01; SE: Deprivation index; F1: Exposure inhalation indicator; F2: Exposure ingestion indicator.
The model explains a moderate proportion of the total variance when the dependent variable y is pleural cancer mortality for different levels of aggregation. The adjusted R2 ranges between 0.24 and 0.28, with significant correlation coefficients of up to 0.5 for F1. The results showed the consistency of the association between trace metal inhalation exposure and pleural mortality across the different aggregation levels. This is explained by the fact that the pleural mortality presented a homogeneous distribution within each county and the exposure indicator was not affected by the use of different spatial structures. Pleural mesothelioma is one of four types of mesothelioma, but it accounts for approximately 75 percent of all diagnoses of asbestos-related cancers. The disease starts in the pleura, the soft protective tissue surrounding the lungs, which can be attributed directly to its cause: repeated or heavy occupational exposure to airborne asbestos fibers. However, Peterson suggests that a significant number of cases of this cancer are apparently not asbestos-related and that even in the absence of asbestos, other agents can induce malignant mesothelioma [40]. Some types of nanoparticles, especially those containing nickel, could also be a risk of pleural diseases [41].

3.4. Geographically Weighted Regression (GWR)

In the aspatial analysis, we implicitly assumed that the impact of covariates was constant across the study area. This assumption is likely unrealistic for large areas, which can display substantial geographic variation in socioeconomic and environmental conditions. Therefore, the global statistics reported by this traditional regression model could potentially hide a number of interesting local relationships. The question is then to examine whether there are any meaningful spatial variations in these relationships.
County-level data were used to optimize the bandwidth of the GWR distance decay function. Figure 7 shows the relationship between AICc and bandwidth size for the two types of cancer. The following bandwidths were found to be optimal: lip, oral cavity and pharynx cancer mortality (47 km), and pleural cancer mortality (54 km). The local regression was based on the following number of closest observations, which represented 15%–20% of the data: 22 for the county level, 39 for the grid level, and 426 for the IRIS level. A comparison of the AICc values suggests that all of the GWR models outperformed the global model (Table 6).
These results strongly suggest that the relationships between cancer (lip, oral cavity and pharynx-pleural) mortality and the environmental and deprivation indexes are not stationary but instead vary over the study area. The application of GWR clearly enhances the explanatory power of the covariates for the three spatial levels: the proportion of variance explained (adjusted R2) is almost doubled (Table 6).
Figure 7. Impact of bandwidth size on the AICc of geographically weighted regression for each cancer.
Figure 7. Impact of bandwidth size on the AICc of geographically weighted regression for each cancer.
Ijerph 11 03765 g007
Figure 8 and Figure 9 show the results of the geographically weighted regression applied to the county, grid and IRIS data. Two statistics are displayed for each spatial scale: (a) the proportion of variance explained (R2) and (b) the local correlation coefficient.
The three spatial scales share the same southeast-northwest trend in the explanatory power of the local regression models for lip, oral cavity and pharynx cancer mortality: the lower mortality values in the south are better explained by the SE index than is the higher risk recorded in the northwest. As a recall, the largest R2 observed in the south (Figure 8 a) corresponds to areas with low SE index variation. Similar to the global model (Table 5), GWR led to local correlation coefficients and R2 values that were higher at the county level than at the IRIS level.
The lower pleural cancer mortality values are better explained in areas of low F1 variation (Figure 9a). The maps of the local correlation coefficients and the R2 values also present the same spatial structure over the different spatial coverages of aggregation: lower in the west оf tһе Aisne department and higher in the north of the Oise department. Very similar results were obtained for the different spatial scales.
Table 6. Comparison of local and global regression models at the three spatial scales.
Table 6. Comparison of local and global regression models at the three spatial scales.
Lip. Oral Cavity and Pharynx Cancer Mortality
Regression modelBandwidth SizeAdjusted R2AICc
County-level correlation
Global model47 km0.22567.00
Local model0.52513.47
Grid-level correlation (64 km)
Global model47 km0.191,530.76
Local model0.481,280.26
IRIS-level correlation
Global model47 km0.1110,932.00
Local model0.2110,112.32
Pleural Cancer Mortality
Regression ModelBandwidth SizeAdjusted R2AICc
County-level correlation
Global model0.28374.35
Local model54 km0.48348.09
Grid-level correlation (64 km)
Global model0.24931.65
Local model54 km0.49803.08
IRIS-level correlation
Global model0.286,219.21
Local model54 km0.465,852.26
Figure 8. Results of the geographically weighted regression applied to the lip mortality kriged rates: (a) maps of the proportions of variance explained by deprivation index (R2); (b) maps of the local correlation coefficients with the deprivation index.
Figure 8. Results of the geographically weighted regression applied to the lip mortality kriged rates: (a) maps of the proportions of variance explained by deprivation index (R2); (b) maps of the local correlation coefficients with the deprivation index.
Ijerph 11 03765 g008
Figure 9. Results of the geographically weighted regression applied to the pleural cancer mortality kriged rates: (a) maps of the proportion of variance explained by the exposure inhalation indicator (R2); (b) maps of the local correlation coefficients with the exposure inhalation indicator.
Figure 9. Results of the geographically weighted regression applied to the pleural cancer mortality kriged rates: (a) maps of the proportion of variance explained by the exposure inhalation indicator (R2); (b) maps of the local correlation coefficients with the exposure inhalation indicator.
Ijerph 11 03765 g009

4. Discussion

Our results substantiate the work on noise filtering described in the introduction section from Oliver et al. [11] and Goovaerts et al. [12,18]. Indeed, we found the following: (1) the mean risk values estimated under each of the three spatial structures were similar; (2) the IRIS risk estimates were non-negative; (3) and their sums were equal to the original county risk. Although the disaggregation of cancer data on a small scale is somewhat arbitrary, in particular when it does not take into account secondary information to guide this downscaling, the approach should facilitate the analysis of the relationships between health data and the putative covariates (i.e., environmental, socioeconomic, behavioral or demographic factors) that are typically estimated for different spatial scales. These covariates can potentially be subsequently used as secondary information in the kriging, leading to more detailed risk maps at finer scales [42].
The other issue was the so-called modifiable areal unit problem (MAUP), for which different geographic scales can lead to inconsistent results for relationships analysis. For example, the mortality rate reported at (1) the county level requires an aggregated deprivation index at the same resolution, and this aggregation obscures the intra-county variation and thus the relationship and (2) the IRIS level, at which the disaggregation leads to a large variance in estimated risk. Exploratory methods, such as the univariate Moran’s can serve as indications of the potential effect of the MAUP in the study how relationships based on the homogeneity and heterogeneity of spatial data are affected by the study level and may affect the ability of the study to detect a relationship [43].
In this study, very similar results were obtained for the different spatial scales between:
  • pleural cancer mortality and the exposure indicator F1.
  • lip, oral cavity, pharynx cancer mortality and the SE index.
Whereas other studies on the relationship between heath and deprivation showed that the use of spatial representations other than the census tract produced different analytical results [35,44], the significant association between lip, oral cavity, pharynx cancer mortality and the SE index estimated using the county structure were stronger than they were under the IRIS structure.
This difference between the significant correlation coefficients is the result of aggregation because the aggregation level is known to increase the strength of correlation [20]. Future work could provide tools to exhaust all possible aggregations and generate empirical frequency distributions of the statistical estimates that could be used to evaluate the sensitivity of results to aggregation effects.
Based on the results obtained, we can confirm that the presence of this significant statistical association was likely not induced by the use of a particular geography. At the three spatial scales, the strongest correlation coefficients were found where low deprivation was associated with low lip, oral cavity and pharynx cancer mortality and where low environmental pollution was associated with low pleural cancer mortality.

5. Conclusions

This paper presents an approach for evaluating spatial relationships between health outcomes (mortality attributable to cancer) initially aggregated at the county level, district socioeconomic covariates, and exposure data modeled on a regular grid. The approach was illustrated using age-adjusted lip, oral cavity and pharynx, and pleural cancer mortality rates measured over the period 2000–2009 for the Picardy region. The deprivation index and trace metal exposure indicators were used as putative risk factors. For the different spatial scales, the strongest associations were found where low deprivation was associated with low lip, oral cavity and pharynx cancer mortality and where low environmental pollution was associated with low pleural cancer mortality. However, applying this approach to other areas, for other causes of death, or with other indicators always requires exploratory analysis to assess the role of the MAUP and downscaling health data in the study of the relationships that will allow decision-makers to develop interventions where they are the most needed.

Acknowledgements

The authors wish to acknowledge the financial support by the French Environment and Energy Management Agency ADEME and the French Picardy Region provided within the framework of the CIRCE project.

Author Contributions

Work presented here was conceived of, carried out and analyzed by Mahdi-Salim Saib, Julien Caudeville and Florence Carre. Olivier Ganry, Alain Trugeon and Andre Cicolella gave important suggestions, and supervised the study.All authors read, revised the manuscript and approved the final version.

Conflicts of Interest

The authors declare no conflict of interest

References

  1. Leclerc, A.; Chastang, J.F.; Menvielle, G.; Luce, D. Socioeconomic inequalities in premature mortality in France: Have they widened in recent decades? Soc. Sci. Med. 2006, 62, 2035–2045. [Google Scholar] [CrossRef]
  2. Melchior, M.; Goldberg, M.; Krieger, N.; Kawachi, I.; Menvielle, G.; Zins, M.; Berkman, L.F. Occupational class, occupational mobility and cancer incidence among middle-aged men and women: A prospective study of the French GAZEL cohort*. Cancer Causes Control 2005, 16, 515–524. [Google Scholar] [CrossRef]
  3. Challier, B.; Viel, J.F. Relevance and validity of a new French composite index to measure poverty on a geographical level. Rev. Epidemiol. Sante Publi. 2001, 49, 41–50. [Google Scholar]
  4. 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–2016. [Google Scholar] [CrossRef] [Green Version]
  5. Salmond, C.; Crampton, P.; Sutton, F. NZDep91: A New Zealand index of deprivation. Aust. N. Z. J. Public Health 1998, 22, 835–837. [Google Scholar] [CrossRef]
  6. Waller, L.A.; Gotway, C.A. Applied Spatial Statistics for Public Health Data; John Wiley and Son: New Jersey, NY, USA, 2004. [Google Scholar]
  7. Lawson, A.B. Tutorial in biostatistics: Disease map reconstruction. Stat. Med. 2001, 20, 2183–2204. [Google Scholar] [CrossRef]
  8. Kafadar, K. Choosing among two-dimensional smoothers in practice. Comput. Stat. Data Anal. 1994, 18, 419–439. [Google Scholar] [CrossRef]
  9. Besag, J.; York, J.; Mollie, A. Bayesian image restoration with two applications in spatial statistics. Ann. Inst. Stat. Math. 1991, 43, 1–59. [Google Scholar] [CrossRef]
  10. Goovaerts, P.; Gebreab, S. How does Poisson kriging compare to the popular BYM model for mapping disease risks? Int. J. Health Geogr. 2008, 7. [Google Scholar] [CrossRef]
  11. Oliver, M.R.; Webster, R.; Lajaunie, C.; Muir, K.R.; Parkes, S.E.; Cameron, A.H.; Stevens, M.C.G. Binomial cokriging for estimating and mapping the risk of childhood cancer. IMA J. Math. Appl. Med. Biol. 1998, 15, 279–297. [Google Scholar]
  12. Goovaerts, P.; Jacquez, G. Accounting for regional background and population size in the detection of spatial clusters and outliers using geostatistical filtering and spatial neutral models: The case of lung cancer in Long Island, New York. Int. J. Health Geogr. 2004, 3. [Google Scholar] [CrossRef] [Green Version]
  13. Goovaerts, P. Simulation-based assessment of a geostatistical approach for estimation and mapping of the risk of cancer. Geostatistics Banff 2005, 2, 787–796. [Google Scholar]
  14. Goovaerts, P. Detection of spatial clusters and outliers in cancer rates using geostatistics filters and spatial neutral models. Geostatistics Environ. Appl. 2005. [Google Scholar] [CrossRef]
  15. Goovaerts, P. Geostatistical analysis of disease data: Estimation of cancer mortality risk from empirical frequencies using poisson kriging. Int. J. Health Geogr. 2005. [Google Scholar] [CrossRef]
  16. Goovaerts, P. Analysis and Detection of Health Disparities Using Geostatistics and a Space-time Information System. The Case of Prostate Cancer Mortality in the United States, 1970–1994. Proceedings of GIS Planet 2005, Estoril, 30 May–2 June; 2005. Available online: http://www.biomedware.com/publications/51paper148pierregoovaerts.pdf (accessed on 31 March 2014).
  17. Goovaerts, P. Geostatistical analysis of disease data: accounting for spatial support and population density in the isopleth mapping of cancer mortality risk using area-to-point Poisson kriging. Int. J. Health Geogr. 2006. [Google Scholar] [CrossRef]
  18. Goovaerts, P. Geostatistical analysis of disease data: Visualization and propagation of spatial uncertainty in cancer mortality risk using poisson kriging and p-field simulation. Int. J. Health Geogr. 2006, 5. [Google Scholar] [CrossRef]
  19. Flowerdew, R.; Geddes, A.; Green, M. Behaviour of Regression Models under Random Aggregation. In Modelling Scale in Geographical Information Science; Tate, N.J., Peter, M.A., Eds.; Wiley: Chichester, UK, 2001; pp. 89–104. [Google Scholar]
  20. Openshaw, S.; Taylor, P.J. A Million or So Correlation Coefficients: Three Experiments on the Modifiable areal Unit Problem. In Statistical Methods in the Spatial Sciences; Pion: London, UK, 1979. [Google Scholar]
  21. Openshaw, S. The Modifiable Areal Unit Problem Concepts and Techniques in Modern Geography; Geo Books: Norwich, UK, 1984. Available online: http://qmrg.org.uk/files/2008/11/38-maup-openshaw.pdf (accessed on 31 March 2014).
  22. Riva, M.; Gauvin, L.; Apparicio, P.; Brodeur, J.M. Disentangling the relative influence of built and socioeconomic environments on walking: The contribution of areas homogeneous along exposures of interest. Soc. Sci. Med. 2009, 69, 1296–1305. [Google Scholar] [CrossRef]
  23. Briant, A.; Combes, P.P.; Lafourcade, M. Dots to boxes: Do the size and shape of spatial units jeopardize economic geography estimations? J. Urban Econ. 2010, 67, 287–302. [Google Scholar] [CrossRef]
  24. Caudeville, J.; Bonnard, R.; Boudet, C.; Denys, S.; Govaert, G.; Cicolella, A. Development of a spatial stochastic multimedia model to assess population exposure at a regional scale. J. Total Environ. 2012, 432, 297–308. [Google Scholar] [CrossRef]
  25. Caudeville, J.; Boudet, C.; Denys, S.; Bonnard, R.; Govaert, G.; Cicolella, A. Caractérisation des inégalités environnementales en Picardie fondée sur l’utilisation couplée d’un modèle multimédia et d’un système d’information géographique (in Frech). Environ. Risque. Sante 2011, 10. [Google Scholar] [CrossRef]
  26. Rey, G.; Jougla, E.; Fouillet, A.; Hémon, D. Ecological association between a deprivation index and mortality in France over the period 1997–2001: Variations with spatial scale, degree of urbanicity, age, gender and cause of death. BMC Public Health 2009. [Google Scholar] [CrossRef]
  27. Atlas de la Mortalité par Cancer en Picardie. (in French). OR2S Picardie: Amiens, France, 2007. Available online: http://www.or2s.fr/Portals/0/Autres%20sanitaire/Rap%20CIRCEnew.pdf (accessed on 31 March 2014).
  28. Monestiez, P.; Dubroca, L.; Bonnin, E.; Durbec, J.P.; Guinet, C. Comparison of Model Based Geostatistical Methods in Ecology: Application to Fin Whale Spatial Distribution in Northwestern Mediterranean Sea. In Geostatistics Banff; Leuangthong, O., Deutsch, C.V., Eds.; Kluwer: Dordrecht, The Netherlands, 2005; Volume 2, pp. 777–786. [Google Scholar]
  29. Anselin, L. Spatial Statistical Modeling in a GIS Environment. In GIS, Spatial Analysis, and Modelling; Maguire, D., Batty, M., Goodchild, M.F., Eds.; ESRI Press: Redlands, CA, USA, 2005; pp. 93–111. [Google Scholar]
  30. Brunsdon, C.; Fotheringham, A.S.; Charlton, M.E. Geographically weighted regression: A method for exploring spatial nonstationarity. Geogr. Anal. 1996, 28, 281–298. [Google Scholar] [CrossRef]
  31. Brunsdon, C.; Fotheringham, A.S.; Charlton, M. Spatial nonstationarity and autoregressive models. Environ. Plan. A 1998, 30, 957–973. [Google Scholar] [CrossRef]
  32. Brunsdon, C.; Fotheringham, A.S.; Charlton, M. Geographically weighted regression: Modelling spatial non-stationarity. J. R. Stat. Soc. Ser. D 1998, 47, 431–443. [Google Scholar] [CrossRef]
  33. Goovaerts, P. Medical geography: Apromising field of application for geostatistics. Math. Geosci. 2009, 41, 243–264. [Google Scholar] [CrossRef]
  34. Fotheringham, A.S.; Brunsdon, C.; Charlton, M.E. Geographically Weighted Regression: The Analysis of Spatially Varying Relationship; Wiley: New York, NY, USA, 2002. [Google Scholar]
  35. Cockings, S.; Martin, D. Zone design for environment and health studies using pre-aggregated data. Soc. Sci. Med. 2005, 60, 2729–2742. [Google Scholar] [CrossRef]
  36. Benach, J.; Yasui, Y. Geographical patterns of excess mortality in Spain explained by two indices of deprivation. J. Epidemiol. Community Health 1999, 53, 423–431. [Google Scholar] [CrossRef]
  37. Lorant, V. Inégalités socio-économiques de la mortalité dans les communes belges. Rev. Epidemiol. Sante Publ. 2000, 48, 239–247. [Google Scholar]
  38. Carstairs, V. Socio-economic Factors at Areal Level and Their Relationship with Health. In Spatial Epidemiology, Methods and Applications; Elliott, P., Wakefield, J., Best, N., Briggs, D., Eds.; Oxford University Press: Oxford, UK, 2000; pp. 51–67. [Google Scholar]
  39. Declercq, C.; labbe, E.; Poirier, G.; Lacoste, O. Inégalités Socio-spatiales de Mortalité Dans la Région Nord—Pas-de-Calais, France. 2004. Available online: http://orsnpdc.gapi.fr/etudes/129868.html (accessed on 31 March 2014).
  40. Peterson, J.T.; Greenberg, S.D.; Buffler, P.A. Non-asbestos-related malignant mesothelioma: A review. Cancer 1984, 54, 951–960. [Google Scholar]
  41. Glista-Baker, E.E.; Taylor, A.J.; Sayers, B.C.; Thompson, E.A.; Bonner, J.C. Nickel nanoparticles enhance platelet-derived growth factor-induced chemokine expression by mesothelial cells via prolonged mitogen-activated protein kinase activation. Am. J. Respir. Cell Mol. Biol. 2012, 47, 552–556. [Google Scholar]
  42. Gotway, C.A.; Young, L.J. A geostatistical approach to linking geographically aggregated data from different sources. J. Comp. Graph. Stat. 2007, 16, 115–135. [Google Scholar] [CrossRef]
  43. Nuckols, J.; Ward, M.H.; Jarup, L. Using geographic information systems for exposure assessment in environmental epidemiology studies. Environ. Health Perspect. 2004, 112, 1007–1015. [Google Scholar] [CrossRef]
  44. Parenteau, M.P.; Sawada, M.C. The modifiable areal unit problem (MAUP) in the relationship between exposure to NO2 and respiratory health. Int. J. Health Geogr. 2011, 10. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Saib, M.-S.; Caudeville, J.; Carre, F.; Ganry, O.; Trugeon, A.; Cicolella, A. Spatial Relationship Quantification between Environmental, Socioeconomic and Health Data at Different Geographic Levels. Int. J. Environ. Res. Public Health 2014, 11, 3765-3786. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110403765

AMA Style

Saib M-S, Caudeville J, Carre F, Ganry O, Trugeon A, Cicolella A. Spatial Relationship Quantification between Environmental, Socioeconomic and Health Data at Different Geographic Levels. International Journal of Environmental Research and Public Health. 2014; 11(4):3765-3786. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110403765

Chicago/Turabian Style

Saib, Mahdi-Salim, Julien Caudeville, Florence Carre, Olivier Ganry, Alain Trugeon, and Andre Cicolella. 2014. "Spatial Relationship Quantification between Environmental, Socioeconomic and Health Data at Different Geographic Levels" International Journal of Environmental Research and Public Health 11, no. 4: 3765-3786. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph110403765

Article Metrics

Back to TopTop