|
|
||||||||||||||
|
|
|||||||||||||||
RESEARCH REPORT |
1 Community Medicine and Public Health, Department of Clinical Sciences, Malmö University Hospital, Lund University, Malmö, Sweden
2 Inserm UMR-S 707 (Inserm-Université Pierre et Marie Curie-Paris6), Paris, France
3 Department of Health Sciences, Malmö University Hospital, Lund University, Malmö, Sweden
Correspondence to:
Correspondence to:
Dr B Chaix
UMR-S 707 Inserm-Université Pierre et Marie Curie-Parise Faculté de Médecine Saint-Antoine, 27 rue Chaligny, 75012 Paris6, France;chaix{at}u707.jussieu.fr
Accepted for publication 29 May 2006
| ABSTRACT |
|---|
Objective: To investigate whether area effects vary according to the individual and contextual characteristics of the population, socioeconomic contextual influences were assessed in different age groups and within territories of differing population densities.
Design: Multilevel survival analysis of a 28-year longitudinal database.
Participants: 341 048 residents of the Scania region in Sweden, reaching age 5079 years in 1996, followed up over 7 years.
Results: After adjustment for several individual socioeconomic indicators over the adult age, Cox multilevel models indicated geographical variations in IHD mortality and socioeconomic contextual effects on the mortality risk. However, the magnitude of geographical variations and strength of contextual effects were modified by the age of individuals and the population density of their residential area: socioeconomic contextual effects were much stronger among non-elderly than among elderly adults, and much larger within urban territories than within rural ones. As a consequence, among non-elderly residents of urban territories, the socioeconomic contextual effect was almost as large as the effect of individual 20-year cumulated income.
Conclusions: Non-elderly residents of deprived urban neighbourhoods constitute a major target for both contextual epidemiology of coronary disease and public health interventions aimed at reducing the detrimental effects of the social environment on IHD.
Abbreviations: ICD, International Classification of Diseases; IDR, interdecile range; IHD, ischaemic heart disease; IHR, interquartile hazard ratio; IQR, interquartile range; MHR, median hazard ratio
Worldwide, ischaemic heart disease (IHD) is the leading cause of death, and will still be so by far in 2020.1,2 In previous literature, an inverse association between socioeconomic position and IHD has been consistently shown.3 In an attempt to explain it, a growing body of literature suggests that effects of the social characteristics of the residential context may independently contribute to the socioeconomic gradient in IHD.4
To orientate future aetiological research and in a public health perspective, it is of importance to assess the magnitude of geographical variability in IHD mortality,57 and examine whether socioeconomic contextual effects are significant as compared with the effect of individual socioeconomic position. Some studies have observed geographical variations in IHD or other cardiovascular diseases between cities8,9,10 or neighbourhoods.1113 Apart from ecological studies,14,15 some authors have found effects of the neighbourhood socioeconomic environment9,12,1625 or other contextual dimensions2630 on IHD or other circulatory diseases, after adjustment for individual socioeconomic factors.
A review of the literature indicates that these effects were of very moderate to quite important magnitude. It is therefore important to identify circumstances that may create differences in the magnitude of geographical variations and contextual effects on IHD. Part of the variations in the strength of contextual effects between studies may be due to dissimilarities in the definition of residential areas,7 in the individual variables used for adjustment of the models, and in the cardiovascular outcomes consideredthat is, incidence,12,17,21,24,27 prevalence8,9,16,18,31,32 or mortality.11,19,20,22,23,26,30,33 Moreover, differences between studies in the individual and contextual characteristics of the populations may also create fluctuations in the magnitude of geographical variations and contextual effects.
In this study, to contribute to the assessment of the strength of contextual effects on IHD, we first examined the extent to which contextual influences varied by age of the population considered. The steep increase in IHD risk over the life course may create differences in the magnitude of contextual effects between adults and elderly people. On the basis of previous literature,11,26,34 we assumed that detrimental social conditions have a stronger role in premature mortality than in older-age mortality, and that the healthy survivor effect may flatten the socioeconomic gradient at older ages. We, therefore, expected a weaker socioeconomic contextual effect among elderly than among non-elderly adults.
On the other hand, our investigation of socioeconomic contextual effects on IHD mortality was conducted within the whole region of Scania, Swedenthat is, within territories of varying population densities. Most previous contextual studies were conducted within urban territories4; it is therefore of interest to compare the magnitude of socioeconomic contextual effects between territories of contrasted population densities. We hypothesised a stronger effect of the socioeconomic environment in densely than in sparsely populated territories, as the causal agents of the socioeconomic contextual effect (whatever they might be) may be more concentrated around individuals in high-density areas.
To deal with these issues, we considered a large cohort comprising all individuals living in the Scania region reaching age 5079 years in 1996. All those individuals were followed up over 7 years, to assess occurrence of IHD mortality. There were two successive steps in our strategy to assess the importance the residential context may have for IHD mortality. Firstly, we used multilevel survival models to quantify the magnitude of variations between areas in mortality risk, using different approaches to express geographical variability on the commonly used hazard ratio (HR) scale.35,36 Secondly, we estimated the effect of the socioeconomic environment on IHD mortality adjusted for individual factors, and compared its magnitude with the effect of individual income. Overall, our study expands on previous literature by using complementary approaches to quantify the importance the residential context may have for IHD mortality, and assessing whether contextual influences differ in magnitude by age of individuals and population density of the area of residence.
| METHODS |
|---|
We constituted a cohort with 1 January 1996 as baseline, comprising all Scanian inhabitants reaching age 5079 years in 1996. The database was split into two parts to investigate separately individuals aged 5064 years and 6579 years at baseline. The cohort comprised 192 840 individuals aged 5064, and 148 208 individuals aged 6579 years at baseline.
Measures
We defined IHD deaths using the International Classification of Diseases (ICD) as ICD-9 codes 410414 or ICD-10 codes I20I25 for the underlying or contributing causes of death.
Besides age and gender, we considered marital status, as an indicator of familial support and help, and information on individual socioeconomic position over the adult lifethat is, education, occupation 25 years before baseline and individual cumulated income over the past 20 years. Age was coded in 5-year categories. Marital status was dichotomised into married/cohabiting and living alone (ie single, divorced or widowed). Education (from the 1970 census) was divided into three classes (
7, 89 and >9 years). Occupation 25 years before baseline was defined from 1970 census data. We distinguished between non-manual, manual, self-employed, farmers and others.
Rather than household income, only individual income was available. To adjust for the cumulated effects of wealth/poverty over time,9 we considered annual income in 1975, 1980, 1985, 1990 and 1995, converted each income variable into a rank between 1 and 100 to obtain values comparable across years, and computed the average of the five variables for each individual. The cumulated income variable was divided into four categories using the quartiles of the distribution.
To investigate contextual effects on a more local scale than heretofore, we used the smallest existing geographical units. However, some of those units contain no or very few inhabitants. To define the socioeconomic contextual factor in areas containing a minimum number of inhabitants, we used an algorithm to group non-inhabited parcels with adjacent ones. This algorithm aggregates each parcel with <100 inhabitants aged 5089 years with the one of the adjacent parcels that
In 1996, the median number of inhabitants in the resulting 1672 local areas was 559 (interquartile range (IQR) 391801; interdecile range (IDR) 3131116), and the median number of individuals aged 5079 years was 169 (IQR 120246; IDR 98348). The median area size was 0.4 km2. Its IQR (0.12.2) and IDR (0.124.7) reflect the larger size of areas in rural settings.
Two contextual factors were considered: population density and socioeconomic environment. Area population density was computed as the number of 5079-year-old residents/km2. To define population density strata, we divided the variable into three categories comprising a similar number of individuals. Socioeconomic environment was defined as area mean income of the 5079-year-old residents at 1 year before baseline. It was divided into four categories that included a similar number of individuals.
Statistical analysis
To quantify area-level variations in IHD mortality, we estimated multilevel Cox proportional hazards regression models,37 with individuals nested within local areas. Individuals were followed up until death, or until the end of a 7-year period (31 December 2002). Such a model extends the usual Cox model by incorporating a random interceptthat is, by allowing the intercept value to vary from one area to another. Including an area-level random intercept in the linear predictor, this model is able to account for the fact that individuals residing in the same local area have more similar survival outcomes than individuals from different areas. Assuming that the area-level random intercept is normally distributed, the model estimates its variance, which allows one to both control for within-area correlation and assess the magnitude of between-area variations. As we found no evidence of variations between areas for the effects of individual variables, only the random intercept was included in the model. Using the R software,38 the random-effect Cox regression models were estimated with a penalised likelihood method.39
To facilitate assessment of area differences,5,6 we first expressed variations on the hazard ratio (HR) scale with an interquartile hazard ratio (IHR)7 and a median hazard ratio (MHR).35,40 Simulating area-level residuals from their distribution, the IHR expresses the difference in IHD mortality between the 25% of all individuals in areas with the lowest risk and the 25% of all individuals in areas with the highest risk.7 The MHR is the median value of the HR between the individual in the area with the lowest risk and the individual in the area with the highest risk when randomly picking out two individuals in different areas.35,40 Secondly, following previous literature,41,42 we computed the intra-area correlation of the logarithmic survival times as an indicator of the magnitude of correlation of the outcome within areas.
After estimating an empty multilevel survival model (ie comprising no explanatory variables) for the whole sample (aged 5079 years), we estimated an empty model in each age group database separately (5064-year-old and 6579-year-old individuals). In each age group, separate empty models were also fitted in the three population density strata. We then introduced the characteristics of individuals into the models, and in a final step the socioeconomic contextual variable. The proportionality assumption of the Cox models was tested and found to be valid for each independent variable. For a correct assessment of the magnitude of the individual income effect and socioeconomic contextual effect, we used specific cut-offs in each sample (ie in each age-group database and in each population-density stratum) to divide the variables into four categories.
| RESULTS |
|---|
An empty Cox multilevel model fitted to the whole study sample (aged 5079 years) indicated significant variations in IHD mortality between areas (area-level variance 0.14, p<0.001; IHR 2.36; MHR 1.43; intra-area correlation 0.08; not reported in a table). However, as shown by empty models estimated separately in each age group (5064 and 6579 years at baseline), a much larger geographical variability was seen among 5064-year-old people than among older people, as assessed with the area-level variance, the IHR, the MHR, and the intra-area correlation coefficient (tables 1
and 2
). The area-level variance was more than four times larger in the youngest group; the IHR was 2.91 among people aged 5064 years, but only 1.71 among those aged 6579 years.
|
|
Individual variables were then introduced into the models. In both age groups, models that were not stratified by population density indicated a higher IHD mortality among men, older individuals, non-cohabiting persons, individuals with low education and manual workers (table 3
). We also found a graded increase in mortality risk with decreasing 20-year cumulated income of individuals. All individual effects tended to be larger among 5064 than among 6579-year-old individuals. Inclusion of individual factors led to a 70% and 57% decrease in the area-level variance in the two age groups (not reported in a table).
|
|
|
| DISCUSSION |
|---|
As a result, IHD mortality of non-elderly residents of urban territories depended almost as much on the socioeconomic environment as on the 20-year cumulated individual income.
The main strengths of our study include the definition of a large cohort of individuals, the geocoding of participants to highly local areas and the adjustment of models for several socioeconomic indicators over the adult age. However, there were limitations to our study. Firstly, the low number of IHD deaths among individuals aged 5064 years at baseline (about 99% of observations were right censored) may have resulted in a weak statistical power to detect between-area variations and socioeconomic contextual effects on IHD. Secondly, and most importantly, we had no information on health-damaging behaviour, healthcare utilisation or social support, preventing us from investigating the mediating mechanisms of the contextIHD association. Even if pioneer studies have started exploring those intermediate processes,16,17,24,28 it is a common view that the appropriate way to identify contextual effects is to estimate correlations between contextual factors and health adjusted for classic individual socioeconomic characteristics. We believe that this approach cannot be considered to be explanatory in itself, as a certain part of the associations between contextual characteristics and diseases probably results from the confounding effect of factors that are systematically neglected. For example, a low control on work organisation combined with a high job demand43,44 or a high demand combined with a low reward45 may increase incidence of IHD; however, no published study has adjusted contextual effects for this factor, which is unevenly distributed among neighbourhoods but not part of the causal pathway between residential context and IHD.16 Accordingly, we believe that it is unwise to conclude that truly contextual influences on IHD exist before some of the main mediating processes of the socioeconomic contextual effect are clearly identified.
Geographical variations in IHD were almost completely captured by individual and contextual socioeconomic factors, suggesting that IHD has a strong social dimension. Our study expands on previous literature in showing that the magnitude of geographical variations in IHD mortality and the strength of socioeconomic contextual effects depend on the age of individuals and on the population density of the residential area.
Firstly, we found that socioeconomic contextual effects on IHD mortality were larger among 5064-year-old individuals than among older individuals.11,22,26 This may be a consequence of the development pattern of IHD over the life course, characterised by an accelerated deterioration of the cardiovascular system and a premature onset of the disease for individuals who are cumulating detrimental exposures related to their individual and contextual socioeconomic characteristics. The more frequent and evenly distributed incidence of the disease at older ages, combined with the better health of survivors in deprived areas,22,23 may explain the weaker socioeconomic contextual effect on IHD mortality observed among elderly people. Therefore, following a perspective presented in a previous work,5,46 it is particularly the public health interventions targeted at premature IHD that need to be contextualisedthat is, adapted to the conditions prevailing in each area.
What is already known
|
Secondly, stronger effects of the socioeconomic environment were found in densely than in sparsely populated territories. It indicates that socioeconomic contextual disparities in IHD mortality only appear within urban territories. Possibly, the clustering of deprived individuals in an area may create detrimental conditions (in terms of resources available and shared values) affecting health-damaging and healthcare-seeking behaviour only when a certain degree of concentration of social disadvantage over a given surface is reached. The burden of contextual poverty may be more deeply felt in densely populated territories in which observable signs of material or social deprivation are more visible.47,48 Complementarily, environmental hazards such as noise or air pollution may be more common in deprived urban neighbourhoods than in any other place.4951
What this study adds
|
Combined, the modifying effects of age and population density lead to a particularly marked geographical variability and strong socioeconomic contextual effect on IHD mortality for non-elderly residents of urban territories. In this population group, the socioeconomic contextual effect was almost as high as the effect of individual 20-year cumulated income. It may be particularly within urban territories that the socioeconomic deprivation of the context confronts individuals with life conditions that result in a premature onset of coronary disease. Therefore, non-elderly residents of deprived urban neighbourhoods constitute a major target for both contextual epidemiology of coronary disease and public health interventions aimed at reducing the detrimental effects of the social environment on IHD risk.
Besides the obvious relevance of global policies dealing with the issue of social disparities, the knowledge available on the relationship between the residential environment and IHD may be too scarce to allow the implementation of definite interventions targeting specific contextual influences on IHD. In any case, socioeconomic contextual effects on premature IHD mortality are strong enough within urban territories to deserve high priority research to elucidate their origin.
| ACKNOWLEDGEMENTS |
|---|
| FOOTNOTES |
|---|
| REFERENCES |
|---|
Relevant Article
This article has been cited by other articles:
![]() |
B Chaix, M Lindstrom, M Rosvall, and J Merlo Neighbourhood social interactions and risk of acute myocardial infarction J Epidemiol Community Health, January 1, 2008; 62(1): 62 - 68. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS | REGISTER |