BAYESIAN MODEL AVERAGING IN THE CONTEXT OF SPATIAL HEDONIC PRICING: AN APPLICATION TO FARMLAND VALUES

ABSTRACT Specification uncertainty arises in spatial hedonic pricing models because economic theory provides no guide in choosing the spatial weighting matrix and explanatory variables. Our objective in this paper is to investigate whether we can resolve uncertainty in the application of a spatial hedonic pricing model. We employ Bayesian Model Averaging in combination with Markov Chain, Monte Carlo Model Composition. The proposed methodology provides inclusion probabilities for explanatory variables and weighting matrices. These probabilities provide a clear indication of which explanatory variables and weighting matrices are most relevant, but they are case specific.     

INTERPRETING SPATIAL ECONOMETRIC ORIGIN‐DESTINATION FLOW MODELS

Spatial interaction or gravity models have been used to model flows that take many forms, for example population migration, commodity flows, traffic flows, all of which reflect movements between origin and destination regions. We focus on how to interpret estimates from spatial autoregressive extensions to the conventional regression‐based gravity models that relax the assumption of independence between flows. These models proposed by LeSage and Pace ( 2008 , 2009 ) define spatial dependence involving flows between regions. We show how to calculate partial derivative expressions for these models that can be used to quantify these various types of effect that arise from changes in the characteristics/explanatory variables of the model.     

The Determinants of Regional Freight Transport: A Spatial,Semiparametric Approach

In the context of modeling regional freight the four‐stage model is a popular choice. The first stage of the model, freight generation and attraction, however, suffers from three shortcomings: first of all, it does not take spatial dependencies among regions into account, thus potentially yielding biased estimates. Second, there is no clear consensus in the literature as to the choice of explanatory variables. Second, sectoral employment and gross value added are used to explain freight generation, whereas some recent publications emphasize the importance of variables which measure the amount of logistical activity in a region. Third, there is a lack of consensus regarding the functional form of the explanatory variables. Multiple recent studies emphasize nonlinear influences of selected variables. This article addresses these shortcomings by using a spatial variant of the classic freight generation and attraction models combined with a penalized spline framework to model the explanatory variables in a semiparametric fashion. Moreover, a Bayesian estimation approach is used, coupled with a penalized Normal inverse‐Gamma prior structure, to introduce uncertainty regarding the choice and functional form of explanatory variables. The performance of the model is assessed on a real‐world example of freight generation and attraction of 258 European NUTS‐2 level regions, covering 25 European countries.     

Bayesian Model Averaging for Spatial Econometric Models

We extend the literature on Bayesian model comparison for ordinary least-squares regression models to include spatial autoregressive and spatial error models. Our focus is on comparing models that consist of different matrices of explanatory variables. A Markov Chain Monte Carlo model composition methodology labeled MC ³ by Madigan and York is developed for two types of spatial econometric models that are frequently used in the literature. The methodology deals with cases where the number of possible models based on different combinations of candidate explanatory variables is large enough such that calculation of posterior probabilities for all models is difficult or infeasible. Estimates and inferences are produced by averaging over models using the posterior model probabilities as weights, a procedure known as Bayesian model averaging. We illustrate the methods using a spatial econometric model of origin–destination population migration flows between the 48 U.S. states and the District of Columbia during the 1990–2000 period.     

Potential Use of Radiometric Data for Wildlife Habitat Modelling

