Email updates

Keep up to date with the latest news and content from IJHG and BioMed Central.

Open Access Highly Accessed Research

Urban sprawl, obesity, and cancer mortality in the United States: cross-sectional analysis and methodological challenges

David Berrigan1*, Zaria Tatalovich2, Linda W Pickle3, Reid Ewing4 and Rachel Ballard-Barbash1

Author Affiliations

1 Applied Research Program, Division of Cancer Control and Population Sciences, National Cancer Institute, Bethesda, MD 20892, USA

2 Surveillance Research Program, Division of Cancer Control and Population Sciences, National Cancer Institute, Bethesda, MD 20892, USA

3 StatNet Consulting, Gaithersburg, MD 20879, USA

4 The University of Utah, College of Architecture and Planning, Salt Lake City, UT 84112, USA

For all author emails, please log on.

International Journal of Health Geographics 2014, 13:3  doi:10.1186/1476-072X-13-3

The electronic version of this article is the complete one and can be found online at: http://www.ij-healthgeographics.com/content/13/1/3


Received:17 October 2013
Accepted:25 December 2013
Published:6 January 2014

© 2014 Berrigan et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Abstract

Background

Urban sprawl has the potential to influence cancer mortality via direct and indirect effects on obesity, access to health services, physical activity, transportation choices and other correlates of sprawl and urbanization.

Methods

This paper presents a cross-sectional analysis of associations between urban sprawl and cancer mortality in urban and suburban counties of the United States. This ecological analysis was designed to examine whether urban sprawl is associated with total and obesity-related cancer mortality and to what extent these associations differed in different regions of the US. A major focus of our analyses was to adequately account for spatial heterogeneity in mortality. Therefore, we fit a series of regression models, stratified by gender, successively testing for the presence of spatial heterogeneity. Our resulting models included county level variables related to race, smoking, obesity, access to health services, insurance status, socioeconomic position, and broad geographic region as well as a measure of urban sprawl and several interactions. Our most complex models also included random effects to account for any county-level spatial autocorrelation that remained unexplained by these variables.

Results

Total cancer mortality rates were higher in less sprawling areas and contrary to our initial hypothesis; this was also true of obesity related cancers in six of seven U.S. regions (census divisions) where there were statistically significant associations between the sprawl index and mortality. We also found significant interactions (p < 0.05) between region and urban sprawl for total and obesity related cancer mortality in both sexes. Thus, the association between urban sprawl and cancer mortality differs in different regions of the US.

Conclusions

Despite higher levels of obesity in more sprawling counties in the US, mortality from obesity related cancer was not greater in such counties. Identification of disparities in cancer mortality within and between geographic regions is an ongoing public health challenge and an opportunity for further analytical work identifying potential causes of these disparities. Future analyses of urban sprawl and health outcomes should consider exploring regional and international variation in associations between sprawl and health.

Keywords:
Urban sprawl; Cancer mortality; Obesity; Spatial heterogeneity; Census division; County; Health disparities; Ecological analysis

Background

Urban areas are growing rapidly in population and land area throughout the world [1]. Much of this growth falls into a pattern characterized as ‘urban sprawl’ [2], an economic and social process associated with low residential density, segregation of land use, and automobile dependence [3]. A number of studies have reported that urban sprawl is associated with health outcomes, health behaviors, and obesity [3-7]. In cross-sectional studies, greater levels of urban sprawl are associated with increased obesity, poorer mental health, more traffic accidents, and lower levels of physical activity including walking [4,5,8-11]. These adverse associations with health could be due to increased exposure to pollution, changes in behaviors that directly influence disease risk such as smoking, sedentary time, or social behavior, or through the association of sprawl with use of health care systems [12]. Alternatively, selection of residential environments because of preference or economic constraints could also produce some of these associations. Two studies, one of childhood obesity [13] and one of walking in adults [9] report statistically significant associations between sprawl and energy balance outcomes in cross-sectional analyses but not in longitudinal ones.

Research on urban sprawl has been particularly prominent in the United States, in part because this kind of development has been very common, notably in the past 30 years [14]. However, international interest in the environmental and health consequences of sprawl is also increasing [15,16]. Independent of the analysis of urban sprawl, there is also a long and global history of examining the effect of cities and residence in urban areas on health outcomes, including cancer [17-19]. In the 19th century, scrotal cancer was associated with soot exposure [20] and more recently, breast and lung cancer rates have sometimes but been reported to be higher in urban areas than rural or suburban areas [21-24]. Wilson and Chakraborty [25] highlight the lack of comparative studies across national boundaries concerning the environmental consequences of sprawl and such studies of health outcomes and health behaviors are also rare [25]. Comparative studies of urban sprawl and health could lead to insights arising from how sprawl develops and manifests itself in different countries.

Comparative studies of geographic areas can be used to examine the contributions of different risk factors to specific cancers. For example, such studies have been used to examine the relative contribution of environmental exposures and reproductive characteristics to breast cancer risk and the contributions of air pollution versus smoking prevalence to lung cancer [22,23]. There is also interest in urban–rural differences and spatial variation in stage at diagnosis. For example a study using SEER data reports that urban as opposed to rural residence is associated with stage at diagnosis for lung and colorectal cancer [26] and there is also substantial geographic variation late stage breast and colorectal cancer diagnosis [27,28]. Access to care and the distribution of race-ethnic groups are associated with county level variation in stage at diagnosis identified in some of these studies, but considerable unaccounted for variation between counties remains. Despite their limitations, ecological analyses or cross case comparisons are sometimes one of the best and most feasible ways to examine the contributions of environmental and contextual variables to health outcomes and to identify targets for focused intervention studies [28-31].

The main goal of this paper is to examine cross-sectional associations between urban sprawl and cancer mortality, accounting for spatial variation in mortality and the effects of obesity, smoking, measures of health care services and access, and socioeconomic status. Because obesity is a risk factor for cancer at a number of sites [32] and because the prevalence of obesity is greater in more sprawling regions of the US [4], we hypothesized that mortality for obesity related cancers would be higher in more sprawling areas. In addition to its association with obesity prevalence, sprawl may also have an effect on access to care. For example; increased travel time could influence access to screening and treatment facilities [33]. Thus our overall conceptual model involves a socio-ecological approach and includes urban form, modifiable risk factors, and demographic variables, as well as measures related to care. In this paper we chose to examine associations between urban sprawl and obesity-related cancer mortality. Future studies could address specific cancers, cancer related to other risk factors such as smoking, and cancer incidence.

