Next Article in Journal
On the Determinants of Bitcoin Returns and Volatility: What We Get from Gets?
Next Article in Special Issue
Transit Travel Community Detection and Evolutionary Analysis: A Case Study of Shenzhen
Previous Article in Journal
One-Dimensional Computational Model of Gyttja Clay for Settlement Prediction
Previous Article in Special Issue
Toward a Healthy Urban Living Environment: Assessing 15-Minute Green-Blue Space Accessibility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Non-Stationarity Effects of Unhealthy Food Environments and Green Spaces for Type-2 Diabetes in Toronto

1
Department of Geography, Geomatics and Environment, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, ON L5L 1C6, Canada
2
Department of Geography and Planning, University of Toronto, 100 St. George St., Toronto, ON M5S 3G3, Canada
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(3), 1762; https://doi.org/10.3390/su15031762
Submission received: 2 December 2022 / Revised: 6 January 2023 / Accepted: 16 January 2023 / Published: 17 January 2023

Abstract

:
Environmental factors may operate differently when relations are measured across different geographical locations, a phenomenon known as spatial non-stationarity. This study investigates the spatial non-stationarity effect of unhealthy food environments and green spaces on the T2DM prevalence rate at the neighborhood level in Toronto. This study also compares how the results vary between age groups, classified as all adults (20 and above), young adults (from 20 to 44), middle adulthood (from 45 to 64), and seniors (65 and above). The geographically weighted regression model is utilized to explore the impacts of spatial non-stationarity effects on the research results, which may lead to biased conclusions, which have often been ignored in past studies. The results from this study reveal that environmental variables dissimilarly affect T2DM prevalence rates among different age groups and neighborhoods in Toronto after controlling for socioeconomic factors. For example, the green space density yields positive associations with diabetes prevalence rates for elder generations but negative relationships for younger age groups in twenty-two and four neighborhoods, respectively, around Toronto East. The observed associations will provide beneficial suggestions to support government and public health authorities in designing education, prevention, and intervention programs targeting different neighborhoods to control the burden of diabetes.

1. Introduction

The rising number of diabetes cases worldwide is emerging as one of the fastest-growing global healthcare emergencies of the 21st century [1,2]. More than half a billion people are living with diabetes, which is predicted to increase by 46% by 2045 [1]. Type-2 diabetes mellitus (T2DM) accounts for more than 90% of diabetes cases [1,3]. It is characterized by excess blood sugar levels caused by insulin deficiency, whereby the pancreas cannot produce enough insulin, and insulin resistance, whereby target organs respond poorly to the insulin that is produced and take in less glucose from the blood [1,3]. Diabetes can cause severe damage to body systems such as the eyes, kidneys, heart, and vascular system [4]. Although T2DM is recognized to be more common in older adults, there has been a recent trend of increasing rates among young adults due to obesity, physical inactivity, and energy-dense diets [5]. Using Toronto, Ontario, Canada, as an example, more than one in ten adults aged 20 years and above have been diagnosed with T2DM, and the prevalence rate has doubled in the past two decades [6]. As a result, governments and public health organizations are working to limit and reduce the prevalence of T2DM by understanding the environmental factors that contribute to its development in the hope of designing effective intervention programs.
Public health and medical researchers have been studying the effects of various environmental factors, such as air pollution, food, greenness, noise, physical activity, and walkability, on diabetes for decades [4]. These factors can affect the risk of T2DM by altering behavioral, psychological, and physical stresses and influencing decision-making [4]. Numerous studies have been performed to further advanced our knowledge of how physical environmental factors [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24] and social environmental factors [25,26,27,28,29] influence T2DM. It has been well established that physical inactivity and unhealthy food environments, which promote the consumption of energy-dense foods, are significant predictors of T2DM [5]. People are likely to exercise more in a well-designed environment. Recreational resources, green spaces, and walkways that encourage physical activity have been shown to reduce the risk of diabetes [7,8,14,15,16,19]. Having supermarkets with fresh food in the community can promote a healthy diet, while unhealthy food outlets can increase the risk of diabetes [11,12,22]. Age and social–economic status also play a critical role in T2DM prevalence. Evidence suggests that T2DM prevalence rates are positively associated with age and unemployment [1,26,30], while education and income are negatively associated with T2DM prevalence rates [26,30,31].
With the increased use of geographic information science (GIS) and spatial analytical techniques, there have been growing concerns regarding how the environment affects human health worldwide [32]. Many studies have discovered associations between environmental factors and T2DM. For instance, some studies have found significant positive associations between the food environment and T2DM [12,13,22,23], while others have found no significant associations [10,25,28]. These examples highlight the uncertain geographic context problem (UGCoP) [33] and the modifiable area unit problem (MAUP) [34], which refer to the fact that associations found at different geographical scales (e.g., neighborhood, city, province, and country levels) may differ. The results may also be impacted by the various methods used to determine the contextual location when measuring exposure effects, such as using road networks or buffers to measure neighborhood health outcomes. However, environmental factors can affect different geographic areas dissimilarly when relationships are measured across space, a phenomenon known as spatial non-stationarity [35,36]. For example, an unhealthy food environment may positively correlate with the T2DM prevalence rate in one city but negatively correlate in another. Nevertheless, research often assumes that the relationships between environmental factors and diabetes effects are stable or stationary over space. Such presumptions may lead to incorrect conclusions about how the environment affects health outcomes [35]. While very few studies have focused on the spatial non-stationarity effect of the environment on T2DM, most have only examined this effect at the county level [37,38,39]. To the best of our knowledge, little research has been performed on the spatial non-stationarity effects between the environment and T2DM at the citywide neighborhood level; thus, this is worthy of investigation. Understanding the geographic variation in T2DM rates at the city level can help when planning and determining where management and prevention resources should be allocated [37].
The goal of this paper is to address the research gap by (1) examining the spatial non-stationarity associations between T2DM and contextual variables (focusing on green space densities and unhealthy food environments) and (2) determining to what extent these associations differ between age groups at the neighborhood level in Toronto, Canada. It is worth noting that spatial autocorrelation may exist for the prevalence of T2DM. Spatial autocorrelation measures the extent to which geographic features and associated values are clustered together or dispersed in space. For example, a higher incidence of diabetes in the city center may indicate that some processes, such as clusters of fast food restaurants, contribute to the increased prevalence in that location. Therefore, this study will use global and local Moran’s indexes to assess the existence of spatial autocorrelation and explore the spatial non-stationarity effect of the environment on T2DM through geographically weighted regression (GWR). The GWR model addresses the issue of spatial autocorrelation and reduces inaccurate spatial data estimations due to non-randomness errors. This research further intends to advance our understanding of how these relationships vary based on Toronto’s geographical locations and age groups. The findings from this study could provide valuable evidence and suggestions to support the government and public health authorities in designing customized education, prevention, and intervention policies. The intervention policies will target different neighborhoods, age groups, and minority populations to control and reduce the increasing rate of diabetes in Toronto.

2. Materials and Methods

This study investigated the spatial non-stationarity relationships of T2DM prevalence rates with both green space and unhealthy food environments in the 158 neighborhoods in Toronto. The neighborhood-level geographic unit was used as a contextual area for all analyses. Four geographical weighted regression (GWR) models were employed to investigate the spatial non-stationarity associations between environmental factors and T2DM prevalence rates in different age groups. The four dependent variables are the T2DM prevalence rates for individuals aged (1) 20 and above (all adults), (2) between 20 and 44 (young adults), (3) from 45 to 64 (middle adulthood), and (4) 65 and above (elder generations). The independent variables included the green space density, unhealthy food outlet density, and social–economic status (as control variables), which were kept the same for all GWR models. Figure 1 illustrates the workflow of this study using the GWR models.

2.1. Study Area

Toronto, the capital city of Ontario, is located in southern Ontario and has a land area of 630.2 km2. It is the most populated municipality in Canada, with an estimated three million culturally diverse inhabitants living at a density of 4720 persons per km2 in 2019 [40]. According to the city of Toronto and the 2016 Canadian census, 47% of the population is made up of immigrants [41]. The city has 158 officially registered neighborhoods, which serve as a microcosm of the population. Unlike census tracts or dissemination areas, these social planning neighborhood boundaries change infrequently over time, making them suitable for statistical reporting [42]. Given the diverse environments and demographic characteristics of Toronto’s neighborhoods, it is likely that spatial non-stationarity effects exist between environmental factors and T2DM. This research, therefore, conducts a large population-based study, covering all of Toronto’s neighborhoods, to investigate whether environmental exposure is associated with the prevalence of T2DM. The 158 neighborhoods in Toronto have an average area of 4.07 km2, with a median of 3.14 km2 and a standard deviation of 3.87 km2. The smallest neighborhood, North Toronto, has an area of 0.40 km2, while the largest neighborhood, West Humber-Clairville, covers 30.16 km2.

