The Saddlepoint Approximation of Moran's I's and Local Moran's Ii's Reference Distributions and Their Numerical Evaluation

Global Moran's I and local Moran's I_i are the most commonly used test statistics for spatial autocorrelation in univariate map patterns or in regression residuals. They belong to the general class of ratios of quadratic forms for whom a whole array of approximation techniques has been proposed in the statistical literature, such as the prominent saddlepoint approximation by Offer Lieberman (1994). The saddlepoint approximation outperforms other approximation methods with respect to its accuracy and computational costs. In addition, only the saddlepoint approximation is capable of handling, in analytical terms, reference distributions of Moran's I that are subject to significant underlying spatial processes. The accuracy and computational benefits of the saddlepoint approximation are demonstrated for a set of local Moran's Ii statistics under either the assumption of global spatial independence or subject to an underlying global spatial process. Local Moran's I_i is known to have an excessive kurtosis and thus void the use of the simple approximation methods of its reference distribution. The results demonstrate how well the saddlepoint approximation fits the reference distribution of local Moran's I_i. Furthermore, for local Moran's I_i under the assumption of global spatial independence several algebraic simplifications lead to substantial gains in numerical efficiency. This makes it possible to evaluate local Moran's I_i's significance in large spatial tessellations.     

Modified Moran's I for Small Samples

The most common indicator used to measure spatial dependence is Moran's I proposed by statistician Patrick A. P. Moran in 1950. The index is simple to use and applies the principle of the Pearson correlation coefficient, although it incorporates a proximity measure between elements. However, Moran's I tends to underestimate real spatial autocorrelation when the number of locations are few. This study aims to present a modified version of Moran's I that can measure real spatial autocorrelation even with small samples and check for spatial dependence.     

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

Map Resolution and Spatial Autocorrelation

Moran's I, a measure of spatial autocorrelation, is affected by map resolution and map scale. This study uses a geographic information system (GIS) to examine the resolution effects. Empirical distribution of wildland fires in Idyllwild, California, and hypothetical distributions of ordered patterns are analyzed. The results indicate that Moran's I increases systematically with the resolution level. The resolution effects can be summarized by a log-linear function relating the I coefficients to resolution levels. Empirical tests that compare the distribution of fire activity in a vegetation map and in a topographic map confirm the resolution effects observed.     

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.     

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

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.     

A Note on the Extremities of Local Moran's Iis and Their Impact on Global Moran's I

By defining local Moran's I_i as a ratio of quadratic forms and making use of its overall additivity to match global Moran's I, we can identify spatial objects with a strong impact on global Moran's I. First, we concentrate on the spatial properties of local Moran's I_i expressed by the local linkage degree. Depending on whether we use the W- or C-coding of the spatial connectivity matrix, the variance of local Moran's I_i for a small local linkage degree will be either large or small. Note that spatial objects associated with a local Moran's I_i with a large variance affect the global statistic much more than spatial objects associated with a local Moran's I_i with a small variance. Counterintuitively, global Moran's I defined in the W-coding is most influenced by spatial objects with a small number of spatial neighbors. In contrast, spatial objects with a large number of spatial neighbors exert more impact on global Moran's I setup in the C-coding. Second, we investigate the impact of the empirical data on local Moran's I_i and show that local Moran's I_i will only be significant for extreme absolute residuals at and around the reference location. Clusters of average regression residuals cannot be detected by local Moran's I_i. Consequently, spatial cliques of extreme residuals contribute more to significance tests on global autocorrelation.     

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.     

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.     

Schooling the Quantum Generations: Textbooks and Quantum Cultures from the 1910s to the 1930s

Ever since Thomas Kuhn's influential The Structure of Scientific Revolutions (1962), textbooks have suffered a bad reputation. They have been accused of distorting—at times purportedly—history and of feeding students with an unacceptably simplified and optimistic view of science. This attitude started to change only in recent times. With the increase of attention paid not only to how theories are conceived, but also how they are practiced, disseminated, and appropriated, historians have rehabilitated textbooks as a legitimate site of knowledge production. In this paper, I adopt textbooks as an instrument to unfold multiple facets of the culture that allowed quantum physics to flourish between 1900 and the early 1930s. I organize the article around two stories about two major textbooks, i.e., Sommerfeld′s Atombau und Spektrallinien and Dirac′s Principles of Quantum Mechanics. I explore the complex pedagogical cultures underlying these two masterpieces and how they intersect local agendas.     

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.     

Network Autocorrelation in Transport Network and Flow Systems

The use of Moran's I to assess the existence of network autocorrelation in flows over the arcs of real (tangible) and abstract (intangible) networks is examined. Residuals of a migration model developed here reveal the presence of such autocorrelation or dependence. Two approaches for removing the observed dependence are examined.     

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.     

Comparison of Moran's I and Geary's c in Multivariate Spatial Pattern Analysis

This article compares multivariate spatial analysis methods that include not only multivariate covariance, but also spatial dependence of the data explicitly and simultaneously in model design by extending two univariate autocorrelation measures, namely Moran's I and Geary's c. The results derived from the simulation datasets indicate that the standard Moran component analysis is preferable to Geary component analysis as a tool for summarizing multivariate spatial structures. However, the generalized Geary principal component analysis developed in this study by adding variance into the optimization criterion and solved as a trace ratio optimization problem performs as well as, if not better than its counterpart the Moran principal component analysis does. With respect to the sensitivity in detecting subtle spatial structures, the choice of the appropriate tool is dependent on the correlation and variance of the spatial multivariate data. Finally, the four techniques are applied to the Social Determinants of Health dataset to analyze its multivariate spatial pattern. The two generalized methods detect more urban areas and higher autocorrelation structures than the other two standard methods, and provide more obvious contrast between urban and rural areas due to the large variance of the spatial component.     

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.     

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.     

A Geospatial Analysis of China's Exports, 1991-2008

Three Taiwan-based economists employ a range of exploratory spatial data analysis tools (e.g., Moran's I and LISA statistics) to investigate trends in the growth of China's exports over the period 1991-2008. A particular focus is on the detection of spatial correlations between China and 40 export destination countries in major world regions. Emphasis in the paper on the key years of 1991, 2001, 2006, and 2008 has enabled the authors to analyze the impacts on China's trade of such major events as the country's accession to the World Trade Organization and the global economic crisis of 2008-2009. The results of the spatial analysis reveal the continuing importance of the U.S. and Asian countries in China's export trade (despite changes in the character of trade relations) and identify the spatial outliers (e.g., in Latin America) that may serve as the basis for new export markets for China in the future.     

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.