A Comparison of Spatially Varying Regression Coefficient Estimates Using Geographically Weighted and Spatial‐Filter‐Based Techniques

Geographically weighted regression (GWR) is a technique that explores spatial nonstationarity in data‐generating processes by allowing regression coefficients to vary spatially. It is a widely applied technique across domains because it is intuitive and conforms to the well‐understood framework of regression. An alternative method to GWR that has been suggested is spatial filtering, which it has been argued provides a superior alternative to GWR by producing spatially varying regression coefficients that are not correlated with each other and which display less spatial autocorrelation. It is, therefore, worthwhile to examine these claims by comparing the output from both methods. We do this by using simulated data that represent two sets of spatially varying processes and examining how well both techniques replicate the known local parameter values. The article finds no support that spatial filtering produces local parameter estimates with superior properties. The results indicate that the original spatial filtering specification is prone to overfitting and is generally inferior to GWR, while an alternative specification that minimizes the mean square error (MSE) of coefficient estimates produces results that are similar to GWR. However, since we generally do not know the true coefficients, the MSE minimizing specification is impractical for applied research.     

The Right to Rule by Thumb: A Comment on Epistemology in “A Route Map for Successful Applications of Geographically-Weighted Regression”

Comber et al. provide an important contribution to the future of quantitative geography and Geographical Analysis. The contribution is chiefly in their development of a “GWR Route Map,” a diagram showing the sequence of analytical steps that “successful” specification searches in local modeling tend to follow. Geographically weighted techniques have been rapidly expanding, both in terms of complexity, users, and disciplinary reach. With geographically weighted methods now in so many more analysts' hands, any new rule of thumb will have a major imprint. But, by what right does the thumb rule the analysts? That is, what “counts” as valid knowledge about local models in general? In the following comment, I argue that we probably should use theory, not route maps to decide specifications. But, if we are pressed to build route maps, we sorely need better epistemological foundations for them. I discuss a few previous examples of strongly grounded route maps and offer a few paths to these better grounds as well as two ways to the exit.     

基于GWR模型的合肥居住地价影响因素研究

以合肥市主城区为例,基于2010-2014年居住用地的出让数据,运用地统计法、GWR模型等方法,对合肥市居住地价的空间异质性及其影响因素进行研究。研究表明：①合肥市居住地价的空间分布呈现出显著的多中心的空间结构,地价的峰值区分别以老城区、政务区天鹅湖及滨湖新区塘西河公园为中心呈现圈层式分布;②不同的地价影响因素表现出不同的空间分布特征,其中容积率对居住地价的贡献度空间差异最大,其次是宗地面积,主干路次之,交通站点对居住地价的贡献度最小;③厘清各影响因素对地价的作用机制,建立动态的数字地价模型,不仅能促进土地资源的集约利用,重塑城市的空间结构,而且能为城市整体价值的发挥提供重要的理论支撑。     

Single and Multiscale Models of Process Spatial Heterogeneity

Recent work in local spatial modeling has affirmed and broadened interest in multivariate local spatial analysis. Two broad approaches have emerged: Geographically Weighted Regression (GWR) which follows a frequentist perspective and Bayesian Spatially Varying Coefficients models. Although several comparisons between the two approaches exist, recent developments, particularly in GWR, mean that these are incomplete and missing some important axes of comparison. Consequently, there is a need for a more thorough comparison of the two families of local estimators, including recent developments in multiscale variants and their relative performance under controlled conditions. We find that while both types of local models generally perform similarly on a series of criteria, some interesting and important differences exist.     

Inference in Multiscale Geographically Weighted Regression