2.2. Data

2.2.1. Type-2 Diabetes Mellitus Prevalence Rates

The type-2 diabetes mellitus prevalence data for Toronto neighborhoods used in this study were obtained from the Ontario Community Health Profiles Partnership (OCHPP) under the Adult Health and Disease section [43]. The diabetes data were sourced from validated disease databases maintained by the Institute for Clinical Evaluative Sciences (ICES) [44]. The database includes records of all citizens and permanent residents aged 20 and above who interact with Canada’s universal healthcare system, regardless of their social status. The population demographics for Toronto’s 158 neighborhoods, classified by different age groups, are summarized in Table 1. These demographics were summarized according to the Ontario Ministry of Health and Long-Term Care Registered Persons Database (RPDB), including individuals who were alive and living in Toronto on 1 April 2019 [43]. If an individual filed an Ontario Health Insurance Plan (OHIP) claim, had two doctor claims, or was admitted to the hospital for diabetes within two years, they were recorded in the database as being diagnosed with diabetes. This study analyzed the most recently published 2018/2019 Adult Health and Disease dataset, focusing on diabetes prevalence rates for individuals aged 20 and above at the neighborhood level in Toronto.
The data were divided into age groups: 20 and above, 20 to 44, 45 to 64, and 65 and above. Each data record was matched and spatially joined with the boundary shapefile for the 158 neighborhoods in the City of Toronto to carry out spatial analyses. Figure 2 illustrates the spatial distribution of the T2DM prevalence rates across Toronto neighborhoods. Higher rates of T2DM were observed in the eastern and northwestern areas of Toronto, while lower rates were found in the central downtown region. For instance, neighborhoods around Humbermede and Scarborough had higher T2DM prevalence rates, while those in the Young-Bay corridor had lower rates.

2.2.2. Unhealthy Food Outlets

Recent studies have shown that unhealthy food environments can impact T2DM [13,22,23]. For instance, fast food restaurants and convenience stores often provide energy-dense and high-carbohydrate food sources that can quickly raise blood sugar levels. In order to analyze the effects of unhealthy food environments on T2DM in Toronto, this study used data on the locations of food outlets in 2021 provided by SafeGraph [45]. SafeGraph is a company that offers points of interest (POI) data containing physical location information at the address level [46,47,48]. The data were classified using the North American Industry Classification System (NAICS). Food outlets were obtained from SafeGraph Core Places data and filtered using the NAICS code for unhealthy food outlets. Previous studies have commonly classified convenience stores, confectionery stores, and limited-outlet restaurants (fast food) as “unhealthy” retail food outlets [49,50,51]. Therefore, this study contextualized unhealthy food outlets using the following NAICS codes: 722513—Limited-Service Restaurants; 722515—Snack and Non-Alcoholic Beverage Bars; 445292—Confectionery Stores; 447110—Gasoline Stations with Convenience Stores; 445120—Convenience Stores; 311811—Retail Bakeries. The locations of unhealthy food outlets are shown in Figure 3. A total of 3964 of these outlets were located around Toronto.

2.2.3. Green Spaces

Studies have found that neighborhood green spaces, such as parks, can reduce the risk of T2DM by promoting physical activity [4,7,8]. Exposure to green spaces and recreational facilities around residential neighborhoods can promote physical activity and benefit individuals’ health [52,53,54]. When there is more green space in a community, residents will have easier access to these spaces for physical activity. The green space data were obtained from the Park Sports Field Region–2019 dataset published by DMTI Spatial, Inc., on Scholars GeoPortal [55]. The dataset includes polygon shapefiles of all parks and recreation areas in Canada, such as sports fields, open spaces, campgrounds, and golf courses. Toronto has a total green space area of 131.9 km2, covering more than 20% of the city, as shown in Figure 3.

2.2.4. Socioeconomic Status

In addition to unhealthy food outlets and green spaces, socioeconomic status also significantly affects T2DM. Previous studies have shown that older adults and the unemployed have a higher risk of being diagnosed with T2DM [1,26,30]. Education and income are often negatively correlated with T2DM [26,30,31]. Regarding ethnicity, recent studies have found that the non-Hispanic white population has a lower prevalence rate of T2DM compared to other ethnic groups [56,57]. Certain immigrant groups were also found to have a higher risk of developing T2DM earlier in their life in Canada [26,58]. For example, South Asian and black immigrants are more likely to develop T2DM earlier than immigrants from the UK. Therefore, this study used socioeconomic characteristics to represent social environmental factors that impact T2DM. The socioeconomic factors used in this analysis included income, unemployment rate, low-education rate, and immigration rate. These variables were extracted from the 2016 Canadian Dissemination Area (DA) Level Census published by Statistics Canada [40]. The low-education rate was calculated by dividing the number of individuals without education certificates by the total population in each DA. The immigration rate was calculated similarly, by dividing the number of individuals identified as immigrants by the total population in each DA. The variables were then aggregated to the neighborhood level using ArcGIS Pro’s (Version 3.0) Spatial Join function by calculating the average value of each variable for each neighborhood, as the boundaries of DA match the neighborhoods [59]. A summary of the variables used in this study is listed in Table 2.

2.2.5. Derivation of Environmental Variables

The problem of an uncertain geographic context (UGCoP) arises from the use of arbitrary areal units while calculating area-based variables. Researchers do not have perfect knowledge of how different spatial and temporal configurations can affect health outcomes [33]. This study quantifies the effects of the environment on T2DM using a small spatial scale. However, it is not accurate to assume that residents in the arbitrary areal units (e.g., neighborhoods in this study) only have access to outlets and locations within their neighborhoods. People living near the borders of these units may have closer access to facilities in adjacent neighborhoods, which can influence the analysis results. This phenomenon is known as the “edge effect” [60]. This study used a container-based (buffer) measurement to address this issue. The analysis unit was given a buffered distance as a boundary to include locations outside each neighborhood as measurements of accessibility. However, there is a lack of agreement on the appropriate buffer threshold, as it depends on the areal units used in different studies [61]. The buffer distances in previous studies ranged from 500 m to 5 km. One study found that most buffers used to measure the density and proximity of retail food outlets varied between 2 and 3 km [62]. Therefore, this study created a buffer zone around the border of each neighborhood to include environmental variables (unhealthy food outlets and green space) up to 2.5 km away in the analysis, as shown in Figure 3.
This study calculates the density of unhealthy food outlets and green spaces within each neighborhood and their surrounding areas (buffer zones) to standardize the measurements of accessibility and environmental exposure. The density of unhealthy food outlets was calculated by dividing the number of locations within a 2.5 km buffer of each neighborhood by the area of the neighborhood with the container buffer. The density of green spaces was calculated in the same way, using the total area of the green spaces within a 2.5 km buffer of each neighborhood, and divided by the area of each neighborhood with the container buffer. A higher density indicates that citizens will have greater access to their surrounding environment. The variables included in the regression models are shown in Table 2.

2.3. Statistical Analysis

2.3.1. Spatial Autocorrelation

Spatial autocorrelation in the prevalence of T2DM in Toronto neighborhoods was examined using global Moran’s index statistics in ArcGIS Pro [59]. Using a default setting with an inverse Euclidian distance, the calculation yielded a Moran’s I of 0.7995, with a z-score of 18.05 (p < 0.001). This indicated that a strong positive spatial autocorrelation exists in Toronto neighborhoods, meaning that similar T2DM prevalence rates tend to be clustered spatially. This can also be seen in Figure 2. Anselin’s local Moran’s I was applied at the neighborhood level to further identify local clusters or spatial outliers in the T2DM prevalence rates. Figure 4 illustrates the results of this analysis, using ArcGIS Pro’s Cluster and Outlier Analysis program with a default setting of an inverse Euclidean distance. Statistically significant (at the 95% confidence level) clusters of high T2DM prevalence rates (HH) were found in northeastern and eastern Toronto. In contrast, clusters of low (LL) prevalence rates were located in central Toronto. The Regent Park neighborhood was an outlier, as it had a high T2DM prevalence rate but was surrounded by neighborhoods with low rates (HL).
Both the global and local Moran’s I analyses indicated that spatial autocorrelation exists for the T2DM prevalence rates in Toronto neighborhoods. However, many studies that examine the associations between T2DM prevalence and environmental factors frequently ignore how spatial autocorrelation may impact the research results. To address this, this study used a spatial regression model called geographically weighted regression to account for the effects brought on by spatial autocorrelation.

2.3.2. Geographically Weighted Regression (GWR)