Our analysis emphasizes the potential role of urban sprawl as an influence on obesity related cancer because of the extensive literature exploring the association between the built environment and energy balance that has appeared over the past decade [4,5,13]. However, recent reviews also suggest that the relationship between obesity and environmental variables is heterogeneous, with a lack of consistent associations across studies [7,34,35]. Few past studies of urban sprawl have examined geographic differences in the association between sprawl and health or health behaviors [12,36,37] and therefore residual spatial heterogeneity associated with such differences could account for some of the inconsistency seen in studies of the environment, health and energy balance.

A further goal of this paper is to improve upon previous approaches to the analysis of ecological data so that the effects of spatial autocorrelation are accounted for properly. Spatial autocorrelation, a form of spatial dependency, occurs when observations that are closer together in space are more similar to each other. Most studies of regional and national variation in cancer incidence and mortality do not explicitly address this issue [12,38], despite the consensus amongst spatial statisticians that it is critical to account for spatial dependency to obtain valid analytical results [39]. Instead most past studies emphasize the standard statistical assumptions of regression models such as independence, normality, and homoscedasticity of error.

In light of these goals, the specific aims of this paper are to 1) analyze associations between urban sprawl and cancer mortality in the US and to determine to what extent these associations differ in different regions of the country and 2) more comprehensively assess and address the presence of spatial autocorrelation in these mortality data. We focus on cancer mortality rather than incidence because urban sprawl has the potential to influence determinants of both cancer incidence and mortality via its relationships to risk factors such as obesity and physical activity and via its influence on factors related to cancer prevention and treatment such as transportation and the distribution of facilities.

Methods

This cross-sectional analysis is based on 2002–2006 county level estimates of cancer mortality in the United States of America from the US Surveillance, Epidemiology and End Results (SEER) program [40], along with data originating from the National Center for Health Statistics [41]. We report results for men and women separately for cancer mortality from the 19 and 21 most common cancers in men and women respectively, and for six cancers for which obesity is a major risk factor.

The geographic coverage of this analysis included 935 urban and suburban counties in the US (Figure 1). These counties represent the complete set of continental US counties of interest because the concept of urban sprawl is not relevant in rural counties. They include all metropolitan statistical areas in the continental United States as defined by the US Census bureau and adjacent (‘suburban’) counties. Variables used in our analysis are summarized below and in Table 1. We also used Census Division as a variable in our analyses to explore regional variation in associations between sprawl and mortality. The US census divides the nation into nine contiguous divisions in four regions (Northeast (New England and Middle Atlantic), Midwest (East and West North Central), South (South Atlantic, East and West South Central), and West (Mountain and Pacific)). The divisions contain from three to eight states [42]. Our sample included 34 to 222 counties per division.

thumbnailFigure 1. County level distributions of A) total cancer mortality (2002–2006), B) obesity (2007) and C) urban sprawl (2000) by tertiles for the continental United States.

Table 1. Distributions of covariates and cancer mortality rates

Cancer mortality

County level mortality rates based on data from 2002–2006 for the 19 most common cancers in adult males (oral cavity and pharynx, esophagus, stomach, colon, liver, pancreas, larynx, lung, melanoma, bladder, kidney, brain, thyroid, Hodgkin’s lymphoma, Non-Hodgkin’s lymphoma, myeloma, leukemia, prostate and testis) and the 21 most common cancers in adult females (similar to the list for males except prostate and testis cancer are replaced with cancer of the cervix, uterus, ovary and breast) were obtained from the SEER database [43]. These cancers were selected because they were sufficiently common to allow stable mortality estimates for all the counties considered. After exploratory analyses of mortality from specific cancers, we chose to analyze mortality rates in two categories: overall cancer mortality and obesity-related cancer (colon, endometrium, kidney, esophagus, pancreas, post-menopausal breast). Obesity related cancers were categorized based on the World Cancer Research Fund [32] report on risk factors. Breast cancer mortality was divided into pre and post-menopausal on the basis of age, with breast cancer mortality in women > 50 years old assumed postmenopausal. No attempt was made to distinguish risk factors for squamous cell versus adenocarcinomas. All rates were directly age-adjusted using the 2000 U.S. population standard.

Results for one county, Madison County, MS, were excluded from regression analyses but included in the summary tables. In Madison County, both male and female cancer mortality rates for all cancers combined were > 6 SD from their respective means of the 935 counties considered. A non-profit group, the Susan G. Komen Foundation, previously identified Madison county as having very high cancer mortality rates and the Centers for Disease Control and Prevention (CDC) has attributed these high rates to the presence of a hospice serving all of central Mississippi [44]. Discussion with SEER staff yielded no further clues as to alternative explanations as to why this county had such high rates.

Sprawl index

The urban sprawl index used here is a composite measure of features of urban form related to population density and street accessibility [4]. In brief, the index, which arrays counties from least to most sprawling, has its highest values in urban counties such as the four urban boroughs of New York City or San Francisco County, and its lowest values in the suburbs of sprawling metropolitan regions throughout the country [4,45]. Because of this inverse scaling, some readers may find it easier to think of this as a ‘compactness’ index [46]. Jackson County, adjacent to Topeka KS, had the greatest level of sprawl in this data set. The index is based on six variables extracted from US census data - 1) Gross population density in persons per square mile, 2) Percentage of population living at densities less than 1,500 persons per square mile (low suburban density), 3) Percentage of population living at densities greater than 12,500 persons per square mile (transit-oriented urban density, 4) net urban density (the ratio of population to urbanized land area from US census data), 5) Average block size in square miles, and 6) Percentage of small blocks, roughly 500 feet on a side. The first four of these variables are related to population density and the remaining two variables are related to block characteristics and street connectivity. For this paper, the sprawl index was calculated based on US Census data from 2000.

Obesity

County level obesity prevalence estimates are based on self reported height and weight summarized by the CDC. Estimates for each county are based on Bayesian analysis combining information from the US Behavioral Risk Factor Surveillance System (BRFSS) and the US Census [47,48].

Smoking

Estimates of county level current smoking prevalence for men and women during 2000–2003 were downloaded from the National Cancer Institute (NCI) Division of Cancer Control and Population Sciences (DCCPS) small area estimates web page [49]. These prevalence estimates were calculated using information from the National Health Interview Survey (NHIS) and the BRFSS using a model-based estimation procedure [50]. Note that this estimation procedure uses covariates such as race/ethnicity, age, and education levels. Some of these covariates are also used in our analysis.

