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.     

The removal estimator: A ‘probable numbers’ statistic that requires no matching

The probable number of individuals (PNI) estimators all require that we be able to match left-right paired bones reliably. The problems of matching are outlined and ways of addressing them are reviewed. It is shown that archaeozoological studies in which reliable matching is difficult raise analytical problems similar to those faced by population biologists unable to use classical release-recapture methods. The removal estimator has been used to estimate the size of whale populations from incomplete census data. It can be generalized into archaeozoology to produce either a simple PNI statistic of the matching type, or a statistic that does not require reliable matching.     

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.     

Development and Evaluation of an Optimal Composite Estimator in Spatial Microsimulation Small Area Estimation

A range of data is of geographic interest but is not available at a small area level from existing data sources. Small area estimation (SAE) offers techniques to estimate population parameters of target variables to detailed scales based on relationships between those target variables and relevant auxiliary variables. The resulting indirect small area estimate can deliver a lower mean squared error compared to its direct survey estimate, given that variance can be reduced markedly even if bias increases. Spatial microsimulation SAE approaches are widely utilized but only beginning to engage with the potential of composite estimators that use a weighted combination of indirect and direct estimators to reduce further the mean squared error of the small area estimate compared to an indirect SAE estimator alone. This article advances these approaches by constructing for the first time in the microsimulation literature an optimal composite estimator for such SAE approaches in which the combining weight is calculated from the mean squared errors of the two estimators; thus, optimizing the reduction in MSE of the resulting small area estimates. This optimal composite estimator is demonstrated and evaluated in a model-based simulation study and application based on the real data.     

INFERENCE BASED ON ALTERNATIVE BOOTSTRAPPING METHODS IN SPATIAL MODELS WITH AN APPLICATION TO COUNTY INCOME GROWTH IN THE UNITED STATES*

ABSTRACT This study examines aggregate county income growth across the 48 contiguous states from 1990 to 2005. To control for endogeneity, we estimate a two‐stage spatial error model and implement a number of spatial bootstrap routines to infer parameter significance. Among the results, we find that outdoor recreation and natural amenities favor positive growth in rural counties and property taxes correlate negatively with rural growth. Comparing bootstrap inference with other models, including the recent General Moment heteroskedastic‐robust spatial error estimator, we find similar conclusions suggesting bootstrapping can be effective in spatial models where asymptotic results are not well established.     

Crime Risk Maps: A Multivariate Spatial Analysis of Crime Data

In crime analyses, maps showing the degree of risk help police departments to make decisions on operational matters, such as where to patrol or how to deploy police officers. This study statistically models spatial crime data for multiple crime types in order to produce joint crime risk maps. To effectively model and map the spatial crime data, we consider two important characteristics of crime occurrences: the spatial dependence between sites, and the dependence between multiple crime types. We reflect both characteristics in the model simultaneously using a generalized multivariate conditional autoregressive model. As a real‐data application, we examine the number of incidents of vehicle theft, larceny, and burglary in 83 census tracts of San Francisco in 2010. Then, we employ a Bayesian approach using a Markov chain Monte Carlo method to estimate the model parameters. Based on the results, we detect the crime hotspots, thus demonstrating the advantage of using a multivariate spatial analysis for crime data.     

DISCRETE SPATIAL CHOICE MODELS FOR AGGREGATE DESTINATIONS*

ABSTRACT. In problems of spatial choice, the choice set is typically more aggregated than the one considered by decision-makers, often because choice data are available only at the aggregate level. These aggregate choice units will exhibit heterogeneity in utility and in size. To be consistent with utility maximization, a choice model must estimate choice probabilities on the basis of the maximum utility within heterogeneous aggregates. The ordinary multinomial logit model applied to aggregate choice units fails this criterion as it is estimated on the basis of average utility. In this paper, we derive and discuss a model which utilizes the theory underlying the nested logit model to estimate the appropriate maximum utilities of aggregates. We also demonstrate that the aggregate alternative error terms are asymptotically Gumbel, thereby relaxing the assumption of extreme value distributed error terms. This is accomplished with help from the asymptotic theory of extremes.     

SELECTION BIAS IN SPATIAL ECONOMETRIC MODELS

