Next Article in Journal
The Study of Utility Valuation of Single-Name Credit Derivatives with the Fast-Scale Stochastic Volatility Correction
Next Article in Special Issue
Viability of Green Roofs as a Flood Mitigation Element in the Central Region of Chile
Previous Article in Journal
Dynamic Change in the Water-Level Fluctuation Zone of the Danjiangkou Reservoir and Its Influence on Water Quality
Previous Article in Special Issue
Natural Stabilized Earth Panels versus Conventional Façade Systems. Economic and Environmental Impact Assessment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geographically Weighted Regression Models in Estimating Median Home Prices in Towns of Massachusetts Based on an Urban Sustainability Framework

Department of Earth and Environment, Boston University, Boston, MA 02215, USA
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(4), 1026; https://0-doi-org.brum.beds.ac.uk/10.3390/su10041026
Submission received: 21 February 2018 / Revised: 25 March 2018 / Accepted: 26 March 2018 / Published: 30 March 2018
(This article belongs to the Collection Sustainable Built Environment)

Abstract

:
Housing is a key component of urban sustainability. The objective of this study was to assess the significance of key spatial determinants of median home price in towns in Massachusetts that impact sustainable growth. Our analysis investigates the presence or absence of spatial non-stationarity in the relationship between sustainable growth, measured in terms of the relationship between home values and various parameters including the amount of unprotected forest land, residential land, unemployment, education, vehicle ownership, accessibility to commuter rail stations, school district performance, and senior population. We use the standard geographically weighted regression (GWR) and Mixed GWR models to analyze the effects of spatial non-stationarity. Mixed GWR performed better than GWR in terms of Akaike Information Criterion (AIC) values. Our findings highlight the nature and spatial extent of the non-stationary vs. stationary qualities of key environmental and social determinants of median home price. Understanding the key determinants of housing values, such as valuation of green spaces, public school performance metrics, and proximity to public transport, enable towns to use different strategies of sustainable urban planning, while understanding urban housing determinants—such as unemployment and senior population—can help modify urban sustainable housing policies.

1. Introduction

1.1. Background

More than 83 percent of the US population now lives in cities, up thirty percent from 50 years ago. By 2050, the US urban population is projected to increase to more than 90 percent (of 423 million in 2050). Urbanization is transforming farmland, wetlands, forests, and other natural ecosystems into urban landscapes at an unprecedented rate resulting in urban sprawl. Urban landscape patterns and dynamics are the physical manifestation of complex interactions between environment, society, and economy [1,2,3,4]. Thus, urban areas are highly relevant, if not central, to any discussion on sustainable development.
The central goal of urban sustainability is efficient use of natural resources within a city region, while simultaneously improving its livability, through social amenities, economic opportunity, and health [5]. We offer a conceptual framework to understanding urban sustainability through the lens of urban housing that is at the intersection of economic, ecological and social dimensions. Housing plays a pivotal role in determining the financial (economic) security and well-being of individuals, neighborhoods, and cities. Owning a home is part of the “American Dream”. Research on housing has focused on determining home prices, traditionally modeled in economics using hedonic models [6]. Some factors that are typically included in hedonic price analyses are socio-economic and demographic factors such as housing age/tenure, vacancy rates, racial and ethnic demographics [7]. Boston metro has been well studied in the context of hedonic pricing models for the last 3 decades. An early paper [8] studied patterns of home price appreciation in the Boston area from 1982 to 1994 and showed that changes in the cross-sectional pattern of home prices are related to differences in manufacturing employment, demographics, new construction, proximity to the downtown, and to aggregate school enrollments. Hedonic pricing models use parcel level housing data to estimate model parameters. Social factors including income, population density, race, and education impact urban housing and school districts. This has led to a renewed interest in social justice addressing accessibility to parks [9], food [10], and environmental health [11] in the context of urban sustainability. Ecological studies in the urban sustainability context have examined challenges and opportunities in addressing imbalances between green and built landscapes [12] and more recently in mitigating urban carbon footprint [13]. Urban sustainability research traditionally has addressed only a single ecological issue in individual neighborhoods [14]. We need to expand this research to include regional or state level scales in order to address sustainability in a larger context.
In the last two decades, spatial econometric studies have been introduced that demonstrate the fundamental role of spatial autocorrelation and spatial heterogeneity in the analysis of housing in varied regional contexts [15,16,17]. Spatial models use a variety of spatial scales, ranging from parcel level data to census town level data, to address different questions in relation to urban patterns and dynamics. Prior studies have focused either on individual housing units [18,19,20], or slightly larger scale, i.e., districts or metro areas [21,22,23]. At a regional or continental scale, remote sensing instruments are the most viable option for data-driven mapping, monitoring, and assessment of the urban systems [24,25,26]. Given the ready availability of remote sensing data, many studies have incorporated the ecosystem service valuation of green space (forest, woodland, or parks) into pricing models [9], thus opening new avenues of research in urban sustainability in the coming decades.
The objective of this study was to assess the significance of key economic, environmental and social determinants of median home prices in Massachusetts towns using spatial models at the level of the state, including many urban agglomerations. We used a sustainability lens in setting our theoretical framework and incorporated the effects of environmental, economic and social factors over a time period. In a sustainable development context, the dynamics of housing growth need to be examined in order to formulate relevant policy or development decisions. Real estate market dynamics are inseparable from the sustainability concept as cities are designing green buildings, public housing, and affordable housing [26]. This analysis uses a select set of determinants based on prior empirical research. Population density [27] is often correlated to income levels; wealthier affluent communities are characterized by higher value homes. Proximity to highways or public transport impacts home values [28]. Residential development is also constrained by specific town conservation laws and regulations, therefore woodlands and forests not protected are open to development [9]. Economic factors, such as unemployment, impacts home prices especially in periods of economic downturn [29]. Social factors such as age demographics [30], public schools [31] and education [32] may also impact home prices. The benefit of investigating the change in housing prices on a town scale is that demographic and socio-economic data are available over multiple years of observation. The data can be used to investigate the changing pattern of spatial non-stationarity that can be modeled and mapped using GIS. Such a mapping shows the non-stationarity in urban housing and sustainability.
The main contribution of our approach is the introduction of determinants related to urban sustainability at a town scale rather than the traditional hedonic home pricing related to individual houses. Our approach incorporates a plethora of data including remote sensing, census, and other published state reports. This “big data” approach to sustainability science is needed to develop public policies on housing, conservation, transportation and employment. We believe that our approach provides a broader perspective on urban home values and enables the characterization of patterns of spatial non-stationarity across Massachusetts. There are policy implications related to development, demographics, housing and transportation. What forest or woodland to cut depends on value placed on open lands in various towns across the state? Suburban towns north and west of Boston place a higher premium on forest compared with rural western Massachusetts. Senior populations living in suburban Boston homes have seen considerable increase in their home prices and seniors choose to continue to live in these towns with greater access to senior services (free rides to malls or hospitals). A town may adopt stricter regulations to protect unprotected forests open from development, or to improve the quality of education by imposing more taxes on its citizens, knowing that school performance positively impacts home prices in these towns.
The structure of the paper is as follows: the next section examines the modeling framework and the long history of hedonic and GWR modeling. Section 2 outlines the data sources and discusses the methodology including spatial modeling considerations. Section 3 provides the results of analysis relating to spatio-temporal patterns of median home prices derived from the ordinary least squares (OLS), and two types of GWR models. We highlight the differences in the degree of spatial non-stationarity in the determinants of the median home price in Massachusetts, as well as characterize the temporal differences in median home prices in the period of bust and boom. Section 4 provides conclusions related to the theory and practice of GWR in this field.

1.2. Modeling Framework