Race/ethnicity data

County level prevalence of self-identified Non-Hispanic White, Non-Hispanic Black, Hispanic and American Indian, Asian, and Pacific Islander (AIAPI) race/ethnic groups were obtained from the 2000 US Census. Prevalence of American Indian, Asian and Pacific Islander were combined into a single category. Analyses that excluded American Indians gave qualitatively similar results (not shown).

Socioeconomic position index

We calculated a version of the Krieger 2002 Socioeconomic Position Index (SEP) index [51] for each county. The SEP represents a summary deprivation measure consisting of a standardized z score combining data on percentage working class, unemployment, percentage below US poverty line, low education, expensive homes, and median household income. Our version of the index was based on a single principal component from the Principal Components Analysis (PCA) of the six variables considered. Higher values of the index indicate higher socioeconomic position. For example, the five counties with the lowest values of the SEP are Tunica MS, Hidalgo, Cameron, and Webb Counties (TX), and Terrell County (GA) whereas the five counties with the highest values of SEP are Douglas CO, Loudon VA, Hamilton IN, Howard MD, and Johnson County KS.

Access and availability of care

Information concerning physician density, hospital density, and health insurance coverage were used as proxies for access and availability to health care resources [33]. These variables were extracted from the 2008 release of the Area Resource File (ARF) for the year 2002 [52]. Total active non-federal and federal physicians and number of hospitals per 1000 people in the county were selected for analysis. Prevalence of uninsured aged 0–64 was also calculated from ARF health care and Census population data. Because prevalence of health insurance coverage was positively associated with cancer mortality, contrary to many results at the individual level [53], we also extracted median age of adults from census data in order to determine if residual variation in age structure accounted for the ecological association of insurance prevalence and cancer mortality (It did not, see below).

Statistical analysis

We used SAS JMP and R to initially examine the association between cancer mortality and the independent variables listed above and to check that the data met the assumptions required of the proposed linear regression models [54,55]. The age-adjusted mortality rates were approximately normally distributed so no data transformation of the rates was necessary. Variables relating to race/ethnicity, physician (MD) density, hospital density, and the sprawl index were natural log transformed to improve the linearity of their associations with mortality. Additionally, cancer mortality rates were weighted by county population size in order to stabilize their variances for large and small counties, an assumption also required by the subsequent models. Separate models were fit for males and females and for both of the cancer groupings.

Our overall approach involved fitting an increasingly complex series of normal linear regression models to the age-adjusted mortality rates and then examining to what extent adding more explanatory covariates, interactions and spatial autocorrelation accounted for the observed variation in mortality rates. Important two way interactions between the main effect covariates were identified for inclusion in the models using the Elastic Net approach, a weighted average of the least absolute shrinkage and selection operator (LASSO) and Ridge regression variable selection methods [56]. This method has been shown to be more robust than conventional stepwise regression variable selection methods and overcomes problems of selection from among highly correlated variables by either the LASSO or Ridge methods alone [56,57]. Interactions from this initial variable selection process were retained in the subsequent regression models only if they were significant (p value < 0.05) Because the Elastic Net method evaluates groups of the covariates and interactions, unlike stepwise regression methods, no further adjustment for multiple comparisons was needed. Note that final models included different interactions for the different gender and cancer groupings.

As noted above, our main study aim involved understanding the relationship between urban sprawl and cancer mortality. Therefore, we began by examining simple models estimating the relationship between mortality and urban sprawl without any other covariates. For each of the four groups (Overall and obesity related cancer mortality by gender) six models were fit and examined for the presence or absence of spatial auto-correlation. These models included successively more variables related to cancer causes and spatial structure: Model 1: ln sprawl, Model 2: ln Sprawl + demographic and other covariates, Model 3: ln Sprawl + demographic and other covariates + Census Division, Model 4: ln Sprawl + demographic and other covariates + interactions + Census Division, Model 5: ln Sprawl + demographic and other covariates + interactions + Census Division + random effects to account for county-level spatial autocorrelation, Model 6: ln Sprawl + demographic and other covariates + interactions + Census Division + County random effects + ln Sprawl*Census Division interaction.

All models were implemented using SAS PROC GLIMMIX [58] by pseudo-likelihood estimation with population weights. For the spatial autocorrelation models, a spatial variance component was included as a random effect on the linear predictor (in SAS terms, a G-side effect). This parameter was estimated by a spline smoothing algorithm over a regular grid that spanned the continental US, i.e., with spatial autocorrelation based on distances between all pair-wise county locations. Location was defined by latitude and longitude of the county geographic centroids. Broad regional spatial patterns were assessed by inclusion of indicator variables for Census Divisions as fixed effects in the models.

Because there is no single statistic, such as a likelihood ratio statistic, that adequately measures the goodness of fit of spatial mixed effects models, we examined fit using several different techniques [39,59]. A standard set of residual plots was generated for each model and examined for model fit and consistency with distributional assumptions. We also examined overall fit of the models by tabulating pseudo R2 and Akaike’s Information Criterion (AIC) values. The pseudo R2 statistic measures the proportion of total variance explained by the model (1-(sum of squared errors)/(total sum of squares)). AIC is a measure of the regression sum of squares that includes a penalty for every added covariate [60]. The “best” of the 6 models is the one with the highest R2 and/or the lowest AIC value.

Because one of our objectives was to explore the extent to which these models accounted for spatial patterns and autocorrelation in mortality rates, we also calculated Moran’s I statistic for the standardized residuals of the six models described above using a one-sided permutation test [39]. This statistic does not account for variation in rates due to population differences, so we used it to rank the degree to which each model reduced the spatial patterns apparent in the original data. Successive efforts to account for spatial effects involved models with a) our complete set of covariates to explain spatial patterns at the county level, b) inclusion of US Census Division to explain broader regional patterns, and c) spatial random effects to account for any remaining spatial autocorrelation, i.e., the tendency of mortality rates for neighboring counties to be similar. All of these three elements were required to adequately account for the spatial patterns and autocorrelation in the data. We included the interaction term between our measure of sprawl and Census Division in our final models because it represents an explicit test of regional heterogeneity in the association between sprawl and mortality.

Results

Descriptive statistics

Characteristics of the 935 counties studied here and mortality rates for the cancers examined are presented in Table 1. Note that the mean and standard deviation of the sprawl (‘compactness’) index (100, 25) reflect the standardization procedure used to develop the index [4]. Obesity prevalence and the sprawl index used here were negatively associated (r = −0.36), this indicates lower prevalence of obesity in less sprawling areas because a higher sprawl index corresponds to a denser more urbanized county. Overall, US counties exhibit a wide range of values for all the variables used in our analyses.

