|
|
||||||||
RESEARCH AND PRACTICE |
The authors are with the Department of Society, Human Development and Health, Harvard School of Public Health, Boston, Mass.
Correspondence: Requests for reprints should be sent to S. V. Subramanian, PhD, Department of Society, Human Development and Health, Harvard School of Public Health, 677 Huntington Ave, Kresge 7th Floor, Boston MA 02115-6096 (e-mail: svsubram{at}hsph.harvard.edu).
| ABSTRACT |
|---|
|
|
|---|
Objectives. We analyzed neighborhood heterogeneity in associations among mortality, race/ethnicity, and area poverty.
Methods. We performed a multilevel statistical analysis of Massachusetts all-cause mortality data for the period 1989 through 1991 (n=142836 deaths), modeled as 79813 cells (deaths and denominators cross-tabulated by age, gender, and race/ethnicity) at level 1 nested within 5532 block groups at level 2 within 1307 census tracts (CTs) at level 3. We also characterized CTs by percentage of the population living below poverty level.
Results. Neighborhood variation in mortality across CTs and block groups was not accounted for by these areas age, gender, and racial/ethnic composition. Neighborhood variation in mortality was much greater for the Black population than for the White population, largely because of CT-level variation in poverty rates.
Conclusions. Neighborhood heterogeneity in the relationship between mortality and race/ethnicity in Massachusetts is statistically significant and is closely related to CT-level variation in poverty.
| INTRODUCTION |
|---|
|
|
|---|
Why would an empirical evaluation of contextual heterogeneity by population subgroups (e.g., White and Black) be important? Patterns of all-cause mortality are shaped by a complex constellation of individual as well as contextual factors that may conceivably vary for Whites and Blacks, as suggested, for example, by the fact that leading causes of death differ for different racial/ethnic groups.11 An investigation of race/ethnicityspecific contextual heterogeneity in mortality at the neighborhood level can give insight into the relative importance to mortality of individual and contextual factors for different racial/ethnic populations. From a methodological standpoint, addressing this question requires that context be included as an intrinsic part of analytical strategies, as achieved by multilevel statistical models.12,13
To date, however, multilevel methods in public health research have principally been applied to estimate the "average" effect of a predictor measured at an area level on individual-level health outcomes, including mortality.14,15 While this approach is important, it can potentially obscure the contextual heterogeneities underlying such average effects. We examined contextual heterogeneity in mortality in relation to race/ethnicity and area poverty by analyzing 19891991 Massachusetts mortality data, prepared for the US Public Health Disparities Geocoding Project,16 in conjunction with 1990 census data. Specifically, we addressed 3 questions: (1) What is the neighborhood variation in mortality rates for different racial/ethnic groups? (2) What is the magnitude of neighborhood variation in the average BlackWhite disparity? (3) Does neighborhood variation in poverty account for the racial/ethnic contextual variation and racial/ethnic disparities in mortality?
To characterize individuals neighborhood context, we employed 2 levels of census geography: the census tract (CT), which on average contains 4000 persons, and the census block group (BG), a subdivision of the CT, which on average contains 1500 persons.17 The 2-fold appeal of CTs is that "when first delineated, [they] are designed to be homogeneous with respect to population characteristics, economic status, and living conditions,"17(ppG10G11) and that, once created, they constitute administrative units used by federal, state, and local governments, including public health departments, to characterize jurisdictions, determine eligibility for diverse programs, and allocate resources.1820 Therefore, CTs have real-life implications for their residents.
| METHODS |
|---|
|
|
|---|
Analytical Structure
The data had a hierarchical multilevel structure of 79813 cells at level 1, consisting of individuals in numerators and denominators cross-tabulated by age, gender, and race/ ethnicity, which were nested within 5532 BGs at level 2, nested within 1307 CTs at level 3. Each BG had between 1 and 30 cells, with the numerator specifying the number of persons who died and the denominator providing the total population used to calculate the proportion of deaths in each cell. In the analytical model we considered both BG and CT levels; however, given the greater policy appeal for health monitoring of populations, our substantive interest was at the CT level. Structurally, our models were identical to models with individuals at level 1.22
Response and Predictors
The response, mortality, is defined as a proportion: the number of deaths as a proportion of the total population for each cell. The cell predictor variables relate to 5 age categories (014 years, 1524 years, 2544 years, 4564 years, and 65 years and older), 2 gender categories (male and female); and 3 racial/ ethnic categories (White, Black, and all others, each including non-Hispanics and Hispanics; these 3 groups made up, respectively, 89.9%, 4.9%, and 5.1% of the population enumerated in Massachusetts in the 1990 census). Because of the heterogeneity of the population in "all others," we focused our interpretation mainly on the WhiteBlack comparisons. On the basis of the cross-tabulation of age, gender, and racial/ethnic categories, we obtained 30 unique population groups, or cells, as shown in Table 1
.
|
Statistical Analysis
We employed multilevel statistical procedures2729 because of their ability to model complex variance structures at multiple levels. The principles underlying multilevel modeling procedures are now well known,13 and in the context of the present analysis they allow estimation of the relationship between mortality and race/ethnicity, conditional on individual age and gender variations ("fixed parameters") and CT- and BG-level variations ("random parameters"). They also enable an estimation of the extent to which the relationship between mortality and race/ethnicity varies across CTs (random parameters) and the degree to which CT poverty explains this variation (fixed parameters).
The response variable, proportion of deaths in each cell, was modeled with allowances made for the varying denominator in each cell.27 The fixed and random parameter estimates (along with their standard errors) for the 3-level binomial logit link model were calibrated using predictive/penalized quasi-likelihood procedures with second order Taylor series expansion,30 as implemented within the MLwiN program.31 We allowed for extrabinomial variations at level 1, because proportions may exhibit more or less variation than a binomial distribution.32 We calibrated 3 models.
Model 1. Model 1 was a 3-level model of cells (level 1) within BGs (level 2) within CTs (level 3), with cell characteristics related to age, gender, and race/ethnicity specified in the fixed part of the model and a residual variation estimated at the CT and BG levels in the random part. The fixed-part specification included the main effects for the age, gender, and racial/ethnic categories, along with a second-order interaction between agegender, agerace/ethnicity, and genderrace/ ethnicity. We did not find any empirical support for third-order interaction terms. The estimates from model 1 allowed an evaluation of the racial/ethnic differences in mortality and the magnitude of variation in mortality at the BG and CT levels, conditional on the relationship between mortality and age, gender, and race/ethnicity within each BG and CT.
Model 2. Model 2 was similar to model 1, but it allowed the fixed racial/ethnic differential on mortality to vary across CTs in the random part to obtain differential CT-level variation in mortality for Whites, Blacks, and others. CT-level random parameters from model 2 were used to test the hypothesis pertaining to race/ethnicitybased contextual heterogeneity, that is, whether neighborhood-level variation in mortality was different for different racial/ethnic groups.
Model 3. Model 3 was similar to model 2, but it included a fixed cross-level interaction effect between CT-level poverty and individual race/ethnicity. In this way, we ascertained the relationship between neighborhood poverty, individual race/ethnicity, and mortality, as well as the extent to which CT-level poverty accounted for the CT-level racial/ethnic variation in mortality.
| RESULTS |
|---|
|
|
|---|
2
0 = 0.09; Table 3
|
|
|
The odds ratio for mortality increased with neighborhood poverty, and the relationship was substantially stronger for Blacks than for Whites (Table 2
, model 3). When Whites living in the 3 higher CT-poverty strata were compared with the reference group (Whites living in CTs with lower than 5% poverty), their odds ratios were 1.05 (CTs with 5%9.9% poverty), 1.23 (CTs with 10%19.9% poverty), and 1.42 (CTs with 20% or higher poverty). For Blacks, however, compared with the same reference group (Whites living in CTs with lower than 5% poverty), the odds ratios for the unique interactions were, respectively, 1.67 (CTs with 5%9.9% poverty), 2.34 (CTs with 10%19.9% poverty), and 3.00 (CTs with 20% or higher poverty).
| DISCUSSION |
|---|
|
|
|---|
Our findings suggest a need to conceptualize contextual effects (e.g., neighborhood poverty) in more complex ways. Typically, the examination of contextual effects has proceeded with an assumption of "main contextual effects" (e.g., that neighborhood poverty affects individual mortality in a similar manner for all racial/ethnic groups). However, it is entirely reasonable (and perhaps more realistic) to anticipate that contextual differences as well as contextual effects inherently interact with individual characteristics. Our finding that neighborhood-level poverty directly contributes to the greater geographic heterogeneity in mortality rates for Blacks suggests that the consequences of neighborhood deprivation may be particularly exacerbated for Blacks, compared with Whites. That is, there is a joint and synergistic shaping of population patterns of mortality by individual (e.g., race/ethnicity) and contextual (e.g., CT poverty) factors.
The following caveats should be considered in interpreting the empirical findings of our study. First, there is a distinct possibility of a potential misspecification of the individual-level factors that are associated with mortality. Our specification of individual variables was constrained by the degree of cross-tabulation (i.e., by age, race/ethnicity, and gender) that was possible on the numerator and denominator populations simultaneously at the BG level. This constraint precluded estimation of the impact of other important (unmeasured) socioeconomic markers on clustering in mortality.
Had we constructed our cells by county, it would have been possible to include educational data at the individual level (thus increasing the "social resolution" at the individual level), but this would have come at the cost of ignoring the critical level of BGs and CTs (thus decreasing the "spatial resolution"). While it is difficult to evaluate the sensitivity of our findings, given the lack of data, it is likely that including additional socioeconomic markers might have attenuated both the magnitude of racial/ethnic disparities in mortality and the effect of CT poverty. It is unlikely, however, that such data would have accounted for the contextual heterogeneity in racial/ethnic disparities in mortality. For this to have happened, the omitted individual socioeconomic markers not only would have to be to perfectly collinear with individual race/ ethnicity but would also need to have exactly the same spatial clustering as race/ethnicity.
It is, however, worth noting that neighborhood effects continue to be thought of as an extension of individual effects on health, as is reflected in the overbearing concern for "controlling" for individual confounders. Such an approach, arguably, obscures the need to reconceptualize the very notion of "individual effects" that results from incorporating neighborhoods in our conceptual model to elucidate health disparities. The differential response to neighborhood environments by individual Blacks and Whites, demonstrated in our analysis, is a simple illustration of this point.
A second limitation is our focus on fixed, discrete, and hierarchical census-defined contexts to define neighborhoods. While these are important, there may be other nonspatial contexts (e.g., households, nongeographic communities) and nonhierarchical contexts (e.g., workplaces and subjectively defined neighborhoods) that are also important to understanding the patterning of mortality. Thus, in addition to the issue of missing covariates, the issue of missing "levels" remains, especially in multilevel models.
Third, our characterization of neighborhoods was based on a single dimensionpoverty. Future researchers might consider systematically developing a "typology" of neighborhoods to take into account the complex interaction of multiple characteristics.
A fourth concern is that it is problematic to demonstrate a "neighborhood effect" on mortality on the basis of cross-sectional observational data, since the magnitude of both neighborhood heterogeneity and neighborhood poverty is for individuals at the time of their death; therefore, the issue of potential mobility (social as well as spatial) over the life course is ignored.
Fifth, the partitioning of variation by different levels (e.g., individual, BG, and CT) in logistic models is not straightforward.35 Specifically, the magnitude of variance at the BG and CT levels cannot be evaluated for its size in relation to the individual-level variance, since the latter is a known function in logistic models. Thus, the magnitude of CT-level contextual heterogeneity in mortality by race/ ethnicity reported here is simply a conditional (not absolute or relative) estimate, that is, conditional on the fixed part and the random part of the model at the BG level.
Finally, our findings can be generalized only to populations (of neighborhoods and individuals) that are similar in their characteristics to those in Massachusetts. The replicability of our results in other states or geographic settings is an empirical question warranting further research.
Despite these challenges, quantifying and studying the distribution of mortality across different local geographic units disaggregated by key sociodemographic and economic markers may provide important input into social and health policy as well as sharpen etiological analysis of variations in mortality. For example, our results suggest that while the routine monitoring of how population groups are doing needs to be area-specific, the monitoring of how areas are doing (a prerequisite for developing area-based public health interventions) needs to be population group specific. Such multilevel thinking has yet to permeate the research geared toward public health practice. Such considerations, we believe, may lead to a fairer, and clearly more sensitive and realistic, means of target-setting and evaluation for population groups and areas. If we rely only on across-the-board, average targets for areas, it is difficult to distinguish improvements consisting of changes only in groups that are already "better-off" in those areas. The analytical approach developed here offers a framework for understanding such public health issues.
In conclusion, our analysis points to the importance of generating a quantitative description of the varying impact of neighborhood on health, in relation to diverse population subgroups, rather than focusing solely on average differences in health disparities across these groups. Extending this approach to investigating health patterns over time could provide opportunities to quantify the impact of various social policies on health inequalities between population groups and areas. Indeed, the complex geographic and social variation in important health outcomes such as mortality is an intrinsically important attribute of society that needs to be routinely described and understood.
| Acknowledgments |
|---|
We thank Mah-Jabeen Soobader for her contribution to earlier versions of these analyses; Daniel Friedman (former assistant commissioner, Bureau of Health Statistics, Research and Evaluation, Massachusetts Department of Health) for facilitating the conduct of this study; and Alice Mroszczyk (Massachusetts Department of Public Health) and Elaine Trudeau and Charlene Zion (Registry of Vital Records and Statistics) for helping us access the mortality data.
Human Participant Protection
Use of the data employed in this study was approved by all relevant institutional review boards/human subjects committees at the Harvard School of Public Health and the Massachusetts Department of Public Health.
| Footnotes |
|---|
Contributors
S. V. Subramanian originated the study, analyzed and interpreted the data, and wrote the article. J. T. Chen contributed to the data preparation, interpretation of the results, and editing of the article. D. H. Rehkopf assisted with data analysis and presentation. P.D. Waterman obtained and arranged geocoding of the mortality data. N. Krieger contributed to the interpretation of the results and editing of the article.
Accepted for publication January 13, 2004.
| References |
|---|
|
|
|---|
2. Casper ML, Barnett E, Halverson JA, et al. Women and Heart Disease: An Atlas of Racial and Ethnic Disparities in Mortality. Morgantown: Office for Social Environment and Health Research, West Virginia University; 2000.
3. Casper ML, Barnett E, Williams GI Jr, Halverson JA, Braham VE, Greenlund KJ. Atlas of Stroke Mortality: Racial, Ethnic and Geographic Disparities in the United States. Atlanta, Ga: Centers for Disease Control and Prevention; 2003.
4. Devesa SS, Grauman DJ, Blot WJ, Pennello GA, Hoover RN, Fraumeni JF Jr. Atlas of Cancer Mortality in the United States, 195094. Bethesda, Md: National Cancer Institute; 1999.
5. Geronimus AT, Bound J, Waidmann TA, Hillemeier MM, Burns PB. Excess mortality among Blacks and Whites in the United States. N Engl J Med. 1996; 335:15521558.
6. Murray CJL, Michaud CM, McKenna MT, Marks JS. US Patterns of Mortality by County and Race: 19651994. Cambridge, Mass: Harvard Center for Population and Development Studies; 1998. Also available at: http://www.hsph.harvard.edu/organizations/bdu/images/usbodi. Accessed September 10, 2004.
7. US Decennial Life Tables for 1989 to 1991. Hyattsville, Md: National Center for Health Statistics; 1997. Also available at: http://www.cdc.gov/nchs/products/pubs/pubd/lftbls/decenn/1991-89.htm. Accessed September 10, 2004.
8. Pickle LW, Mungiole M, Jones GK, White AA. Atlas of United States Mortality. Hyattsville, Md: National Center for Health Statistics; 1996.
9. Singh KG, Yu SM. Trends and differentials in adolescent and young adult mortality in the United States, 1950 through 1993. Am J Public Health. 1996;86: 560564.
10. Health, United States, 2002 With Chartbook on Trends in the Health of Americans. Hyattsville, Md: National Center for Health Statistics; 2002.
11. Kington RS, Nickens HW. Racial and ethnic differences in health: recent trends, current patterns, future directions. In: Smelser NJ, Wilson WJ, Mitchell F, eds. America Becoming: Racial Trends and Their Consequences. Washington, DC: National Academy Press; 2001.
12. OCampo P. Invited commentary: advancing theory and methods for multilevel models of residential neighborhoods and health. Am J Epidemiol. 2003; 157(1):913.
13. Subramanian SV, Jones K, Duncan C. Multilevel methods for public health research. In: Kawachi I, Berkman LF, eds. Neighborhoods and Health. New York, NY: Oxford University Press; 2003.
14. Yen IH, Kaplan GA. Neighborhood social environment and risk of death: multilevel evidence from the Alameda County Study. Am J Epidemiol. 1999;149: 898907.
15. Cubbin C, LeClere F, Smith G. Socioeconomic status and injury mortality: individual and neighbourhood determinants. J Epidemiol Community Health. 2000;54: 517524.
16. Krieger N, Chen J, Waterman PD, Soobader MJ, Subramanian SV, Carson R. Geocoding and monitoring US socioeconomic inequalities in mortality and cancer incidence: does choice of area-based measure and geographic level matter? The Public Health Disparities Geocoding Project. Am J Epidemiol. 2002;156:471482.
17. US Census Bureau. Geographical Areas Reference Manual. Available at: http://www.census.gov/geo/ www/garm.html. Accessed February 5, 2004.
18. Race and the Elimination of Health Disparities in the City of Boston: Promoting the Health of the Underserved. Boston, Mass: Boston Public Health Commission; 2002.
19. Health Resources and Services Administration. Health Professional Shortage Areas. 2002. Available at: http://bhpr.hrsa.gov/shortage. Accessed September 17, 2004.
20. Monmonier M. Cartographies of Danger: Mapping Hazards in America. Chicago, Ill: University of Chicago Press; 1997.
21. Krieger N, Waterman P, Lemieux K, Zierler S, Hogan JW. On the wrong side of the tracts? Evaluating the accuracy of geocoding in public health research. Am J Public Health. 2001;91:11141116.[Abstract]
22. Subramanian SV, Duncan C, Jones K. Multilevel perspectives on modeling census data. Environ Plan A. 2001;33:399417.[CrossRef]
23. Krieger N, Chen J, Waterman PD, Soobader MJ, Subramanian SV, Carson R. Choosing area-based socioeconomic measures to monitor social inequalities in low birthweight and childhood lead poisoning: the Public Health Disparities Geocoding Project. J Epidemiol Community Health. 2003;57(3):186199.
24. Krieger N, Waterman PD, Chen J, Soobader MJ, Subramanian SV. Monitoring socioeconomic inequalities in sexually transmitted infections, tuberculosis, and violence: geocoding and choice of area-based socioeconomic measuresThe Public Health Disparities Geocoding Project (US). Public Health Rep. 2003; 118(3):240260.[CrossRef][ISI][Medline]
25. US Census Bureau. Poverty areas. Available at: http://www.census.gov/population/socdemo/statbriefs/povarea.html. Accessed February 5, 2004.
26. Census of Population and Housing, 1990: Summary Tape File 3 on CD-ROM Technical Documentation. Washington, DC: US Census Bureau; 1991.
27. Goldstein H. Multilevel Statistical Models. 3rd ed. London, UK: Arnold Publishers; 2003.
28. Longford N. Random Coefficient Models. Oxford, UK: Clarendon Press; 1993.
29. Raudenbush S, Bryk A. Hierarchical Linear Models: Applications and Data Analysis Methods. Thousand Oaks, Calif: Sage Publications; 2002.
30. Goldstein H, Rasbash J. Improved approximations for multilevel models with binary responses. J R Stat Soc Ser A Stat Soc. 1996;159:505513.[CrossRef]
31. Rasbash J, Browne W, Goldstein H, et al. A Users Guide to MLwiN, Version 2.1. London, UK: Multilevel Models Project, Institute of Education, University of London; 2000.
32. Collet D. Modelling Binary Data. London, UK: Chapman & Hall; 1991.
33. Krieger N, Rowley DL, Herman AA, Avery B, Phillips MT. Racism, sexism, and social class: implications for studies of health, disease, and well-being. Am J Prev Med. 1993;9(suppl):82122.[ISI][Medline]
34. Williams DR, Collins C. US socioeconomic and racial differences in health: patterns and explanations. Annu Rev Sociol. 1995;21:349386.[CrossRef][ISI]
35. Goldstein H, Browne WJ, Rasbash J. Partitioning variation in multilevel models. Understanding Statistics. 2002;1:223232.
This article has been cited by other articles:
![]() |
N. Volkova, W. McClellan, M. Klein, D. Flanders, D. Kleinbaum, J. M. Soucie, and R. Presley Neighborhood Poverty and Racial Differences in ESRD Incidence J. Am. Soc. Nephrol., February 1, 2008; 19(2): 356 - 364. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. O'Campo, J. G. Burke, J. Culhane, I. T. Elo, J. Eyster, C. Holzman, L. C. Messer, J. S. Kaufman, and B. A. Laraia Neighborhood Deprivation and Preterm Birth among Non-Hispanic Black and White Women in Eight Geographic Areas in the United States Am. J. Epidemiol., January 15, 2008; 167(2): 155 - 163. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Riva, L. Gauvin, and T. A Barnett Toward the next generation of research into small area effects on health: a synthesis of multilevel investigations published since July 1998 J. Epidemiol. Community Health, October 1, 2007; 61(10): 853 - 861. [Abstract] [Full Text] [PDF] |
||||
![]() |
M.-H. Kim, S V Subramanian, I. Kawachi, and C.-Y. Kim Association between childhood fatal injuries and socioeconomic position at individual and area levels: a multilevel study J. Epidemiol. Community Health, February 1, 2007; 61(2): 135 - 140. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. H. Rehkopf, L. T. Haughton, J. T. Chen, P. D. Waterman, S.V. Subramanian, and N. Krieger Monitoring Socioeconomic Disparities in Death: Comparing Individual-Level Education and Area-Based Socioeconomic Measures Am J Public Health, December 1, 2006; 96(12): 2135 - 2138. [Abstract] [Full Text] [PDF] |
||||
![]() |
A M A Smith and S V Subramanian Population contextual associations with heterosexual partner numbers: a multilevel analysis. Sex. Transm. Inf., June 1, 2006; 82(3): 250 - 254. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Chaix, M. Rosvall, J. Lynch, and J. Merlo Disentangling contextual effects on cause-specific mortality in a longitudinal 23-year follow-up study: impact of population density or socioeconomic environment? Int. J. Epidemiol., June 1, 2006; 35(3): 633 - 643. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |