A Geostatistical Linear Regression Model for Small Area Data. 一种适用于小区域数据的地统计线性回归模型

We present a new linear regression model for use with aggregated, small area data that are spatially autocorrelated. Because these data are aggregates of individual‐level data, we choose to model the spatial autocorrelation using a geostatistical model specified at the scale of the individual. The autocovariance of observed small area data is determined via the natural aggregation over the population. Unlike lattice‐based autoregressive approaches, the geostatistical approach is invariant to the scale of data aggregation. We establish that this geostatistical approach also is a valid autoregressive model; thus, we call this approach the geostatistical autoregressive (GAR) model. An asymptotically consistent and efficient maximum likelihood estimator is derived for the GAR model. Finite sample evidence from simulation experiments demonstrates the relative efficiency properties of the GAR model. Furthermore, while aggregation results in less efficient estimates than disaggregated data, the GAR model provides the most efficient estimates from the data that are available. These results suggest that the GAR model should be considered as part of a spatial analyst's toolbox when aggregated, small area data are analyzed. More important, we believe that the GAR model's attention to the individual‐level scale allows for a more flexible and theory‐informed specification than the existing autoregressive approaches based on an area‐level spatial weights matrix. Because many spatial process models, both in geography and in other disciplines, are specified at the individual level, we hope that the GAR covariance specification will provide a vehicle for a better informed and more interdisciplinary use of spatial regression models with area‐aggregated data.     

Inferences for Selected Location Quotients with Applications to Health Outcomes. 选择性区位商推断及其在医疗成效中的应用

The location quotient (LQ) is an index frequently used in geography and economics to measure the relative concentration of activities. This quotient is calculated in a variety of ways depending on which group is used as a reference. Here, we focus on a simultaneous inference for the ratios of the individual proportions to the overall proportion based on binomial data. This is a multiple comparison problem and inferences for LQs with adjustments for multiplicity have not been addressed before. The comparisons are negatively correlated. The quotients can be simultaneously tested against unity, and simultaneous confidence intervals can be constructed for the LQs based on existing probability inequalities and by directly using the asymptotic joint distribution of the associated test statistics. The proposed inferences are appropriate for analysis based on sample surveys. Two real data sets are used to demonstrate the application of multiplicity‐adjusted LQs. A simulation study is also carried out to assess the performance of the proposed methods to achieve a nominal coverage probability. For the LQs considered, the coverage of the simple Bonferroni‐adjusted Fieller intervals for LQs is observed to be almost as good as the coverage of the method that directly takes the correlations into account.     

An Empirically Based Implementation and Evaluation of a Hierarchical Model for Commuting Flows. 一种基于经验的层状通勤流模型实现与评估

This article provides an empirical evaluation of a hierarchical approach to modeling commuting flows. As the gravity family of spatial interaction models represents a benchmark for empirical evaluation, we begin by reviewing basic aspects of these models. The hierarchical modeling framework is the same that Thorsen, Ubøe, and Nævdal (1999) used. However, because some modifications are required to construct a more workable model, we undertake a relatively detailed presentation of the model, rather than merely referring to the presentation in Thorsen, Ubøe, and Nævdal (1999) . The model uses a hierarchical specification of a transportation network and the individual search procedure. Journeys to work are determined by the effects of distance deterrence and of intervening opportunities, and by the location of potential destinations relative to alternatives at subsequent levels in a transportation network. The model calibration uses commuting data from a region in western Norway. The estimated parameter values are reasonable, and the explanatory power is very satisfactory when compared with the results of a competing destinations approach.     

Determinants of Variation in Population–Employment Interaction Findings: A Quasi‐Experimental Meta‐Analysis. 人口‐就业相互作用的决定因素研究：一种准实验元分析法

This article reports about a metaregression analysis of empirical results generated using data for the northern Netherlands (1988–2002) in order to investigate the ambiguity in results in the population–employment interaction literature. Specifically, the analysis deals with the issue whether “jobs follow people” or “people follow jobs.” The article starts with introducing the basics of quasi‐experimental meta‐analysis and with identifying some advantages of using quasi‐experimental meta‐analysis as compared with the standard meta‐analysis approach. Two subsequent sections document the selection of the population–employment interaction model and salient characteristics of the data set as well as the setup of the primary analyses. A total of 4,050 quasi‐experimental empirical results for the jobs–people direction of causality are generated using different specifications and estimators for a spatial econometric interaction model. The subsequent metaregression analysis reveals that the empirical results are largely shaped by the spatial, temporal, and employment characteristics of the data sampling. The results also appear much more sensitive to different measurements of the model's key variables when compared with alternative specifications of the spatial weights matrix. The main determinant driving empirical results about jobs–people causality are differences in model specification and estimation, as revealed by an inherent bias in parameter estimates and misguided inferences for some of the commonly used specifications. Finally, suggestions for future research are identified.     