The twelve variables considered in this analysis are mostly weakly correlated with one another. The strongest correlation (0.79), between male and female smoking prevalence, represents two variables that are not included in the same models. The second strongest correlation (−0.71) was between proportions of non-Hispanic whites (Treated as the reference group in our models) and blacks. There were statistically significant positive correlations from 0.52- 0.62 between the prevalence of obesity and smoking in men and women, between the log transformed sprawl index and the proportion of American Indian, Asian and Pacific Islanders (AIAPI) and physician density and between the proportion Hispanic and AIAPI. Thus, the proportion of AIAPI was greater in more urban areas. There were significant negative correlations (−0.52—0.56) between prevalence of smoking in females and the proportion Hispanic and AIAPI and between current smoking in males and the SEP. The remaining 66 correlations were not statistically significant. The modest correlations among most of the variables suggest that we have identified a robust set of county level characteristics.

Cancer mortality rates for the four groupings we examined showed substantial variation at the county level (Table 1). The range of mortality rates for cancer related to obesity varied by four-fold in females and by almost 20-fold in males among counties and rates for overall cancer mortality varied by 2–3 fold. Geographic distributions of the urban sprawl index, overall cancer mortality for both sexes combined and the prevalence of obesity are illustrated in Figure 1A-C. Note that multiple urban and suburban counties are found in every state in the continental U.S.

Spatial autocorrelation

To determine if our simple models accounted for spatial autocorrelation in cancer mortality rates across these 935 urban and suburban counties, we calculated Moran’s I statistic for a series of models for both genders and the two cancer groups of interest (Table 2). Analysis of residuals from the simplest models (ln Sprawl only) for all four of our cancer mortality groupings and two of the more complex models for females showed highly significant unexplained spatial autocorrelation (Table 2). The addition of random effects to capture spatial autocorrelation of rates at the county level reduced all Moran’s I values to non-significant levels (Table 2, models 5 & 6). Furthermore, the spatial variance component parameters measuring autocorrelation across counties were significant in the most complex models for overall cancer among both men and women and obesity-related cancers among men. Together these results suggest that the simple fixed effects models did not adequately account for spatial patterns in the residuals. In subsequent models we included census divisions to measure broad regional effects and a random effect for spatial autocorrelation between counties. After inclusion of these variables in the models there is little spatial variation in mortality rates that is not explained by the specific demographic, risk, and health care related factors and location effects included in the more complex analyses. This allows us to perform a robust analysis of associations between urban sprawl and cancer mortality.

Table 2. Moran’s I test for spatial autocorrelation for alternative models

Model results

Preliminary analyses indicated that associations between mortality and the sprawl index varied significantly by Census Division. Therefore, we chose to include this interaction term (ln sprawl*Census Division) in the final random effects models presented here. Coefficients for the main effects of the nine covariates in addition to urban sprawl included in each model are given in Table 3 and coefficients for Census Division and ln sprawl*Census Division are presented in Table 4. Because the interactions between census division and natural log transformed sprawl index were significant, we do not report overall coefficients for the sprawl index.

Table 3. Random effects model results

Table 4. Model results for census divisions and sprawl, controlling for factors listed in Table3

We found significant interactions between the natural log (ln) transformed sprawl index and Census Division for all four groupings (Table 4), in other words, the slope of the association between this urban sprawl index and cancer mortality rate differed in different parts of the US. Predicted rates relating mortality from obesity related cancers and the urban sprawl index are given in Figure 2. For females, the interaction terms for census division and ln sprawl for obesity related cancer mortality were statistically significant for the East North Central States, South Atlantic States, and East and West South Central States. For males, the interaction term was significant in the Mid-Atlantic, East North Central, and East South Central Census Divisions and marginally significant in the Pacific Division. Six of the seven statistically significant interactions were positive, with a negative association in males in the Mid-Atlantic division. In other words, obesity related cancer mortality rates were generally higher in less sprawling counties (Figure 2A,B).

thumbnailFigure 2. Model based estimates of the association between natural log-transformed sprawl (‘compactness’) index and obesity-related cancer mortality for males (A) and females (B) from data for 934 urban and suburban counties in the United States. States included in the census divisions are available at the US Census [42]. The length of each regression line reflects the range of observed values in the Divisions (Table 1).

Current smoking was significantly and positively associated with mortality for all four groups examined (Table 4). The associations were strongest for total cancer mortality. Total cancer mortality for both genders was higher in counties with greater prevalence of non-Hispanic blacks, higher prevalence of obesity, and lower where prevalence of American Indians and Asian and Pacific Islanders was greater. A higher SEP was associated with total cancer mortality in females, but not in males. For obesity related cancers, physician and hospital density were negatively associated with cancer mortality in males but not in females. Prevalence of uninsured under age 64 was negatively correlated with all cancer mortality in both sexes. We discuss this counterintuitive result below.

Interactions included in our final models (Table 3) were all statistically significant (p ≤ 0.01). For males, a total of six interactions were included across the four cancer groupings, all including Census Division, physician density or both. For females, seven interactions were included; all seven of these included Census Division, SEP or both. Because our main interest here is in urban sprawl, we do not present detailed results for these interactions; nevertheless, they further highlight the presence of regional heterogeneity in cancer mortality rates at the Census Division level.

To assess goodness of fit, we report AIC and pseudo R2 values in Table 5. AIC values were lowest in our final models in all cases and R2 values were highest or nearly highest (Table 5). Thus, all measures of model fit, including Moran’s I statistic, support our choice of models with Census Division, its interaction with sprawl and spatial autocorrelation random effects included in the models.

Table 5. Model comparisons by cancer group

Discussion

This cross-sectional ecological analysis of urban sprawl and cancer mortality was designed to examine associations between cancer mortality and urban sprawl and to determine if such associations differed in different regions of the United States. The study had two main results. First, obesity related cancer mortality rates are either higher in counties with less sprawling urban form or unrelated to sprawl. This result extends a recent analysis by Fan and Song [12] that reported positive associations between total cancer mortality and sprawl in a subset of the suburban counties examined here and also reported higher mortality associated with cancer in urban compared to suburban counties. Considerable literature has also examined rural–urban differences in mortality, with mixed results [24,61]. Some studies have also reported risk of diagnosis at later stages of cancer is higher in urban compared to rural or suburban areas [19,26]. More specific measures of factors varying along urban-suburban and urban–rural gradients such as the sprawl (‘compactness’) index examined in this paper could help clarify how mortality and stage at diagnosis are associated with features of such environments.