The geographic weighted regression (GWR) model is an extension of the traditional OLS regression model that generates local regression models for each spatial unit (e.g., each neighborhood in this study). The ordinary least squares (OLS) regression model is a widely used statistical model used to determine the relationships between variables [63,64]. The OLS approach models a continuous dependent variable as a linear function of one or more independent variables, allowing for us to understand how the independent variables contribute to the outcome. However, OLS regression models do not consider spatial variability, nor do they handle spatial autocorrelation in georeferenced data. GWR is a spatial regression model that controls the bias created using spatial autocorrelation and can be used to explore and address the issue of spatial non-stationarity [65]. The equation of the GWR model is shown below:
y i = β i 0 + k = 1 L β i k x i k + ϵ i i = 1 ,   2 ,     ,   n
The dependent variable y i represents the neighborhoods, i, in this study. The local regression coefficients β i k and values for the kth independent variable x i k at each location i are used to predict y i ; β i 0 is the intercept and ε i is the random error, where ε   ~   N 0 ,   σ 2 at location i . Unlike the “global” OLS model where the coefficients β k are fixed, the regression coefficients β i k in the GWR model vary by location. Four continuous (Gaussian) GWR models were created to examine the relationships between environmental factors and T2DM prevalence rates for ages 20 and above, from 20 to 44, from 45 to 64, and 65 and above. ArcGIS Pro (version 3.0) was used to run the GWR models for this study [59].
In ArcGIS Pro, a weighting procedure was applied for each location i. Bandwidths were calculated at each location to determine which other neighborhoods would be included in the estimation of each local model. The observations were weighted based on their proximity to location i [66]. The distance bandwidth was determined using the golden selection procedure, selecting the regression model with the lowest model Akaike information criterion (AIC) that best fits the data. A local t-value was also calculated to estimate the significance of the regression coefficients in this study. This was calculated by dividing the local regression coefficient by the local standard error for each observation (in this case, neighborhoods in Toronto). These t-values act as pseudo-t-statistics, testing the null hypothesis that the regression coefficients are equal to zero. Since GWR tends to overfit the model, the t-values can be used to identify the regional areas where relationships occur [67].

2.3.3. Multicollinearity

Multicollinearity can occur when two or more independent variables are highly correlated in a regression model. The existence of multicollinearity in a regression model inflates the standard errors of regression coefficients and hinders the interpretation of regression results. To ensure that multicollinearity did not exist between the independent variables, this study used the variable inflation factors (VIF) measured from the companion to applied regression (car) package in R [68,69]. A higher VIF value indicates severe multicollinearity. Independent variables included in the regression analyses should have VIFs that are less than four [70]. The selection of variables started with all variables included in the regression model. If any variables were found to have VIF values greater than four, indicating multicollinearity, the variable with the highest VIF value was removed from the model. This process was repeated until all the model variables had VIF values of less than four. There was no evidence of multicollinearity among the selected independent variables (all VIF values < 4); hence, all of them were included in the regression analyses, as shown in Table 2.

3. Results

3.1. Descriptive Statistics

The summary statistics of the T2DM prevalence rates, environment, and socioeconomic variables are shown in Table 3. The prevalence rates of T2DM among citizens aged 20 and above ranged from 2.3% to 22.5% in the neighborhoods of Toronto. The average prevalence rate of T2DM in Toronto neighborhoods was 12.2%, with a standard deviation of 4.4%. The average prevalence rate for the age group of 20–44 was 2.6%, while this increased to 14.6% for the age group of 45–64 and 31.6% for the age group of 65 and above. The standard deviation also increased from 1.14 for the 20–44 age group to 5.81 for the 45–64 age group and 7.8 for the 65 and above age group. On average, 19% of the Toronto neighborhoods were covered with green space. There were an average of 7.6 unhealthy food outlets per km2 in Toronto neighborhoods. Toronto’s neighborhoods also had diverse socioeconomic characteristics, as the percentages of immigrant populations varied from 20.73% and 70.23%, including neighborhoods with both local residences and immigrants. The percentages of residents without certificates, diplomas, or degrees varied from 2.47% to 30.46%. There were neighborhoods with annual incomes ranging from CAD 19,797 to CAD 65,639 and unemployment rates ranging from 5.43% to 12.28%. The variations in environment and socioeconomic status among Toronto’s neighborhoods may contribute to spatial non-stationarity and environmental influences on T2DM prevalence rates.

3.2. GWR Regression

This study created four geographically weighted regression models to examine spatial non-stationarity and environmental effects on the prevalence rates of T2DM in different age groups: “20 and above” (model 1), “20 to 44” (model 2), “45 to 64” (model 3), and “65 and above” (model 4). The results of these models are shown in Table 4, Table 5, Table 6 and Table 7. The global R2 values that measure the proportions of variation in the T2DM prevalence rates explained by the relationships with independent variables were 0.9173, 0.9018, 0.8769, and 0.9012 for the four models, respectively. The local R2 values at the neighborhood level ranged from 0.74 to 0.90, 0.50 to 0.94, 0.67 to 0.93, and 0.76 to 0.95 for the four models, respectively, as shown in Figure 5. Most neighborhoods had local R2 values of 0.80 or above (blue areas in Figure 5e–h), indicating that the GWR modes accurately predicted the T2DM prevalence rates in these neighborhoods. The residuals for the GWR models, shown in Figure 5a–d, were also examined. Positive residuals (mapped as green) represented underestimated T2DM prevalence rates, while negative residuals (mapped as purple) represented overestimated rates. No significant spatial autocorrelations existed in the residuals of any of the four GWR models, as indicated by the Moran’s I values of 0.1227, −0.0318, 0.047, and 0.0529 for the residuals in Figure 5a–d, respectively. The positive and negative coefficients, which estimate the percentages of the corresponding statistically significant coefficients (with |t-values| ≥ 2.0), are presented in Table 4, Table 5, Table 6 and Table 7.
Environmental effects on the T2DM prevalence rates existed after controlling for socioeconomic factors, since most coefficients differed from zero. However, some regression coefficients were not statistically significant (|t-values| ≤ 2.0). This may have been due to the small sample size or high degree of random variation in the regression variables. It is not recommended to interpret the non-significant t-values because we cannot be sure that the values of the associated parameters in the regression model have an effect. Therefore, the regression results were interpreted cautiously by only reporting the geographical areas with significant coefficients.
According to the results of the GWR model 1 shown in Table 4, 18% of the neighborhoods in east and west Toronto showed significant positive associations between the green space density (β1) and the T2DM prevalence rate among those aged 20 and above. In contrast, 3% of the neighborhoods, such as those in the central north, had significant negative associations between the green space density (β1) and T2DM prevalence. A total of 87% of neighborhoods had significant negative associations between the unhealthy food outlet density (β2) and T2DM prevalence. Additionally, the immigration rate (β3), low-education rate (β4), and unemployment rate (β6) had statistically significant positive associations (70%, 92%, and 50%, respectively) with T2DM prevalence in neighborhoods in central Toronto. A total of 20% of the neighborhoods in the northwest or northeast of Toronto had negative relationships between the medium annual income (β5) and T2DM prevalence rates. Figure 6 shows the spatial variation in the local coefficients and t-values for GWR model 1. These results revealed that inconsistent associations existed between the green spaces and T2DM prevalence rates among those aged 20 and above across neighborhoods in Toronto.
The results also showed that spatial non-stationarity effects existed for environmental and socioeconomic factors on the prevalence of T2DM among those aged from 20 to 44 across Toronto neighborhoods. The results are reported in Table 5. A total of 29% of Toronto neighborhoods across the central–west area showed significant positive associations between the green space density (β1) and T2DM prevalence rate. However, 3% of the neighborhoods in the east of Toronto showed significant negative associations between the green space density (β1) and T2DM prevalence rate. Furthermore, an unhealthy food outlet density (β2) showed significant negative associations with T2DM prevalence in 26% of the neighborhoods, primarily in eastern Toronto. The immigration rate (β3), low-education rate (β4), and unemployment rate (β6) showed statistically significant positive associations (45%, 45%, and 55%, respectively) with T2DM prevalence rates in most of the Toronto neighborhoods. However, inconsistent associations also existed with the immigration rate (β3) and low-education rate (β4), where 10% and 5% of Toronto neighborhoods had significant negative associations with T2DM prevalence rates, respectively. The spatial non-stationarity effects can also be seen in Figure 7, where local coefficients change from blue (negative) to red (positive) across Toronto neighborhoods.
Similar to model 2, the green space density (β1) had a significant positive association with the prevalence of T2DM in 24% of the neighborhoods across Toronto for the age group from 45 to 64 years. An unhealthy food outlet density (β2) negatively correlated with T2DM prevalence in 28% of neighborhoods in this age group. Immigration (β3), low education (β4), and unemployment (β6) rates were positively correlated with T2DM prevalence rates in 49%, 55%, and 63% of neighborhoods, respectively, for the age group from 45 to 64 years. In contrast, the prevalence of T2DM negatively correlated with the immigration rate (β3), low-education rate (β4), and median annual income (β5) in 11%, 3%, and 85% of neighborhoods, respectively. The results for model 3 are shown in Table 6 and Figure 8.
Spatial non-stationarity effects were also discovered for both contextual and control variables on the T2DM prevalence rate with the fourth GWR model. The T2DM prevalence rate among those aged 65 and above was included as the dependent variable in GWR model 4. The spatial variations in the local coefficients and t-values of the fourth GWR model can be found in Table 7 and Figure 9. The green space density (β1) yielded significant positive associations with the T2DM prevalence rates among those aged 65 and above in 18% of the neighborhoods across Toronto. These neighborhoods were mainly located in the east and southwest of Toronto. However, unhealthy food outlet locations (β2) had positive associations with T2DM prevalence rates in 3% of neighborhoods (located southwest of Toronto), but negative associations in 17% of neighborhoods (in the east of Toronto). Among the socioeconomic factors, the medium annual income (β5) was negatively correlated with the T2DM prevalence rate among those aged 65 and above in 99% of the neighborhoods. The immigration rate (β3), low-education rate (β4), and unemployment rate (β6) had positive correlations with the T2DM prevalence rates among those aged 65 and above in 30%, 92%, and 14% of the neighborhoods, respectively. Nevertheless, the immigration rate (β3), low-education rate (β4), and unemployment rate (β6) had negative correlations with the T2DM prevalence rates among those aged 65 and above in 11%, 3%, and 10% of the neighborhoods, respectively.