From the methodological point of view, the hedonic price function f, typically describes the property price P as a function of three categories of independent variables: structural, locational, and environmental characteristics. Traditional hedonic approaches adopt a model structure which reduces heteroscedasticity and nonlinearity to produce a single solution for the intercept term, along with the coefficients that determine the significance of independent locational, structural, or environmental characteristics, and the overall model’s goodness-of-fit. Hedonic modeling has been executed at a variety of spatial scales ranging from a block or neighborhood [18,19,20], to the district or metro scale [21,22]. However, these models cannot account for spatial autocorrelation resulting from spatially correlated omitted variables or spatial externalities and spatial heterogeneity [15].
Spatial autocorrelation indicates that homes in a neighborhood tend to be more similar. Real estate companies (including the popular Zillow), in effect, use spatial autocorrelation to determine the price of a home at a certain location based on the prices of nearby (similar) homes. Furthermore, many homes in a neighborhood tend to be built around the same time, and proximity to both positive and negative externalities has similar effects on the market values of nearby properties [33]. However, spatial heterogeneity or spatial non-stationarity results when the relationship between two or more variables determining the median home price is not constant across space, resulting in locally varying submarkets. Ignoring spatial non-stationarity leads to misspecification in the model including missing local effects that can have profound implications for understanding the temporal and spatial relationship between housing prices, location and housing attributes [12,18].
The spatial modeling literature provides a variety of local and global models to deal with spatial dependence [34,35], as well as models that explicitly incorporate spatial heterogeneity such as the geographically weighted regression (GWR) methodology [36]. GWR calibrates a series of local regression models separately at each location and offers the ability to map local estimates of the intercept, variable coefficients, and other regression diagnostics and a check for spatial variations in the relationships between dependent and independent variables at each location [19,37]. The keys to GWR modeling are the spatial kernel functions, with fixed or varying bandwidth that impact the shape and size of local neighborhood at each location (e.g., a circular neighborhood of fixed radius or a fixed number of neighbors for each location), and the weighting functions that determine the significance of neighbors with fixed or varying bandwidth. The definition of spatial neighborhood is a critical consideration in the analysis. GWR has been applied in the contexts of determining home prices in many cities [20,38,39,40,41].
The basic GWR model assumes the same degree of spatial smoothness for each coefficient, which may not hold true in all contexts. GWR therefore overfits the data and produces a bias. Hence the basic GWR has undergone the following significant revisions: First, traditional GWR models define distances as straight line or Euclidean, while more recent modifications of the distance function adopt non-Euclidean distance metrics [39] to improve the model fit. Second, traditional GWR models use a fixed bandwidth for all variables to estimate the spatial relationship between variables while a revised GWR can use a flexible bandwidth [42] to estimate spatially varying relationships at various geographical scales within one model. The resulting model estimates coefficient surfaces that may vary at different spatial scales for different variables leading to better model fits. Third, Wheeler [43,44] proposed regularized GWR models, by combining ridge and/or lasso regression with GWR that have shown robustness in addressing the multicollinearity problem. Fourth, there has been a focus on diagnostics to check the model fit such as cross-validation (CV) score [45] to derive an optimal kernel bandwidth for GWR regression to reduce model bias. Another measure is the Akaike Information Criterion (AIC) [46] that is traditionally used to account for model parsimony dealing with the trade-off between prediction accuracy and complexity. In GWR, a corrected version of the AIC is used that accounts for sample size [47] and entails fitting bandwidths with different penalty functions. Fifth, the incorporation of temporal non-stationarity into GWR model is providing more insights on market trends and depreciation of home prices through time [18]. Finally, not all variables in GWR models exhibit non-stationarity in all contexts; hence assuming that all of the independent variables in the GWR exert a spatial influence on the dependent variable can lead to biased estimations [48]. For example, real estate markets may be economically connected through common federal policies such as governmental subsidies while some price-determining effects vary across space resulting in spatial heterogeneity. It would be wrong to assume that both factors exert a spatial influence. Research addressing this issue has led to the formulation of the Mixed Geographically Weighted Regression model (MGWR) that incorporates both linear regression and the GWR [36]. MGWR is a regression model in which the first step involves differentiating non-stationary and stationary variables by explicitly testing for spatial variability. This testing results in some independent variable coefficients being held constant, which are considered global parameters, while some others spatially vary, denoted as local parameters [49]. The second step in MGWR involves in mapping spatial variability of local parameters while global parameters have no spatial variability. Thus, MGWR differs from GWR, where all independent variables are assumed to have spatial variability resulting in different spatial distributions of local parameters.

2. Materials and Methods

2.1. Study Area