Second, the association between urban sprawl and obesity related cancer mortality differs in different regions of the United States, with more positive associations in the southern part of the country and the Great Lakes area for females and the central part of the US for males (Table 3; Figure 2). Obesity related cancer mortality was higher in more sprawling counties in the Mid-Atlantic States for males. The least sprawling county in the US (New York County NY) has a sprawl index of ~480 (natural log = 6.17) whereas Montgomery County, MD, adjacent to Washington DC, a county with an intermediate level of sprawl, has a sprawl index of ~120 (natural log = 4.79). The most sprawling counties in the US, have indices of ~ 70 (natural log = 4.25). Thus, large differences in urban form are associated with modest differences in cancer mortality.

Nevertheless, these associations between sprawl and cancer mortality could have public health significance because of the large number of people involved. Greater than 70% of the US population lives in the urban and suburban counties analyzed here. The association between obesity related cancers and sprawl was not qualitatively altered by the inclusion of county level obesity prevalence in the model despite the presence of cross-sectional associations between obesity and urban sprawl found in several past studies [4,5,7]. This suggests that the elevated mortality rates observed in less sprawling areas in the south and Great Lakes regions could be caused by social, medical care related, or risk factors not adequately accounted for in our models rather than obesity prevalence.

Comparative studies of diverse countries concerning the health and environmental consequences of sprawl are lacking [25]. Recent studies have reported that active transportation is more common amongst Canadian 12–15 year olds residing in more sprawling areas [15] and in a large study of Portuguese adults population density was not associated with physical activity or BMI [62]. In both the US and in other countries, analyses of rural–urban differences in cancer mortality have had mixed results [21-24,61,63]. Careful comparative studies of urban sprawl and cancer incidence and sprawl across diverse regions and countries could provide further insight into associations of specific features of urban form and health outcomes.

Strengths of this analysis include that we make a rigorous attempt to account for spatial variation in cancer mortality and in the distribution of cancer risk factors at the county level across the United States. This effort involves successive elaboration of our regression models plus allowing for possible spatial autocorrelation among neighboring counties. We argue that this is important, because of the proliferation of ecological models that include geographic variables, but fail to test whether or not they account for spatial variation in the dependent variables of interest [12,38].

Our finding that the association between sprawl and cancer mortality varies by region suggests that future analyses of the health effects of sprawl should formally include analysis of regional differences. Most past studies of sprawl and energy balance have not focused on regional variation [3-5,9,64]. Recently, Troped et al. [36] reported regional difference in associations between perceived built environment and physical activity but associations between perceived crime and activity did not appear to differ by region [36]. The analyses in our paper adjust for the effects of demographic, economic and health care system related variables. Additional explanations for the regional heterogeneity observed in relations between sprawl, mortality and behavior could also be sought through consideration of regional differences in transportation systems, built environment aspects of sprawl, climate related factors, or inadequate control for residual confounding by risk factors included in our models. For example, counties classified as urban or suburban in the western United States included some very large counties such as Coconino, Yuma and Pima Counties in Arizona, San Bernardino County in California, and San Juan County in New Mexico. Additionally, suburban counties in some areas of the US are almost entirely residential, whereas in the Midwest and Western United States, suburban counties may have a large amount of agricultural or undeveloped land. Nevertheless, our results suggest that future geospatial analyses of cancer mortality and urbanicity across the US should attempt to account for regional variation in such associations and further explore sources of such variation.

Insurance coverage by county

One ecological puzzle emerged from this analysis. We found that the prevalence of lack of insurance coverage in subjects under age 65 was negatively associated with total cancer mortality rate in both males and females. In other words, cancer mortality rates were lower in counties where prevalence of insurance coverage was lower. Two recent analyses using SEER mortality data concerning urological cancers (Prostate, Bladder, and Kidney) also report negative associations between prevalence of uninsured and cancer mortality [38,65]. These results contrast with a large body of evidence indicating that insurance coverage decreases mortality rate at the individual level [53]. Incomplete adjustment for age does not seem to account for this observation; inclusion of median age in our regression models did not alter these findings, although more detailed analysis of the potential effects of age structure may be necessary since everyone over 65 in the US is covered by Medicare health insurance.

Heterogeneity in age, health and insurance coverage among Hispanics may also account for the observed association between lack of insurance coverage and lower cancer mortality at the county level despite the inclusion of prevalence of Hispanic race/ethnicity in our regression models. Lack of insurance coverage and Hispanic race ethnicity are correlated at the county level in this data set and the counties with high prevalence of uninsured and lower cancer mortality are among the counties with the highest proportion Hispanic. The Hispanic population also has a younger age distribution that non-Hispanic Whites and Blacks. However, the interactions between prevalence of insurance and Hispanic race/ethnicity were not significant, thus the associations observed here occurred independent of each other. Furthermore, these counties span the whole range of socioeconomic conditions as measured by the SEP index. More attention to this ecological puzzle may be warranted in light of the changes occurring in US health care policy.

Study weaknesses

A significant weakness of this study is the reliance on cross sectional data. A cross sectional analysis, such as this one, cannot be used to draw causal inferences concerning dependent variables. Nevertheless, we think that our analysis is useful in that it draws attention to regional heterogeneity in the association between sprawl and cancer mortality and supports a renewed focus on the higher cancer mortality rates observed in urban areas. Additionally, urban form changes slowly, therefore cross sectional studies do have the unique ability to allow examination of large differences in urban form.

Our analysis lacks information concerning two important cancer related exposures. First, environmental exposures to pollutants; such data are difficult to obtain in a uniform format for this complete sample of urban and suburban US counties. Exposure data might also be most relevant to analyses of specific cancers and we hope the analytic framework developed here will foster interest in pursuing such analyses. Second, data concerning physical activity and sedentary time. Sedentary behavior and physical activity are also associated with cancer incidence and mortality at several sites [32]. Physical inactivity data based on self-report are available at the county level [66]. However, we chose not to include this potential covariate in part because it is highly correlated with obesity prevalence, and in part because of the limits of self reported PA data [67]. Self report of obesity and smoking behavior are also subject to measurement error, but these errors are better understood than measurement error in physical activity [68,69]. Lack of information about physical activity and environmental pollutants could result in omitted variable bias if they are part of a causal connection between urban sprawl and cancer mortality. Such omitted variable bias or unmeasured confounders can influence regression results [70].

