Engineering » Environmental Engineering » "Air Pollution", book edited by Vanda Villanyi, ISBN 978-953-307-143-5, Published: August 17, 2010 under CC BY-NC-SA 3.0 license. © The Author(s).

Chapter 15

Ambient Air Pollution, Social Inequalities and Asthma Exacerbation in Greater Strasbourg (France) Metropolitan Area: the PAISA Study

By Denis Bard, Olivier Laurent, Sabrina Havard, Severine Deguen, Gaelle Pedrono, Laurent Filleul, Claire Segala, Agnes Lefranc, Charles Schillinger and Emmanuel Riviere

Article top

Overview

Mapping deprivation in the Strasbourg (France) Metropolitan Area.
Figure 1. Mapping deprivation in the Strasbourg (France) Metropolitan Area.
Maps of age-standardized incidence ratios of emergency calls for asthma and of SAMU-confirmed asthma cases, across Strasbourg Metropolitan Area’s census blocks
Figure 2. Maps of age-standardized incidence ratios of emergency calls for asthma and of SAMU-confirmed asthma cases, across Strasbourg Metropolitan Area’s census blocks
Maps of ambient air pollutant concentrations (µg/m3) averaged by census block across the Strasbourg Metropolitan Area (2000–2005).
Figure 3. Maps of ambient air pollutant concentrations (µg/m3) averaged by census block across the Strasbourg Metropolitan Area (2000–2005).
Scatterplots of the rates of emergency calls for asthma according to deprivation of census blocks, for different age groups (on the left, the most privileged census blocks, on the right the most disadvantaged ones). SIR: Standardized Incidence Ratio.
Figure 4. Scatterplots of the rates of emergency calls for asthma according to deprivation of census blocks, for different age groups (on the left, the most privileged census blocks, on the right the most disadvantaged ones). SIR: Standardized Incidence Ratio.
Odds ratios between asthma calls and pollutant concentrations*, for various lag times (Strasbourg Metropolitan Area, 2000-2005).
Figure 5. Odds ratios between asthma calls and pollutant concentrations*, for various lag times (Strasbourg Metropolitan Area, 2000-2005).
Log odds ratio for emergency asthma calls per 10 µg/m3 increase in ambient pollutant concentrations, according to stratum of socioeconomic deprivation, Strasbourg, France, 2000–2005. Stratum 1 was the least deprived, and stratum 5 was the most deprived. PM10: particulate matter less than 10 µm in aerodynamic diameter; NO2: nitrogen dioxide; O3: ozone. Bars, 95% confidence interval.
Figure 6. Log odds ratio for emergency asthma calls per 10 µg/m3 increase in ambient pollutant concentrations, according to stratum of socioeconomic deprivation, Strasbourg, France, 2000–2005. Stratum 1 was the least deprived, and stratum 5 was the most deprived. PM10: particulate matter less than 10 µm in aerodynamic diameter; NO2: nitrogen dioxide; O3: ozone. Bars, 95% confidence interval.
Odds ratios between asthma calls and pollutants*, for five strata of deprivation levels† and different age groups (Strasbourg Metropolitan Area, 2000-2005).
Figure 7. Odds ratios between asthma calls and pollutants*, for five strata of deprivation levels† and different age groups (Strasbourg Metropolitan Area, 2000-2005).
Odds ratios between asthma calls and pollutants in the 136 statistical units retained, ranked according to deprivation (from the least deprived on the bottom to the most deprived at the top). Strasbourg Metropolitan Area, 2000-2005.
Figure 8. Odds ratios between asthma calls and pollutants in the 136 statistical units retained, ranked according to deprivation (from the least deprived on the bottom to the most deprived at the top). Strasbourg Metropolitan Area, 2000-2005.
Odds ratios for an increase of 10 µg/m3 in the pollutants concentration, subjects <40 years old; for various lags.
Figure 9. Odds ratios for an increase of 10 µg/m3 in the pollutants concentration, subjects <40 years old; for various lags.
Interactions by socioeconomic level in those aged < 40 years: associations for the five different strata of socioeconomic levels.
Figure 10. Interactions by socioeconomic level in those aged < 40 years: associations for the five different strata of socioeconomic levels.

Ambient air pollution, social inequalities and asthma exacerbation in Greater Strasbourg (France) metropolitan area: The PAISA study

D. Bard1, O. Laurent1, 2, S. Havard1, 3, S. Deguen1, G. Pedrono4, L. Filleul5, C. Segala4, A. Lefranc6, C. Schillinger7 and E. Rivière7

1. Introduction

The socio-economic status (SES) of populations has an influence on the incidence or mortality rates of numerous health outcomes, among which respiratory diseases (Prescott et al., 2003; Ellison-Loschmann et al., 2007). Considering asthma, the possible contribution of SES to overall prevalence –regardless of asthma severity-, remains controversial in industrialized countries. Several studies indicate that allergic asthma is more prevalent in more well-off populations whereas the non-allergic forms of asthma are more common in the deprived ones (Cesaroni et al., 2003; Blanc et al., 2006). On the other hand, severe asthma whatever its etiology appears to be more frequent in the latter populations, as compared to the more affluent (Basagana et al., 2004). Risk factors for exacerbations (e.g., passive smoking (Wright Subramanian, 2007), psychosocial stress (Gold & Wright, 2005), cockroach allergens (Kitch et al., 2000), and suboptimal compliance with anti-inflammatory medication (Gottlieb et al., 1995)) are generally more common among people with asthma and low SES than their better-off counterparts. These observations support the hypothesis that some factors more present in deprived populations contribute to asthma exacerbation (Mielck et al., 1996).

One candidate factor is ambient air pollution, a well-demonstrated risk factor for asthma (Samet & Krewski, 2007). Epidemiologic studies report short-term associations (latency of several hours to several days) between exposure to air pollutants (particulate matter, nitrogen dioxide [NO2], sulfur dioxide, and ozone [O3]) and various asthma attack indicators routinely available: hospitalization (Sunyer et al., 1997), emergency department (ED) visits (Romeo et al., 2006), calls to mobile emergency medical service networks (Medina et al., 1997) and visits to doctors’ offices for asthma attacks (Hajat et al., 1999). Associations are also reported with indicators measured directly in individuals with asthma: respiratory symptoms and consumption of short-acting -agonist (SABA) drugs (von Klot et al., 2002; Schildcrout et al., 2006).

It should be stressed that in such studies, be they time series (Sunyer et al., 1997), case-crossover (Jalaludin et al., 2008), or panel (Ward & Ayres, 2004; Romeo et al., 2006) studies, investigators average citywide ambient pollutant concentrations for estimating exposure, although these concentrations often vary spatially and strongly within cities (Jerrett et al., 2005a), with few exceptions (Tolbert et al., 2000; Neidell, 2004; Erbas et al., 2005). Failure to consider these spatial variations at the sub-city level can lead to exposure misclassification and resulting bias (Jerrett et al., 2005a). Conversely, taking these spatial variations into account in studying the relations between air pollution and health would allow more rigorous measurement of the magnitude of these relations (Jerrett et al., 2005a). Another problem inherent in ecologic studies is spatial autocorrelation, which expresses the non-independence of geographic observations (Haynes et al., 2001). With very rare exceptions (Jerrett et al., 2001; Buzzelli et al., 2003), studies of environmental equity do not consider this phenomenon, although it is essential to avoid violating the hypotheses that underlie the application of statistical models, and may thus lead to erroneous conclusions about associations.

Another key issue is identification of the populations most susceptible to air pollution and the precise measurement of exposure-response relations in these subsets (Maynard et al., 2003; Samet & Krewski, 2007). Several studies show that some specific populations, such as children, the elderly, and people with chronic diseases, including people with diabetes and chronic obstructive pulmonary disease, are more susceptible to air pollution than the general population (Annesi-Maesano et al., 2003; Samet & Krewski, 2007). Some suspect that socioeconomic deprivation influences the relations between air pollution and health (O'Neill et al., 2003). Many investigators have tested this hypothesis with regard to mortality, and their overall results suggest the existence of larger relative risks for more deprived populations (Laurent et al., 2007). Substantially fewer researchers have tested this hypothesis with regard to asthma attacks (Norris et al., 1999; Lin et al., 2004; Neidell, 2004), and the few that have done so have considered socioeconomic indicators measured at resolutions ranging from the individual level (Neidell, 2004; Kim et al., 2007) to the level of geographic residence areas of various sizes (Nauenberg & Basu, 1999; Norris et al., 1999; Lin et al., 2004; Son et al., 2006; Kim et al., 2007). Whether or not a socioeconomic indicator is an effect modifier, however, may well depend on the resolution at which it is measured (Laurent et al., 2007). To our knowledge, only two studies have tested the influence of socioeconomic indicators measured at the level of small areas (Canadian enumeration areas (Lin et al., 2004)and US zip codes (Neidell, 2004) on the short-term relations between air pollution and asthma attacks. Disadvantaged people may be more susceptible to asthma attacks when they face additional environmental insults, such as ambient air pollution (Maynard et al., 2003). Nonetheless, the few studies (Nauenberg & Basu, 1999; Annesi-Maesano et al., 2003; Maynard et al., 2003; O'Neill et al., 2003; Lin et al., 2004; Laurent et al., 2007) that have tested this hypothesis have produced divergent results. This may be because the distribution according to SES of factors modulating the relations between air pollution and asthma morbidity differs between the study settings or because of lack of statistical power.

In the absence of individual data, which are not generally routinely available, ecological (or contextual) measures of SES are frequently used to describe health inequalities. Although some epidemiological studies are based on only one socioeconomic indicator (income, educational level, or occupation), e.g. (Kunst et al., 1998; Finkelstein et al., 2003), SES is usually recognized as complex and multidimensional, integrating different components that may be either material (e.g., housing conditions, income, or occupation), social (e.g., social position or isolation, or family support) or both (Folwell, 1995; Braveman et al., 2005).

Area-based deprivation indices for the measurement of the economic or social disadvantages of urban areas were proposed in the 1980s (Townsend, 1987). Initially designed for health care planning and resource allocation, they have been recently used to evaluate and analyse health inequalities (Carstairs, 1995; Niggebrugge et al., 2005; Eibner Sturm, 2006). These measures, including Townsend’s, Carstairs’, and Jarman’s indices, as well as the more recent Index of Multiple Deprivation, combine contextual indicators such as unemployment rate or proportions of overcrowded or of non-owner-occupied households (Jarman, 1983; Townsend et al., 1988; Carstairs & Morris, 1991 ; DETR, 2000). Since the end of the 1990s, numerous other area-based deprivation indices have emerged – in the United States (Singh, 2003; Eibner & Sturm, 2006; Messer et al., 2006), Canada (Pampalon & Raymond, 2000), New Zealand (Salmond et al., 1998), Japan (Fukuda et al., 2007), Italy (Tello et al., 2005), Spain (Benach & Yasui, 1999), and Belgium (Lorant, 2000).