4. Discussion

This study investigated the relationship between both unhealthy food and green space environments with the prevalence of T2DM at the neighborhood level in Toronto. The study used geographically weighted regression (GWR) models to analyze the associations between contextual factors and T2DM prevalence rates in four different age groups: 20 and above, from 20 to 44, from 45 to 64, and 65 and above. The results of all four models indicated that spatial non-stationarity effects existed between environmental characteristics and T2DM prevalence rates. The use of the GWR model successfully addressed the issue of spatial autocorrelations in the T2DM prevalence rate (Moran’s I statistics of 0.7995) across Toronto and produced uncorrelated residuals. Previous studies have often ignored spatial non-stationarity and spatial autocorrelation, which can lead to biased results in environmental health research. Additionally, the inconsistent spatial relationships found between the prevalence rates of T2DM and environmental and socioeconomic factors create uncertainties in environmental health research, in addition to the modifiable areal unit problem (MAUP) and the uncertain geographic context problem (UGCoP). The conflicting findings from previous studies on the relationship between diabetes and the environment may be due to the spatial non-stationarity effects at different geographical locations. Hence, understanding these spatial non-stationarity associations can help inform the development of targeted education, prevention, and intervention policies for specific geographic areas.
The relationships between environmental factors and the prevalence rates of T2DM from four GWR models showed spatial non-stationarity across Toronto neighborhoods. The inconsistent coefficient of the estimates in the GWR models (shown in Table 4, Table 5, Table 6 and Table 7) indicated spatial variation in the effects of the environment on the T2DM prevalence rates. Statistical tests were conducted on the estimates’ coefficients to identify statistically significant ones with | t-values| ≥ 2.00. These significant coefficients of the estimates, representing the results of local regression models, can help us understand the impacts of contextual variables on the T2DM prevalence rates at various geographic locations. For example, the green space density (β1) was significantly positively associated with the prevalence rate of T2DM for 27% of the neighborhoods among young people aged from 20 to 44, but was negatively associated in 3% of the neighborhoods. Figure 6a shows that the associations between the green space density and T2DM prevalence rates varied from negative (blue polygons) to positive (red polygons) across neighborhoods in Toronto. The figure further illustrates that the significant positive coefficients were clustered around the central–west area of Toronto, but the significant negative coefficients were located in the east of Toronto. This showed that when the green space density increased in the central–west area of Toronto, the prevalence rates of T2DM also increased among citizens aged between 20 and 44; in contrast, when the green space density decreased around the east of Toronto, the prevalence rates of T2DM decreased among citizens aged between 20 and 44. On the other hand, when the green space density increased in the east and southwest of Toronto (green polygons in Figure 9g), the prevalence rates of T2DM also increased among the elder generation aged 65 and above. The positive correlation between green spaces and T2DM prevalence rates contradicts the previous research, which found that green spaces promote physical activity and reduce the risks of diabetes [7,8]. It is possible that the quality of green spaces and the presence of noise and unsafe environments in the central–west area of Toronto may have contributed to this contradiction. Poorly designed environments and unsafe neighborhoods may discourage physical activity and increase social isolation due to fear [71]. Elder generations could also have less access to green spaces as activity sites due to their limited mobility in the community compared to the youth.
Unhealthy food outlets (β2) were found to have inconsistent effects on the prevalence rates of T2DM between younger and elder generations. The unhealthy food outlet density was negatively associated with T2DM prevalence among people aged from 20 to 64 in around 25% of the neighborhoods, primarily in eastern Toronto (shown in blue in Figure 7b and Figure 8b). This meant that on average, as the number of unhealthy food outlets decreased, the T2DM prevalence rate for ages between 20 and 64 tended to increase. However, the unhealthy food outlet densities were found to have significant positive relationships with T2DM prevalence in four neighborhoods in southwest Toronto (labeled in green in Figure 9h) among the elder generation aged 65 and above. At the same time, some were negatively associated with T2DM prevalence in eastern Toronto for this age group. These results showed that unhealthy food outlets were associated with low diabetes prevalence rates in some neighborhoods when controlling for socioeconomic status. The results again contradict previous research showing that unhealthy food environments increase the risk of diabetes in neighborhoods [12,22,23,72]. This discrepancy may have been due to the complex interactions between T2DM and socioeconomic factors. For instance, the communities across Toronto were well-educated, where only 13% of the population held no diploma or certificate, as indicated in Table 3. The citizens recognized the adverse health outcomes of consuming unhealthy food and chose healthy food choices, such as cooking at home. On the other hand, the unhealthy fast food restaurants were mainly clustered around central–downtown Toronto due to the agglomeration effects around commercial districts, but were sparsely located around neighborhoods in which negative associations were discovered, as illustrated in Figure 3a. However, the citizens living around communities in downtown Toronto were usually wealthy and well-educated, with greater accessibility and more choices of various food sources. As a result, the socioeconomic and demographic characteristic variations among the neighborhoods may explain the inconsistent and contradictory results.
It is worth noting that the immigration rates (β3) were also found to have spatial non-stationarity effects on the T2DM prevalence rates across Toronto’s neighborhoods. When the immigrant population increased in central–downtown Toronto, the T2DM prevalence rate also increased among all age groups. In contrast, the T2DM prevalence rate decreased when the immigrant population increased north of Toronto. The positive correlation between the immigrant population and T2DM may have been caused by the presence of unhealthy food outlets located around downtown Toronto, which could lead to changes in dietary structures and habits for immigrants, increasing their risk of being diagnosed with T2DM. Additionally, the low-education rate (β4), medium annual income (β5), and unemployment rate (β6) were found to increase, decrease, and increase the T2DM prevalence rates, respectively, across Toronto’s neighborhoods.
This study’s findings indicate the presence of spatial non-stationarity in environmental health studies of T2DM. This is consistent with the expectations but differs from the majority of the previous studies. As a result, it is possible and not surprising that research on the effects of environmental factors on T2DM conducted in different geographic locations would find inverse associations. This is a possible reason why the research findings on the impacts of environmental factors on diabetes health behaviors and outcomes are frequently inconsistent. Recognizing the spatial non-stationarity effects of the environment on T2DM across various geographical locations could help us to better understand the ignored spatial phenomenon of T2DM.
The findings of this study will help to inform the development of customized diabetes intervention and prevention policies. For instance, communities around central Toronto may consider improving the design and attractiveness of green spaces to promote physical activities. Intervention policies and education programs targeting specific social groups and communities with a higher risk of T2DM (such as immigrants and low-income populations) should be developed, as the negative impact of diabetes on health continue to increase. Other potential intervention policies that could be implemented include promoting active transportation (such as walking), providing diabetes education, and establishing clinics in areas with high rates of T2DM prevalence.
This study has several limitations that should be addressed in future research. Firstly, while six contextual variables were selected for the regression analyses, other factors may also influence T2DM. Hence, it would be useful to examine and investigate additional environmental characteristics in future studies, such as the density of sidewalks and the length of cycling trails. Furthermore, the quality of green spaces and neighborhood safety, which can impact physical activity, should also be considered as proxies for contextual variables. A safe neighborhood with attractive landscapes can promote physical activity and lower the risk of diabetes [9]. Secondly, this study did not consider environmental exposures outside the neighborhood, such as workplaces. Although residents spend most of their time in their communities (the most relevant areas that affect health behaviors and outcomes), people’s exposure to social and physical environments is determined by the locations they visit and the time spent moving around for daily activities [73]. As a result, it would be useful to consider people’s movement in space and time when estimating environmental exposures and their effects on health behaviors and outcomes. Thirdly, this cross-sectional study investigated the environmental associations with T2DM in 2019. Longitudinal data analyses using multiple years of data may provide additional insights into how the environment affects the T2DM prevalence rate. Lastly, this study explored the spatial non-stationarity associations between T2DM and environmental factors in various geographical locations at the city level, which has not been well-studied in previous research. However, further research is needed to investigate the causality, particularly the non-stationarity impacts of environmental factors on T2DM.