The potential use of existing radiometric data sets, previously collected for prospecting purposes, has very rarely been used as a variable predictor in wildlife habitat modelling. The utility of radiometric data for predicting vegetation community patterns and wildlife habitat was investigated in the Australian arid zone using the Burt Plain bioregion as a case study. Using spatial datasets and a Species Distribution Modelling Toolkit, arid zone vertebrate species were modelled with Generalised Linear Modelling (GLM) regression modelling techniques. These models were used to predict the probability of occurrence of a species at any given location, defined in terms of its environmental attributes. A statistical correlation between the radioactive elements uranium, thorium and potassium, and terrain aspect was found. No statistical correlations were established between the radioactive elements and vegetation patterns; although we suspect these exist at finer scales of mapping. Radiometric data were identified as explanatory variables in the habitat models of all of the 32 vertebrate species examined, and used as illustration in the development of probabilistic spatial predictions of three species (Red Kangaroo, Macropus rufus; Lesser Hairy‐footed Dunnart, Sminthopsis youngsoni; and Rabbit, Oryctolagus cuniculus) in the bioregion. Our analyses suggest that radiometric data sets involving the radioactive elements: (uranium, thorium, and potassium), and vegetation could be used as predictors of biodiversity patterns at the bioregional and landscape level. This is an important finding given the challenges posed in undertaking broad‐scale biological surveys in the arid zone of Australia.     

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

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

Scale,Factor Analyses,and Neighborhood Effects

Studies of potential neighborhood effects have been constrained in most situations by the absence of small‐area data generated to characterize the local contexts within which individuals operate. Using small‐area data from the U.K. Census, this paper Illustrates the creation of bespoke neighborhoods—local areas defined separately for each individual in a sample survey—at a variety of scales, and their characterization using factor analysis techniques. Theories of neighborhood effects are uncertain as to the spatial scale at which the relevant processes operate, hence the value of exploring patterns consistent with those processes at a range of spatial scales. One problem with such comparative study is the incommensurability of regression coefficients derived from analyses using factor scores as the independent variables. The work reported here adapts a procedure introduced for reconstituting partial regression coefficients to circumvent that problem, and illustrates that patterns of voting at a recent British general election showed neighborhood‐effect‐like patterns at two separate scales simultaneously—with individual voter characteristics held constant.     

Une histoire de résidus : à propos des facteurs généraux et locaux de croissance régionale au Canada,de 1971 à 2001

Starting from an econometric model of local employment growth, applied to Canada (1971–2001), residuals—relative to model predictions—are analyzed over time and over space, in turn allowing us to draw a distinction between general explanatory variables and factors of a more local, cyclical or accidental nature. The model's explanatory power grows over time, founded on variables such as urban size, market access and industrial structure, allowing us to conclude that local employment growth in Canada follows an increasingly geographically predictable pattern. However, an examination of the residuals reveals more localized processes. Growth volatility is most manifest in Alberta and British Columbia, home to the most erratic local economies. Emerging patterns are visible in the last period, most notably the underperformance of Northern Ontario and of non‐metropolitan communities between Windsor and Québec City, lying along the Great Lakes and the Saint Lawrence. The over‐performance—compared to model predictions — of small and mid‐sized towns in south‐eastern Québec can, on the other hand, be interpreted as a sign of truly local social processes, generally associated with a particularly dynamic local entrepreneurial class.     

Model Uncertainty in Matrix Exponential Spatial Growth Regression Models

This article considers the most important aspects of model uncertainty for spatial regression models, namely, the appropriate spatial weight matrix to be employed and the appropriate explanatory variables. We focus on the spatial Durbin model (SDM) specification in this study that nests most models used in the regional growth literature, and develop a simple Bayesian model‐averaging approach that provides a unified and formal treatment of these aspects of model uncertainty for SDM growth models. The approach expands on previous work by reducing the computational costs through the use of Bayesian information criterion model weights and a matrix exponential specification of the SDM model. The spatial Durbin matrix exponential model has theoretical and computational advantages over the spatial autoregressive specification due to the ease of inversion, differentiation, and integration of the matrix exponential. In particular, the matrix exponential has a simple matrix determinant that vanishes for the case of a spatial weight matrix with a trace of zero. This allows for a larger domain of spatial growth regression models to be analyzed with this approach, including models based on different classes of spatial weight matrices. The working of the approach is illustrated for the case of 32 potential determinants and three classes of spatial weight matrices (contiguity‐based, k‐nearest neighbor, and distance‐based spatial weight matrices), using a data set of income per capita growth for 273 European regions.     

