Local Spatial Autocorrelation in a Biological Model

We review the recently developed local spatial autocorrelation statistics I_i, c_i, G_i, and G_i*. We discuss two alternative randomization assumptions, total and conditional, and then newly derive expectations and variances under conditional randomization for I_i and c_i, as well as under total randomization for c_i. The four statistics are tested by a biological simulation model from population genetics in which a population lives on a 21 × 21 lattice of stepping stones (sixty-four individuals per stone) and reproduces and disperses over a number of generations. Some designs model global spatial autocorrelation, others spatially random surfaces. We find that spatially random designs give reliable test results by permutational methods of testing significance. Globally autocorrelated designs do not fit expectations by any of the three tests we employed. Asymptotic methods of testing significance failed consistently, regardless of design. Because most biological data sets are autocorrelated, significance testing for local spatial autocorrelation is problematic. However, the statistics are informative when employed in an exploratory manner. We found that hotspots (positive local autocorrelation) and coldspots (negative local autocorrelation) are successfully distinguished in spatially autocorrelated, biologically plausible data sets.     

What Were We Thinking?

This article outlines the context in geography and statistics in the mid 1960s, at the height of geography's so-called "quantitative revolution," that led us into a long-term collaboration about spatial statistics, which has continued in surges and lulls for some 40 years. We focus upon problems in spatial autocorrelation, including the measurement of autocorrelation, distribution theory, and variable geographical lattices. This narrative may not describe how it was, but it does describe how we remember the events of the time.     

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.     

Properties of Tests for Spatial Dependence in Linear Regression Models

Based on a large number of Monte Carlo simulation experiments on a regular lattice, we compare the properties of Moran's I and Lagrange multiplier tests for spatial dependence, that is, for both spatial error autocorrelation and for a spatially lagged dependent variable. We consider both bias and power of the tests for six sample sizes, ranging from twenty-five to 225 observations, for different structures of the spatial weights matrix, for several underlying error distributions, for misspecified weights matrices, and for the situation where boundary effects are present. The results provide an indication of the sample sizes for which the asymptotic properties of the tests can be considered to hold. They also illustrate the power of the Lagrange multiplier tests to distinguish between substantive spatial dependence (spatial lag) and spatial dependence as a nuisance (error autocorrelation).     

Application of LASSO to the Eigenvector Selection Problem in Eigenvector‐based Spatial Filtering

Eigenvector‐based spatial filtering is one of the often used approaches to model spatial autocorrelation among the observations or errors in a regression model. In this approach, a subset of eigenvectors extracted from a modified spatial weight matrix is added to the model as explanatory variables. The subset is typically specified via the selection procedure of the forward stepwise model, but it is disappointingly slow when the observations n take a large number. Hence, as a complement or alternative, the present article proposes the use of the least absolute shrinkage and selection operator (LASSO) to select the eigenvectors. The LASSO model selection procedure was applied to the well‐known Boston housing data set and simulation data set, and its performance was compared with the stepwise procedure. The obtained results suggest that the LASSO procedure is fairly fast compared with the stepwise procedure, and can select eigenvectors effectively even if the data set is relatively large (n = 10⁴), to which the forward stepwise procedure is not easy to apply.     

Evaluating the Temporal and Spatial Heterogeneity of the European Convergence Process, 1980–1999*

Abstract. In this paper, we suggest a framework that allows testing simultaneously for temporal heterogeneity, spatial heterogeneity, and spatial autocorrelation in β‐convergence models. Based on a sample of 145 European regions over the 1980–1999 period, we estimate a Seemingly Unrelated Regression Model with spatial regimes and spatial autocorrelation for two sub‐periods: 1980–1989 and 1989–1999. The assumption of temporal independence between the two periods is rejected, and the estimation results point to the presence of spatial error autocorrelation in both sub‐periods and spatial instability in the second sub‐period, indicating the formation of a convergence club between the peripheral regions of the European Union.     

A Wavelet‐Based Extension of Generalized Linear Models to Remove the Effect of Spatial Autocorrelation. 基于小波扩展广义线性模型消除空间自相关的影响

Biogeographical studies are often based on a statistical analysis of data sampled in a spatial context. However, in many cases standard analyses such as regression models violate the assumption of independently and identically distributed errors. In this article, we show that the theory of wavelets provides a method to remove autocorrelation in generalized linear models (GLMs). Autocorrelation can be described by smooth wavelet coefficients at small scales. Therefore, data can be decomposed into uncorrelated and correlated parts. Using an appropriate linear transformation, we are able to extend GLMs to autocorrelated data. We illustrate our new method, called the wavelet‐revised model (WRM), by applying it to multiple regression with response variables conforming to various distributions. Results are presented for simulated data and real biogeographical data (species counts of the plant genus Utricularia [bladderworts] in grid cells throughout Germany). The results of our WRM are compared with those of GLMs and models based on generalized estimating equations. We recommend WRMs, especially as a method that allows for spatial nonstationarity. The technique developed for lattice data is applicable without any prior knowledge of the real autocorrelation structure.     

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 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.     

A New Perspective about Moran's Coefficient: Spatial Autocorrelation as a Linear Regression Problem. Moran系数的新视角：空间自相关视为线性回归问题