A recent paper expands the well-known geographically weighted regression (GWR) framework significantly by allowing the bandwidth or smoothing factor in GWR to be derived separately for each covariate in the model—a framework referred to as multiscale GWR (MGWR). However, one limitation of the MGWR framework is that, until now, no inference about the local parameter estimates was possible. Formally, the so-called “hat matrix,” which projects the observed response vector into the predicted response vector, was available in GWR but not in MGWR. This paper addresses this limitation by reframing GWR as a Generalized Additive Model, extending this framework to MGWR and then deriving standard errors for the local parameters in MGWR. In addition, we also demonstrate how the effective number of parameters can be obtained for the overall fit of an MGWR model and for each of the covariates within the model. This statistic is essential for comparing model fit between MGWR, GWR, and traditional global models, as well as for adjusting multiple hypothesis tests. We demonstrate these advances to the MGWR framework with both a simulated data set and a real-world data set and provide a link to new software for MGWR (MGWR1.0) which includes the novel inferential framework for MGWR described here.     

Geographical and Temporal Weighted Regression (GTWR)

Both space and time are fundamental in human activities as well as in various physical processes. Spatiotemporal analysis and modeling has long been a major concern of geographical information science (GIScience), environmental science, hydrology, epidemiology, and other research areas. Although the importance of incorporating the temporal dimension into spatial analysis and modeling has been well recognized, challenges still exist given the complexity of spatiotemporal models. Of particular interest in this article is the spatiotemporal modeling of local nonstationary processes. Specifically, an extension of geographically weighted regression (GWR), geographical and temporal weighted regression (GTWR), is developed in order to account for local effects in both space and time. An efficient model calibration approach is proposed for this statistical technique. Using a 19‐year set of house price data in London from 1980 to 1998, empirical results from the application of GTWR to hedonic house price modeling demonstrate the effectiveness of the proposed method and its superiority to the traditional GWR approach, highlighting the importance of temporally explicit spatial modeling.     

A Note on the Mixed Geographically Weighted Regression Model*

Abstract. A mixed, geographically weighted regression (GWR) model is useful in the situation where certain explanatory variables influencing the response are global while others are local. Undoubtedly, how to identify these two types of the explanatory variables is essential for building such a model. Nevertheless, It seems that there has not been a formal way to achieve this task. Based on some work on the GWR technique and the distribution theory of quadratic forms in normal variables, a statistical test approach is suggested here to identify a mixed GWR model. Then, this note mainly focuses on simulation studies to examine the performance of the test and to provide some guidelines for performing the test in practice. The simulation studies demonstrate that the test works quite well and provides a feasible way to choose an appropriate mixed GWR model for a given data set.     

Using Spatial Filters and Exploratory Data Analysis to Enhance Regression Models of Spatial Data

Residual spatial autocorrelation is a situation frequently encountered in regression analysis of spatial data. The statistical problems arising due to this phenomenon are well‐understood. Original developments in the field of statistical analysis of spatial data were meant to detect spatial pattern, in order to assess whether corrective measures were required. An early development was the use of residual autocorrelation as an exploratory tool to improve regression analysis of spatial data. In this note, we propose the use of spatial filtering and exploratory data analysis as a way to identify omitted but potentially relevant independent variables. We use an example of blood donation patterns in Toronto, Canada, to demonstrate the proposed approach. In particular, we show how an initial filter used to rectify autocorrelation problems can be progressively replaced by substantive variables. In the present case, the variables so retrieved reveal the impact of urban form, travel habits, and demographic and socio‐economic attributes on donation rates. The approach is particularly appealing for model formulations that do not easily accommodate positive spatial autocorrelation, but should be of interest as well for the case of continuous variables in linear regression.     

Designs for Detecting Spatial Dependence

The aim of this article is to find optimal or nearly optimal designs for experiments to detect spatial dependence that might be in the data. The questions to be answered are: how to optimally select predictor values to detect the spatial structure (if it is existent) and how to avoid to spuriously detect spatial dependence if there is no such structure. The starting point of this analysis involves two different linear regression models: (1) an ordinary linear regression model with i.i.d. error terms—the nonspatial case and (2) a regression model with a spatially autocorrelated error term, a so-called simultaneous spatial autoregressive error model. The procedure can be divided into two main parts: The first is use of an exchange algorithm to find the optimal design for the respective data collection process; for its evaluation an artificial data set was generated and used. The second is estimation of the parameters of the regression model and calculation of Moran's I , which is used as an indicator for spatial dependence in the data set. The method is illustrated by applying it to a well-known case study in spatial analysis.     