The Spatial Association of Demographic and Population Health Characteristics with COVID-19 Prevalence Across Districts in India

In less-developed countries, the lack of granular data limits the researcher's ability to study the spatial interaction of different factors on the COVID-19 pandemic. This study designs a novel database to examine the spatial effects of demographic and population health factors on COVID-19 prevalence across 640 districts in India. The goal is to provide a robust understanding of how spatial associations and the interconnections between places influence disease spread. In addition to the linear Ordinary Least Square regression model, three spatial regression models—Spatial Lag Model, Spatial Error Model, and Geographically Weighted Regression are employed to study and compare the variables explanatory power in shaping geographic variations in the COVID-19 prevalence. We found that the local GWR model is more robust and effective at predicting spatial relationships. The findings indicate that among the demographic factors, a high share of the population living in slums is positively associated with a higher incidence of COVID-19 across districts. The spatial variations in COVID-19 deaths were explained by obesity and high blood sugar, indicating a strong association between pre-existing health conditions and COVID-19 fatalities. The study brings forth the critical factors that expose the poor and vulnerable populations to severe public health risks and highlight the application of geographical analysis vis-a-vis spatial regression models to help explain those associations.     

THE SLX MODEL

We provide a comprehensive overview of the strengths and weaknesses of different spatial econometric model specifications in terms of spillover effects. Based on this overview, we advocate taking the SLX model as point of departure in case a well‐founded theory indicating which model is most appropriate is lacking. In contrast to other spatial econometric models, the SLX model also allows for the spatial weights matrix W to be parameterized and the application of standard econometric techniques to test for endogenous explanatory variables. This starkly contrasts commonly used spatial econometric specification strategies and is a complement to the critique of spatial econometrics raised in a special theme issue of the Journal of Regional Science (Volume 52, Issue 2). To illustrate the pitfalls of the standard spatial econometrics approach and the benefits of our proposed alternative approach in an empirical setting, the Baltagi and Li (2004) cigarette demand model is estimated.     

Using Spatial Filters and Exploratory Data Analysis to Enhance Regression Models of Spatial Data

Residual spatial autocorrelation is a situation frequently encountered in regression analysis of spatial data. The statistical problems arising due to this phenomenon are well‐understood. Original developments in the field of statistical analysis of spatial data were meant to detect spatial pattern, in order to assess whether corrective measures were required. An early development was the use of residual autocorrelation as an exploratory tool to improve regression analysis of spatial data. In this note, we propose the use of spatial filtering and exploratory data analysis as a way to identify omitted but potentially relevant independent variables. We use an example of blood donation patterns in Toronto, Canada, to demonstrate the proposed approach. In particular, we show how an initial filter used to rectify autocorrelation problems can be progressively replaced by substantive variables. In the present case, the variables so retrieved reveal the impact of urban form, travel habits, and demographic and socio‐economic attributes on donation rates. The approach is particularly appealing for model formulations that do not easily accommodate positive spatial autocorrelation, but should be of interest as well for the case of continuous variables in linear regression.     

A Simulation Study on Specifying a Regression Model for Spatial Data: Choosing between Autocorrelation and Heterogeneity Effects

In this simulation study, regressions specified with autocorrelation effects are compared against those with relationship heterogeneity effects, and in doing so, provides guidance on their use. Regressions investigated are: (1) multiple linear regression, (2) a simultaneous autoregressive error model, and (3) geographically weighted regression. The first is nonspatial and acts as a control, the second accounts for stationary spatial autocorrelation via the error term, while the third captures spatial heterogeneity through the modeling of nonstationary relationships between the response and predictor variables. The geostatistical‐based simulation experiment generates data and coefficients with known multivariate spatial properties, all within an area‐unit spatial setting. Spatial autocorrelation and spatial heterogeneity effects are varied and accounted for. On fitting the regressions, that each have different assumptions and objectives, to very different geographical processes, valuable insights to their likely performance are uncovered. Results objectively confirm an inherent interrelationship between autocorrelation and heterogeneity, that results in an identification problem when choosing one regression over another. Given this, recommendations on the use and implementation of these spatial regressions are suggested, where knowledge of the properties of real study data and the analytical questions being posed are paramount.     