The computation of Moran's index of spatial autocorrelation requires the definition of a spatial weighting matrix. The eigendecomposition of this doubly centered matrix (i.e., one that forces the sums of all rows and columns to equal zero) has interesting properties that have been exploited in various contexts: distribution properties of the Moran coefficient (MC), spatial filtering in linear models, generalized linear models, and multivariate analysis. In this article, this eigendecomposition is used to propose a new view of MC based on its interpretation in the simple context of linear regression. I use this interpretation to demonstrate the different properties of MC and also the inefficiency of this index in some situations involving simultaneous positive and negative spatial autocorrelation. I propose some new statistics and procedures for testing spatial autocorrelation, and conduct a simulation study to evaluate these new approaches.     

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.     

Symptoms of Instability in Models of Spatial Dependence

In recent years, there has been a growing interest in the problems caused by the existence of instability in cross-sectional regressions. The results about local autocorrelation measures are part of this debate, as are the proposals concerning the concept of geographically weighted regressions. This article also deals with the problem of stability (or the lack thereof), but focusing the discussion on the supposition of constancy in the parameter of spatial dependence. In most cases, this assumption is treated, with the risks that this involves, as a maintained hypothesis, which should be ascertained before continuing with the modeling exercise. In the article, we present a simple heterogeneity test for this type of parameters, based on the Lagrange Multiplier principle. To illustrate its use, we take the distribution of per capita income among the European regions as our discussion case. According to our results, there are clear signs of structural breaks in the spatial distribution of this variable and the scale factor and the autocorrelation coefficient appear to be principal actors.     

Beyond Moran's I: Testing for Spatial Dependence Based on the Spatial Autoregressive Model

The statistic known as Moran's I is widely used to test for the presence of spatial dependence in observations taken on a lattice. Under the null hypothesis that the data are independent and identically distributed normal random variates, the distribution of Moran's I is known, and hypothesis tests based on this statistic have been shown in the literature to have various optimality properties. Given its simplicity, Moran's I is also frequently used outside of the formal hypothesis-testing setting in exploratory analyses of spatially referenced data; however, its limitations are not very well understood. To illustrate these limitations, we show that, for data generated according to the spatial autoregressive (SAR) model, Moran's I is only a good estimator of the SAR model's spatial-dependence parameter when the parameter is close to 0. In this research, we develop an alternative closed-form measure of spatial autocorrelation, which we call APLE , because it is an approximate profile-likelihood estimator (APLE) of the SAR model's spatial-dependence parameter. We show that APLE can be used as a test statistic for, and an estimator of, the strength of spatial autocorrelation. We include both theoretical and simulation-based motivations (including comparison with the maximum-likelihood estimator), for using APLE as an estimator. In conjunction, we propose the APLE scatterplot, an exploratory graphical tool that is analogous to the Moran scatterplot, and we demonstrate that the APLE scatterplot is a better visual tool for assessing the strength of spatial autocorrelation in the data than the Moran scatterplot. In addition, Monte Carlo tests based on both APLE and Moran's I are introduced and compared. Finally, we include an analysis of the well-known Mercer and Hall wheat-yield data to illustrate the difference between APLE and Moran's I when they are used in exploratory spatial data analysis.     

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.     

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.     

A SPATIAL CLIFF‐ORD‐TYPE MODEL WITH HETEROSKEDASTIC INNOVATIONS: SMALL AND LARGE SAMPLE RESULTS*

ABSTRACT In this paper, we specify a linear Cliff‐and‐Ord‐type spatial model. The model allows for spatial lags in the dependent variable, the exogenous variables, and disturbances. The innovations in the disturbance process are assumed to be heteroskedastic with an unknown form. We formulate multistep GMM/IV‐type estimation procedures for the parameters of the model. We also give the limiting distributions for our suggested estimators and consistent estimators for their asymptotic variance‐covariance matrices. We conduct a Monte Carlo study to show that the derived large‐sample distribution provides a good approximation to the actual small‐sample distribution of our estimators.     

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.     

Extending Moran's Index for Measuring Spatiotemporal Clustering of Geographic Events

Moran's Index for spatial autocorrelation and localized index for spatial association have been widely applied in many research fields as the first step to explore and assess the spatial dependency in a set of geographic events. This article presents extensions to the equations for calculating global and localized spatial autocorrelation so to include the temporal attribute values of the geographic events being analyzed. The extended equations were successfully implemented and tested with a real world data set. In addition, simulated data sets were used to reveal how the extended equations performed. Beyond the usefulness of the extended equations, we suggest that care be taken with regard to assessing spatiotemporal patterns under the normality and randomization assumptions as different outcomes from different assumptions would require different approaches for interpretation.     

Maximum Getis–Ord Statistic Adjusted for Spatially Autocorrelated Data

Local statistics test the null hypothesis of no spatial association or clustering around the vicinity of a location. To carry out statistical tests, it is assumed that the observations are independent and that they exhibit no global spatial autocorrelation. In this article, approaches to account for global spatial autocorrelation are described and illustrated for the case of the Getis–Ord statistic with binary weights. Although the majority of current applications of local statistics assume that the spatial scale of the local spatial association (as specified via weights) is known, it is more often the case that it is unknown. The approaches described here cover the cases of testing local statistics for the cases of both known and unknown weights, and they are based upon methods that have been used with aspatial data, where the objective is to find changepoints in temporal data. After a review of the Getis–Ord statistic, the article provides a review of its extension to the case where the objective is to choose the best set of binary weights to estimate the spatial scale of the local association and assess statistical significance. Modified approaches that account for spatially autocorrelated data are then introduced and discussed. Finally, the method is illustrated using data on leukemia in central New York, and some concluding comments are made.