Local Spatial Autocorrelation Statistics: Distributional Issues and an Application

The statistics G_i(d) and G_i*(d), introduced in Getis and Ord (1992) for the study of local pattern in spatial data, are extended and their properties further explored. In particular, nonbinary weights are allowed and the statistics are related to Moran's autocorrelation statistic, I. The correlations between nearby values of the statistics are derived and verified by simulation. A Bonferroni criterion is used to approximate significance levels when testing extreme values from the set of statistics. An example of the use of the statistics is given using spatial-temporal data on the AIDS epidemic centering on San Francisco. Results indicate that in recent years the disease is intensifying in the counties surrounding the city.     

Local Indicators of Spatial Association—LISA

The capabilities for visualization, rapid data retrieval, and manipulation in geographic information systems (GIS) have created the need for new techniques of exploratory data analysis that focus on the “spatial” aspects of the data. The identification of local patterns of spatial association is an important concern in this respect. In this paper, I outline a new general class of local indicators of spatial association (LISA) and show how they allow for the decomposition of global indicators, such as Moran's I, into the contribution of each observation. The LISA statistics serve two purposes. On one hand, they may be interpreted as indicators of local pockets of nonstationarity, or hot spots, similar to the G_i and G*_i statistics of Getis and Ord (1992). On the other hand, they may be used to assess the influence of individual locations on the magnitude of the global statistic and to identify “outliers,” as in Anselin's Moran scatterplot (1993a). An initial evaluation of the properties of a LISA statistic is carried out for the local Moran, which is applied in a study of the spatial pattern of conflict for African countries and in a number of Monte Carlo simulations.     

Comparative Spatial Filtering in Regression Analysis

One approach to dealing with spatial autocorrelation in regression analysis involves the filtering of variables in order to separate spatial effects from the variables’ total effects. In this paper we compare two filtering approaches, both of which allow spatial statistical analysts to use conventional linear regression models. Getis’ filtering approach is based on the autocorrelation observed with the use of the G_i local statistic. Griffith's approach uses an eigenfunction decomposition based on the geographic connectivity matrix used to compute a Moran's I statistic. Economic data are used to compare the workings of the two approaches. A final comparison with an autoregressive model strengthens the conclusion that both techniques are effective filtering devices, and that they yield similar regression models. We do note, however, that each technique should be used in its appropriate context.     

Testing for Local Spatial Autocorrelation in the Presence of Global Autocorrelation

A fundamental concern of spatial analysts is to find patterns in spatial data that lead to the identification of spatial autocorrelation or association. Further, they seek to identify peculiarities in the data set that signify that something out of the ordinary has occurred in one or more regions. In this paper we provide a statistic that tests for local spatial autocorrelation in the presence of the global autocorrelation that is characteristic of heterogeneous spatial data. After identifying the structure of global autocorrelation, we introduce a new measure that may be used to test for local structure. This new statistic Oi is asymptotically normally distributed and allows for straightforward tests of hypotheses. We provide several numerical examples that illustrate the performance of this statistic and compare it with another measure that does not account for global structure.     

USING THE GETIS STATISTIC TO EXPLORE AGGREGATION EFFECTS IN METROPOLITAN TORONTO CENSUS DATA

In an earlier paper by the same authors (Amrhein and Reynolds 1996), the variation of the G statistic (Getis 1991; Getis and Ord 1992) was demonstrated to correlate very highly with calculated aggregation effects. In the conclusions to this earlier paper, we called for experiments with different data sets, different aggregation rules, and different definitions of the connectivity matrix. In this paper, we explore the behaviour of census data from enumeration areas of the Toronto Census Metropolitan Area. Experiments are conducted that describe the aggregation effects that creep into the data at different scales, and the effect of the aggregation algorithm on the results. Finally, we test the ability of a modified G statistic, our Getis statistic, to predict aggregation effects     

A Descriptive Analysis of Discrete U.S. Industrial Complexes*

Abstract. We use the Getis/Ord local G statistic and detailed county‐level industry employment data from the U.S. Bureau of Labor Statistics to isolate discrete industrial complexes—or groups of nominally linked industries clustered in particular locations—for two recent years: 1989 and 1997. We describe the characteristics of the complexes in terms of their number, spatial extent, broad regional distribution, and other factors. Data from the two periods help illustrate key shifts in industrial locations, including the continuing concentration of the apparel industry in the Southeast and the ongoing southern shift in U.S. vehicle production.     

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.     

