PROBIT WITH SPATIAL AUTOCORRELATION

ABSTRACT. Commonly-employed spatial autocorrelation models imply heteroskedastic errors, but heteroskedasticity causes probit to be inconsistent. This paper proposes and illustrates the use of two categories of estimators for probit models with spatial autocorrelation. One category is based on the EM algorithm, and requires repeated application of a maximum-likelihood estimator. The other category, which can be applied to models derived using the spatial expansion method, only requires weighted least squares.     

Lagrange Multiplier Test Diagnostics for Spatial Dependence and Spatial Heterogeneity

Several diagnostics for the assessment of model misspecification due to spatial dependence and spatial heterogeneity are developed as an application of the Lagrange Multiplier principle. The starting point is a general model which incorporates spatially lagged dependent variables, spatial residual autocorrelation and heteroskedasticity. Particular attention is given to tests for spatial residual autocorrelation in the presence of spatially lagged dependent variables and in the presence of heteroskedasticity. The tests are formally derived and illustrated in a number of simple empirical examples.     

基于空间自相关和时空扫描统计量的聚集比较分析

聚集是区域经济研究中的重点问题之一,而聚集定位问题又是深入分析聚集需要解决的首要问题。由于聚集具有高度的尺度敏感性,采用空间自相关方法分析时,尺度选择容易受研究者主观判断的影响,而且空间自相关方法也未考虑聚集的时间特征。与之相比,Kulldorff等学者提出的扫描统计量方法表现出了明显的优势。研究探索性地选用浙江省各市、县工业从业人口的聚集问题,从尺度选择、尺度转换和时空融合三个方面,比较了空间自相关和时空扫描统计量方法在探测聚集问题上的差异性,进而证实了时空扫描统计量方法不仅有效解决了人为选择尺度的偏倚问题,实现了尺度推绎、转换的自动化,而且更加有利地融合了立体、动态、多尺度的时空分析优势。     

BAYSIAN INFERENCE FOR ORDERED RESPONSE DATA WITH A DYNAMIC SPATIAL‐ORDERED PROBIT MODEL

ABSTRACT Many databases involve ordered discrete responses in a temporal and spatial context, including, for example, land development intensity levels, vehicle ownership, and pavement conditions. An appreciation of such behaviors requires rigorous statistical methods, recognizing spatial effects and dynamic processes. This study develops a dynamic spatial‐ordered probit (DSOP) model in order to capture patterns of spatial and temporal autocorrelation in ordered categorical response data. This model is estimated in a Bayesian framework using Gibbs sampling and data augmentation, in order to generate all autocorrelated latent variables. It incorporates spatial effects in an ordered probit model by allowing for interregional spatial interactions and heteroskedasticity, along with random effects across regions or any clusters of observational units. The model assumes an autoregressive, AR(1), process across latent response values, thereby recognizing time‐series dynamics in panel data sets. The model code and estimation approach is tested on simulated data sets, in order to reproduce known parameter values and provide insights into estimation performance, yielding much more accurate estimates than standard, nonspatial techniques. The proposed and tested DSOP model is felt to be a significant contribution to the field of spatial econometrics, where binary applications (for discrete response data) have been seen as the cutting edge. The Bayesian framework and Gibbs sampling techniques used here permit such complexity, in world of two‐dimensional autocorrelation.     

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.     

Toward a Theory of Spatial Statistics: Another Step Forward

An update is presented to an earlier piece outlining some steps that needed to be taken toward the establishment of a theory of spatial statistics. Findings that have appeared since this first paper are summarized and interpreted, and extensions and suggestions are offered for the further establishment of a basis for a theory of spatial statistics. Topics include boundary considerations, the role of latent spatial dependencies, and small-sample-size issues. These topics embrace the problems of data transformations, edge effect bias, reference sampling distributions, multivariate autocorrelation models, conditional expectations, and higher-order autoregressive structures. In part, a course is charted for the next step to be taken.     

The Majority Theorem for the Single (p = 1) Median Problem and Local Spatial Autocorrelation

Except for about a half dozen papers, virtually all (co)authored by Griffith, the existing literature lacks much content about the interface between spatial optimization, a popular form of geographic analysis, and spatial autocorrelation, a fundamental property of georeferenced data. The popular p-median location-allocation problem highlights this situation: the empirical geographic distribution of demand virtually always exhibits positive spatial autocorrelation. This property of geospatial data offers additional overlooked information for solving such spatial optimization problems when it actually relates to their solutions. With a proof-of-concept outlook, this paper articulates connections between the well-known Majority Theorem of the 1-median minisum problem and local indices of spatial autocorrelation; the LISA statistics appear to be the more useful of these later statistics because they better embrace negative spatial autocorrelation. The relationship articulation outlined here results in the positing of a new proposition labeled the egalitarian theorem.     

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.     

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.     

Robustness of Spatial Autocorrelation Specifications: Some Monte Carlo Evidence

Abstract This paper examines the robustness of various models of spatial autocorrelation through a series of Monte Carlo experiments in which each model takes a turn at the data generator. The generated data are then used to estimate all of the models. The estimated models are evaluated primarily on their predictive power.     

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.     

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.     

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.     

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.     

A General Misspecification Test for Spatial Regression Models: Dependence, Heterogeneity, and Nonlinearity

There is an increasing awareness of the potentials of nonlinear modeling in regional science. This can be explained partly by the recognition of the limitations of conventional equilibrium models in complex situations, and also by the easy availability and accessibility of sophisticated computational techniques. Among the class of nonlinear models, dynamic variants based on, for example, chaos theory stand out as an interesting approach. However, the operational significance of such approaches is still rather limited and a rigorous statistical-econometric treatment of nonlinear dynamic modeling experiments is lacking. Against this background this paper is concerned with a methodological and empirical analysis of a general misspecification test for spatial regression models that is expected to have power against nonlinearity, spatial dependence, and heteroskedasticity. The paper seeks to break new research ground by linking the classical diagnostic tools developed in spatial econometrics to a misspecification test derived directly from chaos theory—the BDS test, developed by Brock, Dechert, and Scheinkman (1987). A spatial variant of the BDS test is introduced and applied in the context of two examples of spatial process models, one of which is concerned with the spatial distribution of regional investments in The Netherlands, the other with spatial crime patterns in Columbus, Ohio.     

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.     

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

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

What Problem? Spatial Autocorrelation and Geographic Information Science

In retrospect it is the word "problem" in the title that seems most remarkable about the Cliff and Ord article. Spatial autocorrelation is indeed a problem in standard inferential statistics, which was developed to handle controlled experiments, when these methods are used to generalize from natural experiments. From the perspective of geographic information science, however, spatial dependence is a defining characteristic of geographic data that makes many of the functions of geographic information systems possible. The almost universal presence of spatial heterogeneity in such data also argues against generalization and is made explicit in the recent development of place-based analytic techniques. The final section argues for a new approach to the teaching of quantitative methods in the environmental and social sciences that treats natural experiments, spatial dependence, and spatial heterogeneity as the norm.