A Structural Equation Approach to Models with Spatial Dependence

We introduce the class of structural equation models (SEMs) and corresponding estimation procedures into a spatial dependence framework. SEM allows both latent and observed variables within one and the same (causal) model. Compared with models with observed variables only, this feature makes it possible to obtain a closer correspondence between theory and empirics, to explicitly account for measurement errors, and to reduce multicollinearity. We extend the standard SEM maximum likelihood estimator to allow for spatial dependence and propose easily accessible SEM software like LISREL 8 and Mx. We present an illustration based on Anselin's Columbus, OH, crime data set. Furthermore, we combine the spatial lag model with the latent multiple-indicators–multiple-causes model and discuss estimation of this latent spatial lag model. We present an illustration based on the Anselin crime data set again.     

A Generalized Randomization Approach to Local Measures of Spatial Association

This article establishes a unified randomization significance testing framework upon which various local measures of spatial association are commonly predicated. The generalized randomization approach presented is composed of two testing procedures, the extended Mantel test and the generalized vector randomization test. These two procedures employ different randomization assumptions, namely total and conditional randomization, according to the way in which they incorporate local measures. By properly specifying necessary matrices and vectors for a particular local measure of spatial association under a particular randomization assumption, the generalized randomization approach as a whole yields a reliable set of equations for expected values and variances, which then is confirmed by a Monte Carlo simulation utilizing random permutations.     

A Spatial Logit Association Model for Cluster Detection

In this paper, I propose to set out a logit spatial association model for binary spatial events and develop a scan algorithm to search for spatial associations. I extend the traditional logit model with a spatial autocorrelated component so that the model includes not only known risk factors, but also spatially autocorrelated regions as control or explanatory factors. The case study of West Virginia lung cancer shows that the model effectively captures cool and hot spots in lung cancer mortality.     

Assessing Spatial Dependence in Count Data: Winsorized and Spatial Filter Specification Alternatives to the Auto-Poisson Model

The auto-Poisson probability model furnishes an obvious tool for modeling counts of geographically distributed rare events. Unfortunately, its original specification can accommodate only negative spatial autocorrelation, which itself is a rare event. More recent alternative reformulations, namely, the Winsorized and spatial filter specifications, circumvent this drawback. A comparison of their performances presented in this article reveals some of their relative advantages and disadvantages.     

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.     

A Surface-Based Approach to Measuring Spatial Segregation

Quantitative indices of residential segregation have been with us for half a century, but suffer significant limitations. While useful for comparison among regions, summary indices fail to reveal spatial aspects of segregation. Such measures generally consider only the population mix within zones, not between them. Zone boundaries are treated as impenetrable barriers to interaction between population subgroups, so that measurement of segregation is constrained by the zoning system, which bears no necessary relation to interaction among population subgroups. A segregation measurement approach less constrained by the chosen zoning system, which enables visualization of segregation levels at the local scale and accounts for the spatial dimension of segregation, is required. We propose a kernel density estimation approach to model spatial aspects of segregation. This provides an explicitly geographical framework for modeling and visualizing local spatial segregation. The density estimation approach lends itself to development of an index of spatial segregation with the advantage of functional compatibility with the most widely used index of segregation (the dissimilarity index D ). We provide a short review of the literature on measuring segregation, briefly describe the kernel density estimation method, and illustrate how the method can be used for measuring segregation. Examples using a simulated landscape and two empirical cases in Washington, DC and Philadelphia, PA are presented.     

A Programming Approach to Minimizing and Maximizing Spatial Autocorrelation Statistics

A programming approach is presented for identifying the form of the weights matrix W which either minimizes or maximizes the value of Moran's spatial autocorrelation statistic for a specified vector of data values. Both nonlinear and linear programming solutions are presented. The former are necessary when the sum of the links in W is unspecified while the latter can be used if this sum is fixed. The approach is illustrated using data examined in previous studies for two variables measured for the counties of Eire. While programming solutions involving different sets of constraints are derived, all yield solutions in which the number of nonzero elements in W is considerably smaller than that in W defined using the contiguity relationships between the counties. In graph theory terms, all of the Ws derived define multicomponent graphs. Other characteristics of the derived Ws are also presented.     

Exploring Nonlinear Spatial Dependence in the Tails

In many instances it is of interest to measure the degree of similarity between neighboring regions. Spatial autocorrelation measures are the most popular means of doing it. However, such measures only capture a global linear relationship between regions, whereas in many circumstances a more general instrument is required. For instance, in economic poverty analysis or environmental applications (and in other cases where we are interested in extreme events and threshold exceedances) we should be more interested in the spatial pattern in the tails of the joint distributions. In this article we introduce some exploratory tool that focuses on the bivariate joint tails behavior to detect a pattern of spatial regularities. The method will be illustrated with reference to simulated environmental data.     