A potential weakness, but a difficult one to evaluate, involves the fact that our county level estimates of BMI and smoking are not actual measurements but the results of statistical models that include county demographic information also included in our models. This approach has the potential to increase the correlations between some of the explanatory variables in our models. However, because the covariates used to improve estimates of county level obesity and smoking prevalence are not strongly correlated with the other demographic variables used in our models, this problem may be a small one. Nevertheless, our overall analysis would be more robust if we could use estimates from data collected with sufficient accuracy in the field or, if that is not possible, values that are not contingent on these demographic variables [71,72].

Additional challenges involve the level of spatial aggregation examined here and our use of a random effects regression framework as opposed to a hierarchical linear modeling approach [73] with levels of variance estimated for several nested geographic areas. It is certainly not obvious that county is the appropriate level of spatial aggregation in which to examine cancer mortality rates in relation to urban sprawl and the various covariates considered in this analysis. Multilevel analysis that accounted for individual as well as regional characteristics could also be informative, for example by helping to account for the unexpected direction of the association between insurance coverage and mortality observed in this analysis. Such analysis might also shed further light on the observation that mortality rates for cancer were higher in more urban areas. Analyses at different levels of geographic aggregation and analyses incorporating individual, as well as aggregate, variables are possible. However, at present, US cancer mortality statistics are publicly available only at the county level and death certificates do not include information on diverse cancer related covariates at the individual level [43]. Furthermore, access to the geocoded mortality data needed for such analyses requires negotiation with each SEER registry site and with the CDC. Thus significant barriers remain to more complete and multi-level analysis of mortality data.

Conclusions

Urban sprawl appears to have diverse negative health and environmental consequences. However, this study demonstrates that the cross-sectional associations between sprawl and mortality differ in different regions of the Unites States. Our results and those of other recent papers on this topic [12,36] should encourage a more nuanced view of the association between sprawl and health as well as future studies that explore regional differences in the effects of urban sprawl. More longitudinal studies of sprawl and health [9,13] as well as multilevel analyses, regardless of study design [11,37,74] will continue to clarify the potential for causal associations between urban form and health. There are many existing cohort studies that would allow such analyses.

Renewed attention to the causes of elevated cancer mortality rates in urban areas could help address disparities in cancer mortality. Preliminary results suggest that there are also regional differences in the association between urban sprawl and tobacco related cancer mortality. Thus, future studies could examine tobacco-related (oral cavity and pharynx, esophagus, larynx, lung, bladder) cancers using a similar approach to the one outlined here. Analyses of associations between sprawl and cancer incidence are also a logical next step in an effort to identify specific factors increasing cancer mortality in less sprawling urban and suburban counties both in the US and potentially in other countries where sprawl has emerged as a common pattern of development.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DB, ZT, LP, RE, and RBB conceived of the project. RE developed the sprawl index. DB, ZT, and LP carried out the statistical analysis. DB wrote the first draft of the paper. All authors reviewed and approved the final version of the paper.

Acknowledgements

We thank David Stinchcomb, Eric Augustson, Susan Devesa and Anne Hartman for helpful discussions of this topic, Patty Griffin for excellent programming assistance, and Penny Randall-Levy for bibliographic formatting.