A Local Model for Multivariate Analysis: Extending Wartenberg’s Multivariate Spatial Correlation

Local analysis can provide specific information about individual observations that is often useful in understanding nonstationary interactions among variables. This paper extends the application of Wartenberg’s Multivariate Spatial Correlation (MSC) method to a local setting. The original MSC can be considered as an adaptation of Principal Component Analysis for spatial effects with respect to spatial autocorrelation. The extended MSC method described in this paper, however, further incorporates another spatial effect, spatial heterogeneity, by the addition of geographic weights in standardizing the data and in calculating the spatial association weight matrix. The extension allows more local analysis and facilitates additional visualization of the results. The geographically weighted MSC is illustrated and justified using the classic dataset collected by André-Michel Guerry on moral statistics in 1830s France.     

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.     

Power of the Rank Adjacency Statistic to Detect Spatial Clustering in a Small Number of Regions

The rank adjacency statistic D is a statistical method for assessing spatial autocorrelation or clustering of geographical data. It was originally proposed for summarizing the geographical patterns of cancer data in Scotland (IARC 1985). In this paper, we investigate the power of the rank adjacency statistic to detect spatial clustering when a small number of regions is involved. The investigations were carried out using Monte Carlo simulations, which involved generating patterned/clustered values and computing the power with which the D statistic would detect it. To investigate the effects of region shapes, structure of the regions, and definition of weights, simulations were carried out using two different region shapes, binary and nonhinary weights, and three different lattice structures. The results indicate that in the typical example of considering Canadian total mortality at the electoral district level, the D statistic had adequate power to detect general spatial autocorrelation in twenty‐five or more regions. There was an inverse relationship between power and the level of connectedness of the regions, which depends on the weighting function, shape, and arrangement of the regions. The power of the D statistic was also found to compare favorably with that of Moran's I statistic.     

"The Problem of Spatial Autocorrelation" and Local Spatial Statistics

This article examines the relationship between spatial dependency and spatial heterogeneity, two properties unique to spatial data. The property of spatial dependence has led to a large body of research into spatial autocorrelation and also, largely independently, into geostatistics. The property of spatial heterogeneity has led to a growing awareness of the limitation of global statistics and the value of local statistics and local statistical models. The article concludes with a discussion of how the two properties can be accommodated within the same modelling framework.     

A Statistical Method for the Detection of Geographic Clustering

Kernel‐based, smoothed estimates of spatial variables are useful in exploratory analyses because they yield a clear visual image of geographic variability in the underlying variable. In this paper I suggest an approach for assessing the significance of peaks in the surface that result from the application of the smoothing kernel. The approach may also be thought of as a method for assessing the maximum among a set of suitably defined local statistics. Local statistics for data on a regular grid of cells are first defined by using a Gaussian kernel. Results from integral geometry are then used to find the probability that the maximum local statistic (M) exceeds a given critical value (M). Approximations are provided that make implementation of the approach straightforward. Future work will address several other issues associated with local statistics that have been defined in this way, including edge effects, and the effects of global spatial autocorrelation on the choice of critical value.     

Using AMOEBA to Create a Spatial Weights Matrix and Identify Spatial Clusters

The creation of a spatial weights matrix by a procedure called AMOEBA, A Multidirectional Optimum Ecotope-Based Algorithm , is dependent on the use of a local spatial autocorrelation statistic. The result is (1) a vector that identifies those spatial units that are related and unrelated to contiguous spatial units and (2) a matrix of weights whose values are a function of the relationship of the ith spatial unit with all other nearby spatial units for which there is a spatial association. In addition, the AMOEBA procedure aids in the demarcation of clusters, called ecotopes, of related spatial units. Experimentation reveals that AMOEBA is an effective tool for the identification of clusters. A comparison with a scan statistic procedure (SaTScan) gives evidence of the value of AMOEBA. Total fertility rates in enumeration districts in Amman, Jordan, are used to show a real-world example of the use of AMOEBA for the construction of a spatial weights matrix and for the identification of clusters. Again, comparisons reveal the effectiveness of the AMOEBA procedure.     

The Colocation Quotient: A New Measure of Spatial Association Between Categorical Subsets of Points. 协同区位商：点集分类子集间空间关联性的新度量标准