Weighted‐Average Least Squares Applied to Spatial Econometric Models: A Monte Carlo Investigation

Recently, model averaging techniques have been employed widely in empirical investigations as an alternative to the conventional model selection procedure, a procedure criticized because it disregards a major component of uncertainty, namely, uncertainty regarding the model itself, and, thus, it leads to the underestimation of uncertainty regarding the quantities of interest. Bayesian model averaging (BMA) is one of the most popular model averaging techniques. Some studies indicate that BMA has cumbersome aspects. One of the major practical issues of using BMA is its substantial computational burden, which obstructs the process of obtaining exact estimates. A simulation method, such as Markov chain Monte Carlo (MCMC), is required to resolve this problem. Weighted‐average least squares (WALS) estimation has been proposed as an alternative to BMA. The computational burden of WALS estimation is negligible; therefore, it does not require the MCMC method. Furthermore, WALS estimation has theoretical advantages over BMA estimation. This article presents two contributions to the WALS literature. First, it applies WALS to spatial lag/error models in order to consider spatial dependence. Second, it extends WALS in order to consider explicitly the problem of multicollinearity by employing the technique of principal component regression. The small sample properties of the estimators of the proposed models are examined using Monte Carlo experiments; the results of these experiments suggest that the standard WALS may produce biased estimates when the underlying data‐generating process is a spatial lag process. Results also indicate that when the correlation among the regressors is high, the standard WALS estimators may suffer from large variances and root mean squared errors. Both of these problems are significantly mitigated by using the proposed models. Las técnicas de promediado de modelos (model averaging) vienen siendo empleadas con creciente frecuencia en las investigaciones empíricas como una alternativa a los procedimientos convencionales de selección de modelos estadísticos. Dichos procedimientos convencionales han sido criticados por no tomar en cuenta un componente clave de la incertidumbre: la incertidumbre del modelo en sí, y por lo tanto, conducen a la subestimación de la incertidumbre en la cuantificación de las valores estimados. El promediado bayesiano de modelos (Bayesian Model Averaging‐BMA) es una de las técnicas de promediado más usadas. Algunos estudios indican que BMA tiene aspectos engorrosos: uno de los principales aspectos prácticos a considerar en su uso es su pesada carga computacional, la cual obstruye el proceso de obtención de estimaciones exactas. Esta limitación hace necesario el uso de métodos de simulación, como el de la cadena de Markov de Monte Carlo (Markov Chain Monte Carlo‐MCMC). La estimación de mínimos cuadrados usando un promediado ponderado (Weighted‐Average Least Squares‐WALS) ha sido propuesta como alternativa a BMA. La carga computacional de la estimación WALS es mínima y por lo tanto no requiere del uso de MCMC. Más aun, la estimación WALS posee ventajas teóricas sobre BMA. Este artículo presenta dos contribuciones a la literatura especializada de WALS. En primer lugar, aplica WALS a modelos espaciales tipo lag/error con el fin de incorporar la dependencia espacial. En segundo lugar, modifica el método WALS, a fin de considerar explícitamente el problema de la multicolinealidad entre variables mediante el empleo de la técnica de regresión de componentes principales (Principal Component Regression‐PCR). Luego los autores utilizan experimentos Monte Carlo para examinar las propiedades de tipo “muestra pequeña” (small simple) de los estimadores de los modelos propuestos. Los resultados de los experimentos sugieren que el método WALS estándar puede producir estimaciones sesgadas cuando el proceso generador de datos subyacente (Data Generating Process‐DGP) es un proceso de retardo espacial (Spatial Lag Process‐SLP). Los resultados también indican que cuando la correlación entre las variables es alta, los estimadores estándar de WALS pueden padecer de varianzas y errores cuadráticos medios (root mean squared errors‐RMSEs) atípicamente grandes. Ambos problemas son mitigados significativamente mediante el uso de los modelos propuestos en el presente artículo. 近来,模型平均技术作为与传统模型选择流程可替换的方法,在经验调查中得到广泛应用。传统的模型选择流程忽视了模型本身的不确定性,进而低估了感兴趣样本数量的不确定性而受到批评。贝叶斯模型平均技术(BMA)是最为流行的模型平均技术之一。但已有研究表明,BMA在某些方面较为繁琐复杂,一个最主要的问题是其巨大的计算负荷阻碍模型了精确估计的过程,因此需要利用马尔可夫‐蒙特卡洛(MCMC) 之类的模拟方法进行解决。加权平均最小二乘(WALS)估计可作为BMA的可替换方法,其优点在于计算负荷可以忽略不计,因此不需要采用MCMC方法解决计算负荷问题。此外,WALS估计相比于BMA估计在理论上有一定的优势。本文针对WALS的贡献有两点：将WALS应用于空间滞后/空间误差模型以考虑空间依赖性,并利用主成分回归(PCR)拓展WALS以明确考虑多重共线性问题。本文利用蒙特卡洛实验对所提模型估计的小样本特征进行测试,结果显示当潜在数据生成过程(DGP)是一个空间滞后过程时,标准WALS可能产生有偏估计;此外,当回归量的相关性较高时,标准WALS估计量可能有较大的方差和根均方差(RMSEs).而本文提出的加权平均最小二乘估计模型能很好地缓解这两个问题。     