Identifying the Impacts of Land-Use Spatial Patterns on Street-Network Accessibility Using Geospatial Methods

While the land use-street network nexus is well acknowledged, evidence for the one-way impacts of land-use patterns on street accessibility is still inadequate. The measurements of land-use patterns and street accessibility lack systematic knowledge. Their empirical correlations also lack geographical variability, constraining site-specific land-use practices. Therefore, this study overcame the aforementioned limitations by examining the two-level spatial models to formulate accessibility-oriented land plans, using a well-developed Chinese city as an example. Firstly, two landscape metrics—Euclidean Nearest-Neighbor Distance (ENN) and Similarity Index (SIMI)—were used to quantify the intra- and inter-land-use configurations, respectively. Both city-level and local accessibility were measured using spatial design network analysis. Performing both ordinary least squares (OLS) and geographically weighted regression (GWR) models, results identified the statistically significant effects of inter-land-use patterns on two-level street accessibility. An exception was that land-use configurations within residential and industrial regions were irrelevant to street accessibility. We also found GWR was a better-fitting model than OLS when estimating locally-varied accessibility, suggesting hierarchical multiscale land-use planning. Overall, locally heterogeneous evidence in this study can substantialize land use-street network interactions and support the decision-making and implementation of place-specific accessibility-oriented land use.     

The Spatial Association of Demographic and Population Health Characteristics with COVID-19 Prevalence Across Districts in India

In less-developed countries, the lack of granular data limits the researcher's ability to study the spatial interaction of different factors on the COVID-19 pandemic. This study designs a novel database to examine the spatial effects of demographic and population health factors on COVID-19 prevalence across 640 districts in India. The goal is to provide a robust understanding of how spatial associations and the interconnections between places influence disease spread. In addition to the linear Ordinary Least Square regression model, three spatial regression models—Spatial Lag Model, Spatial Error Model, and Geographically Weighted Regression are employed to study and compare the variables explanatory power in shaping geographic variations in the COVID-19 prevalence. We found that the local GWR model is more robust and effective at predicting spatial relationships. The findings indicate that among the demographic factors, a high share of the population living in slums is positively associated with a higher incidence of COVID-19 across districts. The spatial variations in COVID-19 deaths were explained by obesity and high blood sugar, indicating a strong association between pre-existing health conditions and COVID-19 fatalities. The study brings forth the critical factors that expose the poor and vulnerable populations to severe public health risks and highlight the application of geographical analysis vis-a-vis spatial regression models to help explain those associations.     

The unequal geographies of climate finance: Climate injustice and dependency in the world system

Central to climate justice is the question of who will pay for the mitigation and adaptation efforts needed as the climate crisis worsens, particularly in countries that bear little responsibility for global greenhouse gas emissions. Climate finance is a complex set of mechanisms intended to address this concern. World-systems theory has long understood international development assistance as a tool that reproduces spatial dependency between states. In this paper, we ask whether climate finance follows the expectations of world-systems theory and reproduces relationships of dependency, or if it instead advances climate justice and challenges spatial dependency in the world-system. Through this analysis, we consider the implications of climate finance for world-systems theory. We use recent empirical data to ask whether climate finance follows or challenges world-systems theory expectations, focusing on five areas: (1) spatial flows of climate finance between the core, semi-periphery, and periphery; (2) the governance of climate finance institutions; (3) the types of projects supported by climate finance; (4) the relationship of projects to dominant systems of extraction, production, and consumption; and (5) the agency of peripheral state and non-state actors in shaping climate finance in relation to their interests. Taken together, we argue that climate finance in many ways reproduces relationships of dependency, though potential avenues exist for contesting this unequal balance of power and for advocating for climate justice. This case illustrates the need to approach analyses of dependency in a nuanced way, interrogating specific processes through which dependency is produced and contested across scales.     