Gaussian Process Regression and Bayesian Model Averaging: An Alternative Approach to Modeling Spatial Phenomena

Gaussian Process Regression (GPR) is a nonparametric technique that is capable of yielding reliable out‐of‐sample predictions in the presence of highly nonlinear unknown relationships between dependent and explanatory variables. But in terms of identifying relevant explanatory variables, this method is far less explicit about questions of statistical significance. In contrast, more traditional spatial econometric models, such as spatial autoregressive models or spatial error models, place rather strong prior restrictions on the functional form of relationships, but allow direct inference with respect to explanatory variables. In this article, we attempt to combine the best of both techniques by augmenting GPR with a Bayesian Model Averaging (BMA) component that allows for the identification of statistically relevant explanatory variables while retaining the predictive performance of GPR. In particular, GPR‐BMA yields a posterior probability interpretation of model‐inclusion frequencies that provides a natural measure of the statistical relevance of each variable. Moreover, while such frequencies offer no direct information about the signs of local marginal effects, it is shown that partial derivatives based on the mean GPR predictions do provide such information. We illustrate the additional insights made possible by this approach by applying GPR‐BMA to a benchmark BMA data set involving potential determinants of cross‐country economic growth. It is shown that localized marginal effects based on partial derivatives of mean GPR predictions yield additional insights into comparative growth effects across countries.     

Semiparametric Maximum Likelihood Estimates of Spatial Dependence

We semiparametrically model spatial dependence via a combination of simpler weight matrices (termed spatial basis matrices) and fit this model via maximum likelihood. Estimation of the model relies on the intuition that bounds to the log‐determinant term in the log‐likelihood can provide penalties to overfitting both the level and pattern of spatial dependence. By relying on symmetric and doubly stochastic spatial basis matrices that reflect different weight specifications assigned to neighboring observations, we are able to derive a mathematical expression for bounds on the log‐determinant term that appears in the likelihood function. These bounds can be conveniently calculated allowing us to solve for maximum likelihood estimates at the bounds using a simple optimization over two quadratic forms that involve small matrices. An intuitively pleasing aspect of our approach is that the objective function for the bounded log‐likelihoods contains one quadratic form equal to the sum‐of‐squared errors measuring the quality of fit, and another quadratic form reflecting a penalty to overfitting spatial dependence. We apply our semiparametric estimation method to a housing model using 57,647 U.S. census tracts.     

Re-investing in the Built Environment: A Spatial Network Approach

Before launching ambitious and expensive development programmes to induce new regional technology corridors and clusters, it is critical to appreciate existing spatial economic patterns in a region. Initial economic conditions drive location decisions of firms and a labour force such that any changes must intercede onto an existing landscape built for current economic conditions. This work adopts a simple regional economic model to integrate and review traditional and modern urban location theories in order to illustrate the power of initial conditions to determine a final result. A simple spatial dynamic simulation model captures many of the pertinent effects of real estate pricing patterns to frame both opportunities and constraints to re-shape an urban landscape. Attention to 'ground up' spatially correlated location patterns revealed in price data that suggests close attention to strategic zoning can have profound impacts on the success or failure of economic development. Relatively modest policy interventions that carefully utilize existing preferences for urban amenities and concurrent real property investments involve fewer policy risks with potentially more powerful stimulative economic consequences than promised in more ambitious programmes.     

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.     

A Spatial Autoregressive Poisson Gravity Model