Multi‐Type,Multi‐Zone Facility Location

The placement of facilities according to spatial and/or geographic requirements is a popular problem within the domain of location science. Objectives that are typically considered in this class of problems include dispersion, median, center, and covering objectives—and are generally defined in terms of distance or service‐related criteria. With few exceptions, the existing models in the literature for these problems only accommodate one type of facility. Furthermore, the literature on these problems does not allow for the possibility of multiple placement zones within which facilities may be placed. Due to the unique placement requirements of different facility types—such as suitable terrain that may be considered for placement and specific placement objectives for each facility type—it is expected that different suitable placement zones for each facility type, or groups of facility types, may differ. In this article, we introduce a novel mathematical treatment for multi‐type, multi‐zone facility location problems. We derive multi‐type, multi‐zone extensions to the classical integer‐linear programming formulations involving dispersion, centering and maximal covering. The complexity of these formulations leads us to follow a heuristic solution approach, for which a novel multi‐type, multi‐zone variation of the non‐dominated sorting genetic algorithm‐II algorithm is proposed and employed to solve practical examples of multi‐type, multi‐zone facility location problems.     

Optimal Sampling Design for Variables with Varying Spatial Importance

It is often desirable to sample in those locations where uncertainty associated with a variable is highest. However, the importance of knowing the variable's value may vary across space. We are interested in the spatial distribution of Received Signal Strength Indicator (RSSI), a measure of the signal strength from a cell tower received at a particular location. It is crucial to estimate RSSI values accurately in order to evaluate the effectiveness of mayday systems designed for rapid emergency notification following vehicle crashes. RSSI estimation is less important for locations where the probability of a crash is low and where the likelihood of call completion is either close to zero or one. We develop a method for augmenting an initial spatial sample of RSSI values to achieve a high‐precision estimate of the probability of call completion following a crash. We illustrate the approach using data on RSSI and vehicle crashes in Erie County, NY.     

Growth and Convergence in a Multiregional Model with Space–Time Dynamics. 多区域时空动态模型的增长与收敛

The goal of this article is to test four distinct hypotheses about whether the relative location of an economy affects economic growth and economic well‐being using an extended Solow–Swan neoclassical growth model that incorporates both space and time dynamics. We show that the econometric specification takes the form of an unconstrained spatial Durbin model, and we investigate whether the results depend on some methodological issues, such as the choice of the time span and the inclusion of fixed effects. To estimate the fixed effects spatial Solow–Swan model, we adjust the Arrelano and Bond (1991) generalized method‐of‐moments (GMM) estimator to deal with endogeneity not only arising from the initial income level, as in the basic model, but also from the initial income levels and economic growth rates observed in neighboring economies.     

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.     

Gibrat Revisited: An Urban Growth Model Incorporating Spatial Interaction and Innovation Cycles. 再评吉尔布瑞特：顾及空间相互作用和创新周期的城市增长模型

We designed a geographical model for simulating the distribution of urban growth in systems of cities. The model incorporates the hierarchical and spatial diffusion of innovation cycles through gravitational interactions within a set of cities. Using theoretical simulations, we demonstrate that this model is able to reproduce the observed properties of urban systems for the log‐normal distribution of city sizes as well as the observed distribution of growth rates. Our experimentation was performed on a large harmonized historical database that includes a few hundred French urban agglomerations between 1831 and 1999 (Pumain‐INED database). Both spatial interaction and innovation cycles are necessary ingredients to explain the evolution of urban hierarchies. We suggest that Gibrat's generic stochastic growth model based on independent entities should be replaced by a more relevant model of spatially and temporally interdependent geographical entities.     

Aggregation of Regional Economic Time Series with Different Spatial Correlation Structures. 具有不同空间相关结构的区域经济时间序列集成

In this article, we compare the relative efficiency of different forecasting methods of space‐time series when variables are spatially and temporally correlated. We consider two cases: (1) univariate forecasting (i.e., a space‐time series aggregated into a single time series) and (2) the more general instance of multivariate forecasting (i.e., a space‐time series aggregated into a coarser spatial partition). We extend the results in the literature by including the consideration of larger datasets and the treatment of edge effects and of negative spatial correlation. We first introduce a statistical framework based on the space‐time autoregressive class of random field models, which constitutes the basis of our simulation study, and we present the various alternative forecasting methods considered in the simulation. We then present the results of a Monte Carlo study related to univariate forecasting. In order to allow a comparison with the findings of Giacomini and Granger (2004), we consider the same forecasting strategies and the same combinations of the parameter values used there, but with a larger parametric set. Finally, we extend our analysis to the case of multivariate forecasting. The outcomes obtained provide operational suggestions about how to choose between alternative forecasting methods in empirical circumstances.     