Navigating the Methodological Landscape in Spatial Analysis: A Comment on “A Route Map for Successful Applications of Geographically-Weighted Regression”

The development of “route maps” for spatial analytical methods is a pursuit with important ramifications. Comber et al. propose a route map to guide applications of geographically weighted regression consisting of a three-step primary pathway and a series of secondary arterials. This comment first highlights some concerns about the underlying “map” (i.e., experimental setup and assumptions) and then with the proposed “route” (i.e., core decisions and evaluation criteria). It closes by suggesting a more general focus on identifying modeling issues with the highest impact and facilitating consensus-building, which could improve the future production of route maps for navigating the methodological landscape in spatial analysis.     

Some Notes on Parametric Significance Tests for Geographically Weighted Regression

The technique of geographically weighted regression (GWR) is used to model spatial 'drift' in linear model coefficients. In this paper we extend the ideas of GWR in a number of ways. First, we introduce a set of analytically derived significance tests allowing a null hypothesis of no spatial parameter drift to be investigated. Second, we discuss 'mixed' GWR models where some parameters are fixed globally but others vary geographically. Again, models of this type may be assessed using significance tests. Finally, we consider a means of deciding the degree of parameter smoothing used in GWR based on the Mallows C_p statistic. To complete the paper, we analyze an example data set based on house prices in Kent in the U.K. using the techniques introduced.     

Geographically Weighted Regression: A Method for Exploring Spatial Nonstationarity

Spatial nonstationarity is a condition in which a simple “global” model cannot explain the relationships between some sets of variables. The nature of the model must alter over space to reflect the structure within the data. In this paper, a technique is developed, termed geographically weighted regression, which attempts to capture this variation by calibrating a multiple regression model which allows different relationships to exist at different points in space. This technique is loosely based on kernel regression. The method itself is introduced and related issues such as the choice of a spatial weighting function are discussed. Following this, a series of related statistical tests are considered which can be described generally as tests for spatial nonstationarity. Using Monte Carlo methods, techniques are proposed for investigating the null hypothesis that the data may be described by a global model rather than a non-stationary one and also for testing whether individual regression coefficients are stable over geographic space. These techniques are demonstrated on a data set from the 1991 U.K. census relating car ownership rates to social class and male unemployment. The paper concludes by discussing ways in which the technique can be extended.     

A Note on the Extremities of Local Moran's Iis and Their Impact on Global Moran's I

By defining local Moran's I_i as a ratio of quadratic forms and making use of its overall additivity to match global Moran's I, we can identify spatial objects with a strong impact on global Moran's I. First, we concentrate on the spatial properties of local Moran's I_i expressed by the local linkage degree. Depending on whether we use the W- or C-coding of the spatial connectivity matrix, the variance of local Moran's I_i for a small local linkage degree will be either large or small. Note that spatial objects associated with a local Moran's I_i with a large variance affect the global statistic much more than spatial objects associated with a local Moran's I_i with a small variance. Counterintuitively, global Moran's I defined in the W-coding is most influenced by spatial objects with a small number of spatial neighbors. In contrast, spatial objects with a large number of spatial neighbors exert more impact on global Moran's I setup in the C-coding. Second, we investigate the impact of the empirical data on local Moran's I_i and show that local Moran's I_i will only be significant for extreme absolute residuals at and around the reference location. Clusters of average regression residuals cannot be detected by local Moran's I_i. Consequently, spatial cliques of extreme residuals contribute more to significance tests on global autocorrelation.     

Set-determined Process and the Growth of Spatial Structure

