Exploring Spatial Patterns Using an Expanded Spatial Autocorrelation Framework

Spatial autocorrelation (SA) is regarded as an important dimension of spatial pattern. SA measures usually consist of two components: measuring the similarity of attribute values and defining the spatial relationships among observations. The latter component is often represented by a spatial weights matrix that predefines spatial relationship between observations in most measures. Therefore, SA measures, in essence, are measures of attribute similarity, conditioned by spatial relationship. Another dimension of spatial pattern can be explored by controlling observations to be compared based upon the degree of attribute similarity. The resulting measures are spatial proximity measures of observations, meeting predefined attribute similarity criteria. Proposed measures reflect degrees of clustering or dispersion for observations meeting certain levels of attribute similarity. An existing spatial autocorrelation framework is expanded to a general framework to evaluate spatial patterns and can accommodate the proposed approach measuring proximity. Analogous to the concept of variogram, clustergram is proposed to show the levels of spatial clustering over a range of attribute similarity, or attribute lags. Specific measures based on the proposed approach are formulated and applied to a hypothetical landscape and an empirical example, showing that these new measures capture spatial pattern information not reflected by traditional spatial autocorrelation measures.     

SPATIAL DEPENDENCY OF SEGREGATION INDICES

A few researchers have mentioned the scale sensitivity of segregation index, D. In this paper, I discuss analytically and empirically why using large enumeration areal units usually results in low segregation measures, and using small areal units produces relatively high segregation measures. The discussion is also applicable to the multi-group variant of D. A major finding is that if people of the same ethnic groups are positively spatially auto-correlated, increasing the size of areal units of analysis may not lower D initially, because only people of the same group are added. But enlarging the areal units subsequently may include population of other ethnic groups, and therefore could lower D. However, if the boundaries of the larger enumeration units are drawn to include only population of the same group, then D will not change significantly. Both the spatial autocorrelation of ethnic group population and zonal pattern are critical factors in determining the scale sensitivity of D.     

Measures of Feature-Based and Areal Differentiations and their use in City Services

Two types of geographical differentiations are distinguished: feature-based differentiation, in which a selected property assumes values distributed over a particular area, and areal differentiation, in which an area is broken down into weighted segments. The concept of entropy, a measure of the amount of information transmitted, is borrowed from information theory to provide measures of feature-based and areal differentiations. A possible application of these measures of differentiation is suggested in the study and planning of city services.     

ACCOUNTING FOR NEIGHBORING EFFECTS IN MEASURES OF SPATIAL CONCENTRATION*

ABSTRACT A common problem with spatial economic concentration measures based on areal data (e.g., Gini, Herfindhal, entropy, and Ellison‐Glaeser indices) is accounting for the position of regions in space. While they purport to measure spatial clustering, these statistics are confined to calculations within individual areal units. They are insensitive to the proximity of regions or to neighboring effects. Clearly, since spillovers do not recognize areal units, economic clusters may cross regional boundaries. Yet with current measures, any industrial agglomeration that traverses boundaries will be chopped into two or more pieces. Activity in adjacent spatial units is treated in exactly the same way as activity in far‐flung, nonadjacent areas. This paper shows how popular measures of spatial concentration relying on areal data can be modified to account for neighboring effects. With a U.S. application, we also demonstrate that the new instruments we propose are easy to implement and can be valuable in regional analysis.     

Aggregation of Sampling Units: An Analytical Solution to Predict Variance

Geographical variables generally show spatially structured patterns corresponding to intrinsic characteristics of the environment. The size of the sampling unit has a critical effect on our perception of phenomena and is closely related to the variance and correlation structure of the data. Geostatistical theory uses analytical relationships for change of support (change of sampling unit size), allowing prediction of the variance and autocorrelation structure that would be observed if a survey was conducted using different sampling unit sizes. To check the geostatistical predictions, we use a test case about tree density in the tropical rain forest of the Pasoh Reserve, Malaysia. This data set contains exhaustive information about individual tree locations, so it allows us to simulate and compare various sampling designs. The original data set was reorganized to compute tree densities for 5 times 5-, 10 times 10-, and 20 times 20-meter quadrat sizes. Based upon the 5 times 5-meter data set, the spatial structure is modeled using a nugget effect (white noise) plus an exponential model. The change of support relationships, using within-quadrat variances inferred from the variogram model, predict the spatial autocorrelation structure and new variances corresponding to 10 times 10-meter and 20 times 20-meter quadrats. The theoretical and empirical results agreed closely, whereas neglecting the autocorrelation structure would have led to largely underestimating the variance. As the quadrat size increases, the range of autocorrelation increases, while the variance and the proportion of noise in the data decrease.     

Geostatistical Prediction and Simulation of Point Values from Areal Data