A Fractal Approach to Identifying Urban Boundaries. 城市边界识别的分形方法

Fractal geometry can be used for determining the morphological boundaries of metropolitan areas. A two‐step method is proposed here: (1) Minkowski's dilation is applied to detect any multiscale spatial discontinuity and (2) a distance threshold is located on the dilation curve corresponding to a major change in its behavior. We therefore measure the maximum curvature of the dilation curve. The method is tested on theoretical urban patterns and on several European cities to identify their morphological boundaries and to track boundary changes over space and time. Results obtained show that cities characterized by comparable global densities may exhibit different distance thresholds. The less the distances separating buildings differ between an urban agglomeration and its surrounding built landscape, the greater the distance threshold. The fewer the buildings that are connected across scales, the greater the distance threshold.     

A Framework for Exploring Urban Retail Discontinuities. 探索城市零售业不连续性的框架

A key area in the analysis of urban structural evolution is identifying discontinuities. Effective analysis could improve long‐term forecasting and provide a better understanding of how to steer an urban system toward a desirable future state. We use a simple aggregate retail model to demonstrate an algorithm for identifying discontinuities in model parameter space. Explorations of retailing in both Greater London and South Yorkshire in the United Kingdom illustrate how understanding a system's potential for discontinuity can provide insights for both policy makers and retail businesses. The Harris and Wilson model, described in the section so‐named, is used as a simple archetype to illustrate the new framework. This model can be developed in a straightforward way to incorporate further refinement. In “ Executing the model and visualizing the results ,” we describe a single model run and in “ Investigating discontinuities ,” we explain our framework for detecting and analyzing discontinuities. “ Identifying discontinuities in the London retail system ” shows the results of applying this methodology to the Greater London retail system, and in “ Practical applications ,” we explore the policy applications for this technique as related to the decline of town centers in the South Yorkshire retail system. Some concluding comments are offered in “ Conclusions .”     

Dealing with Timing and Synchronization in Opportunities for Joint Activity Participation. 公共活动参与机会中实时性与同步性处理

The ability of people to access opportunities offered by the built environment is circumscribed by various sets of space–time constraints, including the requirements to meet other persons at particular times and places to undertake activities together. While models of space–time accessibility recognize that joint activities may constrain the performance of activities in space and time, their specifications do not explicitly acknowledge the opportunities that individuals of a group have for joint activity participation. Therefore, this article focuses on joint activity participation and argues that collective activity decisions are the outcome of a complex process involving various aspects of timing, synchronization, and social hierarchy. The utility‐theoretic model proposed here quantifies the extent to which opportunities can be jointly accessed by a particular group of people within a specific time period. Central to the approach are three key variables: the attractiveness of an opportunity, the time available for activity participation, and the travel time to an activity location. Because of the multiperson character of joint activities, the determination of these variables is subject to individual preferences, privileges, and power differentials within a group. Specific attention is given to how time‐of‐day and synchronization effects influence the opportunities accessible to a group of individuals. The impact of these factors on joint accessibility is illustrated by a real‐world example of an everyday rendezvous scenario. The outcomes of a simulation exercise suggest that time‐of‐day and synchronization effects significantly affect the benefits that can be gained from opportunities for joint activities.     

Weighted Bidimensional Regression. 加权二维回归

Shape analysis is useful for a wide variety of disciplines and has many applications. One of the many approaches to shape analysis focuses on shapes that are represented by predefined landmarks on an object. Some landmarks may be measured with greater precision, exhibit more natural variation, or be more important than others to an analysis. This article introduces a method for including this information when estimating mapping relations or assessing the degree of similarity between two objects that are represented by a set of two‐dimensional landmarks. Weighted bidimensional regression combines aspects of weighted least squares regression and bidimensional regression as a way to weight variables that are represented by a set of two‐dimensional spatial coordinates. One possible weighting scheme is suggested, and the effect of weighting is demonstrated through a face‐matching application. Results indicate that appropriate weighting increases the ability to correctly match two faces and that weighting has the largest effect when used with a projective transformation.     

Dispersion of Nodes Added to a Network