We present an analysis of the effects of different information flow patterns on the growth of spatial structure in a broad class of systems. Examples are drawn from plant and message diffusion processes, the growth of herding in animal populations and, in particular, the evolution of cities in human populations. We model the actions of an “individual” as a response to gradients in a transformation of the “real” space, this transformation being the individual's information map of the space. Interesting deductions occur as a result of allowing the individual's information map to vary from extremely local to extremely global. In the city growth model for example, if migration decisions are based upon extremely local information, then the system is characterized by instability and the lack of any well-defined spatial structure. As the localness of this information map decreases, the system is characterized by the development of well defined spatial structure with a characteristic distance separating the “cities.” In a particular case, we show that the wavelength of the dominant spatial structure is directly related to the size of an individual's information field. These effects are a consequence of interactions between individuals that arise from the nonlocalness of the information fields. We interpret the early growth stages in these models in terms of linear filters.     

The Saddlepoint Approximation of Moran's I's and Local Moran's Ii's Reference Distributions and Their Numerical Evaluation

Global Moran's I and local Moran's I_i are the most commonly used test statistics for spatial autocorrelation in univariate map patterns or in regression residuals. They belong to the general class of ratios of quadratic forms for whom a whole array of approximation techniques has been proposed in the statistical literature, such as the prominent saddlepoint approximation by Offer Lieberman (1994). The saddlepoint approximation outperforms other approximation methods with respect to its accuracy and computational costs. In addition, only the saddlepoint approximation is capable of handling, in analytical terms, reference distributions of Moran's I that are subject to significant underlying spatial processes. The accuracy and computational benefits of the saddlepoint approximation are demonstrated for a set of local Moran's Ii statistics under either the assumption of global spatial independence or subject to an underlying global spatial process. Local Moran's I_i is known to have an excessive kurtosis and thus void the use of the simple approximation methods of its reference distribution. The results demonstrate how well the saddlepoint approximation fits the reference distribution of local Moran's I_i. Furthermore, for local Moran's I_i under the assumption of global spatial independence several algebraic simplifications lead to substantial gains in numerical efficiency. This makes it possible to evaluate local Moran's I_i's significance in large spatial tessellations.     

A Hybrid Approach for Mass Valuation of Residential Properties through Geographic Information Systems and Machine Learning Integration

Geographic Information Systems (GIS) and Machine Learning methods are now widely used in mass property valuation using the physical attributes of properties. However, locational criteria, such as as proximity to important places, sea or forest views, flat topography are just some of the spatial factors that affect property values and, to date, these have been insufficiently used as part of the valuation process. In this study, a hybrid approach is developed by integrating GIS and Machine Learning for mass valuation of residential properties. GIS-based Nominal Valuation Method was applied to carry out proximity, terrain, and visibility analyses using Ordnance Survey and OpenStreetMap data, than land value map of Great Britain was produced. Spatial criteria scores obtained from the GIS analyses were included in the price prediction process in which global and spatially clustered local regression models are built for England and Wales using Price Paid Data and Energy Performance Certificates data. Results showed that adding locational factors to the property price data and applying a novel nominally weighted spatial clustering algorithm for creating a local regression increased the prediction accuracy by about 45%. It also demonstrated that Random Forest was the most accurate ensemble model.     

A Simulation Study on Specifying a Regression Model for Spatial Data: Choosing between Autocorrelation and Heterogeneity Effects

In this simulation study, regressions specified with autocorrelation effects are compared against those with relationship heterogeneity effects, and in doing so, provides guidance on their use. Regressions investigated are: (1) multiple linear regression, (2) a simultaneous autoregressive error model, and (3) geographically weighted regression. The first is nonspatial and acts as a control, the second accounts for stationary spatial autocorrelation via the error term, while the third captures spatial heterogeneity through the modeling of nonstationary relationships between the response and predictor variables. The geostatistical‐based simulation experiment generates data and coefficients with known multivariate spatial properties, all within an area‐unit spatial setting. Spatial autocorrelation and spatial heterogeneity effects are varied and accounted for. On fitting the regressions, that each have different assumptions and objectives, to very different geographical processes, valuable insights to their likely performance are uncovered. Results objectively confirm an inherent interrelationship between autocorrelation and heterogeneity, that results in an identification problem when choosing one regression over another. Given this, recommendations on the use and implementation of these spatial regressions are suggested, where knowledge of the properties of real study data and the analytical questions being posed are paramount.