The spatial prediction and simulation of point values from areal data are addressed within the general geostatistical framework of change of support (the term support referring to the domain informed by each measurement or unknown value). It is shown that the geostatistical framework (i) can explicitly and consistently account for the support differences between the available areal data and the sought-after point predictions, (ii) yields coherent (mass-preserving or pycnophylactic) predictions, and (iii) provides a measure of reliability (standard error) associated with each prediction. In the case of stochastic simulation, alternative point-support simulated realizations of a spatial attribute reproduce (i) a point-support histogram (Gaussian in this work), (ii) a point-support semivariogram model (possibly including anisotropic nested structures), and (iii) when upscaled, the available areal data. Such point-support-simulated realizations can be used in a Monte Carlo framework to assess the uncertainty in spatially distributed model outputs operating at a fine spatial resolution because of uncertain input parameters inferred from coarser spatial resolution data. Alternatively, such simulated realizations can be used in a model-based hypothesis-testing context to approximate the sampling distribution of, say, the correlation coefficient between two spatial data sets, when one is available at a point support and the other at an areal support. A case study using synthetic data illustrates the application of the proposed methodology in a remote sensing context, whereby areal data are available on a regular pixel support. It is demonstrated that point-support (sub-pixel scale) predictions and simulated realizations can be readily obtained, and that such predictions and realizations are consistent with the available information at the coarser (pixel-level) spatial resolution.     

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.     

Geographical Differentiation and its Measures in a Discrete System

The preceding preliminary paper is expanded into a full-fledged theory of differentiation. Two types are considered: feature-based differentiation and areal differentiation, both in their geographical and nongeographical contexts. Each type is examined both in a broad sense, without weighting, and in a narrow sense, with the values of a selected feature or other weights attached to the parts of a whole. Measures of differentiation are then borrowed from information theory. Differentiation in the broad sense (without weights) can be measured by the simpler Hartley measure of information transmission, which applies when any one of distinct possible outcomes of an event is equally likely to occur. Differentiation in the narrow sense (weighted) can be measured by the more complicated Shannon measure of entropy, which takes into consideration the probability that a given event will have a particular outcome.     

flowAMOEBA: Identifying Regions of Anomalous Spatial Interactions

This study aims at developing a data‐driven and bottom‐up spatial statistic method for identifying regions of anomalous spatial interactions (clusters of extremely high‐ or low‐value spatial flows), based on which it creates a spatial flow weights matrix. The method, dubbed flowAMOEBA, upgrades a multidirectional optimum ecotope‐based algorithm (AMOEBA) from areal data to spatial flow data through a proper spatial flow neighborhood definition. The method has the potential to dramatically change the way we study spatial interactions. First, it breaks the convention that spatial interaction data are always collected and modeled between spatial entities of the same granularity, as it delineates the OD region of anomalous spatial interactions, regardless of the size, shape, scale, or administrative level. Second, the method creates an empirical spatial flow weights matrix that can handle network autocorrelation embedded in spatial interaction modeling, thus improving related policy‐making or problem‐solving strategies. flowAMOEBA is tested and demonstrated on a synthetic data set as well as a county‐to‐county migration data set.     

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.     

Measuring Spatial Autocorrelation of Vectors

This article introduces measures to quantify spatial autocorrelation for vectors. In contrast to scalar variables, spatial autocorrelation for vectors involves an assessment of both direction and magnitude in space. Extending conventional approaches, measures of global and local spatial associations for vectors are proposed, and the associated statistical properties and significance testing are discussed. The new measures are applied to study the spatial association of taxi movements in the city of Shanghai. Complications due to the edge effect are also examined.     

GeoDa : An Introduction to Spatial Data Analysis

This article presents an overview of GeoDa™, a free software program intended to serve as a user-friendly and graphical introduction to spatial analysis for non-geographic information systems (GIS) specialists. It includes functionality ranging from simple mapping to exploratory data analysis, the visualization of global and local spatial autocorrelation, and spatial regression. A key feature of GeoDa is an interactive environment that combines maps with statistical graphics, using the technology of dynamically linked windows. A brief review of the software design is given, as well as some illustrative examples that highlight distinctive features of the program in applications dealing with public health, economic development, real estate analysis, and criminology.     

中国城市信息服务业发展影响因素研究——基于中国191个城市的空间横截面数据分析

采用空间统计和空间计量经济学原理和方法对2005-2010年间中国191个城市信息服务业发展的影响因素进行了实证研究。结果表明:城市信息化水平和人力资本对信息服务业的发展具有显著的正向作用;对外开放水平对信息服务业的发展影响不明显;城市间的产业关联对信息服务业发展没有通过显著性水平检验,这表明城市间物质层面的相互交换并没有对信息服务业的发展产生应有的正效应,而相反城市间的信息流动和知识溢出对信息服务业的发展具有显著的促进作用。     

Spatial Autocorrelation Statistics of Areal Prevalence Rates under High Uncertainty in Denominator Data