Hospitalizations or emergency department visits reflect severe asthma exacerbation episodes. As to increase sensitivity, alternative indicators integrating less severe events may be used, such as telephone calls to medical services or asthma exacerbation specific drugs.

Naureckas et al. (Naureckas et al., 2005) showed that SABA sales are valuable indicators of asthma morbidity and its temporal variability. These sales are especially good predictors of the risk of emergency department visits and hospitalization for asthma attacks in the days immediately afterwards (Naureckas et al., 2005). Moreover, they generally reflect less severe morbidity states than those requiring emergency department visits or hospitalizations and are therefore very sensitive indicators. Finally, they are very common, likely to involve large numbers of individuals, even within small geographic units and short observation periods (Pitard et al., 2004; Naureckas et al., 2005). While SABA drugs are the treatment of choice for asthma attacks (Global Initiative for Asthma, 2007), they are also prescribed for acute exacerbation of chronic obstructive pulmonary disorders (COPD) (Urbano & Pascual, 2005). Nevertheless, this disease is rare in people less than 40 years old (Viegi et al., 2007); within this age group, SABA prescriptions and sales are therefore very specific for asthma attacks.

The PAISA study objectives were to investigate the short-term relations between air pollution estimated at the level of small geographic areas and asthma attacks and to test the influence of socioeconomic deprivation, also measured at this level, on these relations. Two different, although related, indicators were tested, that is i) calls for emergency services for asthma exacerbation and ii) SABA drugs deliveries.

2. Methods

2.1. Setting and statistical unit

The setting of our study was the Strasbourg metropolitan area (SMA), gathering municipalities (316 km2) with a population of about 450,000 inhabitants located in the Bas-Rhin district (or ‘dé́partement’, an administrative subdivision) in the Alsace region of North-Eastern France. The statistical small-area unit used was the French census residential block (Ilots Regroupés pour l’Information Statistique-IRIS), a submunicipal division devised by the National Institute for Statistics and Economic Studies (INSEE). This unit is the smallest geographic area in France for which socioeconomic and demographic information from the national census is available. It is, in terms of population size, intermediate between US census tracts (about 4,000 inhabitants) and US census block groups (about 1,000 inhabitants) with an average number of inhabitants of 2,000. The division of neighborhoods into census blocks takes into account the physical obstacles that may break up urban landscapes (important traffic arteries, bodies of water, green spaces, etc.) and aims to maximize their homogeneity in population size, socioeconomic characteristics, land use and zoning. The SMA is subdivided into 190 blocks, with areas ranging from 0.05 km2 to 19.6 km2 (median values of 0.45 km2 for the entire SMA and 0.29 km2 for the city center). Fifteen blocks covering a very small population (250 people, altogether 0.2 percent of the total population) were removed from the study because of confidentiality rules in force in France.

2.2. Socioeconomic Status Indicator

The demographic and socioeconomic data come from the 1999 national population census, conducted by INSEE. The census database is structured into different domains, including employment, family and household, educational level, housing, immigration status, and income; together these regroup a collection of very diverse quantitative variables, such as proportions of unemployed people, foreigners, blue-collar workers, households without cars, people aged 15 years or older with general or vocational maturity certificates, etc. Accordingly, we selected 52 variables from the available data, endeavoring to apply theoretical concepts of deprivation (Krieger et al., 1997), relying on the indicators most often used in the literature (Morris & Carstairs, 1991; Pampalon & Raymond, 2000; Challier & Viel, 2001; Jordan et al., 2004), and using the same definitions as INSEE most often uses in its studies. All but 2 of the 52 variables were proportions (the exceptions were mean number of people per room and median income per consumption unit). Some variables, intentionally redundant, were introduced into the analysis to determine statistically the most discriminating for characterizing deprivation (e.g., unemployment among total labor force, among men, among women).

Constructing the index

We developed the socioeconomic deprivation index by principal component analysis (PCA) of the 52 variables. All variables were first standardized (i.e. centered and reduced) to remove the influence of different units of measure and thus to give them the same weight in the analysis. This statistical technique permits the synthesis of information contained in a great number of variables by constructing new and independent synthetic variables – the principal components. These new variables are linear combinations of some of the initial variables, defined by a maximum variance. The principal advantage of this method is that it makes it possible to consider relations between the variables, by attributing to each a weight that takes the relations into account (that is, treating them as coefficients of the linear combinations).

In our study, to facilitate the use of the deprivation index in a public health context, we chose to construct a single index for all the blocks to maximizing the variance of the first principal component. Several consecutive PCA were thus performed to conserve only the variables most strongly correlated with the first component and contributing most to its construction. If two redundant variables were conserved, only the one most closely correlated with the first component was finally retained (and in each case, that turned out to be the variable not stratified by sex). Finally, this first component, used as our socioeconomic deprivation index, was calculated for each block, as the linear combination of the final variables retained.

To obtain the most discriminant cartographic representation of socioeconomic disparities, the blocks were grouped into five classes of index values, approximately equivalent to the quintiles (Beguin & Pumain, 1994). The first class comprised the least deprived blocks and the fifth the most deprived blocks. The choropleth mapping of the index was performed with the Geographic Information System ArcView™ v. 9.1 (ESRI Inc., Redlands, CA, USA).

Index validity

According to Coste et al. (1995) (Coste et al., 1995), several criteria make it possible to judge the validity of composite measurement scales (in our case, a socioeconomic deprivation index), in particular, its content validity and its construct validity (Fermanian, 1996).

Content validity is assessed by the extent to which the items composing the scale or index are relevant to and representative of all of the possible items that can describe the phenomenon measured. To select the initial subset of variables for analysis, we carefully ensured that each census domain (income, job, housing, family and household, etc.) was represented by at least one variable, and again at the end we verified this representativeness for all of the final variables.

Construct validity verifies that from a theoretical point of view the instrument is associated with the concept it is supposed to measure. It tests the relations between the index variables, both internally (factorial validity and internal consistency) and externally (convergent validity), and thus identifies more precisely the real meaning of the concept measured by the indicator. The internal validity of our index was estimated with the results of the PCA and Cronbach’s alpha coefficient. This coefficient judges how the variables retained by the PCA measure a one-dimensional concept – here, socioeconomic deprivation. The closer the coefficient value is to one, the better the verification that the index variables are homogenous. The convergent validity was tested using Pearson’s correlation coefficients with the most widely used British indices (those of Carstairs and Townsend) (Townsend et al., 1988; Carstairs & Morris, 1991 ). These indices combine four variables, three common to both: (1) proportion of unemployed, (2) proportion of households without a car, and (3) proportion of over- crowded households. Townsend’s index also includes the concept of rental versus owner-occupied housing (proportion of non-owner-occupied housing) while Carstairs’ index considers a variable related to the household head’s place in the British social scale (proportion of households in which the householder belongs to social classes 4 or 5). These indices were compiled at the census block level by adapting or translating some British variables when they had no direct correspondence in the French census database (for example, the variable ‘‘proportion of households in which the householder belongs to social class 4 or 5’’ was translated for the French census database as ‘‘the proportion of blue-collar workers’’).

Index reliability

To test the robustness of our index, the same methodological procedure as that used for the SMA (the same 52 initial variables, same algorithm decision for the final choice of variables, the same tests of validity) was conducted for another French metropolitan area at the resolution of the census block – the Lille metropolitan area (LMA). The LMA was selected because its demographic and socioeconomic characteristics are quite distinct from those of the SMA (e.g., 489 census blocks for a total population around 1,100,000 inhabitants). The objective of this stage was to demonstrate that similar results could be obtained in a different part of France and to show that the variables selected for the construction of the index were independent of the study area and characterized socioeconomic level in as discriminating a manner as possible.

Application

Census blocks were divided into five deprivation categories according to their index value; the first category comprised the most privileged blocks and the fifth the most deprived ones. More details can be found in (Havard et al., 2008).

Calls to emergency medical services for asthma exacerbation