ABSTRACT. The problem of spatial autocorrelation has been ignored in selection-bias models estimated with spatial data. Spatial autocorrelation is a serious problem in these models because the heteroskedasticity with which it commonly is associated causes inconsistent parameter estimates in models with discrete dependent variables. This paper proposes estimators for commonly-employed spatial models with selection bias. A maximum-likelihood estimator is applied to data on land use and values in 1920s Chicago. Evidence of significant heteroskedasticity and selection bias is found.     

Spatial Analysis of Employment and Population Density: The Case of the Agglomeration of Dijon 1999

The aim of this paper is to analyze the intraurban spatial distributions of population and employment in the agglomeration of Dijon (regional capital of Burgundy, France). We study whether this agglomeration has followed the general tendency of job decentralization observed in most urban areas or whether it is still characterized by a monocentric pattern. To that purpose, we use a sample of 136 observations at the communal and at the IRIS (infraurban statistical area) levels with 1999 census data and the employment database SIRENE (INSEE). First, we study the spatial pattern of total employment and employment density using exploratory spatial data analysis. Apart from the CBD, few IRIS are found to be statistically significant, a result contrasting with those found using standard methods of subcenter identification with employment cut‐offs. Next, in order to examine the spatial distribution of residential population density, we estimate and compare different specifications: exponential negative, spline‐exponential, and multicentric density functions. Moreover, spatial autocorrelation, spatial heterogeneity, and outliers are controlled for by using the appropriate maximum likelihood, generalized method of moments, and Bayesian spatial econometric techniques. Our results highlight again the monocentric character of the agglomeration of Dijon.     

Exploring scale-dependent correlations between cancer mortality rates using factorial kriging and population-weighted semivariograms

This article presents a geostatistical methodology that accounts for spatially varying population size in the processing of cancer mortality data. The approach proceeds in two steps: (1) spatial patterns are first described and modeled using population-weighted semivariogram estimators, (2) spatial components corresponding to nested structures identified on semivariograms are then estimated and mapped using a variant of factorial kriging. The main benefit over traditional spatial smoothers is that the pattern of spatial variability (i.e., direction-dependent variability, range of correlation, presence of nested scales of variability) is directly incorporated into the computation of weights assigned to surrounding observations. Moreover, besides filtering the noise in the data, the procedure allows the decomposition of the structured component into several spatial components (i.e., local versus regional variability) on the basis of semivariogram models. A simulation study demonstrates that maps of spatial components are closer to the underlying risk maps in terms of prediction errors and provide a better visualization of regional patterns than the original maps of mortality rates or the maps smoothed using weighted linear averages. The proposed approach also attenuates the underestimation of the magnitude of the correlation between various cancer rates resulting from noise attached to the data. This methodology has great potential to explore scale-dependent correlation between risks of developing cancers and to detect clusters at various spatial scales, which should lead to a more accurate representation of geographic variation in cancer risk, and ultimately to a better understanding of causative relationships.     

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 Geostatistical Analysis of Soil, Vegetation, and Image Data Characterizing Land Surface Variation

The elucidation of spatial variation in the landscape can indicate potential wildlife habitats or breeding sites for vectors, such as ticks or mosquitoes, which cause a range of diseases. Information from remotely sensed data could aid the delineation of vegetation distribution on the ground in areas where local knowledge is limited. The data from digital images are often difficult to interpret because of pixel-to-pixel variation, that is, noise, and complex variation at more than one spatial scale. Landsat Thematic Mapper Plus (ETM+) and Satellite Pour l'Observation de La Terre (SPOT) image data were analyzed for an area close to Douna in Mali, West Africa. The variograms of the normalized difference vegetation index (NDVI) from both types of image data were nested. The parameters of the nested variogram function from the Landsat ETM+ data were used to design the sampling for a ground survey of soil and vegetation data. Variograms of the soil and vegetation data showed that their variation was anisotropic and their scales of variation were similar to those of NDVI from the SPOT data. The short- and long-range components of variation in the SPOT data were filtered out separately by factorial kriging. The map of the short-range component appears to represent the patterns of vegetation and associated shallow slopes and drainage channels of the tiger bush system. The map of the long-range component also appeared to relate to broader patterns in the tiger bush and to gentle undulations in the topography. The results suggest that the types of image data analyzed in this study could be used to identify areas with more moisture in semiarid regions that could support wildlife and also be potential vector breeding sites.     