For location problems in which optimal locations can be at nodes or along arcs but no finite dominating set has been identified, researchers may desire a method for dispersing p additional discrete candidate sites along the m arcs of a network. This article develops and tests minimax and maximin models for solving this continuous network location problem, which we call the added-node dispersion problem (ANDP). Adding nodes to an arc subdivides it into subarcs. The minimax model minimizes the maximum subarc length, while the maximin model maximizes the minimum subarc length. Like most worst-case objectives, the minimax and maximin objectives are plagued by poorly behaved alternate optima. Therefore, a secondary MinSumMax objective is used to select the best-dispersed alternate optima. We prove that equal spacing of added nodes along arcs is optimal to the MinSumMax objective. Using this fact we develop greedy heuristic algorithms that are simple, optimal, and efficient (O( mp )). Empirical results show how the maximum subarc, minimum subarc, and sum of longest subarcs change as the number of added nodes increases. Further empirical results show how using the ANDP to locate additional nodes can improve the solutions of another location problem. Using the p-dispersion problem as a case study, we show how much adding ANDP sites to the network vertices improves the p-dispersion objective function compared with (a) network vertices only and (b) vertices plus randomly added nodes. The ANDP can also be used by itself to disperse facilities such as stores, refueling stations, cell phone towers, or relay facilities along the arcs of a network, assuming that such facilities already exist at all nodes of the network.     

The Probabilistic Minisum Flow Interception Problem: Minimizing the Expected Travel Distance until Intercept under Probabilistic Interception

We develop a variant of the flow interception problem (FIP) in which it is more desirable for travelers to be intercepted as early as possible in their trips. In addition, we consider flows being intercepted probabilistically instead of the deterministic view of coverage assumed in the FIP literature. We call the proposed model the probabilistic minisum FIP (PMFIP); it involves minimizing the sum of the expected distance that each flow travels until intercepted at a facility among placed facilities. This extension allows us to evaluate the effect of facility location under any given value of the interception probability and to apply the model to a variety of situations. We apply the proposed model to an example network by assuming a hypothetical situation in which people gather at a stadium from various nodes on the network, and receive some goods or services on the way to the stadium. We analyze optimal solutions obtained by varying the number of facilities and interception probability. It is shown that the expected travel distance until intercept is greatly reduced by means of a few optimally located facilities under a moderate interception probability.     

On Covering Problems and the Central Facilities Location Problem

A number of variations of facilities location problems have appeared in the research literature in the past decade. Among these are problems involving the location of multiple new facilities in a discrete solution space, with the new facilities located relative to a set of existing facilities having known locations. In this paper a number of discrete solution space location problems are treated. Specifically, the covering problem and the central facilities location problem are shown to be related. The covering problem involves the location of the minimum number of new facilities among a finite number of sites such that all existing facilities (customers) are covered by at least one new facility. The central facilities location problem consists of the location of a given number of new facilities among a finite number of sites such that the sum of the weighted distances between existing facilities and new facilities is minimized. Computational experience in using the same heuristic solution procedure to solve both problems is provided and compared with other existing solution procedures.     

A Theory‐Based Cellular Automata for the Simulation of Land‐Use Change. 基于理论细胞自动机的土地利用变化模拟

This study proposes a theory‐based cellular automata (CA) in which CA rules are constructed on the basis of a theory of urban land‐use change and other local conditions. Rent gap theory that describes urban redevelopment and exogenous conditions that reflect urban development are embedded into a CA to simulate changes of land use. The theoretical framework provides a strong conceptual background for the CA simulation. A case study demonstrates the flexibility of the integrated framework for simulating land‐use changes in complex settings. This research aims to help decision‐makers formulate appropriate development plans.     

Nonparametric Significance Test for Discovery of Network‐Constrained Spatial Co‐Location Patterns

Spatial co‐location patterns are useful for understanding positive spatial interactions among different geographical phenomena. Existing methods for detecting spatial co‐location patterns are mostly developed based on planar space assumption; however, geographical phenomena related to human activities are strongly constrained by road networks. Although these methods can be simply modified to consider the constraints of networks by using the network distance or network partitioning scheme, user‐specified parameters or priori assumptions for determining prevalent co‐location patterns are still subjective. As a result, some co‐location patterns may be wrongly reported or omitted. Therefore, a nonparametric significance test without priori assumptions about the distributions of the spatial features is proposed in this article. Both point‐dependent and location‐dependent network‐constrained summary statistics are first utilized to model the distribution characteristics of the spatial features. Then, by using these summary statistics, a network‐constrained pattern reconstruction method is developed to construct the null model of the test, and the prevalence degree of co‐location patterns is modeled as the significance level. The significance test is evaluated using the facility points‐of‐interest data sets. Experiments and comparisons show that the significance test can effectively detect network‐constrained spatial co‐location patterns with less priori knowledge and outperforms two state‐of‐the‐art methods in excluding spurious patterns.