Notes of papers read before the society through the Yangtse Gorges

Critics of spatial modeling make several damaging points: (1) results are often trivial; (2) they contradict a priori insights derived from theory or experience; and (3) they preclude an evaluation of the substantive coherence of theories. These critiques arise because analysts have selected criteria of statistical explanation over substantive interpretation. This paper reviews the trade‐off as exemplified by regression models. The paper then presents an alternative methodology, staged regression analysis, and discusses its three‐fold advantage over traditional regression models: (1) avoiding statistical‐substantive trade‐offs; (2) facilitating comparisons among spatial systems; and (3) evaluating and refining existing theory.     

Fitting Constrained Poisson Regression Models to Interurban Migration Flows

This paper demonstrates the effects of fitting singly and doubly constrained spatial interaction models using the Poisson regression approach. A large data set containing migration flows between labor market areas in Great Britain in 1970–71 is used. The results of fitting unconstrained, singly constrained, and doubly constrained models are compared with respect to goodness of fit and the interpretability of parameter estimates. The addition of other explanatory variables to the model is also explored.     

Duration Modeling of Spatial Point Patterns

This paper attempts to further the research by Odland and Ellis (1992) in applying event history methodology to the analysis of spatial point patterns (that is, event patterns). Its empirical focus is the event pattern derived from the adoption of an agricultural innovation, the Harvestore, in southern Ontario, Canada, from 1963 to 1986. Event history analysis involves the use of discrete-state, continuous-time stochastic models to investigate a temporal longitudinal record on discrete variables. Event history models are usually concerned with durations of time between events and the effects of intertemporal time dependencies on future event occurrences. As such, they are often referred to as duration models. Many of the methods used in event history analysis allow the use of other nonnegative interval measurements in place of standard temporal intervals to investigate a series of events. In particular, spatial intervals (or durations) of distances between events may also be accommodated by event history models. Our analysis extends the previous research of Odland and Ellis (1992) by using a wider range of parametric models to explore duration dependence, investigating the role of spatial censoring, and using a more extensive set of explanatory variables. In addition, simulation experiments and graphical tests are used to evaluate the empirical event pattern against one generated from Complete Spatial Randomness. Results indicate that the event pattern formed by the Harvestore adopter farms is clustered (that is, is described by positive duration dependency), the sales agent is a significant factor in the distribution of adopters, and that contrasting results are obtained from the analysis using censored data versus uncensored data.     