This article presents a new metric we label the colocation quotient (CLQ), a measurement designed to quantify (potentially asymmetrical) spatial association between categories of a population that may itself exhibit spatial autocorrelation. We begin by explaining why most metrics of categorical spatial association are inadequate for many common situations. Our focus is on where a single categorical data variable is measured at point locations that constitute a population of interest. We then develop our new metric, the CLQ, as a point‐based association metric most similar to the cross‐k‐function and join count statistic. However, it differs from the former in that it is based on distance ranks rather than on raw distances and differs from the latter in that it is asymmetric. After introducing the statistical calculation and underlying rationale, a random labeling technique is described to test for significance. The new metric is applied to economic and ecological point data to demonstrate its broad utility. The method expands upon explanatory powers present in current point‐based colocation statistics.     

Spatial Autocorrelation in Ecological Studies: A Legacy of Solutions and Myths

A major aim of including the spatial component in ecological studies is to characterize the nature and intensity of spatial relationships between organisms and their environment. The growing awareness by ecologists of the importance of including spatial structure in ecological studies (for hypothesis development, experimental design, statistical analyses, and spatial modeling) is beneficial because it promotes more effective research. Unfortunately, as more researchers perform spatial analysis, some misconceptions about the virtues of spatial statistics have been carried through the process and years. Some of these statistical concepts and challenges were already presented by Cliff and Ord in 1969. Here, we classify the most common misconceptions about spatial autocorrelation into three categories of challenges: (1) those that have no solutions, (2) those where solutions exist but are not well known, and (3) those where solutions have been proposed but are incorrect. We conclude in stressing where new research is needed to address these challenges.     

Constructing the Spatial Weights Matrix Using a Local Statistic

Spatial weights matrices are necessary elements in most regression models where a representation of spatial structure is needed. We construct a spatial weights matrix, W , based on the principle that spatial structure should be considered in a two‐part framework, those units that evoke a distance effect, and those that do not. Our two‐variable local statistics model (LSM) is based on the G_i* local statistic. The local statistic concept depends on the designation of a critical distance, d_c, defined as the distance beyond which no discernible increase in clustering of high or low values exists. In a series of simulation experiments LSM is compared to well‐known spatial weights matrix specifications—two different contiguity configurations, three different inverse distance formulations, and three semi‐variance models. The simulation experiments are carried out on a random spatial pattern and two types of spatial clustering patterns. The LSM performed best according to the Akaike Information Criterion, a spatial autoregressive coefficient evaluation, and Moran's I tests on residuals. The flexibility inherent in the LSM allows for its favorable performance when compared to the rigidity of the global models.     

A History of the Concept of Spatial Autocorrelation: A Geographer's Perspective

Spatial autocorrelation is a concept that helps to define the field of spatial analysis. It is central to studies using spatial statistics and spatial econometrics. In this paper, we trace the early development of the concept and explain the academic links that brought the concept to the fore in the late 1960s. In geography, the importance of the work of Michael F. Dacey, Andrew D. Cliff, and J. Keith Ord is emphasized. Later, with the publication of a volume on spatial econometrics by Luc Anselin, spatial research and the use of the concept of spatial autocorrelation received a considerable boost. These developments are outlined together with comments about recent and possible future trends in spatial autocorrelation-based research.     

Spatial autocorrelation tests and the Classic Maya collapse: Methods and inferences

Spatial autocorrelation, resulting in pattern or structure in geographically distributed data, is discussed in theoretical and practical terms. Tests for spatial autocorrelation are presented, along with an explication of the relationship between autocorrelation models, the product-moment correlation coefficient and the spatial autocorrelation test statistic. Two archaeological examples illustrate the application of the auto-correlation test statistic. The first uses a hypothetical data set, which shows the type of map patterns that appear with various levels of spatial autocorrelation, and the second examines the terminal distribution of long-count-dated monuments at lowland Classic Maya sites. The results of the second example fail to support arguments for simple patterning in the cessation of the erection of these monuments and, by inference, in the Maya collapse itself. Finally, it is argued that while the identification of spatial autocorrelation is often the goal of spatial analyses, the presence of autocorrelation violates the assumptions of certain statistics used in such analyses.     

Impact of Cliff and Ord (1969, 1981) on Spatial Epidemiology

The influence in spatial epidemiology of the seminar work on autocorrelation by Cliff and Ord is discussed. Quantifying the evidence of spatial clustering was an important step in the development of modern statistical methods for analyzing spatial variations of diseases. Autocorrelation is nowadays mostly accounted for at a latent level within a hierarchical framework to small area disease mapping. The importance of accounting for autocorrelation in geographical correlation studies is also reviewed.