5. Conclusions

This study found that the relationship between T2DM and environmental factors varied depending on geography and age in Toronto. The results indicated that spatial non-stationarity effects of unhealthy food and green space environments on T2DM exist in Toronto neighborhoods. Future research should consider regional differences to accurately understand the relationships between T2DM and environmental factors and to fully understand the causes of spatial non-stationarity.

Author Contributions

H.G. conceived, designed, and implemented the experiments; analyzed the results; and wrote the paper. J.W. supervised the project and contributed to reviewing and editing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Center for Urban Environments Undergraduate Research Award at the University of Toronto Mississauga, the University of Toronto Mississauga Undergraduate Research Grant, and the Connaught Fund [509252].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset of Toronto’s neighborhood boundary, diabetes prevalence, and population census, supporting this study’s findings, are openly available at Toronto Open Data Portal [42], Ontario Community Health Profiles Partnership [43], and Statistics Canada [40], respectively. Restrictions applied to the availability of green space data and unhealthy food outlet data from DMTI Spatial Inc. [55] and SafeGraph [45], respectively, which were used under license for this study.

Acknowledgments

The authors would like to thank the anonymous reviewers and editors for their valuable comments to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. International Diabetes Federation. IDF Diabetes Atlas, 10th ed. Available online: www.diabetesatlas.org (accessed on 1 October 2022).
  2. Bommer, C.; Heesemann, E.; Sagalova, V.; Manne-Goehler, J.; Atun, R.; Bärnighausen, T.; Vollmer, S. The Global Economic Burden of Diabetes in Adults Aged 20–79 Years: A Cost-of-Illness Study. Lancet Diabetes Endocrinol. 2017, 5, 423–430. [Google Scholar] [CrossRef] [PubMed]
  3. Chatterjee, S.; Khunti, K.; Davies, M.J. Type 2 Diabetes. Lancet 2017, 389, 2239–2251. [Google Scholar] [CrossRef] [PubMed]
  4. Dendup, T.; Feng, X.; Clingan, S.; Astell-Burt, T. Environmental Risk Factors for Developing Type 2 Diabetes Mellitus: A Systematic Review. Int. J. Environ. Res. Public Health 2018, 15, 78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. International Diabetes Federation Type 2 Diabetes. Available online: https://idf.org/aboutdiabetes/type-2-diabetes.html (accessed on 8 March 2022).
  6. Toronto Public Health. Health Surveillance Indicators: Diabetes. Available online: https://www.toronto.ca/wp-content/uploads/2017/12/8c72-tph-hsi-diabetes-july18f.pdf (accessed on 21 August 2022).
  7. Astell-Burt, T.; Feng, X.; Kolt, G.S. Is Neighborhood Green Space Associated with a Lower Risk of Type 2 Diabetes? Evidence from 267,072 Australians. Diabetes Care 2014, 37, 197–201. [Google Scholar] [CrossRef] [Green Version]
  8. Bodicoat, D.H.; O’Donovan, G.; Dalton, A.M.; Gray, L.J.; Yates, T.; Edwardson, C.; Hill, S.; Webb, D.R.; Khunti, K.; Davies, M.J. The Association between Neighbourhood Greenspace and Type 2 Diabetes in a Large Cross-Sectional Study. BMJ Open 2014, 4, e006076. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ulmer, J.M.; Wolf, K.L.; Backman, D.R.; Tretheway, R.L.; Blain, C.J.A.; O’Neil-Dunne, J.P.M.; Frank, L.D. Multiple Health Benefits of Urban Tree Canopy: The Mounting Evidence for a Green Prescription. Health Place 2016, 42, 54–62. [Google Scholar] [CrossRef] [PubMed]
  10. Paquet, C.; Coffee, N.T.; Haren, M.T.; Howard, N.J.; Adams, R.J.; Taylor, A.W.; Daniel, M. Food Environment, Walkability, and Public Open Spaces Are Associated with Incident Development of Cardio-Metabolic Risk Factors in a Biomedical Cohort. Health Place 2014, 28, 173–176. [Google Scholar] [CrossRef] [PubMed]
  11. Mezuk, B.; Li, X.; Cederin, K.; Rice, K.; Sundquist, J.; Sundquist, K. Beyond Access: Characteristics of the Food Environment and Risk of Diabetes. Am. J. Epidemiol. 2016, 183, 1129–1137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ahern, M.; Brown, C.; Dukas, S. A National Study of the Association between Food Environments and County-level Health Outcomes. J. Rural. Health 2011, 27, 367–379. [Google Scholar] [CrossRef]
  13. Gebreab, S.Y.; Hickson, D.M.A.; Sims, M.; Wyatt, S.B.; Davis, S.K.; Correa, A.; Diez-Roux, A.V. Neighborhood Social and Physical Environments and Type 2 Diabetes Mellitus in African Americans: The Jackson Heart Study. Health Place 2017, 43, 128–137. [Google Scholar] [CrossRef]
  14. Sundquist, K.; Eriksson, U.; Mezuk, B.; Ohlsson, H. Neighborhood Walkability, Deprivation and Incidence of Type 2 Diabetes: A Population-Based Study on 512,061 Swedish Adults. Health Place 2015, 31, 24–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Creatore, M.I.; Glazier, R.H.; Moineddin, R.; Fazli, G.S.; Johns, A.; Gozdyra, P.; Matheson, F.I.; Kaufman-Shriqui, V.; Rosella, L.C.; Manuel, D.G. Association of Neighborhood Walkability with Change in Overweight, Obesity, and Diabetes. JAMA 2016, 315, 2211–2220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Glazier, R.H.; Creatore, M.I.; Weyman, J.T.; Fazli, G.; Matheson, F.I.; Gozdyra, P.; Moineddin, R.; Shriqui, V.K.; Booth, G.L. Density, Destinations or Both? A Comparison of Measures of Walkability in Relation to Transportation Behaviors, Obesity and Diabetes in Toronto, Canada. PLoS ONE 2014, 9, e85295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Booth, G.L.; Creatore, M.I.; Moineddin, R.; Gozdyra, P.; Weyman, J.T.; Matheson, F.I.; Glazier, R.H. Unwalkable Neighborhoods, Poverty, and the Risk of Diabetes among Recent Immigrants to Canada Compared with Long-Term Residents. Diabetes Care 2013, 36, 302–308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Shin, S.; Bai, L.; Oiamo, T.H.; Burnett, R.T.; Weichenthal, S.; Jerrett, M.; Kwong, J.C.; Goldberg, M.S.; Copes, R.; Kopp, A.; et al. Association between Road Traffic Noise and Incidence of Diabetes Mellitus and Hypertension in Toronto, Canada: A Population-Based Cohort Study. J. Am. Heart Assoc. 2020, 9, e013021. [Google Scholar] [CrossRef]
  19. Müller-Riemenschneider, F.; Pereira, G.; Villanueva, K.; Christian, H.; Knuiman, M.; Giles-Corti, B.; Bull, F.C. Neighborhood Walkability and Cardiometabolic Risk Factors in Australian Adults: An Observational Study. BMC Public Health 2013, 13, 1–9. [Google Scholar] [CrossRef] [Green Version]
  20. Requia, W.J.; Adams, M.D.; Koutrakis, P. Association of PM2.5 with Diabetes, Asthma, and High Blood Pressure Incidence in Canada: A Spatiotemporal Analysis of the Impacts of the Energy Generation and Fuel Sales. Sci. Total Environ. 2017, 584–585, 1077–1083. [Google Scholar] [CrossRef]
  21. Heidemann, C.; Niemann, H.; Paprott, R.; Du, Y.; Rathmann, W.; Scheidt-Nave, C. Residential Traffic and Incidence of Type 2 Diabetes: The German Health Interview and Examination Surveys. Diabet. Med. 2014, 31, 1269–1276. [Google Scholar] [CrossRef]
  22. Bodicoat, D.H.; Carter, P.; Comber, A.; Edwardson, C.; Gray, L.J.; Hill, S.; Webb, D.; Yates, T.; Davies, M.J.; Khunti, K. Is the Number of Fast-Food Outlets in the Neighbourhood Related to Screen-Detected Type 2 Diabetes Mellitus and Associated Risk Factors? Public Health Nutr. 2015, 18, 1698–1705. [Google Scholar] [CrossRef] [Green Version]
  23. Ntarladima, A.M.; Karssenberg, D.; Poelman, M.; Grobbee, D.E.; Lu, M.; Schmitz, O.; Strak, M.; Janssen, N.; Hoek, G.; Vaartjes, I. Associations between the Fast-Food Environment and Diabetes Prevalence in the Netherlands: A Cross-Sectional Study. Lancet Planet Health 2022, 6, e29–e39. [Google Scholar] [CrossRef]
  24. Li, R.; Chen, G.; Jiao, A.; Lu, Y.; Guo, Y.; Li, S.; Wang, C.; Xiang, H. Residential Green and Blue Spaces and Type 2 Diabetes Mellitus: A Population-Based Health Study in China. Toxics 2021, 9, 11. [Google Scholar] [CrossRef] [PubMed]
  25. Drewnowski, A.; Rehm, C.D.; Moudon, A.V.; Arterburn, D. The Geography of Diabetes by Census Tract in a Large Sample of Insured Adults in King County, Washington, 2005–2006. Prev. Chronic Dis. 2014, 11, E125. [Google Scholar] [CrossRef] [Green Version]
  26. Kolpak, P.; Wang, L. Exploring the Social and Neighbourhood Predictors of Diabetes: A Comparison between Toronto and Chicago. Prim. Health Care Res. Dev. 2017, 18, 291–299. [Google Scholar] [CrossRef] [PubMed]
  27. Christine, P.J.; Auchincloss, A.H.; Bertoni, A.G.; Carnethon, M.R.; Sánchez, B.N.; Moore, K.; Adar, S.D.; Horwich, T.B.; Watson, K.E.; Roux, A.V.D. Longitudinal Associations between Neighborhood Physical and Social Environments and Incident Type 2 Diabetes Mellitus: The Multi-Ethnic Study of Atherosclerosis (MESA). JAMA Intern. Med. 2015, 175, 1311–1320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Frankenfeld, C.L.; Leslie, T.F.; Makara, M.A. Diabetes, Obesity, and Recommended Fruit and Vegetable Consumption in Relation to Food Environment Sub-Types: A Cross-Sectional Analysis of Behavioral Risk Factor Surveillance System, United States Census, and Food Establishment Data. BMC Public Health 2015, 15, 1–9. [Google Scholar] [CrossRef] [Green Version]
  29. Piccolo, R.S.; Duncan, D.T.; Pearce, N.; McKinlay, J.B. The Role of Neighborhood Characteristics in Racial/Ethnic Disparities in Type 2 Diabetes: Results from the Boston Area Community Health (BACH) Survey. Soc. Sci. Med. 2015, 130, 79–90. [Google Scholar] [CrossRef] [Green Version]
  30. Chowdhury, M.A.B.; Uddin, M.J.; Khan, H.M.R.; Haque, M.R. Type 2 Diabetes and Its Correlates among Adults in Bangladesh: A Population Based Study. BMC Public Health 2015, 15, 1070. [Google Scholar] [CrossRef] [Green Version]
  31. Wu, H.; Meng, X.; Wild, S.H.; Gasevic, D.; Jackson, C.A. Socioeconomic Status and Prevalence of Type 2 Diabetes in Mainland China, Hong Kong and Taiwan: A Systematic Review. J. Glob. Health 2017, 7, 011103. [Google Scholar] [CrossRef] [PubMed]
  32. Wang, F. Why Public Health Needs GIS: A Methodological Overview. Ann. GIS 2019, 26, 1–12. [Google Scholar] [CrossRef] [Green Version]
  33. Kwan, M.P. The Uncertain Geographic Context Problem. Environ. Health 2012, 102, 958–968. [Google Scholar] [CrossRef]
  34. Buzzelli, M. Modifiable Areal Unit Problem. Int. Encycl. Hum. Geogr. 2020, 169. [Google Scholar] [CrossRef]
  35. Kwan, M.-P. The Stationarity Bias in Research on the Environmental Determinants of Health. Health Place 2021, 70, 102609. [Google Scholar] [CrossRef]
  36. Wang, J.; Lee, K.; Kwan, M.-P. Environmental Influences on Leisure-Time Physical Inactivity in the US: An Exploration of Spatial Non-Stationarity. ISPRS Int. J. Geoinf. 2018, 7, 143. [Google Scholar] [CrossRef] [Green Version]
  37. Siordia, C.; Saenz, J.; Tom, S.E. An Introduction to Macro- Level Spatial Nonstationarity: A Geographically Weighted Regression Analysis of Diabetes and Poverty. Hum. Geogr. 2012, 6, 5. [Google Scholar] [CrossRef] [PubMed]
  38. Hipp, J.A.; Chalise, N. Peer Reviewed: Spatial Analysis and Correlates of County-Level Diabetes Prevalence, 2009–2010. Prev. Chronic Dis. 2015, 12, 140404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Margolis, D.J.; Hoffstad, O.; Nafash, J.; Leonard, C.E.; Freeman, C.P.; Hennessy, S.; Wiebe, D.J. Location, Location, Location: Geographic Clustering of Lower-Extremity Amputation Among Medicare Beneficiaries with Diabetes. Diabetes Care 2011, 34, 2363–2367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Statistics Canada 2016 Census of Population—Data Products. Available online: https://www12.statcan.gc.ca/census-recensement/2016/dp-pd/index-eng.cfm (accessed on 16 October 2022).
  41. City of Toronto. Population Demographics of Toronto. Available online: https://www.toronto.ca/wp-content/uploads/2019/11/99b4-TOHealthCheck_2019Chapter1.pdf (accessed on 14 October 2022).
  42. Toronto Open Data Portal. Neighbourhoods—City of Toronto Open Data Portal. Available online: https://open.toronto.ca/dataset/neighbourhoods/ (accessed on 16 January 2023).
  43. Ontario Community Health Profiles Partnership. Available online: http://www.ontariohealthprofiles.ca/dataTablesON.php?varTab=HPDtbl&select1=7 (accessed on 8 March 2022).
  44. ICES. Available online: https://www.ices.on.ca (accessed on 11 September 2022).
  45. SafeGraph. Available online: https://www.safegraph.com/ (accessed on 14 October 2022).
  46. Chang, S.; Pierson, E.; Koh, P.W.; Gerardin, J.; Redbird, B.; Grusky, D.; Leskovec, J. Mobility Network Models of COVID-19 Explain Inequities and Inform Reopening. Nature 2020, 589, 82–87. [Google Scholar] [CrossRef] [PubMed]
  47. Kupfer, J.A.; Li, Z.; Ning, H.; Huang, X. Using Mobile Device Data to Track the Effects of the COVID-19 Pandemic on Spatiotemporal Patterns of National Park Visitation. Sustainability 2021, 13, 9366. [Google Scholar] [CrossRef]
  48. Zhou, H.; Kim, G.; Wang, J.; Wilson, K. Investigating the Association between the Socioeconomic Environment of the Service Area and Fast Food Visitation: A Context-Based Crystal Growth Approach. Health Place 2022, 76, 102855. [Google Scholar] [CrossRef]
  49. Polsky, J.Y.; Moineddin, R.; Glazier, R.H.; Dunn, J.R.; Booth, G.L. Foodscapes of Southern Ontario: Neighbourhood Deprivation and Access to Healthy and Unhealthy Food Retail. Can. J. Public Health 2014, 105, e369–e375. [Google Scholar] [CrossRef] [PubMed]
  50. Mahendra, A.; Polsky, J.Y.; Robitaille, É.; Lefebvre, M.; McBrien, T.; Minaker, L.M. Geographic Retail Food Environment Measures for Use in public Health. Health Promot. Chronic Dis. Prev. Can. 2017, 37, 357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Powell, L.M.; Han, E.; Zenk, S.N.; Khan, T.; Quinn, C.M.; Gibbs, K.P.; Pugach, O.; Barker, D.C.; Resnick, E.A.; Myllyluoma, J.; et al. Field Validation of Secondary Commercial Data Sources on the Retail Food Outlet Environment in the U.S. Health Place 2011, 17, 1122–1131. [Google Scholar] [CrossRef] [PubMed]
  52. De la Fuente, F.; Saldías, M.A.; Cubillos, C.; Mery, G.; Carvajal, D.; Bowen, M.; Bertoglia, M.P. Green Space Exposure Association with Type 2 Diabetes Mellitus, Physical Activity, and Obesity: A Systematic Review. Int. J. Environ. Res. Public Health 2020, 18, 97. [Google Scholar] [CrossRef] [PubMed]
  53. Hayes, C.; Kriska, A. Role of Physical Activity in Diabetes Management and Prevention. J. Am. Diet. Assoc. 2008, 108, S19–S23. [Google Scholar] [CrossRef] [PubMed]
  54. Richardson, E.A.; Pearce, J.; Mitchell, R.; Kingham, S. Role of Physical Activity in the Relationship between Urban Green Space and Health. Public Health 2013, 127, 318–324. [Google Scholar] [CrossRef] [Green Version]
  55. DMTI Spatial Inc. Park Sports Field Region—2019. Available online: https://geo.scholarsportal.info/#r/details/_uri@=11120733$DMTI_2019_CMCS_ParksSportsFieldRegion (accessed on 15 October 2022).
  56. Cheng, Y.J.; Kanaya, A.M.; Araneta, M.R.G.; Saydah, S.H.; Kahn, H.S.; Gregg, E.W.; Fujimoto, W.Y.; Imperatore, G. Prevalence of Diabetes by Race and Ethnicity in the United States, 2011–2016. JAMA 2019, 322, 2389–2398. [Google Scholar] [CrossRef]
  57. Ledford, C.J.W.; Seehusen, D.A.; Crawford, P.F. Geographic and Race/Ethnicity Differences in Patient Perceptions of Diabetes. J. Prim. Care Community Health 2019, 10, 2150132719845819. [Google Scholar] [CrossRef]
  58. Tenkorang, E.Y. Early Onset of Type 2 Diabetes among Visible Minority and Immigrant Populations in Canada. Ethn. Health 2016, 22, 266–284. [Google Scholar] [CrossRef]
  59. Esri Inc. ArcGIS Pro 3.0. Available online: https://www.esri.com/en-us/arcgis/products/arcgis-pro/overview (accessed on 16 January 2023).
  60. Chen, X. Take the Edge off: A Hybrid Geographic Food Access Measure. Appl. Geogr. 2017, 87, 149–159. [Google Scholar] [CrossRef]
  61. Chen, X.; Ye, X.; Widener, M.J.; Delmelle, E.; Kwan, M.-P.; Shannon, J.; Racine, E.F.; Adams, A.; Liang, L.; Jia, P. A Systematic Review of the Modifiable Areal Unit Problem (MAUP) in Community Food Environmental Research. Urban Inform. 2022, 1, 1–11. [Google Scholar] [CrossRef]
  62. Needham, C.; Sacks, G.; Orellana, L.; Robinson, E.; Allender, S.; Strugnell, C. A Systematic Review of the Australian Food Retail Environment: Characteristics, Variation by Geographic Area, Socioeconomic Position and Associations with Diet and Obesity. Obes. Rev. 2020, 21, e12941. [Google Scholar] [CrossRef] [PubMed]
  63. Leung, Y.; Mei, C.L.; Zhang, W.X. Statistical Tests for Spatial Nonstationarity Based on the Geographically Weighted Regression Model. Environ. Plan. A Econ. Space 2016, 32, 9–32. [Google Scholar] [CrossRef] [Green Version]
  64. Wheeler, D.; Tiefelsdorf, M. Multicollinearity and Correlation among Local Regression Coefficients in Geographically Weighted Regression. J. Geogr. Syst. 2005, 7, 161–187. [Google Scholar] [CrossRef]
  65. Brunsdon, C.; Fotheringham, A.S.; Charlton, M.E. Geographically Weighted Regression: A Method for Exploring Spatial Nonstationarity. Geogr. Anal. 1996, 28, 281–298. [Google Scholar] [CrossRef]
  66. Powers, S.L.; Matthews, S.A.; Mowen, A.J. Does the Relationship between Racial, Ethnic, and Income Diversity and Social Capital Vary across the United States? A County-Level Analysis Using Geographically Weighted Regression. Appl. Geogr. 2021, 130, 102446. [Google Scholar] [CrossRef]
  67. Gebreab, S.Y.; Diez Roux, A.V. Exploring Racial Disparities in CHD Mortality between Blacks and Whites across the United States: A Geographically Weighted Regression Approach. Health Place 2012, 18, 1006–1014. [Google Scholar] [CrossRef] [Green Version]
  68. John Fox, S.W. Car: Companion to Applied Regression. Available online: https://cran.r-project.org/package=car (accessed on 16 October 2022).
  69. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.r-project.org/ (accessed on 17 October 2022).
  70. Garson, G.D. Testing Statistical Assumptions; Statistical Associates Publishing: Asheboro, NC, USA, 2012. [Google Scholar]
  71. Bennett, G.G.; McNeill, L.H.; Wolin, K.Y.; Duncan, D.T.; Puleo, E.; Emmons, K.M. Safe to Walk? Neighborhood Safety and Physical Activity among Public Housing Residents. PLoS Med. 2007, 4, e306. [Google Scholar] [CrossRef] [Green Version]
  72. Salois, M.J. Obesity and Diabetes, the Built Environment, and the ‘Local’Food Economy in the United States, 2007. Econ. Hum. Biol. 2012, 10, 35–42. [Google Scholar] [CrossRef] [Green Version]
  73. Kwan, M.P.; Wang, J.; Tyburski, M.; Epstein, D.H.; Kowalczyk, W.J.; Preston, K.L. Uncertainties in the Geographic Context of Health Behaviors: A Study of Substance Users’ Exposure to Psychosocial Stress Using GPS Data. Int. J. Geogr. Inf. Sci. 2019, 33, 1176–1195. [Google Scholar] [CrossRef]