In this article, a Poisson gravity model is introduced that incorporates spatial dependence of the explained variable without relying on restrictive distributional assumptions of the underlying data‐generating process. The model comprises a spatially filtered component—including the origin‐, destination‐, and origin‐destination‐specific variables—and a spatial residual variable that captures origin‐ and destination‐based spatial autocorrelation. We derive a two‐stage nonlinear least‐squares (NLS) estimator (2NLS) that is hetero‐scedasticity‐robust and, thus, controls for the problem of over‐ or underdispersion that often is present in the empirical analysis of discrete data or, in the case of overdispersion, if spatial autocorrelation is present. This estimator can be shown to have desirable properties for different distributional assumptions, like the observed flows or (spatially) filtered component being either Poisson or negative binomial. In our spatial autoregressive (SAR) model specification, the resulting parameter estimates can be interpreted as the implied total impact effects defined as the sum of direct and indirect spatial feedback effects. Monte Carlo results indicate marginal finite sample biases in the mean and standard deviation of the parameter estimates and convergence to the true parameter values as the sample size increases. In addition, this article illustrates the model by analyzing patent citation flows data across European regions. En el presente artículo, se introduce un modelo de gravedad Poisson, que incorpora la dependencia espacial de la variable explicada, sin apoyarse en presunciones de distribución restrictivas del proceso subyacente de generación de datos. El modelo comprende de un componente espacialmente filtrado, que incluye las variables de origen, destino y origen‐destino específico; y una variable espacial residual que captura la auto‐correlación espacial basada en el origen y destino. Se deriva del calculador (2NLS) de dos etapas no lineales de mínimos cuadrados (NLS), el cual es robusto en heterocedasticidad, y por ello controla el problema de sobre‐dispersión o baja‐dispersión (over and under dispersion), que a menudo se presenta en el análisis empírico de datos discretos; o, en el caso de de sobre‐dispersión, cuando se presenta la auto correlación espacial. Este calculador puede demostrar tener propiedades deseables para diferentes supuestos distribucionales, como los flujos observados un componente (espacialmente) filtrado, ya sea Poisson o binomial negativo. En nuestra especificación de modelo espacial auto regresivo (SAR), las estimaciones de los parámetros resultantes se pueden interpretar como los efectos de impacto total implícitos, definidos como la suma de efectos espaciales, directos o indirectos, de retroalimentación (feedback). Los resultados Monte Carlo indican sesgos marginales de muestras finitas en la media y la desviación estándar de los parámetros estimados, y la convergencia de los valores de los parámetros reales, a medida que aumenta el tamaño de muestra. Este artículo ilustra el modelo mediante el análisis de flujos de datos de citas de patentes, a través de las regiones europeas. 本文提出了一种蕴含空间依赖的泊松引力模型,该模型中解释变量无需依赖潜在数据生成过程的限制性分布假设。该模型由包含起点、终点、起点‐终点特定变量的空间滤波组分和空间残差变量组成,能捕捉到基于起点和终点的空间自相关。我们推导出一个二阶非线性最小二乘（NLS）估计（2NLS）,它对异方差具有鲁棒性,从而可控制对于离散或过离散数据经验性分析中经常出现的过离散和低离散问题。如果空间自相关存在,过离散数据分析就是一个例子。对于不同的分布假设,如或泊松分布或是负二项式分布的观测流或（空间）滤波组分,该估计量显示出令人满意的性能。在本文的空间自回归（SAR）模型设定中,参数估计结果可解释为隐含的全局影响效应,并可被定义为直接和间接的空间反馈效应之和。蒙特卡罗结果给出了参数估计中均值、标准差的临界有限样本偏差,且随样本量增大收敛于真正参数值。此外,本文基于欧洲地区专利引用的流数据进行了模型验证。     

A Spatial Search-Location Model of Retail Centers

Central place theory can be modeled using two types of optimization problems. Location-allocation models have been widely applied to operationalize central place theory as an aggregate optimization problem. This paper constructs a spatial search-location model to formulate central place theory as an individual optimization problem.     

A Spatial Mixture Model of Innovation Diffusion

The diffusion of new product or technical innovation over space is here modeled as an event‐based process in which the likelihood of the next adopter being in region r is influenced by two factors: (i) the potential interactions of individuals in r with current adopters in neighboring regions, and (ii) all other attributes of individuals in r that may influence their adoption propensity. The first factor is characterized by a logit model reflecting the likelihood of adoption due to spatial contacts with previous adopters, and the second by a logit model reflecting the likelihood of adoption due to other intrinsic effects. The resulting spatial diffusion process is then assumed to be driven by a probabilistic mixture of the two. A number of formal properties of this model are analyzed, including its asymptotic behavior. But the main analytical focus is on statistical estimation of parameters. Here it is shown that standard maximum‐likelihood estimates require large sample sizes to achieve reasonable results. Two estimation approaches are developed which yield more sensible results for small sample sizes. These results are applied to a small data set involving the adoption of a new Internet grocery‐shopping service by consumers in the Philadelphia metropolitan area.     

A Generalization of the Utility Theory Approach To the Problem of Spatial Interaction

The derivation of the gravity model from utility theory is reformulated and generalized in order to show (1) how both the number of trips to individual destinations and the total travel budget may be determined simultaneously, and (2) how the effect of both distance and destination quality on trip distribution and total budget may be analyzed. Results are compared with the revealed space preference approach and are found to be superior in dealing theoretically with trip frequency and the effects of spatial context. For empirical applications, however, the revealed preference approach is advantageous.