References

  1. Glaeser EL: Cities, productivity, and quality of life.

    Science 2011, 333:592-594. PubMed Abstract | Publisher Full Text OpenURL

  2. Duany A, Plater-Zyberk E, Speck J: Suburban Nation: the Rise of Sprawl and the Decline of the American Dream. New York: North Point Press; 2000. OpenURL

  3. Sturm R, Cohen DA: Suburban sprawl and physical and mental health.

    Public Health 2004, 118:488-496. PubMed Abstract | Publisher Full Text OpenURL

  4. Ewing R, Schmid T, Killingsworth R, Zlot A, Raudenbush S: Relationship between urban sprawl and physical activity, obesity, and morbidity.

    Am J Health Promot 2003, 18:47-57. PubMed Abstract | Publisher Full Text OpenURL

  5. Lopez R: Urban sprawl and risk for being overweight or obese.

    Am J Public Health 2004, 94:1574-1579. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Vandegrift D, Yoked T: Obesity rates, income, and suburban sprawl: an analysis of US states.

    Health Place 2004, 10:221-229. PubMed Abstract | Publisher Full Text OpenURL

  7. Papas MA, Alberg AJ, Ewing R, Helzlsouer KJ, Gary TL, Klassen AC: The built environment and obesity.

    Epidemiol Rev 2007, 29:129-143. PubMed Abstract | Publisher Full Text OpenURL

  8. Stewart CJ: Snake bite in Australia: first aid and envenomation management.

    Accid Emerg Nurs 2003, 11:106-111. PubMed Abstract | Publisher Full Text OpenURL

  9. Lee IM, Ewing R, Sesso HD: The built environment and physical activity levels: the Harvard Alumni Health Study.

    Am J Prev Med 2009, 37:293-298. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Garden FL, Jalaludin BB: Impact of urban sprawl on overweight, obesity, and physical activity in Sydney, Australia.

    J Urban Health 2009, 86:19-30. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. James P, Troped PJ, Hart JE, Joshu CE, Colditz GA, Brownson RC, Ewing R, Laden F: Urban sprawl, physical activity, and body mass index: nurses’ health study and nurses’ health study II.

    Am J Public Health 2013, 103:369-375. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Fan Y, Song Y: Is sprawl associated with a widening urban-suburban mortality gap?

    J Urban Health 2009, 86:708-728. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Ewing R, Brownson RC, Berrigan D: Relationship between urban sprawl and weight of United States youth.

    Am J Prev Med 2006, 31:464-474. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Wheeler SM: The evolution of built landscapes in metropolitan regions.

    J Plan Educ Res 2008, 27:400-416. Publisher Full Text OpenURL

  15. Seliske L, Pickett W, Janssen I: Urban sprawl and its relationship with active transportation, physical activity and obesity in Canadian youth.

    Health Rep 2012, 23:17-25. OpenURL

  16. Nayha S, Lankila T, Rautio A, Koiranen M, Tammelin TH, Taanila A, Rusanen J, Laitinen J: Body mass index and overweight in relation to residence distance and population density: experience from the Northern Finland birth cohort 1966.

    BMC Public Health 2013, 13:938. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  17. Vlahov D, Galea S: Urbanization, urbanicity, and health.

    J Urban Health 2002, 79:S1-S12. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Freudenberg N, Klitzman S, Saegert S: Urban Health and Society: Interdisciplinary Approaches to Research and Practice. San Francisco: John Wiley and Sons; 2009. OpenURL

  19. McLafferty S, Wang F: Rural reversal? Rural–urban disparities in late-stage cancer risk in Illinois.

    Cancer 2009, 115:2755-2764. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Waldron HA: A brief history of scrotal cancer.

    Br J Ind Med 1983, 40:390-401. PubMed Abstract | PubMed Central Full Text OpenURL

  21. Stocks P: On the relations between atmospheric pollution in urban and rural localities and mortality from cancer, bronchitis and pneumonia, with particular reference to 3:4 benzopyrene, beryllium, molybdenum, vanadium and arsenic.

    Br J Cancer 1960, 14:397-418. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Pearce J, Boyle P: Is the urban excess in lung cancer in Scotland explained by patterns of smoking?

    Soc Sci Med 2005, 60:2833-2843. PubMed Abstract | Publisher Full Text OpenURL

  23. Robert SA, Strombom I, Trentham-Dietz A, Hampton JM, McElroy JA, Newcomb PA, Remington PL: Socioeconomic risk factors for breast cancer: distinguishing individual- and community-level effects.

    Epidemiology 2004, 15:442-450. PubMed Abstract | Publisher Full Text OpenURL

  24. Pozet A, Westeel V, Berion P, Danzon A, Debieuvre D, Breton JL, Monnier A, Lahourcade J, Dalphin JC, Mercier M: Rurality and survival differences in lung cancer: a large population-based multivariate analysis.

    Lung Cancer 2008, 59:291-300. PubMed Abstract | Publisher Full Text OpenURL

  25. Wilson B, Chakraborty A: The environmental impacts of sprawl: emergent themes from the past decade of planning research.

    Sustainability 2013, 5:3302-3327. Publisher Full Text OpenURL

  26. Paquette I, Finlayson SR: Rural versus urban colorectal and lung cancer patients: differences in stage at presentation.

    J Am Coll Surg 2007, 205:636-641. PubMed Abstract | Publisher Full Text OpenURL

  27. Kuo TM, Mobley LR, Anselin L: Geographic disparities in late-stage breast cancer diagnosis in California.

    Health Place 2011, 17:327-334. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Mobley LR, Kuo TM, Watson L, Gordon BG: Geographic disparities in late-stage cancer diagnosis: multilevel factors and spatial interactions.

    Health Place 2012, 18:978-990. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Morgenstern H: Ecologic studies in epidemiology: concepts, principles, and methods.

    Annu Rev Public Health 1995, 16:61-81. PubMed Abstract | Publisher Full Text OpenURL

  30. Oakes JM: The (mis) estimation of neighborhood effects: causal inference for a practicable social epidemiology.

    Soc Sci Med 2004, 58:1929-1952. PubMed Abstract | Publisher Full Text OpenURL

  31. Gerring J: Case Study Research: Principles and Practices. New York: Cambridge University Press; 2007. OpenURL

  32. World Cancer Research Fund: American Institute for Cancer Research: Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective. Washington, DC: American Institute for Cancer Research; 2007. OpenURL

  33. Meersman SC, Breen N, Pickle LW, Meissner HI, Simon P: Access to mammography screening in a large urban population: a multi-level analysis.

    Cancer Causes Control 2009, 20:1469-1482. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Lovasi GS, Hutson MA, Guerra M, Neckerman KM: Built environments and obesity in disadvantaged populations.

    Epidemiol Rev 2009, 31:7-20. PubMed Abstract | Publisher Full Text OpenURL

  35. Feng J, Glass TA, Curriero FC, Stewart WF, Schwartz BS: The built environment and obesity: a systematic review of the epidemiologic evidence.

    Health Place 2010, 16:175-190. PubMed Abstract | Publisher Full Text OpenURL

  36. Troped PJ, Tamura K, Whitcomb HA, Laden F: Perceived built environment and physical activity in U.S. women by sprawl and region.

    Am J Prev Med 2011, 41:473-479. PubMed Abstract | Publisher Full Text OpenURL

  37. Griffin BA, Eibner C, Bird CE, Jewell A, Margolis K, Shih R, Ellen SM, Whitsel EA, Allison M, Escarce JJ: The relationship between urban sprawl and coronary heart disease in women.

    Health Place 2013, 20:51-61. PubMed Abstract | Publisher Full Text OpenURL

  38. Colli J, Sartor O, Thomas R, Lee BR: Does urological cancer mortality increase with low population density of physicians?

    J Urol 2011, 186:2342-2346. PubMed Abstract | Publisher Full Text OpenURL

  39. Waller LA, Gotway CA: Applied Spatial Statistics for Public Health Data. Hoboken, NJ: John Wiley & Sons; 2004. OpenURL

  40. Hankey BF, Ries LA, Edwards BK: The surveillance, epidemiology, and end results program: a national resource.

    Cancer Epidemiol Biomarkers Prev 1999, 8:1117-1121. PubMed Abstract | Publisher Full Text OpenURL

  41. Centers for Disease Control and Prevention: National Program of Cancer Registries (NPCR). http://www.cdc.gov/cancer/npcr/ webcite

  42. United States Census Bureau: Geographic terms and concepts--Census divisions and census regions. http://www.census.gov/geo/reference/gtc/gtc_census_divreg.html webcite

  43. Howlader N, Noone AM, Krapcho M, Neyman N, Aminou R, Waldron R, Altekruse SF, Kosary CL, Ruhl J, Tatalovich Z, Cho H, Mariotto A, Eisner MP, Lewis DR, Chen HS, Feurer EJ, Cronin KA, Edwards BK: SEER cancer statistics review, 1975–2008. http://seer.cancer.gov/csr/1975_2008/ webcite

  44. Centers for Disease Control and Prevention: Promoting data use for program planning: Mississippi Cancer Registry. http://www.cdc.gov/cancer/npcr/success/south/ms/ms1.htm webcite

  45. Ewing R: Can the physical environment determine physical activity levels?

    Exerc Sport Sci Rev 2005, 33:69-75. PubMed Abstract | Publisher Full Text OpenURL

  46. Tsai Y: Quantifying urban form: compactness versus ‘sprawl’.

    Urban Stud 2005, 42:141-161. Publisher Full Text OpenURL

  47. Gregg EW, Kirtland KA, Cadwell BL, Rios Burrows N, Barker LE, Thompson TJ, Geiss L, Pan L: Estimated county-level prevalence of diabetes and obesity - United States, 2007.

    MMWR Morb Mortal Wkly Rep 2009, 58:1259-1263. PubMed Abstract | Publisher Full Text OpenURL

  48. Centers for Disease Control and Prevention: Diabetes data & trends. County level estimates of obesity--U.S. maps: 2009 age-adjusted estimates of the percentage of adults who are obese. http://apps.nccd.cdc.gov/DDT_STRS2/NationalDiabetesPrevalenceEstimates.aspx?mode=OBS webcite

  49. National Cancer Institute: Small area estimates for cancer risk factors & screening behaviors. Current smoking prevalence (age 18+): maps and data of model-based small area estimates. http://sae.cancer.gov/estimates/current.html webcite

  50. Raghunathan TE, Xie D, Schenker N, Parsons VL, Davis WW, Dodd KW, Feuer EJ: Combining information from two surveys to estimate county-level prevalence rates of cancer risk factors and screening.

    J Am Stat Assoc 2007, 102:474-486. Publisher Full Text OpenURL

  51. Krieger N, Chen JT, Waterman PD, Soobader M-J, Subramanian SV, Carson R: Geocoding and monitoring of US socioeconomic inequalities in mortality and cancer incidence: does the choice of area-based measure and geographic level matter?: the Public Health Disparities Geocoding Project.

    Am J Epidemiol 2002, 156:471-482. PubMed Abstract | Publisher Full Text OpenURL

  52. U.S. Department of Health and Human Services, Health Resources and Services Administration, Bureau of Health Professions, Office of Workforce Policy and Performance Managenment:

    User Documentation for the Area Resource File (ARF), 2008 Release. 2009. OpenURL

  53. Marlow NM, Pavluck AL, Bian J, Ward EM, Halpern MT: The Relationship Between Insurance Coverage and Cancer Care: a Literature Synthesis. Research Triangle Park, NC: RTI Press Publication; 2009.

    No. RR-0005-0905

    OpenURL

  54. SAS Institute Inc: JMP®, Version 8.0. Cary, NC: SAS Institute, Inc; 1989–2007. OpenURL

  55. R Development Core Team: R: a Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2008. OpenURL

  56. Friedman J, Hastie T, Tibshirani R: Regularization paths for generalized linear models via coordinate descent.

    J Stat Softw 2010, 33:1-22. PubMed Abstract | PubMed Central Full Text OpenURL

  57. Zou H, Hastie T: Regularization and variable selection via the elastic net.

    J R Stat Soc Series B Stat Methodol 2005, 67:301-320. Publisher Full Text OpenURL

  58. SAS Institute Inc: SAS/STAT software: Documentation of GLIMMIX procedure in SAS/STAT 9.1. http://www.sas.com/statistics/doc.html webcite

  59. Brownie C, Gumpertz ML: Validity of spatial analyses for large field trials.

    J Agric Biol Environ Stat 1997, 2:1-23. Publisher Full Text OpenURL

  60. Akaike H: A new look at the statistical model identification.

    IEEE Trans Automat Contr 1974, 19:716-723. Publisher Full Text OpenURL

  61. Obertova Z, Brown C, Holmes M, Lawrenson R: Prostate cancer incidence and mortality in rural men–a systematic review of the literature.

    Rural Remote Health 2039, 2012:12. OpenURL

  62. Santana P, Santos R, Nogueira H: The link between local environment and obesity: a multilevel analysis in the Lisbon Metropolitan Area, Portugal.

    Soc Sci Med 2009, 68:601-609. PubMed Abstract | Publisher Full Text OpenURL

  63. Monroe AC, Ricketts TC, Savitz LA: Cancer in rural versus urban populations: a review.

    J Rural Health 1992, 8:212-220. PubMed Abstract | Publisher Full Text OpenURL

  64. Zhao Z, Kaestner R: Effects of urban sprawl on obesity.

    J Health Econ 2010, 29:779-787. PubMed Abstract | Publisher Full Text OpenURL

  65. Odisho AY, Cooperberg MR, Fradet V, Ahmad AE, Carroll PR: Urologist density and county-level urologic cancer mortality.

    J Clin Oncol 2010, 28:2499-2504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  66. Centers for Disease Control and Prevention: Physical inactivity estimates, by county. Staying active during leisure time can reduce your risk of obesity, type 2 diabetes and related diseases. http://www.cdc.gov/Features/dsPhysicalInactivity/ webcite

  67. Troiano RP, Berrigan D, Dodd KW, Masse LC, Tilert T, McDowell M: Physical activity in the United States measured by accelerometer.

    Med Sci Sports Exerc 2008, 40:181-188. PubMed Abstract | Publisher Full Text OpenURL

  68. Connor Gorber S, Tremblay M, Moher D, Gorber B: A comparison of direct vs. self-report measures for assessing height, weight and body mass index: a systematic review.

    Obes Rev 2007, 8:307-326. PubMed Abstract | Publisher Full Text OpenURL

  69. Patrick DL, Cheadle A, Thompson DC, Diehr P, Koepsell T, Kinne S: The validity of self-reported smoking: a review and meta-analysis.

    Am J Public Health 1994, 84:1086-1093. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  70. Fewell Z, Davey SG, Sterne JA: The impact of residual and unmeasured confounding in epidemiologic studies: a simulation study.

    Am J Epidemiol 2007, 166:646-655. PubMed Abstract | Publisher Full Text OpenURL

  71. Malec D, Sedransk J, Moriarity CL, LeClere FB: Small area inference for binary variables in the national health interview survey.

    J Am Stat Assoc 1997, 92:815-826. Publisher Full Text OpenURL

  72. Knorr-Held L, Best NG: A shared component model for joint and selective clustering of two diseases.

    J R Stat Soc Ser A 2001, 164:73-85. Publisher Full Text OpenURL

  73. Raudenbush SW, Byrk AS: Hierarchical linear models: applications and data analysis methods. 2nd edition. Newbury Park, CA: Sage Publications; 2002. OpenURL

  74. Kelly-Schwartz AC, Stockard J, Doyle S, Schlossberg M: Is sprawl unhealthy? A multilevel analysis of the relationship of metropolitan sprawl to the health of individuals.

    J Plan Educ Res 2004, 24:184-196. Publisher Full Text OpenURL