基于GWR模型的合肥居住地价影响因素研究

以合肥市主城区为例,基于2010-2014年居住用地的出让数据,运用地统计法、GWR模型等方法,对合肥市居住地价的空间异质性及其影响因素进行研究。研究表明：①合肥市居住地价的空间分布呈现出显著的多中心的空间结构,地价的峰值区分别以老城区、政务区天鹅湖及滨湖新区塘西河公园为中心呈现圈层式分布;②不同的地价影响因素表现出不同的空间分布特征,其中容积率对居住地价的贡献度空间差异最大,其次是宗地面积,主干路次之,交通站点对居住地价的贡献度最小;③厘清各影响因素对地价的作用机制,建立动态的数字地价模型,不仅能促进土地资源的集约利用,重塑城市的空间结构,而且能为城市整体价值的发挥提供重要的理论支撑。     

Frequencies of Rare Incisor Variations Reflect Factors Influencing Precontact Population Relationships in Mexico and the American Southwest

Mexican prehistory is characterized by migration and population isolation in its earliest phase, followed by later inter‐population cultural interactions, such as politics and trade. While shared, common morphological variants are often informative about recent population events, rare trait variants have not been widely investigated to see whether they might be informative about earlier events. Here, we consider populations that show several rare variants at frequencies high enough to warrant such an investigation. Examining past population structure can help us understand population interactions across different periods of time, levels of migration, and population isolation/drift. Multiple‐origin populations may have more variation, including more and higher rates of rare traits. We examined maxillary lateral incisors in 1117 dentitions from 76 samples of Preclassic/Classic and Postclassic pre‐European contact Mexico and the American Southwest for the prevalence of rare lateral incisor variants. Variants observed included barrel, congenital absence, interruption groove, mesial marginal bending, peg, reduced, supernumerary, and talon teeth. The most common variant overall is interruption groove. Central Mexico, Huasteca, and the Lowland Yucatan region samples present the highest overall number of incisor variants. Higher frequencies of single unusual traits are commonly found in samples from smaller populations, while samples from larger population centres show a greater range of these rare variants. We investigated to see whether the pattern of rare incisor variant frequencies reflected early or later population history by comparing similarity/distance matrices and factor model matrices using Mantel tests and Generalized Procrustes analysis. Results show a weak relationship with the Postclassic period and shared migration/language and trade interaction models. We suggest that later cultural interactions have acted to mask earlier population history. Research on serial founder effects should be approached multiregionally and across time, to avoid missing inter‐region biological relationships. Copyright © 2015 John Wiley & Sons, Ltd.     

The population ecology of hunter-gatherers and their prey

In this paper we combine foraging theory and population biology models to simulate dynamic relationships between hunter-gatherers and their prey resources. Hunter-gatherer population growth responds to the net marginal rate of foraging; prey population growth responds logistically to exploitation. Thus conceived, the relationship between forager and prey biomass is time-dependent and nonlinear. It changes from stable equilibrium to damped and stable cycles with modest adjustments of input parameters. And, it produces the largest sustainable human population at intermediate levels of individual work effort. At equilibrium the forager takes all prey types with a pursuit and handling rate greater than or equal to its maintenance foraging rate. The structural properties of the model compel us to reject standard anthropological interpretations of the carrying capacity concept; they provide new insights on old issues such as original affluence and intensification. Analysis of the interaction of human population, diet selection, and resource depletion requires microecological models in part because the relevant processes occur on time scales largely invisible to both ethnography and archaeology.