We propose a new estimator of spatial autocorrelation of areal incidence or prevalence rates in small areas, such as crime and health indicators, for correcting spatially heterogeneous sampling errors in denominator data. The approach is dubbed the heteroscedasticity‐consistent empirical Bayes (HC‐EB) method. As American Community Survey (ACS) data have been released to the public for small census geographies, small‐area estimates now form the demographic landscape of neighborhoods. Meanwhile, there is growing awareness of the diminished statistical validity of global and local Moran’s I when such small‐area estimates are used in denominator data. Using teen birth rates by census tracts in Mecklenburg County, North Carolina, we present comparisons of conventional and new HC‐EB estimates of Global and Local Moran’s I statistics created on ACS data, along with estimates on ground truth values from the 2010 decennial census. Results show that the new adjustment method dramatically enhances the statistical validity of global and local spatial autocorrelation statistics.     

PERSISTENCE OF REGIONAL UNEMPLOYMENT: APPLICATION OF A SPATIAL FILTERING APPROACH TO LOCAL LABOR MARKETS IN GERMANY*

ABSTRACT The geographical distribution and persistence of regional/local unemployment rates in heterogeneous economies (such as Germany) have been, in recent years, the subject of various theoretical and empirical studies. Several researchers have shown an interest in analyzing the dynamic adjustment processes of unemployment and the average degree of dependence of the current unemployment rates or gross domestic product from the ones observed in the past. In this paper, we present a new econometric approach to the study of regional unemployment persistence, in order to account for spatial heterogeneity and/or spatial autocorrelation in both the levels and the dynamics of unemployment. First, we propose an econometric procedure suggesting the use of spatial filtering techniques as a substitute for fixed effects in a panel estimation framework. The spatial filter computed here is a proxy for spatially distributed region‐specific information (e.g., the endowment of natural resources, or the size of the “home market”) that is usually incorporated in the fixed effects coefficients. The advantages of our proposed procedure are that the spatial filter, by incorporating region‐specific information that generates spatial autocorrelation, frees up degrees of freedom, simultaneously corrects for time‐stable spatial autocorrelation in the residuals, and provides insights about the spatial patterns in regional adjustment processes. We present several experiments in order to investigate the spatial pattern of the heterogeneous autoregressive coefficients estimated for unemployment data for German NUTS‐3 regions. We find widely heterogeneous but generally high persistence in regional unemployment rates.     

Redundancy as an Index of Change in Point Pattern Analysis

Point pattern analysis based on concepts from information theory can go beyond existing techniques. Direct measurement of spatial form is achieved when Thiessen polygons are constructed around the points; in this scheme, each point's proportion of total area may be treated like a probability. Three information-theoretic indices are available for analysis of a distribution of such probabilities. Entropy is density dependent. Redundancy, defined as the difference between observed and maximum entropy, seems to avoid this problem when the number of individuals in a pattern exceeds twenty. Comparisons of prior and posterior redundancy provide an indication of change in overall pattern form. An information gain expectation reflects changes for each individual in a pattern. Here, point-area redundancy parameters are determined for Poisson-generated patterns, using a gamma distribution of polygon areas and a computer-generated set. An application to an urban crime problem illustrates the use of these parameters in the analysis of pattern change.     

Potential Contributions of Spatial Analysis to Geographic Information Systems for Transportation (GIS-T)

Geographic information systems (GIS) are having tremendous impacts on many scientific and application domains. The traditional subfield of spatial analysis is witnessing a major resurgence and enhancement due to GIS and geographical information science (GISci), an interdisciplinary field focusing on the theory and methodology underlying GIS software. The interdisciplinary field of geographic information systems for transportation (GIS-T) has emerged to focus on the role of GIS in transportaton analysis and planning. This paper suggests the benefits of closer linkages between spatial analysis, GISci, and transportation through a focused review of spatial analytical issues and their potential contributions to GIS-T. Specifically, this paper reviews the following issues: (i) modifiable areal units; (ii) boundary problems and spatial sampling; (iii) spatial dependence and spatial heterogeneity; and (iv) alternative representations of geographic environments. The discussion highlights the general issues as well as identifies their specific relevance to GIS-T. In addition, this paper identifies some emerging tools from GISci that can address these spatial analytical issues in GIS-T.     

The Size-adjusted Critical Region of Moran's I Test Statistics for Spatial Autocorrelation and Its Application to Geographical Areas

In this article, we explore the expression of the asymptotic approximation of the distribution function of Moran's I test statistic for the check of spatial autocorrelation, and we derive a more accurate critical value, which gives the rejection region similar to significant level α to the order of N^?1 (N = sample size). We show that in some cases our procedure effectively finds the significance of positive spatial autocorrelation in the problem testing for the lack of the spatial autocorrelation. Compared with our method, the testing procedure of Cliff and Ord (1971) is clearly ad hoc and should not be applied blindly, as they pointed out. Our procedure is derived from the theory of asymptotic expansion. We numerically analyze four types of county systems with rectangular lattices and three regional areas with irregular lattices.