The two mobile emergency and healthcare networks operating in the SMA (Service d’Aide Médicale d’Urgence (SAMU) and SOS Médecins provided data about emergency telephone calls made to physicians for asthma attacks. The Bas-Rhin emergency medical services unit (SAMU 67) is a public service that coordinates pre- hospital emergency medical services. SOS Médecins of Strasbourg is a private service providing emergency general medical care.

We included each call regarding an asthma attack that reached either SAMU or SOS Médecins from January 1, 2000, to December 31, 2005. The two databases were merged, duplicates being excluded. Each call was geo-coded to the census block where the patient was located at the time of the call, based on the postal address. Only 2 percent of the calls could not be geocoded; they were excluded from the analysis.

SABA sales

The data about SABA sales to people living in the SMA during 2004 came from four health insurance funds: the regional union of health insurance funds (Union Régionale des Caisses d’Assurance Maladie [URCAM], for salaried workers and very low-income families), agricultural social insurance fund (Mutualité Sociale Agricole [MSA], for farmers and similar workers), the fund for self-employed workers (Régime Social des Indépendants [RSI]), and the Lorraine students’ insurance fund (Mutuelle Générale des Etudiants de Lorraine [MGEL]). These funds cover >90% of the local population (Com-Ruelle et al.). In 2004, these funds either paid or reimbursed (fully or partially) all SABA sales prescribed by a physician. All records of SABA sales in pharmacies were extracted from the databases of these four insurance funds by requests for code R3A4 of the European Pharmaceutical Market Research Association nomenclature (Tollier et al., 2005). The following information was furnished for each sale: date, age group (0 to 9, 10 to 19, 20 to 39 years), sex, and census block of residence of the person for which SABA were prescribed. All sales were combined into a single file for analysis.

Ambient air pollution

Hourly ambient concentrations of particulate matter less than 10 m in aerodynamic diameter (PM10), nitrogen dioxide, sulfur dioxide, and ozone were modeled by the local air quality monitoring association (Association pour la Surveillance et l’Etude de la Pollution Atmosphérique en Alsace, Schiltigheim) for each block during the entire study period (Bard et al., 2007). Modeling was conducted with Atmospheric Dispersion Modeling System (ADMS-Urban), a deterministic Gaussian dispersion model (McHugh et al., 1997) that integrates emissions inventories, meteorological data (wind direction and speed, temperature and cloudiness, supplied by Météo France, the French meteorological service), and background pollution measurements as input parameters. A limitation of the model was the estimation of Monin-Obukhov length from cloudiness measurements, but the correlation coefficients for the modeled and effectively measured ambient concentrations were 0.73 for PM10, 0.87 for nitrogen dioxide, 0.84 for ozone. However, this correlation was only 0.06 for sulfur dioxide. Therefore, sulfur dioxide was no longer considered in the analyses.

Potential confounders

Daily meteorological variables (temperature, atmospheric pressure, and relative humidity) were obtained from Météo France, and daily pollen counts were obtained from the National Network of Aerobiological Surveillance (Thibaudon, 2004). Weekly influenza case counts came from the Sentinelles network (Flahault et al., 2006) of the National Institute of Health and Medical Research.

Statistical analyses

As some census blocks featured a relatively small population size, yielding unstable rates, we carried out an empirical Bayesian smoothing using STIS software (Space Time Intelligence System 2007 (Ann Arbor, MI, USA: Terraseer), except for drug sales, where smoothing was not warranted in view of the much greater number of events per census block (Marshall, 1991). Positive and statistically significant spatial autocorrelations were found for our SES index (Moran’s I = 0.54, p<0.01), rates and standardized incidence ratios (SIR) for asthma calls (I from 0.67 to 0.77 according to age group, p<0.01). Associations between SES and health event rates or SIR were quantified by Pearson’s correlation test and their significance assessed using a modified Student’s t-test as proposed by Clifford et al., (Clifford et al., 1989) as to accounting for spatial autocorrelation of neighboring census blocks variables.

Associations between asthma calls and air pollution were assessed with case-crossover models (Maclure & Mittelman, 1991), which are similar to case-control models, except that in the former, the subjects serve as both cases and controls, depending on when they are considered. The subject serves as a case on the day of the health event and a control on days without any health event. Control days were defined according to a monthly time-stratified design (Janes et al., 2005). For an asthma call occurring on a given weekday (e.g., Monday), control days were chosen as the same days of the week throughout the rest of the month (thus, three or four days; here, the other Mondays of the month). Conditional logistic regression was employed for analyses. The statistical significance of the odds ratios was tested by means of a two-sided 2 test at the 5 percent level.

Base models

We first analyzed all calls or SABA sale, without differentiating them according to socioeconomic deprivation. Associations between health events and ambient air pollution concentrations modeled by census block were estimated, adjusting for holidays, meteorological variables (daily maximum temperature, maximum atmospheric pressure, and mean relative humidity), influenza epidemics, and pollen counts.

We tested the influence of air pollution indicators averaged on the day of the call (lag 0) and then averaged on the day of the call and the 1–5 previous days (lag 0–1 to 0–5) or up to 10 days in the case of SABA sales. The daily air pollution indicator considered for PM10, nitrogen dioxide, and sulfur dioxide was the 24-hour average concentration, and for ozone it was the maximum daily value of the 8-hour moving average. The analysis for ozone considered only health events occurring between April 1 and September 30 of each year, because of the very low concentrations of this pollutant in winter.

These associations were assessed for cases of all ages and then for groups of cases aged 0–19, 20–64, and >64 years, except for SABA sales (below 40). These groups were defined according to the guidelines of the French Data Protection Authority, to ensure the confidentiality of the health data geocoded by census blocks.

Testing the influence of deprivation

Interactions with socioeconomic deprivation were tested for the lag times for which the strongest associations were observed in the basic models.

Deprivation was introduced first as a discrete variable. The SMA population was divided into five strata with contrasting deprivation levels, according to the quintiles of the distribution of the deprivation index. An odds ratio for the relation between asthma calls and pollution was estimated for each stratum. The heterogeneity of the odds ratios between these strata was assessed by means of a two-sided 2 test at the 5 percent level (Atkinson et al., 2001). These analyses were conducted for the same age groups as in the base models.

In an alternative method, deprivation was introduced as a continuous variable. For that purpose, a case-crossover model was constructed for each of the 174 census blocks in which asthma calls occurred during the study period. The models did not converge in the census blocks with fewer than 11 asthma calls. These census blocks were treated as follows. Those adjoining neighboring census blocks of a similar deprivation level (differences of deprivation index values less than one 10th of the deprivation index scale) were merged with them. The deprivation index value attributed to each newly created geographic unit was the mean of the deprivation index values from the merged blocks, weighted by their respective population counts. Census blocks (n = 26) that bordered no neighboring census block of a similar deprivation level (differences of deprivation index values one 10th or more of the deprivation index scale) were excluded from these analyses.

The heterogeneity of the odds ratios estimated in the 136 geographic units so defined was assessed as described above (Atkinson et al., 2001). The influence of the deprivation level of the geographic units on the odds ratios was tested by linear regression (weighted by the inverse of the variance of these odds ratios) in both fixed- and random-effects models. These analyses were conducted for PM10, nitrogen dioxide, and sulfur dioxide for cases of all ages but were not feasible for ozone (April–September) and specific age groups, because of the small numbers of calls in most census blocks.

All analyses were carried out using SAS statistical software v 9.1 and 8.2 (SAS Institute, Cary, NC, USA).

3. Results

3.1. Socioeconomic Status Indicator

Our deprivation index was defined by the PCA as the first principal component, which explained 66% of the total variance of the model, while the second explained only 17%. The impossibility of drawing any residual structure from the second component justifies the conclusion that all of the discriminating information useful for characterizing socioeconomic deprivation was in fact captured by the first component. Our index is composed of 19 socio-economic variables that describe different dimensions of socioeconomic deprivation. Seventeen can be qualified as primarily material variables and only two directly refer to social dimension of the deprivation (that is, single-parent family and foreign population). Although several of the variables retained belong to the same domain, we consider that they provide complementary information for characterizing the deprivation. Our index opposes the concepts of material and, to a lesser extent, social deprivation (it is significantly and positively correlated with variables relative to unemployment, immigration, social isolation, household overcrowding, low educational level, and low income) to those of well-being and economic comfort (it is significantly and negatively correlated with variables relative to job stability, high educational level, better living conditions, and high income). Moreover, these variables had differing relative influences in the construction of the index (e.g., proportion of unemployed people contributed to 6.72% of the variance and the percentage of households without cars to only 3.91%).

Table I presents the mean values of the index and of the 19 variables by deprivation class. The first class (C1) is characterized by the lowest mean index value and includes the most privileged blocks, defined by the best living conditions (low percentages of the variables describing socioeconomic deprivation and high percentages of variables illustrating material comfort). At the other end, the fifth class (C5) is characterized by the highest mean index value and includes the most deprived blocks. The blocks of this class are characterized by more socioeconomic deprivation, a greater cumulative lack of material and social resources (the highest percentages of variables describing deprivation and the lowest percentages of variables related to positive living conditions). For example, on average, only 1.28% of homes in the C1 blocks were low income subsidized housing while this percentage exceeded 75% for the C5 blocks.

VariablesCensus domainDeprivation categories
C1C2C3C4C5
Blue collar workers in the labor forceJob17.40 15.80 17.2021.7036.40
Primary residences that are houses or farmsHousing61.4028.607.105.803.60
Primary residences that are multiple dwelling unitsHousing37.1069.5090.2090.9095.10
Households without a car Housing8.8016.0027.6035.8036.50
Households with 2 carsHousing48.9032.4018.9014.9014.00
People with insecure jobs in the labor forceJob 8.3010.4013.7016.9017.40
People aged 15 yrs with general or vocational maturity certificatesEducation13.3012.2011.909.506.60
People aged 15 yrs with at least a lower postsecondary educationEducation10.9010.3010.308.103.50
People aged "/ 15 yrs with no more than a completed elementary educationEducation9.809.8010.8016.2033.50
Non owner-occupied primary residencesHousing28.2048.90 72.1080.3091.30
People with stable jobs in the labor forceJob76.4072.0068.2063.1054.90
Subsidized housing among all primary residencesHousing1.306.809.3027.3075.60
Single-parent familiesFamilies and households7.7010.3013.2015.7023.60
Median income/consumption unit (in Euros/yr) c Income19,97719,81917,12014,2708,402
Primary residences with more than 1 person/roomHousing3.20 4.305.907.9021.40
Foreigners in the populationImmigration status3.206.207.9012.3022.0
Mean no. people/room c Housing0.600.600.64 0.690.86
Unemployed people in the labor forceJob5.307.5010.2013.5024.20
People unemployed "/ 1 yr in the labor force Job2.403.20 4.405.8011.70
Deprivation index--1.091.752.252.673.44
a All variables are percentages (%), unless otherwise stated; p value for trend < 0.01 for all variables. b C1 is the least deprived category, C5 is the most. c Not a percentage .

Table 1.

Deprivation Index and mean values of its 19 variables, according to deprivation categoriesa.

Conversely, nearly 50% of the households in C1 blocks had two cars or more compared with 14% for the C5 blocks (Table I).

Our index thus demonstrates the existence of a socio-economic gradient in the SMA, along which we observe a progressive degradation of level and quality of life between the first and fifth classes (p < 0.01). The different variables used to construct the index sometimes vary significantly in the same direction as the gradient, and sometimes in the opposite direction (Table I), depending on the aspect (positive or negative) of SES that they express (p < 0.01).

This gradient appears clearly in the cartographic representation (Figure 1). The most affluent blocks are located on the outskirts of the study area and constitute a peri-urban ring around Strasbourg. These are principally blocks with a relatively low population, comprising rural municipalities. Conversely, socioeconomic deprivation is accentuated as we approach Strasbourg and reaches its maximum at the centre of the metropolitan area. The most disadvantaged blocks are principally those in Strasbourg and the inner suburbs with a relatively high population and comprising the urban municipalities.

Validity tests

Content validity. The 19 variables retained cover all of the census domains and thus insure that our index has fully integrated the multiple socioeconomic aspects of deprivation (Table I).

Construct validity. The high Cronbach’s coefficient (0.92) supports the hypothesis that these variables measure a one-dimensional concept – socioeconomic deprivation. Convergent validity is verified by very high correlation coefficients with the indices of both Townsend (0.97; p<0.01) and Carstairs (0.96; p < 0.01).

Reliability

Following the same process, 20 variables were retained for Lille metropolitan area, 18 of them were the same as for SMA and they explained the same proportion of variance (66%). Only three variables differed between the two indices: the percentages of the self-employed and of foreign immigrants were retained only for the Lille index and the percentage of blue-collar workers only for Strasbourg. Mapping the index for LMA again illustrates the strong spatial heterogeneity of SES between urban centers and outer suburbs. Finally, as for Strasbourg, all the tests of the index conclusively demonstrated its validity (i.e., the Lille index had a high Cronbach’s coefficient and was very highly correlated with both British indices; data not shown).

Calls to emergency medical services for asthma exacerbation

There were 4,682 calls for asthma (SAMU and SOS Médecins combined) and 1,173 SAMU-confirmed asthma cases (from 1,578 calls given to SAMU - some but not all initially reported to be for asthma). Cases finally diagnosed as asthma by participating physicians were treated as such in our analyses, regardless of the initial reason for call. We observed no trend in the annual number of health events.

Figure 2 shows the results of empirical Bayesian smoothing of calls for asthma exacerbation SIRs.

Maps of ambient air pollutant concentrations averaged by census block across the SMA (for the period 2000–2005) appear on Figure 3.

As expected for a small-area study, high and statistically significant spatial autocorrelation was found for the deprivation index (I =0.54, p<0.01) and for health event indicators (Table II). These are mapped respectively in Figure 4 and 5.

The age-specific rates and SIRs of calls for asthma increased from the least towards the most deprived blocks (Figure 4). Elevated (R≥0.65, p<0.01) correlations were observed between the deprivation index and both the SIRs and rates of calls for the groups aged 20 years or older (Table II). Correlations were lower but still statistically significant for the younger groups.

Results were similar for SAMU-confirmed asthma cases (Table II). Correlations with the deprivation index were elevated (R ≈ 0.62-0.63, p<0.01) for those aged 20 years or older, and lower (but still significant) for the younger groups.

Table III presents the distribution of asthma calls and air pollutant concentrations for the entire SMA and for the five socioeconomic deprivation strata. Overall, there were 4,677 usable asthma calls made during the study period. As we previously reported, the numbers of calls increased from the least deprived stratum to the most deprived, despite their similar population denominators.

media/image1.jpg

Figure 1.

Mapping deprivation in the Strasbourg (France) Metropolitan Area.

All calls for asthma Asthma cases confirmed by SAMU
Age groupnumber of callsMean number of calls per block% of all callsAnnual rate per 1,000 inhabitantsMoran's indexCorrelation with deprivation indexnumber of casesMean number of cases per block% of all SAMU-confirmed casesAnnual rate per 1,000 inhabitantsMoran's indexCorrelation with deprivation index
0-95393.111.51.7600.67 §0.52 †1350.811.50.4370.71 §0.26 †
10-19 4162.48.81.2080.74 §0.47 †1200.710.20.3490.68 §0.20 ‡
20-399515.420.31.0370.75 §0.65 †2491.421.20.2720.64 §0.63 †
40-6411386.524.31.4640.70 §0.70 †3772.1132.10.4850.71 §0.62 †
65 +16379.3354.5890.75 §0.68 †2921.724.90.8190.66 §0.63 †
Total *468226.81001.7520.77 §0.77 †11736.71000.4420.73 §0.66 †
* Moran's index and correlations with SES index calculated for the age-standardised incidence ratio (SIR) § p<0.01, estimated by a standard t -test † p<0.01 estimated by Clifford, Richardson and Hemon's test ‡ p<0.05 estimated by Clifford, Richardson and Hemon's test

Table 2.

Summary of health data for the Strasbourg metropolitan area and their correlations with deprivation index.

media/image2.jpg

Figure 2.

Maps of age-standardized incidence ratios of emergency calls for asthma and of SAMU-confirmed asthma cases, across Strasbourg Metropolitan Area’s census blocks

Deprivation stratum † PopulationCalls, all agesCalls, age 0-19Calls, age 19-64Calls, age over 64PM 10 ‡SO 2 ‡NO 2 ‡O 3 ‡
Stratum 1 80,91731370125118217.630.263.4
Stratum 292,53478415930032522.18.63558.6
Stratum 393,39295916938440623.29.63955
Stratum 4105,3671,243226574443239.338.355.5
Stratum 574 6951,37833070334522.58.835.758
Overall446,9054,6779542,0861,63722.68.936.057.7
† Stratum 1 is the least deprived, and stratum 5 the most deprived ‡ PM 10 , particulate matter less than 10 micrometers in aerodynamic diameter; NO 2 , nitrogen dioxide; SO 2 , sulfur dioxide. Concentrations for these three pollutants are averaged for period 2000-2005, in microgram per cubic meter § O 3 , ozone. Numbers reported are maximum 8-hour daily concentrations, averaged for summer months ( April 1 to September 30) of period 2000-2005, in microgram per cubic meter

Table 3.

Distribution of population, asthma calls and air pollutant concentrations (Strasbourg Metropolitan Area, 2000-2005).

Base models

Figure 5 shows the odds ratios for the relations between asthma calls (all ages) and pollutant concentrations averaged on the day of the call (lag 0) and the previous 1–5 days (lag 0–1 to lag 0–5). The highest positive associations were observed for lag 0–1 for PM10 (for a 10-µg x m-3 increase, odds ratio (OR) = 1.035, 95 percent confidence interval (CI): 0.997, 1.075) and sulfur dioxide (OR = 1.056, 95 percent CI: 0.979, 1.139) and for lag 0 for nitrogen dioxide (OR = 1.025, 95 percent CI: 0.990, 1.062). The association for PM10 approached statistical significance (p= 0.07). For ozone, no association was observed (OR = 0.998, 95 percent CI: 0.965, 1.032).

Table IV presents the results of analyses by age group for these same lags. For PM10 (lag 0–1), associations were stronger for people younger than age 20 years and older than age 64 years as compared with people of all ages. For nitrogen dioxide (lag 0), higher odds ratios were observed in people over age 64 years than in people of all ages. For ozone, positive associations were observed only in persons older than 64 years.

Testing the influence of deprivation

Deprivation as a discrete variable. Figure 6 presents the odds ratios for the five deprivation strata among people of all ages. The odds ratios show no clear trend according to deprivation level for any pollutant. We detected no heterogeneity of the odds ratios between strata. No clear trends by deprivation were observed in specific age groups (Figure 7).

Deprivation as a continuous variable. Significant heterogeneity (p < 0.05) was detected between the odds ratios estimated for the 136 geographic units for each of the three pollutants considered (PM10, and nitrogen dioxide). Heterogeneity remained significant after we discarded the highest and lowest 5 percent of the odds ratios in the distribution. Negative linear regression coefficients (table V) suggested that the values of the odds ratios for these pollutants tended to decrease slightly from the least deprived geographic units to the most deprived (Figure 8). The associations between the odds ratios and deprivation were not significant, however, regardless of the pollutant or the model (fixed or random effects).

SABA sales

Table VI presents the distribution of SABA sales and ambient pollutant concentrations for the entire SMA and for the five socioeconomically different strata. In 2004, SABA drugs were dispensed on 15,121 separate occasions for subjects aged 0 to 39 years. The number of sales per inhabitant slightly increased from the least to the most deprived stratum(Bard et al., 2007). The mean ambient pollutant concentrations over the study period were 20.8g/m3 for PM10 (range, 1.2 to 106.3 g/m3), 35.0 g/m3 for nitrogen dioxide (range, 6.4 to 84.3 g/m3), and 58.7g/m3 for ozone (range, 2.6 to 220.0 g/m3).

media/image3.jpeg

Figure 3.

Maps of ambient air pollutant concentrations (µg/m3) averaged by census block across the Strasbourg Metropolitan Area (2000–2005).

media/image4.jpeg

Figure 4.

Scatterplots of the rates of emergency calls for asthma according to deprivation of census blocks, for different age groups (on the left, the most privileged census blocks, on the right the most disadvantaged ones). SIR: Standardized Incidence Ratio.

media/image5.jpeg

Figure 5.

Odds ratios between asthma calls and pollutant concentrations*, for various lag times (Strasbourg Metropolitan Area, 2000-2005).

* PM10, particulate matter less than 10 micrometers in aerodynamic diameter; NO2, nitrogen dioxide; O3, ozone.

† Odds ratios reported for a 10 µg.m-3 increase in pollutant concentrations, adjusted for temperature, relative humidity, atmospheric pressure, holidays, influenza epidemics and pollen counts. Lag 0 is for pollutant concentrations averaged on the day of the call, lag 01 for pollutant concentrations averaged on the day of the call and the previous one, and so on.

media/image6.jpeg

Figure 6.

Log odds ratio for emergency asthma calls per 10 µg/m3 increase in ambient pollutant concentrations, according to stratum of socioeconomic deprivation, Strasbourg, France, 2000–2005. Stratum 1 was the least deprived, and stratum 5 was the most deprived. PM10: particulate matter less than 10 µm in aerodynamic diameter; NO2: nitrogen dioxide; O3: ozone. Bars, 95% confidence interval.

media/image7.jpeg

Figure 7.

Odds ratios between asthma calls and pollutants*, for five strata of deprivation levels and different age groups (Strasbourg Metropolitan Area, 2000-2005).

* PM10, particulate matter less than 10 micrometers in aerodynamic diameter; NO2, nitrogen dioxide; O3, ozone. Odds ratios reported for a 10 μg/m3 increase in pollutant concentrations, adjusted for temperature, relative humidity, atmospheric pressure, holidays, influenza epidemics and pollen counts

Stratum 1 is the least deprived, and stratum 5 is the most

Base model

Figure 9 presents the odds ratios for the individual lags for subjects aged 0 to 39 years. For all three pollutants, odds ratios greater than one were observed for lags from 4 to 10 days. For PM10 and nitrogen dioxide, most of the odds ratios were statistically significant (p<0.05) for lags from 5 to 10 days. For a lag of 2 days, odds ratios significantly lower than one were observed for all three pollutants. The structure of the lags reported in Figure 9 did not change when the analyses were stratified by smaller age subgroups (0 to 19 years, 20 to 39 years). No association was observed between SABA sales and the mean pollutant concentrations of the day of the sale and the 1 to 10 preceding days (lags 0 to 1 to 0 to 10).

On the basis of the associations reported in Figure 9, the lags kept for the subsequent analyses were lag 4 to 7 (mean concentrations of days 4, 5, 6, and 7 for PM10, lag 4 to 10 for nitrogen dioxide, and lag 4 to 6 for ozone). A 10 g/m3 increase in ambient concentrations of PM10, nitrogen dioxide, and ozone for these lags was associated with increases of 7.5% (95% CI, 4 to 11.2%), 8.4% (95% CI, 5 to 11.9%), and 1 (95% CI, - 0.3 to 2.2%), respectively, in SABA sales.

Testing the influence of deprivation

Figure 10 presents the odds ratios for the five different socioeconomic strata for the optimal lags mentioned above for all subjects aged 0 to 39 years. We observed no trend toward an increase or decrease of these odds ratios according to SES or when we stratified the analyses by smaller age subgroups (0 to 19 years, 20 to 39 years).

4. Discussion

4.1. Socioeconomic Status Indicator

From available census variables, we developed and validated a new French socioeconomic deprivation index at the census block level. Working at this resolution ensures a better homogeneity of the residents’ characteristics. As a result, it is possible to better identify and measure with greater precision their relation with the health events observed. Available French census data combines material and, to a lesser extent, social aspects of deprivation, with no indication of a possible association with the health events to which it may be associated. In this study, it allowed to identify a socioeconomic gradient in the utilization of mobile emergency medical services for asthma attacks in the SMA (Laurent et al., 2008). This scale has also been shown to be robust (validity criterion) and reproducible for another French metropolitan area with different socioeconomic characteristics.

All agesAge 0-19
Lag ‡Odds ratios §95% confidence intervalsp #Odds ratios §95% confidence intervalsp #
PM 10 *0-11.0350.997, 1.0750.071.0470.961, 1.1410.29
NO 2 *01.0250.990, 1.0620.161.0030.926, 1.0860.94
SO 2 *0-11.0560.979, 1.1390.151.1220.945, 1.3340.19
O 3 †0-10.9980.965, 1.0320.900.9660.891, 1.0480.40
Age 20-64Age 64 +
Lag ‡Odds ratios §95% confidence intervalsp #Odds ratios §95% confidence intervalsp #
PM 10 *0-11.0200.964, 1.0790.491.0490.985, 1.1160.14
NO 2 *01.0150.964, 1.0700.571.0500.990, 1.1130.10
SO 2 *0-10.9900.882, 1.1120.871.1020.975, 1.2450.12
O 3 †0-10.9970.949, 1.0470.901.0150.960, 1.0740.59
* PM 10 , particulate matter less than 10 micrometers in aerodynamic diameter; NO 2 , nitrogen dioxide. † O 3 , ozone. Study period extending from April 1 to September 30 of each year ‡ Lag 0 is for pollutant concentrations averaged on the day of the call, lag 01 for pollutant concentrations averaged on the day of the call and the previous one. § Odds ratios reported for a 10 µg.m -3 increase in pollutant concentrations, adjusted for temperature, relative humidity, atmospheric pressure, holidays, influenza epidemics and pollen counts # Estimated by a two-sided 2 test

Table 4.

Odds ratios between asthma calls and pollutant concentrations for different age groups and for selected lag times (Strasbourg Metropolitan Area, 2000-2005).

Dependent variables: odds ratios between pollutants and asthma calls
NO 2 †SO 2 †PM 10 †
Independent variable: deprivation indexFixed effects model-0.00272, p =0.49 ‡-0.01032, p =0.18 ‡-0.00244, p = 0.48 ‡
Random effects model-0.00284, p =0.48 ‡-0.01052, p =0.18 ‡-0.00249, p =0.47 ‡
† Odds ratios reported for a 10 µg.m -3 increase in pollutant concentrations, adjusted for temperature, relative humidity, atmospheric pressure, holidays, influenza epidemics and pollen counts. PM 10 , particulate matter less than 10 micrometers in aerodynamic diameter; NO 2 , nitrogen dioxide. ‡ Beta coefficients for regression of the values of the odds ratios on the deprivation index. Statistical significance was estimated by a two-sided χ 2 test

Table 5.

Influence of deprivation introduced as a continuous variable on the relations between pollutants and asthma calls (Strasbourg Metropolitan Area, 2000-2005).

Nonetheless, the strength of the correlations observed with the British indices by Townsend and Carstairs (correlations greater than 0.95 for both) suggest that we might have simply used them in our study, rather than constructing our own index. Indeed, their indices are more straightforward to use than ours (unweighted sum of standardized variables) and require fewer variables for their construction (4 compared with 19 in ours). However, they rely exclusively on purely material indicators, thus missing the social dimension of deprivation. The direct use of the British indices in a French setting is also limited by several factors. Firstly, these indices use data from the British census that do not correspond directly to French census items and for which an adaptation or translation is required (for example, the variable ‘‘proportion of households in which the householder belongs to social class 4 or 5’’ considered by Carstairs must be translated as ‘‘the proportion of blue-collar workers’’ for use in the French census database). Secondly, they were developed at a specific geographical level (enumeration district or ward) that is not equivalent to the French census block in terms of population size and socioeconomic characteristics. The effects of theses adaptations (i.e., variable definitions and geographical units) on the associations observed with health events have never, as far as we know, been assessed and do not encourage their use in our context. A last limitation of these indices is their inability to take into account differences between rural and urban areas because they include in their composition only variables specific to urban conditions (Gilthorpe & Wilson, 2003). The use in our study of the Index of Multiple Deprivation, constructed especially to overcome this weakness (Jordan et al., 2004) should therefore be a useful option. The French census data do not, however, lend themselves to the use of this index, for its composition includes 33 variables, many of which were either unavailable or could not be transposed from the census data. Our index thus appears to be a good alternative to the existing British indices: it is simpler to implement than the Index of Multiple Deprivation (19 compared with 33 variables, no weighting, etc.) and at the same time it overcomes the principal limitations of Carstairs and Townsend’s indices, that is, it is validated at the resolution of census blocks and able to contrast deprived urban and advantaged rural areas in our study. This study nonetheless has its own limitations, the first of which is the need to acknowledge that one size does not fit all (Braveman et al., 2005) Our index is a comprehensive and composite measure of SES and was constructed independently of the health events to which it might be associated. It includes variables related to education, income, occupation, etc. and in that respect is a better indicator than a single variable taken in isolation to describe SES as a whole.

Characterizing in more detail the social component of deprivation, both at the individual level and accounting for neighborhood influences would allow shedding more light on the complex and intricate relationships between socioeconomic characteristics and the health impact of exposure to environmental stressors. Nonetheless, in the absence of individual data, ecological studies of this type remain necessary and useful for a better understanding of the interactions between socioeconomic factors and health.

Calls to emergency medical services for asthma exacerbation

Emergency calls for asthma attacks were positively, although not significantly, associated with concentrations of PM10, and nitrogen dioxide modeled by census blocks. No association was observed for ozone. Overall, associations were higher among people younger than 20 years and older than 64 years. Socioeconomic deprivation measured by census block did not appear to influence these relations.

The daily exposure estimates we used were modeled for small areas and are geographically more precise than those usually employed to study short-term relations between air pollution and health. Compared with ambient concentrations averaged citywide, our exposure estimate likely reduced ecological biases (Jerrett et al., 2005b). Practically, the exposure measurement attributed to each subject was the concentration estimated for the census block where each patient was when the emergency network was called. However, we do not actually know whether the patients were in the same block in the hours to days preceding the call, and this fact obviously determines the extent to which our exposure measurement adequately reflects subjects’ true exposure.

This point matters mainly for subjects who are frequently away from the neighborhood they live in -principally in the 20-64 year-old age group, globally characterized by mobility and autonomy, and who work for a living (often outside the neighborhood of residence). Conversely, in general the elderly rarely commute out of their neighborhoods of residence and when they do so, they generally go shorter distances (Benlahrech et al., 2001). Children also have more limited mobility than people aged 20-64 years (Agence de l'Environnement et de la Maîtrise de l'Energie/Institut de Radioprotection et de Sûreté Nucléaire, 2003) and generally attend the schools closest to their homes.

These points support the idea that our exposure measurement is more accurate for people aged younger than 20 and older than 64 years. These subjects are more likely to have called the emergency networks from their neighborhood of residence (and thus be geocoded in it). Moreover, for these subjects, the measurement of air pollution in the neighborhood of residence is more likely to provide an adequate reflection of exposure integrated over the days preceding the call than for more mobile subjects.

The ranges of the associations we found for the base models were similar to those reported by other studies of emergency calls (Medina et al., 1997) and visits to hospital emergency departments (Galan et al., 2003) for asthma. Above all, the associations observed for PM10 were very close to those reported by two meta-analyses of the associations between this pollutant and asthma symptoms (Weinmayr et al., 2010; Kunzli et al., 2000).

Small-area deprivation (introduced either as a discrete or as continuous variable) did not appear to influence the relations between ambient air pollution and asthma attacks. Of the six previous studies that investigated the influence of socioeconomic indicators on these relations (Nauenberg & Basu, 1999; Norris et al., 1999; Lin et al., 2004; Neidell, 2004; Son et al., 2006; Kim et al., 2007), five reported higher relative risks for populations with less advantageous socioeconomic characteristics (Nauenberg & Basu, 1999; Lin et al., 2004; Neidell, 2004; Lee et al., 2006; Kim et al., 2007). However, one of these studies found evidence of interaction according to the ecological socioeconomic indicator considered, and no interaction with the individual indicator (Kim et al., 2007). The sixth study reported slightly higher relative risks for the most deprived populations (Norris et al., 1999). Nevertheless, formal comparison of the results of these studies is difficult, as they focused on socioeconomic indicators measured at very heterogeneous resolutions (Laurent et al., 2007).

Three of these studies focused on socioeconomic indicators measured at very coarse geographic resolutions. Lee et al (Lee et al., 2006) and Kim et al (Kim et al., 2007) focused on Gu neighborhoods (of around 400,000 inhabitants) in the city of Seoul (Korea), for children aged 0-15 (Lee et al., 2006) and people of all ages, respectively (Kim et al., 2007). Overall, these studies reported higher relative risks between PM10, SO2, and NO2 concentrations and asthma hospitalizations in the most deprived Gu. Norris et al, in Seattle (Washington, USA) (Norris et al., 1999) studied associations between the same pollutants and visits to emergency departments for asthma in children aged 0-18. They reported slightly lower relative risks for residents of the inner city, that is, the most deprived areas.

Two other studies focused on individual socioeconomic indicators in people of all ages. Kim et al (Lee et al., 2006), in Seoul (Korea) found that associations between PM10 and visits to emergency departments for asthma did not vary according to the annual amount of taxes paid to the national health insurance system. In contrast, Nauenberg and Basu, in Los Angeles (USA), found that effects of PM10 were greater in people with a less favorable health insurance status (Nauenberg & Basu, 1999).

Last, two studies focused on socioeconomic indicators measured by small areas. Neidell (Neidell, 2004) observed that carbon monoxide and ozone had a greater effect on asthma hospitalizations of Californian children aged 3-18 in ZIP codes characterized by lower educational attainment. This study, however, estimated associations between asthma hospitalizations and air pollutants on the basis of monthly indicators, which are inadequate to study short-term relations between these factors. Lin et al studied children aged 6-12 years in Vancouver (Canada) (Lin et al., 2004) and reported higher associations between nitrogen dioxide, sulfur dioxide and hospitalization for asthma in enumeration areas with lower household income levels.

The study by Lin et al (Lin et al., 2004) is the most comparable in design to ours, but reports somewhat different results. The reasons for this are unclear. Pollutant concentrations were very similar in the two settings. Although the studies used different types of exposure measurements (pollutant concentrations averaged citywide for Lin et al (Lin et al., 2004) and modeled by census blocks in ours), this difference does not explain the variation in findings. Indeed, alternative analysis in the SMA with citywide average exposure measurements did not noticeably change our results about interactions with deprivation. Lack of statistical power appears to be a plausible explanation for the difference between our results and those of Lin et al. For comparable age groups (0-20 versus 6-12) we had one quarter the number of health events to analyze.

Another point is that findings of interaction with deprivation are not necessarily transposable from one setting to another. If small-area deprivation does exert an influence on the relations between air pollution and asthma attacks, it would most likely be mediated through “third” factors that in some (but perhaps not in all) settings would be distributed unequally according to deprivation. Previous studies report that several factors thought to strengthen the associations between air pollution and asthma attacks are more present in deprived than in well-off neighborhoods. Among these are the prevalence of (both active and passive) smoking (Diez-Roux et al., 2003), psychosocial stress (Gold & Wright, 2005), unhealthy eating habits (Diez-Roux et al., 1999), amounts of indoor allergens (Kitch et al., 2000), inadequate compliance with anti-inflammatory medication (Gottlieb et al., 1995) and (a plausible result of the factors mentioned above) a higher ratio of severe to moderate forms of asthma among subjects with asthma (Basagana et al., 2004). Nevertheless, the distribution of these factors according to small-area deprivation may differ between study settings, due, for instance, to differences in climate (affecting allergen proliferation), social and cultural characteristics of local populations (influencing, among other things, eating and smoking habits) or effectiveness of health systems (influencing prescription of and compliance with anti-inflammatory medication).

Moreover, although we observed no interactions with small-area deprivation in the SMA, this does not rule out the existence here of interactions by socioeconomic factors measured at other resolutions (individual, household, geographic areas more or less fine than the French census block), as the study by Kim et al clearly illustrates (Kim et al., 2007). The use of multilevel models, which make it possible to assess more precisely the influence of factors (e.g. socioeconomic characteristics) measured at different resolutions, would be useful in studying this question further (O'Neill et al., 2003).

SABA sales

We observed positive associations between ambient concentrations of atmospheric pollutants and SABA sales for subjects < 40 years old. These are expressed with latency periods of 4 to 10 days and do not tend to increase or decrease according to SES.

This study is the first to examine the relations between exposure to urban air pollution and SABA sales. The use of this indicator, obtained from the four primary French health insurance funds, allowed us to cover > 90% of the local population and to capture the entire range of SES in the SMA. People not covered by these funds are mainly employees of various sectors once publicly owned (railway, electricity, gas), with jobs ranging from manual workers to administrators and mainly in the middle classes. The large number of SABA sales allowed us to measure their associations with air pollution modeled by small areas. This resolution is particularly pertinent for studying this risk factor because its spatial distribution varies strongly within urban areas (Laurent et al., 2008). These large numbers also allowed us to test the existence of interactions by neighborhood SES with satisfactory statistical power. The event analyzed is the patient’s purchase of one (or sometimes more) box of drugs, and not a quantity of active ingredient delivered or really inhaled. Naureckas et al. (Naureckas et al., 2005) nonetheless showed that this indicator is a good predictor of the risk of emergency department visits and of hospitalization for asthma attacks in the days immediately afterwards. These purchases generally reflect asthma morbidity less severe than that requiring hospitalization or emergency treatment (Naureckas et al., 2005). A large portion of SABA sales anticipate the respiratory disorders the drugs are intended to treat. These sales, which need not be associated in time with asthma, add some “noise” to the data. Noise is standard in ecologic studies of the short-term health effects of air pollution because of the influence of unmeasured competing factors (other than air pollution) on the temporal distribution of the health outcomes studied. Nevertheless, if the signal-to-noise ratio is sufficiently high, the effect specifically due to air pollution can be observed. In our study, despite the additional noise due to anticipatory sales, this ratio appears to be high enough to detect statistically significant associations. These associations are consistent with those reported by most panel studies (von Klot et al., 2002; Romeo et al., 2006; Schildcrout et al., 2006) that have investigated the relation between air pollution and SABA consumption.

The associations observed involved latency periods of 4 to 10 days, an order of magnitude similar to the delayed responses observed in two earlier temporal ecologic studies of the relation between air pollution and drug sales mucolytic and antitussive agents (Zeghnoun et al., 1999) “cough and cold preparations,” and all types of anti-COPD/antiasthma drugs (Pitard et al., 2004). This is probably because the latency periods between exposure to air pollution and drug purchases result from a mixed process involving both pathophysiological response and management of medicine supplies. The literature (von Klot et al., 2002; Rabinovitch et al., 2006; Romeo et al., 2006; Schildcrout et al., 2006) shows that in people with asthma, air pollution induces respiratory disorders expressed by increased SABA consumption, with low latency periods (several hours 34 to several days (von Klot et al., 2002; Schildcrout et al., 2006)). This increased consumption requires a successive -but not necessarily immediate- replenishment of the SABA supply. Several days of delay, possibly marked by a return to a normal rhythm of use, may pass before the purchase, which may explain the particularly long lags (up to 10 days) observed. For lag 2, we observed odds ratios significantly less than one. Several authors, in other settings and with different methods, report similar findings. Zeghnoun et al. (Zeghnoun et al., 1999) also observed low associations for lag 2, especially for mucolytic and antitussive drugs. Moreover, in a panel study of SABA consumption, Rabinovitch et al. (Rabinovitch et al., 2006) also observed relative risks less than one for the same lag. Von Klot et al. (von Klot et al., 2002) report comparable observations specifically for lag 1. No satisfactory explanation has yet been found for these observations.

SES did not influence the relation between SABA sales and ambient air pollution. This is consistent with the results appearing above for an indicator of more severe asthma morbidity, interventions of mobile emergency medical services for asthma attacks (Laurent et al., 2008).

In conclusion, emergency calls for asthma attacks SABA sales for children, adolescents and young adults were positively (not significantly for the former, but significantly for SABA sales) associated with PM10, NO2 but not O3 concentrations modeled by small areas. Small area deprivation did not influence these associations. Nonetheless, discrepancies between our results on emergency calls and those of the study of Lin et al (Lin et al., 2004) emphasize the need to investigate this question further in other study settings. Similarly, the observations we made regarding SABA sales do not rule out the possibility that SES might be an interaction factor in other settings, for the distribution according to SES of other factors that might modulate the relations between air pollution and asthma morbidity may well differ between countries or even cities. The results on SABA are consistent with those of panels of asthma patients and their SABA consumption, although expressed here with longer time lags. Our results support the usefulness of SABA sales for the analysis of relations between asthma morbidity and air pollution.

media/image8.jpeg

Figure 8.

Odds ratios between asthma calls and pollutants in the 136 statistical units retained, ranked according to deprivation (from the least deprived on the bottom to the most deprived at the top). Strasbourg Metropolitan Area, 2000-2005.

†Odds ratios reported for a 1 µg.m-3 increase in pollutant concentrations (for the sake of figure visibility), adjusted for temperature, relative humidity, atmospheric pressure, holidays, influenza epidemics and pollen counts. PM10, particulate matter less than 10 micrometers in aerodynamic diameter; NO2: nitrogen dioxide

media/image9.jpeg

Figure 9.

Odds ratios for an increase of 10 µg/m3 in the pollutants concentration, subjects <40 years old; for various lags.

media/image10.jpg

Figure 10.

Interactions by socioeconomic level in those aged < 40 years: associations for the five different strata of socioeconomic levels.

Stratum 1 is the most advantaged and stratum 5 the most deprived. Associations were estimated for “optimal” lags defined according to the associations reported for individual lags, see Figure 7. For example, for PM10 among those younger than 40 years, 4-7 corresponds to the mean concentrations for day D4, D5, D6 and D7.

Population, age groupsPollutants (mean, SD )
Deprivation stratum §Total, 0-390-910-1920-39PM 10 †NO 2†O 3 ‡
Stratum 143,6749,34211,48922,84319.4 9.929.9 10.363.1 36.5
Stratum 247,7579,35911,02227,37620.5 10.434.2 10.259.3 36.1
Stratum 354,5279,4989,55035,47921.6 10.637.5 10.756.0 35.7
Stratum 465,99410,40411,91743,67321.4 10.537.0 10.657.0 35.8
Stratum 549,11112,44013,26923,40220.7 10.335.7 10.258.6 35.8
Overall261,06351,04352,247152,77320.8 10.235.0 10.358.7 36.0
SABA sales, age groups
Deprivation stratum §Total, 0-390-910-1920-39
Stratum 12,140733426981
Stratum 22,5387285781,232
Stratum 32,9729314601,581
Stratum 43,7521,1926311,929
Stratum 53,7191,3327631,624
Overall15,1214,9162,8587,347
§ Stratum 1 is the last deprived, and stratum 5 the most deprived † Concentrations averaged for year 2004, in microgram per cubic meter ‡ Maximum 8-hour daily concentrations, averaged for summer months ( April 1 to September 30) of year 2004, in microgram per cubic meter

Table 6.

Distribution of SABA sales and air pollutant concentrations, Strasbourg Metropolitan Area, 2004.

Acknowledgements

The authors thank all the organizations that kindly provided the data used in these analyses: URCAM Alsace (especially Benoît Wollbrett), RSI Alsace (especially Katia Bischoff), MSA Alsace (especially Hervé Hunold), MGEL (especially Gérard Rey), Météo France, National Network of Aerobiologic Surveillance, and the INSERM Sentinelles Network. The authors also thank Professor Frédéric de Blay for providing opinion as a specialist in pulmonary medicine, and Dr. Fabienne Wachet for providing expertise on drug supply and repayment systems. Finally, the authors thank Jo Ann Cahn for editorial assistance.

These studies were made possible by the grant ANR SEST 0057 05 from the French National Research Agency (ANR).

References

1 - Agence de l’Environnement et de la Maîtrise de l’Energie/Institut de Radioprotection et de Sûreté Nucléaire 2003 Ciblex: database of descriptive characteristics of the French population residing near polluted sites (CD-ROM) ref N 4773. Angers, France/Fontenay-aux-Roses, France.
2 - I. Annesi-Maesano, N. Agabiti, R. Pistelli, M. F. Couilliot, F. Forastiere, 2003 Subpopulations at increased risk of adverse health outcomes from air pollution. Eur Respir J Suppl 40: 57s 63s .
3 - R. W. Atkinson, H. R. Anderson, J. Sunyer, J. Ayres, M. Baccini, J. M. Vonk, A. Boumghar, F. Forastiere, B. Forsberg, G. Touloumi, J. Schwartz, K. Katsouyanni, 2001 Acute effects of particulate air pollution on respiratory admissions: results from APHEA 2 project. Air Pollution and Health: a European Approach. Am J Respir Crit Care Med 164(10 Pt 1): 1860 1866.
4 - D. Bard, O. Laurent, L. Filleul, S. Havard, S. Deguen, C. Segala, G. Pedrono, E. Rivière, C. Schillinger, L. Rouïl, D. Arveiler, D. Eilstein, 2007 Exploring the joint effect of atmospheric pollution and socioeconomic status on selected health outcomes: an overview of the PAISARC project. Environ Res Lett 2(045003): 7 pp.
5 - X. Basagana, J. Sunyer, M. Kogevinas, J. P. Zock, E. Duran-Tauleria, D. Jarvis, P. Burney, J. M. Anto, 2004 Socioeconomic status and asthma prevalence in young adults: the European Community Respiratory Health Survey. Am. J. Epidemiol. 160 2 178 88 .
6 - M. Beguin, D. Pumain, 1994 La représentation des données géographiques. Paris, Armand Colin.
7 - J. Benach, Y. Yasui, 1999 Geographical patterns of excess mortality in Spain explained by two indices of deprivation. J. Epidemiol. Community Health 53 7 423 431 .
8 - N. Benlahrech, A. Le Ruyet, C. Livebardon, M. Dejeammes, 2001 The mobility of old people (analysis of the household displacement surveys) [In French]. CERTU. Lyon, France.
9 - P. D. Blanc, I. H. Yen, H. Chen, P. P. Katz, G. Earnest, J. R. Balmes, L. Trupin, N. Friedling, E. H. Yelin, M. D. Eisner, 2006 Area-level socio-economic status and health status among adults with asthma and rhinitis. Eur. Respir. J. 27 1 85 94 .
10 - P. A. Braveman, C. Cubbin, S. Egerter, S. Chideya, K. S. Marchi, M. Metzler, S. Posner, 2005 Socioeconomic status in health research: one size does not fit all. JAMA 294 22 2879 2888 .
11 - M. Buzzelli, M. Jerrett, R. Burnett, N. Finkelstein, 2003 Spatiotemporal perspectives on air pollution and environmental justice in Hamilton, Canada, 1985-1996. Annals of the Association of American Geographers 93 3 557 573 .
12 - V. Carstairs, 1995 Deprivation indices: their interpretation and use in relation to health. J Epidemiol Community Health 49 Suppl 2: S3 8 .
13 - V. Carstairs, R. Morris, 1991 Deprivation and Health in Scotland. Aberdeen, Aberdeen University Press.
14 - G. Cesaroni, S. Farchi, M. Davoli, F. Forastiere, C. A. Perucci, 2003 Individual and area-based indicators of socioeconomic status and childhood asthma. Eur. Respir. J. 22 4 619 624 .
15 - B. Challier, J. F. Viel, 2001 [Relevance and validity of a new French composite index to measure poverty on a geographical level]. Rev Epidemiol Sante Publique 49 1 41 50 .
16 - P. Clifford, S. Richardson, D. Hemon, 1989 Assessing the significance of the correlation between two spatial processes. Biometrics 45 1 123 134 .
17 - L. Com-Ruelle, B. Crestin, S. Dumesnil, 2000 Asthma in France by gravity degrees. IRDES. Paris, CREDES,.
18 - J. Coste, J. Fermanian, A. Venot, 1995 Methodological and statistical problems in the construction of composite measurement scales: a survey of six medical and epidemiological journals. Stat Med 14 4 331 345 .
19 - DETR 2000 Measuring deprivation at the small area level: the indices of deprivation 2000. 31, N.
20 - A. V. Diez-Roux, F. J. Nieto, L. Caulfield, H. A. Tyroler, R. L. Watson, M. Szklo, 1999 Neighbourhood differences in diet: the Atherosclerosis Risk in Communities (ARIC) Study. J Epidemiol Community Health 53 1 55 63 .
21 - A. V. Diez-Roux, S. S. Merkin, P. Hannan, D. R. Jacobs, C. I. Kiefe, 2003 Area characteristics, individual-level socioeconomic indicators, and smoking in young adults: the coronary artery disease risk development in young adults study. Am J Epidemiol 157 4 315 326 .
22 - C. Eibner, R. Sturm, 2006 US-based indices of area-level deprivation: results from HealthCare for Communities. Soc Sci Med 62 2 348 359 .
23 - L. Ellison-Loschmann, J. Sunyer, E. Plana, N. Pearce, J. P. Zock, D. Jarvis, C. Janson, J. M. Anto, M. Kogevinas, 2007 Socioeconomic status, asthma and chronic bronchitis in a large community-based study. Eur Respir J 29 5 897 905 .
24 - B. Erbas, A. M. Kelly, B. Physick, C. Code, M. Edwards, 2005 Air pollution and childhood asthma emergency hospital admissions: estimating intra-city regional variations. Int J Environ Health Res 15 1 11 20 .
25 - J. Fermanian, 1996 [Evaluating correctly the validity of a rating scale: the numerous pitfalls to avoid]. Rev Epidemiol Sante Publique 44 3 278 286 .
26 - M. M. Finkelstein, M. Jerrett, P. De Luca, N. Finkelstein, D. K. Verma, K. Chapman, M. R. Sears, 2003 Relation between income, air pollution and mortality: a cohort study. CMAJ. 169 5 397 402 .
27 - A. Flahault, T. Blanchon, Y. Dorleans, L. Toubiana, J. F. Vibert, A. J. Valleron, 2006 Virtual surveillance of communicable diseases: a 20-year experience in France. Stat Methods Med Res 15 5 413 421 .
28 - K. Folwell, 1995 Single measures of deprivation. J Epidemiol Community Health 49 Suppl 2: S51 56 .
29 - Y. Fukuda, K. Nakamura, T. Takano, 2007 Higher mortality in areas of lower socioeconomic position measured by a single index of deprivation in Japan. Public Health 121 3 163 173 .
30 - I. Galan, A. Tobias, J. R. Banegas, E. Aranguez, 2003 Short-term effects of air pollution on daily asthma emergency room admissions. Eur Respir J 22 5 802 808 .
31 - M. Gilthorpe, R. Wilson, 2003 Rural/urban differences in the association between deprivation and healthcare utilisation. Soc Sci Med 57(2055 2063).
32 - Global Initiative for Asthma 2007 Global strategy for asthma management and prevention.
33 - D. R. Gold, R. Wright, 2005 Population disparities in asthma. Annu. Rev. Public Health 26 89 113 .
34 - D. J. Gottlieb, A. S. Beiser, G. T. O’Connor, 1995 Poverty, race, and medication use are correlates of asthma hospitalization rates. A small area analysis in Boston. Chest 108 1 28 35 .
35 - S. Hajat, A. Haines, S. A. Goubet, R. W. Atkinson, H. R. Anderson, 1999 Association of air pollution with daily GP consultations for asthma and other lower respiratory conditions in London. Thorax 54 7 597 605 .
36 - S. Havard, S. Deguen, J. Bodin, O. Laurent, D. Bard, 2008 A small-area index of socioeconomic deprivation to capture health inequalities in France. Soc Sci Med 2008, 67-2007 2016.
37 - K. Haynes, S. Lall, P. Trice, 2001 Spatial issues in environmental equity. Int. J. Environ. Technol. Manag. 1(1-2): 17 31.
38 - B. Jalaludin, B. Khalaj, V. Sheppeard, G. Morgan, 2008 Air pollution and ED visits for asthma in Australian children: a case-crossover analysis. Int Arch Occup Environ Health 81 8 967 974 .
39 - H. Janes, L. Sheppard, T. Lumley, 2005 Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias. Epidemiology 16 6 717 726 .
40 - B. Jarman, 1983 Identification of underprivileged areas. Br Med J (Clin Res Ed) 286 6379 1705 1709 .
41 - M. Jerrett, A. Arain, P. Kanaroglou, B. Beckerman, D. Potoglou, T. Sahsuvaroglu, J. Morrison, C. Giovis, 2005a A review and evaluation of intraurban air pollution exposure models. J. Expo. Anal. Environ. Epidemiol. 15 2 185 204 .
42 - M. Jerrett, R. Burnett, P. Kanaroglou, J. Eyles, N. Finkelstein, C. Giovis, J. Brook, 2001 A GIS- environmental justice analysis of particulate air pollution in Hamilton, Canada. Environ. Planning A 33 955 973 .
43 - M. Jerrett, R. T. Burnett, R. Ma, C. A. Pope , III.D. Krewski, K. B. Newbold, G. Thurston, Y. Shi, N. Finkelstein, E. E. Calle, M. J. Thun, 2005b Spatial analysis of air pollution and mortality in Los Angeles. Epidemiology 16 6 727 736 .
44 - H. Jordan, P. Roderick, D. Martin, 2004 The Index of Multiple Deprivation 2000 and accessibility effects on health. J Epidemiol Community Health 58 3 250 257 .
45 - S. Y. Kim, M. S. O’Neill, J. T. Lee, Y. Cho, J. Kim, H. Kim, 2007 Air pollution, socioeconomic position, and emergency hospital visits for asthma in Seoul, Korea. Int Arch Occup Environ Health 80 8 701 710 .
46 - B. T. Kitch, G. Chew, H. A. Burge, M. L. Muilenberg, S. T. Weiss, T. A. Platts-Mills, G. O’Connor, D. R. Gold, 2000 Socioeconomic predictors of high allergen levels in homes in the greater Boston area. Environ Health Perspect 108 4 301 307 .
47 - N. Krieger, D. R. Williams, N. E. Moss, 1997 Measuring social class in US public health research: concepts, methodologies, and guidelines. Annu. Rev. Public Health 18 341 378 .
48 - A. E. Kunst, F. Groenhof, J. P. Mackenbach, E. W. Health, 1998 Occupational class and cause specific mortality in middle aged men in 11 European countries: comparison of population based studies. EU Working Group on Socioeconomic Inequalities in Health. BMJ 316 7145 1636 1642 .
49 - N. Kunzli, R. Kaiser, S. Medina, M. Studnicka, O. Chanel, P. Filliger, M. Herry, F. Horak , Jr.V. Puybonnieux-Texier, P. Quenel, J. Schneider, R. Seethaler, J. C. Vergnaud, H. Sommer, 2000 Public-health impact of outdoor and traffic-related air pollution: a European assessment. Lancet 356 9232 795 801 .
50 - O. Laurent, D. Bard, L. Filleul, C. Segala, 2007 Effect of Socioeconomic Status on the Relation between Atmospheric Pollution and Mortality: A Review. J Epidemiol Community Health.
51 - O. Laurent, G. Pedrono, C. Segala, L. Filleul, S. Havard, S. Deguen, C. Schillinger, E. Riviere, D. Bard, 2008 Air pollution, asthma attacks, and socioeconomic deprivation: a small-area case-crossover study. Am J Epidemiol 168 1 58 65 .
52 - J. T. Lee, J. Y. Son, H. Kim, S. Y. Kim, 2006 Effect of air pollution on asthma-related hospital admissions for children by socioeconomic status associated with area of residence. Arch Environ Occup Health 61 3 123 130 .
53 - M. Lin, Y. Chen, P. J. Villeneuve, R. T. Burnett, L. Lemyre, C. Hertzman, K. M. McGrail, D. Krewski, 2004 Gaseous air pollutants and asthma hospitalization of children with low household income in Vancouver, British Columbia, Canada. Am. J. Epidemiol. 159 3 294 303 .
54 - V. Lorant, 2000 [Mortality socio-economic inequalities for small-areas in Belgium: assessing concentration]. Rev Epidemiol Sante Publique 48 3 239 247 .
55 - R. J. Marshall, 1991 Mapping disease and mortality rates using empirical Bayes estimators. J R Stat Soc Ser C Appl Stat 40 2 283 294 .
56 - R. Maynard, D. Krewski, R. T. Burnett, J. Samet, J. R. Brook, G. Granville, L. Craig, 2003 Health and air quality: directions for policy-relevant research. J Toxicol Environ Health A 66(16-19): 1891 1904 .
57 - C. McHugh, D. Carruthers, H. Edmunds, 1997 ADMS-Urban: an air quality management system for traffic, domestic and industrial pollution. Int. J. Environ Pollution 8 437 440 .
58 - M. Mclure, M. Mittelman, 1991 Should we use a case-crossover design? Ann Rev Public Health 21 193 221 .
59 - S. Medina, A. Le Tertre, P. Quenel, Y. Le Moullec, P. Lameloise, J. C. Guzzo, B. Festy, R. Ferry, W. Dab, 1997 Air pollution and doctors’ house calls: results from the ERPURS system for monitoring the effects of air pollution on public health in Greater Paris, France, 1991-1995. Evaluation des Risques de la Pollution Urbaine pour la Sante. Environ Res 75 1 73 84 .
60 - L. C. Messer, B. A. Laraia, J. S. Kaufman, J. Eyster, C. Holzman, J. Culhane, I. Elo, J. G. Burke, P. O’Campo, 2006 The development of a standardized neighborhood deprivation index. J Urban Health 83 6 1041 1062 .
61 - A. Mielck, P. Reitmeir, M. Wjst, 1996 Severity of childhood asthma by socioeconomic status. Int. J. Epidemiol. 25 2 388 393 .
62 - R. Morris, V. Carstairs, 1991 Which deprivation? A comparison of selected deprivation indexes. J Public Health Med 13 4 318 326 .
63 - E. Nauenberg, K. Basu, 1999 Effect of insurance coverage on the relationship between asthma hospitalizations and exposure to air pollution. Public Health Rep 114 2 135 148 .
64 - E. T. Naureckas, V. Dukic, X. Bao, P. Rathouz, 2005 Short-acting beta-agonist prescription fills as a marker for asthma morbidity. Chest 128 2 602 608 .
65 - M. J. Neidell, 2004 Air pollution, health, and socio-economic status: the effect of outdoor air quality on childhood asthma. J. Health Econ. 23 6 1209 1236 .
66 - A. Niggebrugge, R. Haynes, A. Jones, A. Lovett, I. Harvey, 2005 The index of multiple deprivation 2000 access domain: a useful indicator for public health? Soc Sci Med 60 12 2743 2753 .
67 - G. Norris, S. N. YoungPong, J. Q. Koenig, T. V. Larson, L. Sheppard, J. W. Stout, 1999 An association between fine particles and asthma emergency department visits for children in Seattle. Environ. Health Perspect. 107 6 489 493 .
68 - M. S. O’Neill, M. Jerrett, I. Kawachi, J. I. Levy, A. J. Cohen, N. Gouveia, P. Wilkinson, T. Fletcher, L. Cifuentes, J. Schwartz, 2003 Health, wealth, and air pollution: advancing theory and methods. Environ. Health Perspect. 111 16 1861 1870 .
69 - R. Pampalon, G. Raymond, 2000 A deprivation index for health and welfare planning in Quebec. Chronic Dis Can 21 3 104 113 .
70 - A. Pitard, A. Zeghnoun, A. Courseaux, J. Lamberty, V. Delmas, J. L. Fossard, H. Villet, 2004 Short-term associations between air pollution and respiratory drug sales. Environ Res 95 1 43 52 .
71 - E. Prescott, N. Godtfredsen, J. Vestbo, M. Osler, 2003 Social position and mortality from respiratory diseases in males and females. Eur Respir J 21 5 821 826 .
72 - N. Rabinovitch, M. Strand, E. W. Gelfand, 2006 Particulate levels are associated with early asthma worsening in children with persistent disease. Am J Respir Crit Care Med 173 10 1098 1105 .
73 - E. Romeo, M. De Sario, F. Forastiere, P. Compagnucci, M. Stafoggia, A. Bergamaschi, C. A. Perucci, 2006 [PM 10 exposure and asthma exacerbations in pediatric age: a meta-analysis of panel and time-series studies]. Epidemiol Prev 30(4-5): 245-54.
74 - C. Salmond, P. Crampton, F. Sutton, 1998 NZDep91: A New Zealand index of deprivation. Aust N Z J Public Health 22 7 835 837 .
75 - J. Samet, D. Krewski, 2007 Health effects associated with exposure to ambient air pollution. J Toxicol Environ Health A 70(3-4): 227 242 .
76 - J. S. Schildcrout, L. Sheppard, T. Lumley, J. C. Slaughter, J. Q. Koenig, G. G. Shapiro, 2006 Ambient air pollution and asthma exacerbations in children: an eight-city analysis. Am J Epidemiol 164 6 505 517 .
77 - G. K. Singh, 2003 Area deprivation and widening inequalities in US mortality, 1969-1998. Am J Public Health 93 7 1137 1143 .
78 - J. Y. Son, H. Kim, J. T. Lee, S. Y. Kim, 2006 [Relationship between the exposure to ozone in Seoul and the childhood asthma-related hospital admissions according to the socioeconomic status]. J. Prev. Med. Pub. Health 39 1 81 86 .
79 - J. Sunyer, C. Spix, P. Quenel, A. Ponce-de-Leon, A. Ponka, T. Barumandzadeh, G. Touloumi, L. Bacharova, B. Wojtyniak, J. Vonk, L. Bisanti, J. Schwartz, K. Katsouyanni, 1997 Urban air pollution and emergency admissions for asthma in four European cities: the APHEA Project. Thorax 52 9 760 765 .
80 - J. E. Tello, J. Jones, P. Bonizzato, M. Mazzi, F. Amaddeo, M. Tansella, 2005 A census-based socio-economic status (SES) index as a tool to examine the relationship between mental health services use and deprivation. Soc Sci Med 61 10 2096 2105 .
81 - M. Thibaudon, 2004 [Report of the activity of RNSA (National Aero-biological Surveillance Network), 2003]. Eur Ann Allergy Clin Immunol 36 7 252 255 .
82 - P. E. Tolbert, J. A. Mulholland, D. L. MacIntosh, F. Xu, D. Daniels, O. J. Devine, B. P. Carlin, M. Klein, J. Dorley, A. J. Butler, D. F. Nordenberg, H. Frumkin, P. B. Ryan, M. C. White, 2000 Air quality and pediatric emergency room visits for asthma in Atlanta, Georgia, USA. Am J Epidemiol 151 8 798 810 .
83 - C. Tollier, I. Fusier, M. C. Husson, 2005 [ATC and EphMRA classifications: evolution from 1996 to 2003 and comparative analysis]. Therapie 60 1 47 56 .
84 - P. Townsend, 1987 Deprivation. J Soc Pol 16 125 146 .
85 - P. Townsend, P. Phillimore, A. Beattie, 1988 Health and Deprivation: Inequality and the North. London, Croom Helm.
86 - F. L. Urbano, R. M. Pascual, 2005 Contemporary issues in the care of patients with chronic obstructive pulmonary disease. J Manag Care Pharm 11(5 Suppl A): S2 13 ; quiz S4-6.
87 - G. Viegi, F. Pistelli, D. L. Sherrill, S. Maio, S. Baldacci, L. Carrozzi, 2007 Definition, epidemiology and natural history of COPD. Eur Respir J 30 5 993 1013 .
88 - S. von Klot, G. Wolke, T. Tuch, J. Heinrich, D. W. Dockery, J. Schwartz, W. G. Kreyling, H. E. Wichmann, A. Peters, 2002 Increased asthma medication use in association with ambient fine and ultrafine particles. Eur Respir J 20 3 691 702 .
89 - D. J. Ward, J. G. Ayres, 2004 Particulate air pollution and panel studies in children: a systematic review. Occup Environ Med 61(4): e13.
90 - G. Weinmayr, E. Romeo, M. De Sario, S. K. Weiland, F. Forastiere, Short-term effects of PM10 and NO2 on respiratory health among children with asthma or asthma-like symptoms: a systematic review and meta-analysis. Environ Health Perspect 118 4 449 457 .
91 - R. J. Wright, S. V. Subramanian, 2007 Advancing a multilevel framework for epidemiologic research on asthma disparities. Chest 132(5 Suppl): 757S EOF 769S.
92 - A. Zeghnoun, P. Beaudeau, F. Carrat, V. Delmas, O. Boudhabhay, F. Gayon, D. Guincetre, P. Czernichow, 1999 Air pollution and respiratory drug sales in the City of Le Havre, France, 1993-1996. Environ Res 81 3 224 230 .