The Modifiable Areal Unit Problem,in a Regression Model Whose Independent Variable Is a Distance from a Predetermined Point

This paper deals with the modifiable areal unit problem in the context of a regression model where a dependent variable is an attribute value (say, income) of an atomic data unit (say, a household) and an independent variable is a distance from a predetermined point (say, a central business district) to the atomic data unit (a disaggregated model). We apply this disaggregated model to spatially aggregated data in which the dependent variable is the average income over a spatial unit and the independent variable is the average distance from each household in a spatial unit to the predetermined point (an aggregated model). First, estimating the slope coefficient by the least squares method, we prove that the variance of the estimator for the slope coefficient in the aggregated model is larger than that in the disaggregated model. Second, focusing on variations in the variance of the estimator for the slope coefficient in the aggregated model with respect to the number of zones, we obtain the number of zones in which the variance is close to that in the disaggregated model. Third, we obtain the zoning system that has the minimum variance for a fixed number of zones. We also calculate the maximum variance in order to examine the range of the variance.     

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.     

Geostatistical Prediction and Simulation of Point Values from Areal Data

The spatial prediction and simulation of point values from areal data are addressed within the general geostatistical framework of change of support (the term support referring to the domain informed by each measurement or unknown value). It is shown that the geostatistical framework (i) can explicitly and consistently account for the support differences between the available areal data and the sought-after point predictions, (ii) yields coherent (mass-preserving or pycnophylactic) predictions, and (iii) provides a measure of reliability (standard error) associated with each prediction. In the case of stochastic simulation, alternative point-support simulated realizations of a spatial attribute reproduce (i) a point-support histogram (Gaussian in this work), (ii) a point-support semivariogram model (possibly including anisotropic nested structures), and (iii) when upscaled, the available areal data. Such point-support-simulated realizations can be used in a Monte Carlo framework to assess the uncertainty in spatially distributed model outputs operating at a fine spatial resolution because of uncertain input parameters inferred from coarser spatial resolution data. Alternatively, such simulated realizations can be used in a model-based hypothesis-testing context to approximate the sampling distribution of, say, the correlation coefficient between two spatial data sets, when one is available at a point support and the other at an areal support. A case study using synthetic data illustrates the application of the proposed methodology in a remote sensing context, whereby areal data are available on a regular pixel support. It is demonstrated that point-support (sub-pixel scale) predictions and simulated realizations can be readily obtained, and that such predictions and realizations are consistent with the available information at the coarser (pixel-level) spatial resolution.     

Unconditional Maximum Likelihood Estimation of Linear and Log-Linear Dynamic Models for Spatial Panels

This article hammers out the estimation of a fixed effects dynamic panel data model extended to include either spatial error autocorrelation or a spatially lagged dependent variable. To overcome the inconsistencies associated with the traditional least-squares dummy estimator, the models are first-differenced to eliminate the fixed effects and then the unconditional likelihood function is derived taking into account the density function of the first-differenced observations on each spatial unit. When exogenous variables are omitted, the exact likelihood function is found to exist. When exogenous variables are included, the pre-sample values of these variables and thus the likelihood function must be approximated. Two leading cases are considered: the Bhargava and Sargan approximation and the Nerlove and Balestra approximation. As an application, a dynamic demand model for cigarettes is estimated based on panel data from 46 U.S. states over the period from 1963 to 1992.     

Modeling Interaction in a System of Markets

This paper develops a nested family of market interaction models for the general case of spatial or aspatial goods. Market structure effects are identified in the form of intraregional and extraregional competition. The competing destinations and competing central place models are shown to be special cases of the general modeling framework. An empirical example, using survey data collected from Chicago-area Management and Public Relations (SIC 874) firms, illustrates how the family of models may be used to identify market structure.     

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.     

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.     

A Spatial Prior for Bayesian Vector Autoregressive Models

In this paper we develop a Bayesian prior motivated by cross-sectional spatial autoregressive models for use in time-series vector autoregressive forecasting involving spatial variables. We compare forecast accuracy of the proposed spatial prior to that from a vector autoregressive model relying on the Minnesota prior and find a significant improvement. In addition to a spatially motivated prior variance as in LeSage and Pan (1995) we develop a set of prior means based on spatial contiguity. A Theil-Goldberger estimator may be used for the proposed model making it easy to implement.