Massachusetts is located in the northeastern United States and has an area of 27,340 km2, 25.7% of which is water and 61% is covered by forest. It is the 7th-smallest state in the United States and accounts for around 2.75% of the total GDP. The state has a population of 6.8123 M (2016) of which 4.7 M live in the Boston Metro region. Boston is the 10 largest metropolitan area in the US. The population of the state is mostly urban (83%). Figure 1 shows the location of Boston and other towns mentioned throughout this paper. Metro Boston (#7 on Figure 1) is located on the eastern coast of the state.
The GWR analysis is conducted at the town and city scale. There are 351 towns and cities in Massachusetts. Each town’s data across the time period was selected for census years 2000 and 2010 [50] and ACS data for years 2009, 2011–2013 [51]. We only included 336 towns and cities due to data availability. Our data includes the Great Recession, a highly influential period for housing prices in the US economy. The period of the Great Recession, from 2008 to 2009, was characterized by loss of wealth, reduction in consumer spending, and massive job loss that resulted in a decline of home prices [29,52].

2.2. Data Sources—Socio-Economic and Environmental Variables

Our data was collected from various sources including the census, remote sensing (land cover), transportation, tax department, and labor statistics. We describe our data shown in Table 1. Our dependent variable is the Median Home Prices, recorded from the census (2000 and 2010) and ACS (2009, 2011, 2012, and 2013). Our choice of time periods was dictated by the availability of relevant data. We differentiated our temporal analysis into two sets based on the source of data—decennial census collected in 2000 and 2010, and ACS data from 2009 to 2013. (Note ACS data for 2008 is unavailable). The American Community Survey (ACS) data are estimates and therefore cannot be compared with census data directly. We used the ACS data since it represents the period of Great Recession starting in 2009. We used inflation-adjusted median home prices in order to compare today’s real estate prices to their historical norm.
The independent variables used in our models were selected based on past theoretical and empirical works examining their relevance in estimating home values. Population Density and Unemployment Rate have an impact on home prices [27,53,54]. The impact of unemployment on housing during the Great Recession [29,55] shows that employment was key for recovery in some metro areas of the US. Unemployment was a significant factor in the wake of the housing collapse in 2008 in Massachusetts and impacted home prices [56]. Therefore, we assume that unemployment is likely to impact the housing prices in our study context. The impact of education (college degree) and home price highlight that since 1988, those who completed college owned homes at higher rates in the US than those with no college education [32]. Sellers of existing homes provide a major share of the annual supply of homes sold in the US; home sales are driven by the aging of the population since seniors are net home sellers [30]. Older homeowners have emerged as a dominant segment of the housing market following the housing collapse in 2008. The homeownership rate for Americans aged 65 and over has remained at 80 percent while dropping for every other age group. Seniors typically have less mortgage debt than younger homeowners; they typically downsize and sell, or increasingly some stay at the same home [57]. Seniors play a significant role in housing dynamics and hence are included in our analysis.
Residential property tax rates vary across the state related to home values and the mix of residential and commercial holdings in each town. Suburban towns such as Weston and Wellesley have some of the lowest residential tax rates in the state while Longmeadow in western Massachusetts, with few business establishments, relies heavily on residential property taxes to fund town services. Residential taxes may therefore reflect the economic structure of each town and is included in our analysis. Prior studies identified and measured land value uplift (LVU) resulting from rapid transit and other mass transit [28,58]. Thus, we have included distances from town centroids to commuter rail stations as an independent variable.
Vehicle ownership is derived by normalizing the census variable called vehicle availability for each town by the population of each town. This data may provide information on town’s travel accessibility and modal choices. More vehicle ownership may imply more travel by car while less vehicle ownership may imply the use of public transport, walking or biking [59]. Vehicle ownership has direct implications in urban sustainability studies. The city of Boston is planning to go carbon neutral by 2050 and is seeking solutions for reducing air pollution caused by automobile transportation.
There is strong evidence to suggest that school quality substantially affects home prices in the US [60,61]. Earlier studies focused on the relationship between home prices and the quality of local education, using public school expenditures per pupil as the key school variable. However, more recent research highlights that the measure of K-12 student achievement is a more appropriate variable in home value estimations. We use the Composite Performance Index (CPI) scores of school districts in our analysis. The CPI is a measure of the extent to which students are progressing toward proficiency (a CPI of 100) in ELA (English Language Arts) and mathematics on the state’s MCAS (Massachusetts Comprehensive Assessment System) [62]. CPIs are generated separately for ELA and mathematics at all levels in Massachusetts. For this study, we considered CPI for mathematics performance in public school district. In order to estimate the MCAS score for each individual town, we used Area Interpolation Tool in ArcGIS. First, we created an interpolated surface based on CPI from all public school districts. Then we applied Massachusetts Town and City boundary shapefile to estimate the CPI for each district. The State implemented and collected CPI data starting 2003. Hence, we substituted the CPI score of 2003 for 2000 since we did not have data for 2000. Better public schools’ CPI (MCAS) generally correlates with a higher home price [31], exploited by realtors (such as Zillow) in selling homes.
In our study, we obtained the land cover classes designated as residential (low and high) and forest from the Mass Audubon’s publication called Losing Ground [63], derived from processing and classifying Landsat time series of Landsat-5, 7, and 8 data for this time period. Two thematic classes called low and high residential classes described in the Losing Ground report were summed up to define the residential class in our study. A recent study by Cunningham et al., [64] used Landsat archives to examine the change from undeveloped (forest) to developed land-use during the real estate bubble (2000–2006) and the subsequent bust (2006–2013) in Massachusetts. The results in this paper further highlight the significance of the land cover change during this period.
According to US Forest Service, Massachusetts had an estimated 3.0 million acres of forest land in 2014. About 61% of the land area of Massachusetts meets the Forest Inventory and Analysis (FIA) definition of forestland [65]. Forests are not evenly distributed across the state, and are largely influenced by development patterns. The lowest occurrences of forestland are seen in areas surrounding Metro Boston, Springfield, and Worcester, as well as along the coast and the major transportation corridors. Unprotected forest is an important determinant of home price and hence we created this new class of forest that is at the greatest risk for development. We derived this class based on two different data sources—Losing Ground [63] report and MassGIS data, as follows: First, forest areas were extracted from the Losing Ground dataset for the entire study period. Second, MassGIS dataset called Protected and Recreational Open Space layer was extracted; this class includes recreation, conservation, surface water supply protection areas and wellhead protection areas, or scenic sections [66]. Third, we created a new class called unprotected forest by differencing layers created in Step 1 and 2 to derive forest that is at risk for development. This thematic class provides us the areas at risk for conversion from forest to some form of development in the towns of Massachusetts.

2.3. Spatial Model Considerations

In this study, we use conventional ordinary least squares (OLS) as a benchmark, GWR and MGWR to describe spatial heterogeneity in housing across the towns in the state of Massachusetts.
The first step in the analysis was to use a traditional OLS of the form:
y i = a 0 + j = 1 p a j x i j + e i ,   i = 1 , 2 , ,   n
where y i is home price of each town i at a specified time of the study and x i j is a row vector of explanatory variables for town i , a j is a column vector of regression coefficients, and e i is the random error for town i . The first element of the equation is the intercept. The initial model considered the relationship of the median price of a home in each town using independent variables—Population Density, Unprotected Forest, Unemployment rate, Percent Residential Area, Vehicle Ownership, Higher Education, Senior Population, distance to the nearest Commuter Rail station, Residential Property Taxes and CPI (Composite Performance Index) for each town. OLS results were interpreted based on an assessment of multicollinearity, adjusted R2 and Akaike’s information criterion (AIC) [46]. The variance inflation factor (VIF) statistic, which measures redundancy among explanatory variables, was used to assess multicollinearity. Explanatory variables with large VIF values—above a threshold of 7.5—were considered to be multicollinear. This process ensured that the model became unbiased.
The next step in the analysis was to use GWR that explicitly incorporates the spatial structure of the variables into the estimation of the regression and shows how those estimates vary across space. We have selected an adaptive kernel whose bandwidth was found by minimizing the AIC value. The bandwidth is a count of the number of nearest observations to be included under the kernel. Preference is given to lower values of AIC since they indicate a closer fit to the data.
y i = a 0 ( u i , v i ) + j = 1 p a j ( u i , v i ) x i j + e i ,   i = 1 , 2 , ,   n
  • y i is the dependent variable at location i
  • x i j is the j th independent variable at location i
  • p is the number of independent variables
We explore the spatial variability of relationships between median home prices and the explanatory variables by mapping GWR coefficients and local R2 values. We also performed an F3 test to probe whether the GWR estimates are a significant improvement on the conventional globally estimated OLS [67]. The Akaike Information Criterion (AIC) is used in this study as a test diagnostic to select flexible bandwidth b [68,69] from:
A I C c ( b )   = 2 n   l n ( σ ^ ) + n   l n ( 2 π )   + n { n + t r ( S ) n 2 t r ( S ) }
where n is the local sample size (according to bandwidth); σ ^ is the estimated standard deviation of the error term; and t r ( S ) represents the trace of the hat matrix S . The hat matrix denotes the projection matrix from the observed y to the fitted values.
As highlighted before, GWR is not always appropriate if some of the variables do not exhibit spatial non-stationarity and can be held constant. We used a Mixed Geographically Weighted Regression model (MGWR) after testing for spatial variability of all variables [70]. In MGWR, some contributing factors that have no spatial variability will generate a global parameter, while others with spatial variability will produce a local parameter. The MGWR is defined as:
y i = a 0 ( u i , v i ) + j = 1 q a j x i j + j = q + 1 p a j ( u i , v i ) x i j + e i , i = 1 , 2 , ,   n
a j are global coefficients.
We used a Monte Carlo approach to test for significant (spatial) variation in each regression coefficient of the basic GWR against a series of randomized data sets [71]. If the true variance of the coefficient did not fall in the top 5% tail of the ranked results, the corresponding variable was treated as a global variable in the specification of MGWR.

3. Results

We discuss three sets of results. We first examine the spatio-temporal patterns of median home prices derived by comparing and validating the results of the OLS and GWR models. Our second set of results show the application of MGWR to highlight the differences in the degree of spatial non-stationarity in the determinants of the median home price in towns of Massachusetts, while the third set of results describes the urban sustainability from the perspective of economic, social and ecological determinants. The OLS results are presented in Table 2, Table 3 and Table 4; GWR results are presented in Table A1, Table A2, Table A3, Table A4, Table A5 and Table A6 and Table 5; and MGWR are presented in Table A7, Table A8, Table A9, Table A10, Table A11 and Table A12 and Table 6. Summary results for each model are presented in Table 4 (OLS), Table 5 (Basic GWR) and Table 6 (MGWR).

3.1. Spatial-Temporal Patterns of Home Prices

We first explore and model the spatiotemporal variability of median home prices and associated determinants in the state of Massachusetts by benchmarking the performance of the global regression model (OLS) with its GWR counterpart with the same set of variables. Table 2 shows the results from the OLS regression for the median home prices using decennial census data from 2000 and 2010, while Table 3 shows similar results for the ACS years—2009, 2011, 2012, and 2013. Table 4 shows overall summary results for the OLS model.
Of the covariates, only population density, unprotected forest, and residential area are not significant (at significance level 0.05) in the model in any year of observation, while others including unemployment, vehicles owned, residential taxes, educated population above 25 years, senior population, or distance to commuter rail stations, were statistically significant in all years. Unemployment rate was consistently significant from 2009 to 2013 as seen in other studies. Suburban towns exemplify this pattern the best indicating that an increase in unemployment causes a drop in median home prices in suburban towns such as Natick, Framingham, Wayland, Wellesley and Cambridge. Senior population impacted median home prices more significantly in 2013 (at significance level 0.001) and to a lesser degree from 2010 to 2012 (at significance level 0.05). Perhaps this segment of population held on to their homes that depreciated in value (in 2008) and sold their home starting 2010. The ownership of homes by seniors and its impacts on median home prices deserve further empirical scrutiny. Coefficients, associated with unemployment rate and residential taxes, are negative indicating that a decrease in home price is associated with an increase in both unemployment and residential taxes. Unprotected forest may not have overall significance in the state, but it may have more significance in the eastern part of the state around Boston where development of housing in unprotected forest has resulted in urban sprawl [9] and the building of “McMansions”. Such non-stationarity patterns in this determinant have to be explored using GWR since traditional OLS is unable to account for spatial heterogeneity in local submarkets.
Figure 2 shows the OLS residuals in the median price of homes in various towns in different time periods. In general, OLS underestimates median home price in eastern Massachusetts, which is more populated with greater housing density, and overestimates the median home price in the rural western part of the state. This pattern is consistent over the entire time period. The ratio of unprotected forest to residential area in 2013 is computed to highlight such price variations. Eastern towns around Boston have ratios ranging from 0 to 0.626, while in western and central Massachusetts, the ratios range from 5.82 to 75. Higher ranges of ratio in the west suggest that more expansive forest is not influencing median home prices in these towns, while the opposite may be true in the eastern towns. Therefore, unprotected forest may be characterized by non-stationarity and needs to be explored using GWR.
Table 4 shows the overall OLS results. The coefficient of determination (R2) across the time period ranges between 0.71 (2000) and 0.73 (2010) in the census years and from 0.64 (2009) to 0.65 (2011), 0.67 (2012), and 0.66 (2013). The standard error for the intercept is not significant in any year, except 2000. We expect that median home prices exhibit spatial heterogeneity and vary in the state from east to west reflecting proximity to Boston. As highlighted before, most of western Massachusetts is rural with large areas of forest. Hence unprotected forest may not impact overall home prices in the west but may have an impact for residential development in the east.
GWR results are presented in Table A1, Table A2, Table A3, Table A4, Table A5 and Table A6; each table pertains to one year of observation and shows the medians and ranges in the values of each coefficient across all towns (Columns 1–3). These coefficient numbers, in general, have large ranges. Hence, the results are next summarized using percentages of coefficient estimates that were positive and negative (columns 4 and 5) for each variable across all towns. GWR results for census years are shown in Table A1 and Table A2 along with p-value (F3 test) and significance. Population density is not a significant factor in both census years. The public-school CPI scores are less significant in all years. The remaining determining variables are highly significant across both census years.
The ACS survey years are shown in Table A3, Table A4, Table A5 and Table A6. Examining the impact of determinants, CPI score is significant in all ACS years (at significance level 0.001). Population density is only significant in recent years 2012–2013, probably driven by increasing urban growth in towns such as Everett, Lawrence, Malden, Arlington and Brookline (see Figure 1). Overall population increased the most in Boston, Cambridge, Somerville, Chelsea and Brookline 2011–2013. Other factors are significant across all years. To summarize, GWR results highlight that most determinants are characterized by non-stationarity leading to spatial variation in median home prices.
We next interpret columns 4 and 5 in Table A6 for 2013. Most towns (65.77%) seem to value unprotected forest (positive coefficient) while some towns (34.2%) in the rural west (negative coefficient) may not. Senior population accounts for 91.96% in the positive coefficient column, indicating that most senior citizens may have higher incomes, and better homes (longer tenure of ownership). The coefficient of this determinant is negative only in 8.04% of towns. This pattern is exemplified in towns outside of Boston Metro, such as Wales, Richmond, and Bernardston, where the total population of the town may be decreasing due to migration of all segments of population except seniors. Overall, the results across census and ACS years, show that unemployment rate, vehicles ownership, residential property taxes, residential area and distance to commuter rail stations have a larger percentage of towns with negative coefficients (column 5 in all tables). The model predicts a negative relationship between these determinants and median home prices across most towns. On the other hand, CPI, educated population above 25 years, population density, unprotected forest and senior population have a range between 60% and 100% of towns with positive coefficients (column 4 in all tables). The model predicts a positive relationship between these determinants and median home prices across most towns. The model fit in terms of positive and negative coefficient signs is consistent across most years highlighting the strong role of the selected determinants in determining the median home prices in the state.
Figure A1 maps the estimated coefficients of senior population in 2000, and 2009 to 2013 which indicate that proximity to Boston plays a strong role in differentiating eastern (colored red) and western towns (blue). Figure A2 shows coefficients of unprotected forest class for the time period are highly significant in towns northwest of Boston including Lexington, Winchester, Woburn, and Belmont. Finally, the R2 values are examined to measure the fit of GWR in each town displayed in Figure 3 for year 2013. Metro Boston has good GWR R2 values. The towns of Bernardston, East Longmeadow, Hampden, Longmeadow, Ludlow, and Springfield have higher R2 values. These towns are located in the central and western regions of the state. Bernardston (Route 91 south of NH Border) had the highest GWR R2 value. In contrast, Seekonk, Somerset, Swansea, Taunton, and Freetown had lower GWR R2 values. These towns are all located south of Boston closer to Providence RI and therefore may be impacted by the determinants in Rhode Island. Table 5 shows the overall results in terms of bandwidth, RSS, AIC and Adjusted R2. The R2 value was the highest in 2000 and averages around 0.80 over other years. To summarize, most determinants in our study are characterized by spatial non-stationarity.

3.2. Are All Determinants of Median Home Price Non-Stationary?

We next examine the question of non-stationarity of determinants using MGWR that differentiates the spatially non-stationary from stationary determinants. Results are shown in the Table A7, Table A8, Table A9, Table A10, Table A11 and Table A12 where variables with p-value (Monte Carlo simulations) in column 7 greater than 0.05 are spatially stationary and should be treated as fixed global variables (listed in the lower section of the table). Vehicle ownership and residential taxes should be treated as fixed global variables in most years. The importance of vehicle ownership is uniform across the state, as are town taxes. These should be treated as fixed global variables, while MCAS CPI scores were spatially non-stationary (local) in 2009 and 2011. Population density presents contrasting results in GWR and MGWR analysis. It is not a fixed global variable in MGWR suggesting spatial non-stationarity; population density was not significant in the GWR analysis until 2012–2013, coinciding with the increase in population in metro Boston area in this time period.
A closer examination of unprotected forest indicates that towns such as Belmont, Billerica, Carlisle, Lexington and Burlington, within the Greater Boston region, consistently perform highest on this coefficient indicating the greatest potential for residential development in these towns, if town regulations permit. On the other hand, towns such Provincetown (on the Cape), North Attenborough (on the intersection of Routes 95 and 295), and Norfolk rank the lowest for at least two years on this coefficient, indicating they have less potential for development.
Figure A3 displays unemployment coefficient maps of the GWR and MGWR that differentiate the spatial non-stationarity in this determinant around Worcester (west of Boston) and suburbs south of Boston. Similarly, results are shown for the display of senior population coefficient maps of the GWR and MGWR in Figure A4. Differences emerge in the western suburbs of Boston. Both determinants are non-stationary using GWR and MGWR but are producing different spatial patterns of coefficients around Boston. They produce similar patterns of coefficients in western part of the state.
MGWR results presented in Table 6 show better results compared with GWR results shown in Table A7 across both census and ACS years. As shown in Figure A5, GWR and MGWR have similar AIC in 2009, 2010, but display some differences in the remaining years. MGWR produces consistently the lowest AIC values for all years. GWR is the second-best performer using the AIC measure while OLS is the worst performer. This result highlights the spatial non-stationary of some factors—Population Density, Unemployment Rate, Residential Area, Vehicle Ownership, Senior Population, Distance to Commuter Rail Stations and Property Tax—that influence median home prices in Massachusetts towns. The highest difference between two model results occurs in 2013 indicating that MGWR is a better choice for this year since the nature of spatial non-stationarity of some determinants has changed.

3.3. Impact of Housing on Urban Sustainability

Economically, our study captures the Great Recession period of bust and its aftermath based on the ACS data 2009, as well as census 2010 and recovery in 2011–2013 ACS data. Determinants such as unemployment rate show differences from 2009 to 2013, while unprotected forest is significant beginning in 2009, suggesting that housing demand is impacted forests in the period of boom that followed the economic recovery. Population density is significant in recent years 2012–2013 and was likely driven by increasing employment in Metro Boston. Our spatial models capture the economic boom and bust in the various towns from 2009 to 2013 displaying the linkages between economic, social, and ecological factors in urban environments during critical periods. Addressing real estate expansion during periods of boom, one has to consider ecological and social impacts, while during periods of bust, urban planning has to consider the social implications of unemployment and drop in housing value. Social factors such as unemployment impacted certain towns in Metro Boston (such as Natick, Framingham, Wayland, Wellesley, and Somerville) during this time period. This perspective on social determinants may inform policy makers about provisioning social services and employment opportunities in these towns. Ecological factors related to unprotected forest are also significant in Metro Boston towns (such as Bedford, Lexington, Burlington, Woburn, and Waltham), suggesting these towns should balance residential expansion and forest cover in the future.

4. Discussion

Our data driven approach to model housing in this study calls for the integration of remote sensing, socio-economic, town and other data. The study covers 336 towns and cities in the state of Massachusetts for the period 2000–2013 that was characterized by periods of boom and bust in the US economy following the housing market collapse in 2008 [52]. Our modeling approach of GWR and MGWR, enables us to analyze how spatial non-stationarity of environmental, economic and social determinants lead to changes in median home prices. The basic OLS model appears least suited in modeling median home prices, while both GWR and MGWR models using adaptive bandwidths perform better.
The analysis presented here can help us address a series of questions focused on social determinants: What is the relationship between an educated workforce and home prices? Boston and Cambridge are home to sectors that attract educated workforce with high-wage jobs. Hence, people with higher education and income move closer to the suburban towns around Boston. This finding is shown in our analysis with marked differences between 2009 and 2013. What is the impact of senior citizens on home prices? Our GWR and MGWR suggests that senior citizens living in close vicinity to metro Boston exerted an influence on median home prices in 2013 and not in the period 2009–2012. Perhaps this segment of population held on to their homes that depreciated in value (in 2008) and sold their home starting 2010. The ownership of homes by seniors and its impacts on home prices deserve further empirical scrutiny. We used public school CPI scores as a determinant on home value. We found that these scores were significant using GWR model but were fixed as global variables using MGWR in 2000, 2010, 2012, and 2013. For future studies, we could incorporate other test scores (such as English Language Arts) to obtain a more nuanced understanding of the correlation between test performance across school districts in each town and their corresponding home prices.
What are the implications of spatial non-stationarity on environmental determinants? Unprotected forest cover in the eastern part of the state is more valuable than in the rural western Massachusetts. This may give towns an opportunity to build resilient communities that can preserve unprotected forest or mobilize citizen activists to transform the unprotected forest into protected forest. This investigation has implications in urban planning as it can guide us in conservation of key species. Work presented here can be extended to address reducing ecological footprints, reducing carbon emissions through reduction in deforestation, and conservation in towns that could benefit from taking actions to be resilient in the future. How does unemployment impact home prices? Unemployment rate is significant in determining home prices in all years using the GWR model. Suburban towns exemplify this pattern, indicating that an increase in unemployment causes a drop in home prices. The model predicts a negative relationship between these determinants and home prices across most towns. Property tax is another significant variable in our GWR analysis in determining home prices.
Housing impacts household behavior, policy, and the environment and therefore directly relates to urban sustainability [72]. Housing construction includes land, energy, and materials. Housing is directly connected to transportation and other externalities since people need to move back and forth from their homes for employment, recreation and other activities [73,74]. Given the complex connectivity between housing, transportation, and other economic sectors, sustainability in housing has typically been examined as a function of resource and location efficiency [75,76]. Home prices are related to environmental determinants, as described in our paper and elsewhere [77,78,79], that will be transformed given the sustainability focus of many cities, including Boston. Attributes such as green homes and solar panels may increasingly impact future home prices. Many new housing units are built with a focus on energy efficiency. Modern urban architecture emphasizes eco-friendly, or “green” homes that use materials and building methods with less energy requirements as well as result in reduced energy bills for the homeowner. The decreasing price in solar panels and their increasing availability will impact home prices in the future. Accessibility to public transport is relevant in this paper today but in the future, accessibility to electric car charging stations may become equally important. We discussed the nature and valuation of open space that is relevant in any discussion of urban sustainability as cities lower energy consumption, reduce emissions, and promote healthier lifestyles. Our approach can incorporate and estimate the non-stationarity of these sustainability driven factors that may impact future home prices. Thus, home prices can inform us of the changing public attitudes towards sustainability in cities across the US in the coming decades.

5. Conclusions

We proposed an urban sustainability framework centered on housing that could address issues involving economic, social and ecological dimensions. We used GWR model to examine spatial non-stationarity of key economic, ecological and social dimensions that are easily available at a larger spatial scale. Traditional approaches both in economics and in spatial econometrics involve hedonic and spatial models that estimate prices of individual homes using block level census data or town parcel data. Such models do not account for spatial non-stationarity at a regional or state scale. Our approach is to address urban sustainability at a broader spatial scale to encapsulate the entire state. We show that a data driven approach to modeling home prices at this broad scale calls for integration of geospatial data from a number of sources, including remote sensing, and census data to address multiple facets of sustainability. The understanding of key determinants to housing enables towns to take different strategies to sustainability such as valuation of green spaces and proximity to public transport, while unemployment rates and public-school performance can help shape urban housing policies. We hope to incorporate other variables such as green homes, residential solar panels, charging stations, air quality and other components to make our work relevant for assessing future urban sustainability.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2071-1050/10/4/1026/s1. Data used for this study (listed in Table 1) have been included as shapefile in a zipped document.

Acknowledgments

This research is based upon work supported by the National Science Foundation under Grant No. CNH #1617053, CNH-S and Boston University Initiative on Cities, 2017.

Author Contributions

Yaxiong Ma and Sucharita Gopal conceived and designed the study; Yaxiong Ma analyzed the data; Yaxiong Ma and Sucharita Gopal wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. GWR Results for Home Prices in MA Towns 2000.
Table A1. GWR Results for Home Prices in MA Towns 2000.
Median Home Price
Basic GWR Model (2000)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−212.70243.241190.7690.779.231.47 × 10−11***
Population Density−11.920.6828.7860.7139.290.5142965-
Unprotected Forest−217.6868.88388.9084.5215.480.3060173-
Unemployment Rate−189.08−17.6616.9211.3188.69<2.2 × 10−16***
Residential Area−756.43−73.39300.6935.7164.290.0004993***
Vehicle ownership−1323.53−162.52229.869.8290.18<2.2 × 10−16***
Higher Education297.391110.621985.2510006.85 × 10−13***
Senior Population−522.05226.411470.8084.8215.181.04 × 10−11***
Dist. to Stations−8.48−0.1816.7239.8860.123.59 × 10−6***
Property Tax−20.63−3.481.3810.1289.886.51 × 10−12***
CPI−1.371.386.7880.6519.350.040276*
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A2. GWR Results for Home Prices in MA Towns 2010.
Table A2. GWR Results for Home Prices in MA Towns 2010.
Median Home Price
Basic GWR Model (2010)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−130.54347.021501.3397.922.083.07 × 10−2*
Population Density−11.222.3025.1578.2721.730.531882-
Unprotected Forest−61.80100.84745.1777.0822.921.90 × 10−13***
Unemployment Rate−114.97−15.333.102.6897.32<2.2 × 10−16***
Residential Area−781.96−107.39246.6433.6366.374.49 × 10−9***
Vehicle ownership−872.92−274.5364.772.0897.921.57 × 10−7***
Higher Education−2067.55626.962327.9988.9911.01<2.2 × 10−16***
Senior Population−409.73660.412574.7183.6316.37<2.2 × 10−16***
Dist. to Stations−18.26−2.003.392575<2.2 × 10−16***
Property Tax−29.32−7.032.514.1795.833.94 × 10−3**
CPI−3.723.5917.9887.512.50.043231*
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A3. GWR Results for Home Prices in MA Towns 2009.
Table A3. GWR Results for Home Prices in MA Towns 2009.
Median Home Price
Basic GWR Model (2009)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−1507.61−34.44433.6944.9455.063.91 × 10−5***
Population Density−10.475.0730.8889.8810.120.2161-
Unprotected Forest−59.97141.74812.4493.456.55<2.2 × 10−16***
Unemployment Rate−29.38−3.1417.1221.4378.577.21 × 10−8***
Residential Area−842.60−139.32242.0515.1884.824.40 × 10−5***
Vehicle ownership−1018.52−300.08126.042.9897.025.52 × 10−7***
Higher Education−545.121083.322578.5397.922.08<2.2 × 10−16***
Senior Population−625.56537.452741.0782.1417.86<2.2 × 10−16***
Dist. to Stations−12.16−2.293.1626.1973.811.90 × 10−5***
Property Tax−37.52−7.638.017.4492.567.09 × 10−8***
CPI0.095.4728.541000<2.2 × 10−16***
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A4. GWR Results for Home Prices in MA Towns 2011.
Table A4. GWR Results for Home Prices in MA Towns 2011.
Median Home Price
Basic GWR Model (2011)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−1213.561.45339.4950.649.45.81 × 10−1-
Population Density−7.346.9331.5992.867.140.2586052-
Unprotected Forest−154.5697.60857.4679.7620.24<2.2 × 10−16***
Unemployment Rate−35.99−7.292.881.7998.211.72 × 10−8***
Residential Area−950.93−169.73147.899.2390.771.03 × 10−3**
Vehicle ownership−860.20−253.1599.182.3897.621.93 × 10−3**
Higher Education−1494.37800.342474.5394.945.06<2.2 × 10−16***
Senior Population−567.12539.882330.5287.512.5<2.2 × 10−16***
Dist. to Stations−9.94−1.143.5724.475.61.58 × 10−1-
Property Tax−25.86−5.603.563.8796.139.45 × 10−4***
CPI−0.224.1225.7099.40.61.20 × 10−13***
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A5. GWR Results for Home Prices in MA Towns 2012.
Table A5. GWR Results for Home Prices in MA Towns 2012.
Median Home Price
Basic GWR Model (2012)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−918.93−27.37341.9646.1353.874.86 × 10−1-
Population Density−8.037.7340.8192.567.440.005703**
Unprotected Forest−191.2294.72770.027525< 2.2 × 10−16***
Unemployment Rate−40.78−7.784.254.4695.542.99 × 10−16***
Residential Area−1178.58−203.18144.679.2390.773.64 × 10−6***
Vehicle ownership−742.86−222.6993.1913.3986.616.33 × 10−5***
Higher Education−1489.48837.342527.9594.355.65< 2.2 × 10−16***
Senior Population−477.13481.552309.4787.512.55.66 × 10−15***
Dist. to Stations−11.23−0.9410.8229.1770.839.25 × 10−5***
Property Tax−26.11−5.165.013.5796.438.22 × 10−9***
CPI0.764.6223.771000< 2.2 × 10−16***
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A6. GWR Results for Home Prices in MA Towns 2013.
Table A6. GWR Results for Home Prices in MA Towns 2013.
Median Home Price
Basic GWR Model (2013)
VariablesMinimumMedianMax% of Positive% of Negativep-Value (F3)Sig.
Intercept−748.02−159.96279.5630.9569.059.71 × 10−1-
Population Density−7.817.6041.4893.156.851.43 × 10−6***
Unprotected Forest−63.64120.01768.5365.7734.23<2.2 × 10−16***
Unemployment Rate−32.81−6.942.926.5593.453.50 × 10−12***
Residential Area−937.70−202.99148.4413.186.99.43 × 10−6***
Vehicle ownership−796.82−259.97519.415.6594.351.62 × 10−8***
Higher Education−537.89843.471975.5997.622.38<2.2 × 10−16***
Senior Population−408.29656.822340.8091.968.04<2.2 × 10−16***
Dist. to Stations−10.98−0.9510.0525758.52 × 10−6***
Property Tax−21.75−5.721.918.3391.677.10 × 10−9***
CPI1.056.4220.6610002.01 × 10−5***
The first 3 columns refer to range of coefficients and median values of the intercepts and 10 independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the F3 and the significance refer to the coefficients (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A7. MGWR Results for Home Prices in MA Towns 2000.
Table A7. MGWR Results for Home Prices in MA Towns 2000.
Median Home Price
Mixed GWR Model (2000)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Intercept−47.18212.411440.5694.645.360.04.
Population Density−8.631.2625.8563.3936.610.00***
Unemployment Rate−244.13−24.5513.4923.2176.790.00***
Residential Area−730.09−98.82279.0434.2365.770.00***
Vehicle ownership−1352.88−208.48304.0211.6188.390.00***
Senior Population−528.56179.111395.8577.3822.620.00***
Dist. to Stations−8.04−0.2714.4334.2365.770.00***
Property Tax−17.89−3.250.8410.4289.580.04.
Global Variables-------
Unprotected Forest---72.49-0.48-
Higher Education---1063.50-0.11-
CPI---1.08-0.32-
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A8. MGWR Results for Home Prices in MA Towns 2009.
Table A8. MGWR Results for Home Prices in MA Towns 2009.
Median Home Price
Mixed GWR Model (2009)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Population Density1.686.9214.9810000.01.
Unprotected Forest−18.47136.56877.8197.022.980.01.
Unemployment Rate−32.56−4.525.249.5290.480.04.
Higher Education−46.091027.442162.8599.70.30.01.
Senior Population−717.27543.892525.4884.8215.180.00***
Dist. to Stations−17.40−2.754.6126.7973.210.00***
Property Tax1.295.199.5910000.01.
Global Variables-------
Intercept---−22.98-0.23-
Residential Area---−192.74-0.08-
Vehicle ownership---−263.5063-0.21-
Property Tax---−6.2983-0.19-
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A9. MGWR Results for Home Prices in MA Towns 2010.
Table A9. MGWR Results for Home Prices in MA Towns 2010.
Median Home Price
Mixed GWR Model (2010)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Population Density−10.613.0123.3775.624.40.00***
Unprotected Forest−117.2481.151128.2982.1417.860.02.
Unemployment Rate−78.05−17.19−3.4301000.00***
Residential Area−693.29−89.47238.3827.3872.620.01.
Higher Education−50.24901.481603.1099.70.30.00***
Senior Population−676.64641.033500.1077.6822.320.00***
Dist. to Stations−20.46−1.464.6120.5479.460.00***
Global Variables-------
Intercept---273.3707-0.45-
Vehicle ownership---−243.0027-0.13-
Property Tax---−4.75-0.51-
CPI---273.37-0.53-
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A10. MGWR Results for Home Prices in MA Towns 2011.
Table A10. MGWR Results for Home Prices in MA Towns 2011.
Median Home Price
Mixed GWR Model (2011)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Population Density−3.657.9131.4995.544.460.01.
Unprotected Forest−283.68121.03783.8891.378.630.00***
Residential Area−970.74−192.49132.066.5593.450.00***
Higher Education−150.67936.712298.4899.110.890.00***
Senior Population−456.07496.822273.6788.111.90.00***
Dist. to Stations−13.77−1.913.6618.1581.850.00***
CPI−1.254.939.9597.622.380.00***
Global Variables-------
Intercept---−18.07---
Unemployment Rate---−7.38---
Vehicle ownership---−183.00---
Property Tax---−6.56---
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A11. MGWR Results for Home Prices in MA Towns 2012.
Table A11. MGWR Results for Home Prices in MA Towns 2012.
Median Home Price
Mixed GWR Model (2012)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Population Density−7.009.3642.0097.622.380.00***
Unprotected Forest0.7785.71712.421001000.00***
Residential Area−1101.51−262.61198.375.6594.350.00***
Higher Education−183.99927.831689.2999.70.30.00***
Senior Population−474.79444.502876.7284.2315.770.00***
Dist. to Stations−17.50−1.104.1720.8379.170.00***
Global Variables-------
Intercept---7.78-0.77-
Unemployment Rate---−7.8404-0.06-
Vehicle ownership---−196.24-0.64-
Property Tax---−6.95-0.06-
CPI---4.7264-0.05-
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Table A12. MGWR Results for Home Prices in MA Towns 2013.
Table A12. MGWR Results for Home Prices in MA Towns 2013.
Median Home Price
Mixed GWR Model (2013)
Local Variables
VariablesMinimumMedianMax% of Positive% of Negativep-Value (MC)Sig.
Population Density−5.258.7532.6594.945.060.00***
Unprotected Forest−71.06100.28960.4264.5835.420.00***
Unemployment Rate−32.11−7.660.852.6897.320.02.
Residential Area−867.56−232.16119.826.2593.750.02.
Higher Education73.76800.331532.8310000.01.
Senior Population−606.89584.402809.6586.913.10.00***
Dist. to Stations−11.99−1.688.0317.2682.740.00***
Property Tax−19.62−7.011.375.0694.940.03.
Global Variables-------
Intercept---−105.84-0.90-
Vehicle ownership---−206.13-0.24-
CPI---6.03-0.37-
The first 3 columns refer to range of coefficients and median values of the intercepts and non-stationary independent variables. The % positive and % negative refer to the sign of coefficients. p-Values of the Monte Carlo test and the significance refer to the coefficients (critical value 0.05). Global variables are stationary are shown in the lower panel of the table along with their coefficient estimates. (Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’).
Figure A1. Coefficient Estimates of Basic GWR for Senior Population: (A) 2000; (B) 2009; (C) 2010; (D) 2011; (E) 2012; (F) 2013.
Figure A1. Coefficient Estimates of Basic GWR for Senior Population: (A) 2000; (B) 2009; (C) 2010; (D) 2011; (E) 2012; (F) 2013.
Sustainability 10 01026 g0a1
Figure A2. Coefficient Estimates of Basic GWR for Unprotected Forest: (A) 2000; (B) 2009; (C) 2010; (D) 2011; (E) 2012; (F) 2013.
Figure A2. Coefficient Estimates of Basic GWR for Unprotected Forest: (A) 2000; (B) 2009; (C) 2010; (D) 2011; (E) 2012; (F) 2013.
Sustainability 10 01026 g0a2
Figure A3. Unemployment Coefficient for 2013: (A) Basic GWR; (B) MGWR.
Figure A3. Unemployment Coefficient for 2013: (A) Basic GWR; (B) MGWR.
Sustainability 10 01026 g0a3
Figure A4. Senior Population Coefficient for 2013: (A) Basic GWR; (B) MGWR.
Figure A4. Senior Population Coefficient for 2013: (A) Basic GWR; (B) MGWR.
Sustainability 10 01026 g0a4
Figure A5. Senior Population Coefficient for 2013: OLS (Blue); Basic GWR (Brown); MGWR (Gray).
Figure A5. Senior Population Coefficient for 2013: OLS (Blue); Basic GWR (Brown); MGWR (Gray).
Sustainability 10 01026 g0a5

References

  1. Pickett, S.; Burch, W., Jr.; Dalton, S.; Foresman, T.; Grove, J.; Rowntree, R. A conceptual framework for the study of human ecosystems in urban areas. Urban Ecosyst. 1997, 1, 185–199. [Google Scholar] [CrossRef]
  2. Grimm, N.; Grove, J.M.; Pickett, S.; Redman, C. Integrated Approaches to Long-Term Studies of Urban Ecological Systems. BioScience 2000, 50, 571–584. [Google Scholar] [CrossRef]
  3. Alberti, M. Modeling the urban ecosystem: A conceptual framework. Environ. Plan. B 1999, 26, 605–630. [Google Scholar] [CrossRef]
  4. Ostrom, E. A General Framework for Analyzing Sustainability of Social-Ecological Systems. Science 2009, 325, 419–422. [Google Scholar] [CrossRef] [PubMed]
  5. Newman, P.; Kenworthy, J. Sustainability and Cities; Island Press: Washington, DC, USA, 1999. [Google Scholar]
  6. Shiller, R. Understanding Recent Trends in House Prices and Home Ownership; Cowles Foundation for Research in Economics; Yale University: New Haven, CT, USA, 2007. [Google Scholar]
  7. Harris, D. “Property Values Drop When Blacks Move in, Because...”: Racial and Socioeconomic Determinants of Neighborhood Desirability. Am. Sociol. Rev. 1999, 64, 461–479. [Google Scholar] [CrossRef]
  8. Case, K.; Mayer, C. Housing price dynamics within a metropolitan area. Reg. Sci. Urban Econ. 1996, 26, 387–407. [Google Scholar] [CrossRef]
  9. Wolch, J.; Byrne, J.; Newell, J. Urban green space, public health, and environmental justice: The challenge of making cities ‘just green enough’. Landsc. Urban Plan. 2014, 125, 234–244. [Google Scholar] [CrossRef]
  10. Alkon, A.; Agyeman, J. Cultivating Food Justice; MIT Press: Cambridge, MA, USA, 2014. [Google Scholar]
  11. Cutts, B.; Darby, K.; Boone, C.; Brewis, A. City structure, obesity, and environmental justice: An integrated analysis of physical and social barriers to walkable streets and park access. Soc. Sci. Med. 2009, 69, 1314–1322. [Google Scholar] [CrossRef] [PubMed]
  12. Wu, B.; Li, R.; Huang, B. A geographically and temporally weighted autoregressive model with application to housing prices. Int. J. Geogr. Inf. Sci. 2014, 28, 1186–1204. [Google Scholar] [CrossRef]
  13. Wang, S.; Fang, C.; Wang, Y.; Huang, Y.; Ma, H. Quantifying the relationship between urban development intensity and carbon dioxide emissions using a panel data analysis. Ecol. Indic. 2015, 49, 121–131. [Google Scholar] [CrossRef]
  14. Wachsmuth, D.; Cohen, D.; Angelo, H. Expand the frontiers of urban sustainability. Nature 2016, 536, 391–393. [Google Scholar] [CrossRef] [PubMed]
  15. Goodman, A.; Thibodeau, T. Housing market segmentation and hedonic prediction accuracy. J. Hous. Econ. 2003, 12, 181–201. [Google Scholar] [CrossRef]
  16. LeSage, J. An Introduction to Spatial Econometrics. Revue D’économie Industrielle 2008, 123, 19–44. [Google Scholar] [CrossRef]
  17. Anselin, L. Spatial Econometrics; Springer: Dordrecht, The Netherlands, 2010. [Google Scholar]
  18. Huang, B.; Wu, B.; Barry, M. Geographically and temporally weighted regression for modeling spatio-temporal variation in house prices. Int. J. Geogr. Inf. Sci. 2010, 24, 383–401. [Google Scholar] [CrossRef]
  19. Bitter, C.; Mulligan, G.; Dall’erba, S. Incorporating spatial variation in housing attribute prices: A comparison of geographically weighted regression and the spatial expansion method. J. Geogr. Syst. 2006, 9, 7–27. [Google Scholar] [CrossRef] [Green Version]
  20. Helbich, M.; Brunauer, W.; Vaz, E.; Nijkamp, P. Spatial Heterogeneity in Hedonic House Price Models: The Case of Austria. Urban Stud. 2013, 51, 390–411. [Google Scholar] [CrossRef]
  21. Paulsen, K. The Effects of Growth Management on the Spatial Extent of Urban Development, Revisited. Land Econ. 2013, 89, 193–210. [Google Scholar] [CrossRef]
  22. Hasse, J.; Lathrop, R. A Housing-Unit-Level Approach to Characterizing Residential Sprawl. Photogramm. Eng. Remote Sens. 2003, 69, 1021–1030. [Google Scholar] [CrossRef]
  23. Collins, J.; Woodcock, C. An assessment of several linear change detection techniques for mapping forest mortality using multitemporal landsat TM data. Remote Sens. Environ. 1996, 56, 66–77. [Google Scholar] [CrossRef]
  24. Jeon, S.; Olofsson, P.; Woodcock, C. Land use change in New England: A reversal of the forest transition. J. Land Use Sci. 2013, 9, 105–130. [Google Scholar] [CrossRef]
  25. Kennedy, R.; Yang, Z.; Cohen, W. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 1. LandTrendr—Temporal segmentation algorithms. Remote Sens. Environ. 2010, 114, 2897–2910. [Google Scholar] [CrossRef]
  26. Maliene, V.; Malys, N. High-quality housing—A key issue in delivering sustainable communities. Build. Environ. 2009, 44, 426–430. [Google Scholar] [CrossRef]
  27. Case, K.; Shiller, R. Is There a Bubble in the Housing Market? Brook. Pap. Econ. Act. 2003, 2003, 299–362. [Google Scholar] [CrossRef]
  28. Higgins, C.; Kanaroglou, P. Forty years of modelling rapid transit’s land value uplift in North America: Moving beyond the tip of the iceberg. Transp. Rev. 2016, 36, 610–634. [Google Scholar] [CrossRef]
  29. Goodman, C.; Mance, S. Employment loss and the 2007–09 recession: An overview. Mon. Labor Rev. 2011, 134, 3–12. [Google Scholar]
  30. Myers, D.; Ryu, S. Aging Baby Boomers and the Generational Housing Bubble: Foresight and Mitigation of an Epic Transition. J. Am. Plan. Assoc. 2008, 74, 17–33. [Google Scholar] [CrossRef]
  31. Gibbons, S.; Machin, S. Valuing school quality, better transport, and lower crime: Evidence from house prices. Oxf. Rev. Econ. Policy 2008, 24, 99–119. [Google Scholar] [CrossRef]
  32. Aughinbaugh, A. Patterns of Homeownership, Delinquency, and Foreclosure among Youngest Baby Boomers; Bureau of Labor Statistics: Washington, DC, USA, 2013.
  33. Fotheringham, A.; Crespo, R.; Yao, J. Geographical and Temporal Weighted Regression (GTWR). Geogr. Anal. 2015, 47, 431–452. [Google Scholar] [CrossRef]
  34. Anselin, L. GIS Research Infrastructure for Spatial Analysis of Real Estate Markets. J. Hous. Res. 1998, 9, 113–133. [Google Scholar]
  35. Pace, R.; LeSage, J.; Zhu, S. Impact of Cliff and Ord on the Housing and Real Estate Literature. Geogr. Anal. 2009, 41, 418–424. [Google Scholar] [CrossRef]
  36. Fotheringham, A.; Brundson, C.; Charlton, M. Geographically Weighted Regression; Wiley: Chichester, West Sussex, UK, 2002. [Google Scholar]
  37. Yu, D.; Wei, Y.; Wu, C. Modeling Spatial Dimensions of Housing Prices in Milwaukee, WI. Environ. Plan. B 2007, 34, 1085–1102. [Google Scholar] [CrossRef]
  38. Crespo, R.; Fotheringham, S.; Charlton, M. Application of Geographically Weighted Regression to a 19-Year Set of House Price Data in London to Calibrate Local Hedonic Price Models; National University of Ireland Maynooth: Maynooth, Ireland, 2007. [Google Scholar]
  39. Lu, B.; Charlton, M.; Harris, P.; Fotheringham, A. Geographically weighted regression with a non-Euclidean distance metric: A case study using hedonic house price data. Int. J. Geogr. Inf. Sci. 2014, 28, 660–681. [Google Scholar] [CrossRef]
  40. Demšar, U.; Fotheringham, A.; Charlton, M. Exploring the spatio-temporal dynamics of geographical processes with geographically weighted regression and geovisual analytics. Inf. Vis. 2008, 7, 181–197. [Google Scholar] [CrossRef]
  41. Saphores, J.; Li, W. Estimating the value of urban green areas: A hedonic pricing analysis of the single family housing market in Los Angeles, CA. Landsc. Urban Plan. 2012, 104, 373–387. [Google Scholar] [CrossRef]
  42. Yang, W. An Extension of Geographically Weighted Regression with Flexible Bandwidths. Ph.D. Thesis, University of St Andrews, St Andrews, Fife, Scotland, UK, 2014. [Google Scholar]
  43. Wheeler, D. Diagnostic Tools and a Remedial Method for Collinearity in Geographically Weighted Regression. Environ. Plan. A 2007, 39, 2464–2481. [Google Scholar] [CrossRef]
  44. Wheeler, D. Visualizing and Diagnosing Output from Geographically Weighted Regression Models; Department of Biostatistics, Rollins School of Public Health, Emory University: Atlanta, GA, USA, 2008. [Google Scholar]
  45. Bowman, A. An Alternative Method of Cross-Validation for the Smoothing of Density Estimates. Biometrika 1984, 71, 353–360. [Google Scholar] [CrossRef]
  46. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  47. Hurvich, C.; Simonoff, J.; Tsai, C. Smoothing parameter selection in nonparametric regression using an improved Akaike information criterion. J. R. Stat. Soc. Ser. B 1998, 60, 271–293. [Google Scholar] [CrossRef]
  48. Wei, C.; Qi, F. On the estimation and testing of mixed geographically weighted regression models. Econ. Model. 2012, 29, 2615–2620. [Google Scholar] [CrossRef]
  49. Brunsdon, C.; Fotheringham, S.; Charlton, M. Geographically Weighted Regression as a Statistical Model; University of Newcastle: Callaghan, Australia, 2000. [Google Scholar]
  50. Bureau, U. Decennial Census Data. Available online: https://www.census.gov/programs-surveys/decennial-census/data.html (accessed on 12 November 2017).
  51. Bureau, U. American Community Survey (ACS). Available online: https://www.census.gov/programs-surveys/acs/ (accessed on 12 November 2017).
  52. Elsby, M.; Hobijn, B.; Sahin, A. The Labor Market in the Great Recession; The Brookings Institution: Washington, DC, USA, 2010. [Google Scholar]
  53. Shiller, R. Long-Term Perspectives on the Current Boom in Home Prices. Econ. Voice 2006, 3, 1–11. [Google Scholar] [CrossRef]
  54. Rogers, W.; Winkler, A. The relationship between the housing and labor market crises and doubling up: An MSA-level analysis, 2005–2011. Mon. Labor Rev. 2013, 1–25. [Google Scholar] [CrossRef]
  55. Byun, K. The U.S. housing bubble and bust: Impacts on employment. Mon. Labor Rev. 2010, 3, 3–17. [Google Scholar]
  56. Holly, S.; Pesaran, M.; Yamagata, T. A spatio-temporal model of house prices in the USA. J. Econ. 2010, 158, 160–173. [Google Scholar] [CrossRef]
  57. Projections and Implications for Housing a Growing Population: Older Adults 2015–2035 | Joint Center for Housing Studies, Harvard University. Available online: http://www.jchs.harvard.edu/housing-a-growing-population-older-adults (accessed on 12 November 2017).
  58. Mulley, C. Accessibility and Residential Land Value Uplift: Identifying Spatial Variations in the Accessibility Impacts of a Bus Transitway. Urban Stud. 2013, 51, 1707–1724. [Google Scholar] [CrossRef]
  59. Rodrigue, J. The Geography of Transport Systems, 4th ed.; Routledge: New York, NY, USA, 2006. [Google Scholar]
  60. Clapp, J.; Nanda, A.; Ross, S. Which school attributes matter? The influence of school district performance and demographic composition on property values. J. Urban Econ. 2008, 63, 451–466. [Google Scholar] [CrossRef]
  61. Brasington, D.; Haurin, D. Educational Outcomes and House Values: A Test of the value added Approach. J. Reg. Sci. 2006, 46, 245–268. [Google Scholar] [CrossRef]
  62. Mass. ESE. 2010 Glossary of AYP Reporting Terms. Available online: http://profiles.doe.mass.edu/ayp/ayp_report/glossary2010.html#cpi (accessed on 12 November 2017).
  63. Mass Audubon. Losing Ground. Planning for Resilience, 5th ed.; Mass Audubon: Lincoln, MA, USA, 2014. [Google Scholar]
  64. Cunningham, S.; Rogan, J.; Martin, D.; DeLauer, V.; McCauley, S.; Shatz, A. Mapping land development through periods of economic bubble and bust in Massachusetts using Landsat time series data. GISci. Remote Sens. 2015, 52, 397–415. [Google Scholar] [CrossRef]
  65. Butler, B. Forests of Massachusetts, 2015. Resource Update FS-89; U.S. Department of Agriculture, Forest Service, Northern Research Station: Newtown Square, PA, USA, 2016.
  66. MassGIS (Bureau of Geographic Information). Available online: http://www.mass.gov/anf/research-and-tech/it-serv-and-support/application-serv/office-of-geographic-information-massgis (accessed on 12 November 2017).
  67. Leung, Y.; Mei, C.; Zhang, W. Statistical Tests for Spatial Nonstationarity Based on the Geographically Weighted Regression Model. Environ. Plan. A 2000, 32, 9–32. [Google Scholar] [CrossRef]
  68. Gollini, I.; Lu, B.; Charlton, M.; Brunsdon, C.; Harris, P. GWmodel: AnRPackage for Exploring Spatial Heterogeneity Using Geographically Weighted Models. J. Stat. Softw. 2015, 63. [Google Scholar] [CrossRef]
  69. Lu, B.; Harris, P.; Charlton, M.; Brunsdon, C. The GWmodel R package: Further topics for exploring spatial heterogeneity using geographically weighted models. Geo-Spat. Inf. Sci. 2014, 17, 85–101. [Google Scholar] [CrossRef]
  70. Mei, C.; Wang, N.; Zhang, W. Testing the Importance of the Explanatory Variables in a Mixed Geographically Weighted Regression Model. Environ. Plan. A 2006, 38, 587–598. [Google Scholar] [CrossRef]
  71. Lu, B.; Harris, P.; Gollini, I.; Charlton, M.; Brunsdon, C. GWmodel: An R package for exploring spatial heterogeneity. In Proceedings of the GISRUK 2013, Liverpool, UK, 3–5 April 2013. [Google Scholar]
  72. Saiz, A. The Geographic Determinants of Housing Supply. Q. J. Econ. 2010, 125, 1253–1296. [Google Scholar] [CrossRef]
  73. Mills, E. An Aggregative Model of Resource Allocation in a Metropolitan Area. Am. Econ. Rev. 1967, 57, 197–210. [Google Scholar]
  74. Muth, R. Cities and Housing; University of Chicago Press: Chicago, IL, USA, 1975. [Google Scholar]
  75. Koebel, C.; McCoy, A.; Sanderford, A.; Franck, C.; Keefe, M. Diffusion of green building technologies in new housing construction. Energy Build. 2015, 97, 175–185. [Google Scholar] [CrossRef]
  76. Sanderford, A.; McCoy, A.; Keefe, M. Adoption of Energy Star certifications: Theory and evidence compared. Build. Res. Inf. 2017, 46, 207–219. [Google Scholar] [CrossRef]
  77. Tu, C.; Eppli, M. An Empirical Examination of Traditional Neighborhood Development. Real Estate Econ. 2001, 29, 485–501. [Google Scholar] [CrossRef]
  78. Rauterkus, S.; Thrall, G.; Hangen, E. Location Efficiency and Mortgage Default. J. Sustain. Real Estate 2010, 2, 117–141. [Google Scholar]
  79. Tsatsaronis, K.; Zhu, H. What Drives Housing Price Dynamics: Cross-Country Evidence. BIS Q. Rev. 2004, 65–78. [Google Scholar]
Figure 1. Location of Boston and Other Towns in Massachusetts mentioned in the paper.
Figure 1. Location of Boston and Other Towns in Massachusetts mentioned in the paper.
Sustainability 10 01026 g001
Figure 2. OLS Residuals of 2000 and 2009–2013.
Figure 2. OLS Residuals of 2000 and 2009–2013.
Sustainability 10 01026 g002
Figure 3. Local R2 of Basic GWR for 2013.
Figure 3. Local R2 of Basic GWR for 2013.
Sustainability 10 01026 g003
Table 1. Variables used in the study.
Table 1. Variables used in the study.
VariablesDescriptionSource
Median Home PriceMedian home value in thousand dollars (adjusted for inflation)Census, ACS
Population DensityPopulation density (number of people per hectare)Census, ACS
Unprotected ForestPercent coverage of unprotected forest in each townLandsat
Unemployment RatePercent of unemployed people in each townMass. Labor and Workforce Development
Residential AreaPercent coverage of residential areasLandsat
Vehicle ownershipNumber of vehicles per capita Census, ACS
Higher EducationPercent of people have bachelor’s or higher degree above the age of 25Census, ACS
Senior PopulationPercent of senior population (over age 65)Census, ACS
Distance to Commuter Rail Sta.Distance from town centroid to nearest Commuter Rail StationMassGIS, MBTA
Residential Property TaxAmount per $1000 assessed home priceMass. Department of Revenue
Composite Performance IndexStudents’ performance on Mathematics Mass. Department of Elementary and Secondary Education
Table 2. Ordinary Least Squares (OLS) Results for Home Prices in MA Towns 2000 and 2010.
Table 2. Ordinary Least Squares (OLS) Results for Home Prices in MA Towns 2000 and 2010.
Median Home Price
OLS Model (2000)OLS Model (2010)
VariablesCoeff.t-Valuep-ValueSig.VIFCoeff.t-Valuep-ValueSig.VIF
Intercept289.693.53 4.77 × 10−4***-252.911.660.099.-
Population Density0.790.930.353-4.010.040.040.971-3.40
Unprotected Forest30.660.730.468-4.00−52.01−0.910.363-3.94
Unemployment Rate−7.31−1.940.053.1.53−8.69−3.281.15 × 10−3**1.43
Residential Area−63.12−1.350.179-6.61−65.83−1.100.272-5.95
Vehicle ownership−390.28−5.135.11 × 10−7***2.34−427.09−5.918.45 × 10−9***2.09
Higher Education1634.2214.75<2 × 10−16***2.261527.8710.50<2 × 10−16***1.98
Senior Population162.681.360.176-1.90335.002.200.028*2.09
Dist. to Stations−0.60−4.0197.26 × 10−5***2.01−0.95−4.4959.69 × 10−6***2.26
Property Tax−7.61−5.643.80 × 10−8***1.35−13.03−6.522.73 × 10−10***1.29
CPI1.953.211.48 × 10−3**2.245.143.564.28 × 10−4***2.44
(Signif. Codes: 0 ‘***’, 0.001 ’**’, 0.01 ‘*’, 0.05 ‘.’). The coefficients, intercept and 10 independent variables along with t-values, p-values, significance, VIF values are listed. VIF is a metric for testing multicollinearity (threshold 7.5).
Table 3. OLS Results for Home Prices in MA Towns 2009, 2011–2013.
Table 3. OLS Results for Home Prices in MA Towns 2009, 2011–2013.
Median Home Price
OLS Model (2009)OLS Model (2011)
VariablesCoeff.t-Valuep-ValueSig.VIFCoeff.t-Valuep-ValueSig.VIF
Intercept129.380.930.351--220.081.5440.124--
Population Density0.690.640.520-3.541.271.3350.183-3.51
Unprotected Forest16.390.300.767-3.735.810.1190.905-3.65
Unemployment Rate−7.47−2.864.54 × 10−3**1.35−9.60−4.0067.65 × 10−5***1.59
Residential Area−17.39−0.300.766-5.72−70.84−1.3510.177-5.75
Vehicle ownership−394.14−4.822.21 × 10−6***2.25−363.62−4.9111.44 × 10−6***2.31
Higher Education1520.5910.85< 2 × 10−16***1.881326.719.754< 2 × 10−16***2.21
Senior Population237.891.540.125-2.07355.342.5270.012*2.14
Dist. to Stations−0.90−4.292.33 × 10−5***2.20−1.06−5.6962.75 × 10−8***2.20
Property Tax−13.68−6.319.17 × 10−10***1.35−10.27−5.7791.76 × 10−8***1.41
CPI5.684.272.57 × 10−5***2.304.383.2381.33 × 10−3**2.67
OLS Model (2012)OLS Model (2013)
VariablesCoeff.t-Valuep-ValueSig.VIFCoeff.t-Valuep-ValueSig.VIF
Intercept198.771.660.097.-101.550.750.457--
Population Density1.321.550.122-3.370.981.150.251-3.33
Unprotected Forest9.660.210.831-3.630.710.010.988-3.67
Unemployment Rate−10.95−4.999.66 × 10−7***1.55−10.19−4.851.88 × 10−6***1.57
Residential Area−85.82−1.760.079.5.73−77.12−1.580.116-5.68
Vehicle ownership−391.33−5.811.50 × 10−8***2.35−449.96−6.454.04 × 10−10***2.51
Higher Education1234.9910.05< 2 × 10−16***1.991157.219.32< 2 × 10−16***2.05
Senior Population403.833.121.98 × 10−3**2.11487.123.871.31 × 10−4***2.06
Dist. to Stations−1.02−6.093.23 × 10−9***2.08−1.00−6.054.05 × 10−9***2.01
Property Tax−10.31−6.601.64 × 10−10***1.39−9.88−6.463.81 × 10−10***1.41
CPI4.944.538.45 × 10−6***2.096.244.842.05 × 10−6***2.26
(Signif. Codes: 0 ‘***’, 0.001 ‘**’, 0.01 ‘*’, 0.05 ‘.’). The coefficients, intercept and 10 independent variables along with t-values, p-values, significance, VIF values are listed. VIF is a metric for testing multicollinearity (threshold 7.5).
Table 4. OLS Diagnostics for Home Prices in MA Towns.
Table 4. OLS Diagnostics for Home Prices in MA Towns.
Decennial Census YearsACS Years
Year200020102009201120122013
RSS1,609,7452,946,5442,918,6092,332,9151,993,7942,033,210
AIC3825.924029.054025.853950.593897.813904.38
Adjusted R20.710.730.640.650.670.66
RSS, AIC and adjusted R2, for decennial years (2000 and 2010), ACS years (2009, 2011–2013).
Table 5. Basic Geographically Weighted Regression (GWR) Diagnostics for Home Prices in MA Towns.
Table 5. Basic Geographically Weighted Regression (GWR) Diagnostics for Home Prices in MA Towns.
Decennial Census YearsACS Years
Year200020102009201120122013
Bandwidth849898917790
RSS531,883.1111,82761,166,110951,077.8633,734.2786,608.1
AIC3684.433888.253904.783858.533778.993798.11
Adjusted R20.860.810.800.790.840.81
Bandwidth, RSS, AIC and adjusted R2, for decennial years (2000 and 2010), ACS years (2009, 2011–2013).
Table 6. Mixed GWR Diagnostics for Home Prices in MA Towns.
Table 6. Mixed GWR Diagnostics for Home Prices in MA Towns.
Decennial Census YearsACS Years
Year200020102009201120122013
Bandwidth849898847790
RSS635,7971,402,2131,479,2421,188,077944,864922,170
AIC366638833905384037653778
Bandwidth, RSS and AIC for decennial years (2000 and 2010), ACS years (2009, 2011–2013).

Share and Cite

MDPI and ACS Style

Ma, Y.; Gopal, S. Geographically Weighted Regression Models in Estimating Median Home Prices in Towns of Massachusetts Based on an Urban Sustainability Framework. Sustainability 2018, 10, 1026. https://0-doi-org.brum.beds.ac.uk/10.3390/su10041026

AMA Style

Ma Y, Gopal S. Geographically Weighted Regression Models in Estimating Median Home Prices in Towns of Massachusetts Based on an Urban Sustainability Framework. Sustainability. 2018; 10(4):1026. https://0-doi-org.brum.beds.ac.uk/10.3390/su10041026

Chicago/Turabian Style

Ma, Yaxiong, and Sucharita Gopal. 2018. "Geographically Weighted Regression Models in Estimating Median Home Prices in Towns of Massachusetts Based on an Urban Sustainability Framework" Sustainability 10, no. 4: 1026. https://0-doi-org.brum.beds.ac.uk/10.3390/su10041026

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