Figure 1. Using geographically weighted regression model to examine the spatial non-stationarity effects of green space density (a) and unhealthy food environment (b) with container buffer (d) for type-2 diabetes mellitus prevalence rates (c). With mappings of local coefficients (e) and t-value (f).
Figure 1. Using geographically weighted regression model to examine the spatial non-stationarity effects of green space density (a) and unhealthy food environment (b) with container buffer (d) for type-2 diabetes mellitus prevalence rates (c). With mappings of local coefficients (e) and t-value (f).
Sustainability 15 01762 g001
Figure 2. Neighborhood-level type-2 diabetes mellitus prevalence rates in Toronto, April 2019.
Figure 2. Neighborhood-level type-2 diabetes mellitus prevalence rates in Toronto, April 2019.
Sustainability 15 01762 g002
Figure 3. Distribution of (a) unhealthy food outlets and (b) green spaces around Toronto.
Figure 3. Distribution of (a) unhealthy food outlets and (b) green spaces around Toronto.
Sustainability 15 01762 g003
Figure 4. Analysis of Anselin’s local Moran’s index (a) for type-2 diabetes mellitus prevalence rates in Toronto neighborhoods with results of z-score (b), p-value (c), and statistically significant cluster (d).
Figure 4. Analysis of Anselin’s local Moran’s index (a) for type-2 diabetes mellitus prevalence rates in Toronto neighborhoods with results of z-score (b), p-value (c), and statistically significant cluster (d).
Sustainability 15 01762 g004
Figure 5. GWR model residuals (ad) and R2 values (eh) by neighborhood.
Figure 5. GWR model residuals (ad) and R2 values (eh) by neighborhood.
Sustainability 15 01762 g005
Figure 6. Spatial variation of local coefficients and t-values of GWR model 1 with T2DM prevalence rates in those aged 20 and above. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Figure 6. Spatial variation of local coefficients and t-values of GWR model 1 with T2DM prevalence rates in those aged 20 and above. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Sustainability 15 01762 g006
Figure 7. Spatial variation of local coefficients and t-values of GWR model 2 with T2DM prevalence rates in those aged 20 to 44. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Figure 7. Spatial variation of local coefficients and t-values of GWR model 2 with T2DM prevalence rates in those aged 20 to 44. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Sustainability 15 01762 g007
Figure 8. Spatial variation of local coefficients and t-values of GWR model 3 with T2DM prevalence rates in those aged 45 to 64. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Figure 8. Spatial variation of local coefficients and t-values of GWR model 3 with T2DM prevalence rates in those aged 45 to 64. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Sustainability 15 01762 g008
Figure 9. Spatial variation of local coefficients and t-values of GWR model 4 with T2DM prevalence rates in those aged 65 and above. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Figure 9. Spatial variation of local coefficients and t-values of GWR model 4 with T2DM prevalence rates in those aged 65 and above. (af) Local coefficients of the green space density (% km2), unhealthy food outlet density (count per km2), immigration rate (% population), low-education rate (% population), medium annual income ($), unemployment rate (% population), and (gl) corresponding local t-values of all predictors.
Sustainability 15 01762 g009
Table 1. Toronto’s 158 neighborhoods’ population demographics by age group.
Table 1. Toronto’s 158 neighborhoods’ population demographics by age group.
PopulationMinimumMaximumMeanMedianStandard DeviationTotal
All ages 20+536029,45715,16714,39254192,396,337
Age 20 to 44231819,4637241673830641,144,070
Age 45 to 6418459726496646341837784,704
Age 65 and above6797628295927561340467,563
Table 2. List of variables used in geographically weighted regression analyses.
Table 2. List of variables used in geographically weighted regression analyses.
VariableDescriptionYearSource
Dependent Variable
Type-2 Diabetes Mellitus Prevalence RatesTotal cases of diabetes by population in neighborhood2019Ontario Community Health Profiles Partnership
Independent Variable
Unhealthy Food Outlet Density
(Count per km2)
Number of locations of limited-service (fast food) restaurants, confectionery retailers, bakeries, and convenience stores by neighborhood area2019SafeGraph
Green Space Density
(% km2)
Area of parks and recreation spaces by neighborhood area2019DMTI Spatial Inc.
Medium Total IncomeMedian total income among recipients ($)20152016 Canadian Census
Unemployment RatePercentage of residents who are unemployed2016
Low-Education RatePercentage of residents who have not obtained any certificates, diplomas, or degrees2016
Immigration RatePercentage of the residents who are, or who have ever been, landed immigrants and permanent residents2016
Table 3. Descriptive statistics for T2DM, environmental, and socioeconomic variables.
Table 3. Descriptive statistics for T2DM, environmental, and socioeconomic variables.
VariablesMinimumMaximumMeanMedianStandard Deviation
T2DM Prevalence Rate
Age 20+ (%)
2.3022.5012.2012.004.40
T2DM Prevalence Rate
Age 20–44 (%)
0.905.502.602.301.14
T2DM Prevalence Rate
Age 45–64 (%)
4.7031.0014.6013.705.81
T2DM Prevalence Rate
Age 65+ (%)
14.7050.0031.6030.907.80
Green Space Density (% km2)7.7064.0019.0020.308.60
Unhealthy Food Outlet Density
(Count per km2)
1.6652.07.612.612.4
Immigration Rate (% Population)20.733470.233044.946345.849912.4900
Low-Education Rate (% Population)2.469030.462613.285513.11086.1950
Medium Annual Income ($)19,79765,63934,70132,38710,360
Unemployment Rate
(% Population)
4.527212.28397.90107.49821.5618
Table 4. Summary of GWR model 1.
Table 4. Summary of GWR model 1.
VariablesGWR Model 1
Response Variable: T2DM Prevalence Rate (Age 20 and Above)
(R2 = 0.9173; Adjusted R2 = 0.8982; AIC = 581.1; Distance Band = 13.3870 km)
Positive Coefficient Estimates (%)Significant Positive Coefficient Estimates (%) †Negative Coefficient Estimates
(%)
Significant Negative Coefficient Estimates (%) ††Minimum Coefficient EstimateMedian Coefficient EstimateMean Coefficient EstimateMaximum Coefficient Estimate
Green Space Density (β1)7218283−10.52582.455811.520511.5205
Unhealthy Food Outlet Density
(β2)
0010087−1.1155−0.0547−0.1761−0.0251
Immigration Rate (β3)8970113−17.71478.33767.468330.153
Low-Education Rate (β4)979230−24.315125.901526.950853.8229
Medium Annual Income (β5)809220−0.000355−0.000042−0.0000670.000029
Unemployment Rate (β6)965040−0.04750.37330.34750.6921
MeanMedianStandard DeviationMinimumMaximum
Local R20.85650.86480.037510.74080.9063
Residual0.04250.02561.2647−3.14602.830
† Number of positive significant (t-value ≥ 2.00) coefficient estimates/number of neighborhoods. †† Number of negative significant (t-value ≤ −2.00) coefficient estimates/number of neighborhoods. AIC: Akaike information criterion; sample size: n = 158.
Table 5. Summary of GWR model 2.
Table 5. Summary of GWR model 2.
VariablesGWR Model 2
Response Variable: T2DM Prevalence Rate (Age 20 to 44)
(R2 = 0.9018; Adjusted R2 = 0.8568; AIC = 223.89; Distance Band = 9.61 km)
Positive Coefficient Estimates (%)Significant Positive Coefficient Estimates (%) †Negative Coefficient Estimates
(%)
Significant Negative Coefficient Estimates (%) ††Minimum Coefficient EstimateMedian Coefficient EstimateMean Coefficient EstimateMaximum Coefficient Estimate
Green Space Density (β1)7129283−2.88412.18851.78345.8218
Unhealthy Food Outlet Density
(β2)
109926−0.4595−0.01710.05880.1266
Immigration Rate (β3)78452211−13.03592.94341.643715.2293
Low-Education Rate (β4)8845113−15.75214.32533.985912.9259
Medium Annual Income (β5)430575−0.0001−0.000005−0.000010.000059
Unemployment Rate (β6)8955110−0.15510.16210.19200.5348
MeanMedianStandard DeviationMinimumMaximum
Local R20.78940.79740.09200.49460.9362
Residual0.0289−0.01280.9730−3.15912.3771
† Number of positive significant (t-value ≥ 2.00) coefficient estimates/number of neighborhoods. †† Number of negative significant (t-value ≤ −2.00) coefficient estimates/number of neighborhoods. AIC: Akaike information criterion; sample size: n = 152.
Table 6. Summary of GWR model 3.
Table 6. Summary of GWR model 3.
VariablesGWR Model 3
Response Variable: T2DM Prevalence Rate (Age 45 to 64)
(R2 = 0.8769; Adjusted R2 = 0.8479; AIC = 705.04; Distance Band = 13.5573 km)
Positive Coefficient Estimates (%)Significant Positive Coefficient Estimates (%) †Negative Coefficient Estimates
(%)
Significant Negative Coefficient Estimates (%) ††Minimum Coefficient EstimateMedian Coefficient EstimateMean Coefficient EstimateMaximum Coefficient Estimate
Green Space Density (β1)8224181−10.23525.79406.048021.5513
Unhealthy Food Outlet Density
(β2)
2107928−2.389−0.0204−0.23080.2108
Immigration Rate (β3)71492911−59.29939.69306.339858.8553
Low-Education Rate (β4)905593−75.493518.195715.527349.1686
Medium Annual Income (β5)0010085−0.001003−0.000164−0.0002220.000067
Unemployment Rate (β6)10063000.01340.81200.83131.8059
MeanMedianStandard DeviationMinimumMaximum
Local R20.81750.82190.06070.66860.9268
Residual0.1286−0.00972.0344−8.13367.1685
† Number of positive significant (t-value ≥ 2.00) coefficient estimates/number of neighborhoods. †† Number of negative significant (t-value ≤ −2.00) coefficient estimates/number of neighborhoods. AIC: Akaike information criterion; sample size: n = 152.
Table 7. Summary of GWR model 4.
Table 7. Summary of GWR model 4.
VariablesGWR Model 4
Response Variable: T2DM Prevalence Rate (Age 65 and above)
(R2 = 0.9012; Adjusted R2 = 0.8813; AIC = 784.66; Distance Band = 14.55 km)
Positive Coefficient Estimates (%)Significant Positive Coefficient Estimates (%) †Negative Coefficient Estimates
(%)
Significant Negative Coefficient Estimates (%) ††Minimum Coefficient EstimateMedian Coefficient EstimateMean Coefficient EstimateMaximum Coefficient Estimate
Green Space Density (β1)8118190−9.84955.69455.630126.1653
Unhealthy Food Outlet Density
(β2)
6233817−2.09390.0264−0.12920.1017
Immigration Rate (β3)67303311−43.76935.50092.968057.6589
Low-Education Rate (β4)949263−82.598248.938744.385672.6199
Medium Annual Income (β5)0010099−0.00125−0.000373−0.000419−0.000261
Unemployment Rate (β6)57144310−1.15610.14960.11201.2070
MeanMedianStandard DeviationMinimumMaximum
Local R20.87590.87820.04280.75760.9514
Residual0.12200.09912.4626−7.55157.2896
† Number of positive significant (t-value ≥ 2.00) coefficient estimates/number of neighborhoods. †† Number of negative significant (t-value ≤ −2.00) coefficient estimates/number of neighborhoods. AIC: Akaike information criterion; sample size: n = 158.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ge, H.; Wang, J. Spatial Non-Stationarity Effects of Unhealthy Food Environments and Green Spaces for Type-2 Diabetes in Toronto. Sustainability 2023, 15, 1762. https://doi.org/10.3390/su15031762

AMA Style

Ge H, Wang J. Spatial Non-Stationarity Effects of Unhealthy Food Environments and Green Spaces for Type-2 Diabetes in Toronto. Sustainability. 2023; 15(3):1762. https://doi.org/10.3390/su15031762

Chicago/Turabian Style

Ge, Haoxuan, and Jue Wang. 2023. "Spatial Non-Stationarity Effects of Unhealthy Food Environments and Green Spaces for Type-2 Diabetes in Toronto" Sustainability 15, no. 3: 1762. https://doi.org/10.3390/su15031762

APA Style

Ge, H., & Wang, J. (2023). Spatial Non-Stationarity Effects of Unhealthy Food Environments and Green Spaces for Type-2 Diabetes in Toronto. Sustainability, 15(3), 1762. https://doi.org/10.3390/su15031762

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop