Spatial and Scale‐Dependent Controls on North American Pan‐Arctic Minimum River Discharge

Spatial patterns of minimum monthly river discharge in the North American Pan‐Arctic and its potential controls are explored with geographically weighted regression (GWR). Minimum discharge is indicative of soil water conditions; therefore, understanding spatial variability of its controls may provide insights into patterns of hydrologic change. Here, GWR models are applied to determine a suitable combination of independent variables selected from a set of eight variables. A model specification with annual mean river discharge, temperature at time of minimum discharge, and biome describes well the spatial patterns in minimum discharge. However, minimum discharge in larger watersheds is influenced more by temperature and biome distributions than it is in small basins, suggesting that scale is critical for understanding minimum river discharge. This study is the first to apply GWR to explore spatial variation in Pan‐Arctic hydrology. Factores de control espaciales y dependientes de escala en las descargas fluviales mínimas de ríos Pan‐Articos en Norteamérica. El artículo explora los patrones espaciales de caudales fluviales mínimos mensuales la región pan‐ártica de Norteamérica y sus posibles factores de control haciendo uso de una regresión ponderada geográficamente (geographically weigted regression‐GWR). Los caudales mínimos son indicadores de las condiciones del agua en el suelo, y por lo tanto el entendimiento de la variabilidad espacial de los factores que los controlan puede ayudar a comprender los patrones de cambio hidrológico. En el presente estudio, varios modelos de tipo GWR son aplicados para determinar una combinación adecuada de variables independientes seleccionadas a partir de un conjunto de ocho variables. El modelo que utiliza la media anual media de descarga fluvial, la temperatura en el momento de caudal mínimo, y el bioma, proporciona una buena descripción de los patrones espaciales en la descarga mínima. Sin embargo, en las cuencas hidrográficas grandes, la descarga mínima está más influenciada por la temperatura y la distribución de los biomas que en el caso de cuencas más pequeñas, lo que sugiere que la escala es fundamental para entender la descarga mínima fluvial. Este estudio es el primero en aplicar GWR para comprender la variación espacial en la hidrología de la región pan‐ártica. 基于GWR（地理加权回归模型）对北美泛北极地区月份最小河流流量的空间模式和潜在控制进行研究。最小流量暗示水土条件;因此,理解空间分异及控制可深刻理解水文变化的模式。GWR可从8个变量中提取一组独立变量的适当组合。通过年均河流流量、最小流量时的温度和生物群落,来描述最小下泄流量的空间格局。在大范围流域中,最小流量受到温度和生物群落分布的影响大于在小规模的流域,揭示出在河流最小流量分析中尺度是非常重要的。本文首次将GWR应用于泛北极水文空间异质性分析。     

Nonparametric Estimation of the Spatial Connectivity Matrix Using Spatial Panel Data

We use moments from the covariance matrix for spatial panel data to estimate the parameters of the spatial autoregression model, including the spatial connectivity matrix W. In the unrestricted spatial autoregression model, the parameters are underidentified by one when W is symmetric. We show that a special case exists in which W is asymmetric and its parameters are exactly identified. If the panel data are stationary and ergodic, spatially and temporally, the estimates of W and the spatial autoregression coefficients are consistent. Spatial panel data for house prices in Israel are used to illustrate this methodology. Los autores usan momentos de una matriz de covarianza para datos panel espaciales para estimar los parámetros del modelo de autoregresión espacial (spatial autoregressive model), incluyendo la matriz de conectividad (o de ponderación) espacial W. En el modelo de autoregresión espacial sin restricciones, los parámetros están sub‐identificados por un valor de uno en los casos que la matriz W es simétrica. Los autores demuestran que existe un caso especial en el cual la matriz W es asimétrica y sus parámetros tienen cálculo exacto. Si los datos panel son estacionarios y ergódicos, espacial y temporalmente, los estimados de W y el coeficiente de autoregresión espacial son consistentes. Para ilustrar la metodología propuesta, los autores usan datos‐panel espaciales de precios de vivienda en Israel. 本文通过采用空间面板数据的协方差矩阵对包含空间相关矩阵W的空间自回归模型进行参数的矩估计。在无约束空间自回归模型中,W是对称矩阵时,参数可由其估计得到。本文展示了一种W是对称矩阵且其参数能够被精确估计的特殊情况。如果面板数据在时间与空间特征上具有平稳性和遍历性,那么W和空间自回归参数的估计是一致的。最后,针对以色列住房价格的空间面板数据采用此方法进行实证研究。     

Distribution of the Sum of Distances to the First and Second Nearest Facilities

This article presents a derivation of the distribution of the sum of the distances to the first and second nearest facilities. Facilities are represented as points configured in regular and random patterns, and distance is measured as the Euclidean and rectilinear distances on a continuous plane. The sum of the distances represents the service level of facility location when customers are serviced by the first and second nearest facilities. Thus, the distribution of the sum is useful for facility location problems with nonclosest facility service. The distribution of the sum of road network distances also is calculated to evaluate the efficiency of actual facility locations. Este artículo presenta una derivación de la distribución de la suma de las distancias a la primera y segunda instalación más cercana. Las instalaciones son presentadas como puntos en configuraciones espaciales regulares y aleatorias, y la distancias son medidas usando el método euclidiano y el rectilínea sobre un plano continuo. La suma de las distancias representa el nivel de servicio en cada punto, en los casos en los que los clientes son atendidos en la instalación primera y segunda más cercana. El uso de la distribución de la suma de distancias es útil para resolver problemas de localización para casos en los que los clientes son atendidos por las instalaciones de servicios que nos son las más cercanas. Como demostración se calcula la distribución de la suma de las distancias de una red de carreteras para evaluar la eficiencia de la localización de instalaciones en un caso del mundo real. 本文讨论了一种最邻近距离总和与次邻近设施距离总和分布情况的推导方法。具体设施以规则和随机两种模式的点状要素来表示,距离采用连续平面中的欧氏及直线距离测度。距离之和反映了最邻近及次邻近设施为用户提供该设施区位的服务水平。因此,距离和的分布在无最邻近设施服务的设施选址问题中具有参考价值。路网距离和的分布也可计算用于现有设施区位的效率评估。     

Cognitive Value and Numerical Accuracy in Geocomputation: Choosing a Method to Compute the Log‐Jacobian Term

This article addresses the importance of cognitive value, numerical accuracy, and technical feasibility for geospatial models and methods. Researchers pursue large spatial data in an analysis to improve the statistical efficiency of parameter estimates and the cognitive value of the results. The critical requirement for this strategy is the availability of numerically accurate computational methods, without which parameter estimates lose their cognitive value. This article illustrates the relationship between cognitive value and numerical accuracy by establishing properties of popular polynomial schemes for approximating the log‐Jacobian term in spatial autoregressive models. It shows that numerical inaccuracies might be more misleading than is commonly assumed. En este artículo se aborda la importancia del valor cognitivo, la precisión numérica, y la viabilidad técnica de los modelos y métodos geoespaciales. Los investigadores interesados en análisis espacial tienden a preferir el uso de bases de datos de gran tamaño con el fin de mejorar la eficiencia estadística de sus estimaciones de parámetros y el valor cognitivo de los resultados. El requisito crítico para esta estrategia es la disponibilidad de métodos computacionales numéricamente precisos, sin los cuales las estimaciones de parámetros pierden su valor cognitivo. En este artículo se ilustra la relación entre el valor cognitivo y precisión numérica mediante la determinante jacobiano tipo log en modelos espaciales autorregresivos (spatial autorregresive models‐SAR). El estudio demuestra que las imprecisiones numéricas pueden ser más engañosas de lo que se asume comúnmente. 本文描述了地理空间模型和方法中认知价值、数值精度、技术可行性的重要性。研究者一般倾向于通过追求大规模空间数据,来提升参数估计的统计有效性和结果的认知价值。这一策略的关键是有实用的精确数值计算方法,否则参数估计将失去其认知价值。本文通过建立普遍使用的多项式方案来逼近空间自回归（SAR）模型中的Log‐Jacobian,从而说明认知价值和数值精度间的关系。结果表明数值不确定性可能比通常假设的误导性更大。     

Assessing Activity Pattern Similarity with Multidimensional Sequence Alignment Based on a Multiobjective Optimization Evolutionary Algorithm

Due to the complexity and multidimensional characteristics of human activities, assessing the similarity of human activity patterns and classifying individuals with similar patterns remains highly challenging. This article presents a new and unique methodology for evaluating the similarity among individual activity patterns. It conceptualizes multidimensional sequence alignment as a multiobjective optimization problem and solves this problem with an evolutionary algorithm (EA). The study utilizes sequence alignment to code multiple facets of human activities into multidimensional sequences and to treat similarity assessment as a multiobjective optimization problem that aims to minimize the alignment cost for all dimensions simultaneously. A multiobjective optimization evolutionary algorithm is used to generate a diverse set of optimal or near‐optimal alignment solutions. Evolutionary operators are specifically designed for this problem, and a local search method is also incorporated to improve the search ability of the algorithm. We demonstrate the effectiveness of our method by comparing it with a popular existing method called ClustalG using a set of 50 sequences. The results indicate that our method outperforms the existing method for most of our selected cases. The multiobjective EA presented in this article provides an effective approach for assessing activity pattern similarity and a foundation for identifying distinctive groups of individuals with similar activity patterns. Debido a la complejidad y el carácter multidimensional de las actividades humanas, la evaluación de la similitud de los patrones de la actividad humana y la clasificación de individuos con patrones similares sigue siendo una tarea muy difícil. Este artículo presenta una metodología novedosa y única para evaluar la similitud entre los patrones de actividad individual. Se conceptualiza el método de alineación de secuencias multidimensionales (Multidimensional Sequence Alignment‐MDSA) como un problema de optimización multiobjetivo y se resuelve mediante un algoritmo evolutivo (evolutionary algorithm‐EA). El estudio utiliza la alineación de secuencias para codificar múltiples facetas de las actividades humanas en secuencias multidimensionales y para tratar la evaluación de similitud como un problema de optimización multiobjetivo que tiene como objetivo reducir al mínimo el coste de alineación para todas las dimensiones simultáneamente . Los autores utilizan algoritmo evolutivo de optimización multiobjetivo (multiobjective optimization evolutionary algorithm‐MOEA) para generar un conjunto diverso de soluciones de alineación óptimas o casi óptimas. Los operadores evolutivos del algoritmo están diseñados específicamente para este problema, al cual se le incorpora un método de búsqueda local para mejorar su desempeño. Los autores demuestran la eficacia del método propuesto mediante su comparación al método existente ClustalG haciendo uso de un conjunto de 50 secuencias. Los resultados indican que el método propuesto supera a los métodos existentes para la mayoría de casos seleccionados. La EA multiobjetivo presentado en este artículo proporciona un método eficaz para evaluar similitudes del patrón de actividad y constituye un componente fundamental para la identificación de grupos distintivos de personas con patrones de actividad similares. 由于人类活动的复杂性和多维度特性,人类活动模式的相似性评估和基于相似模式的个体分类仍具有较高的挑战性。本文提出了一种新型的、独特的方法来评估个体活动模式的相似性。该研究利用序列比对将人类活动的多个方面编码成多维序列,并将活动模式相似性评估作为一个多目标优化问题,其目的是使所有维度的比对成本同时最低。多目标最优进化算法被用于生成一组不同的最优或接近最优的比对解决方案。为此特别设计了进化算子,并集成一种局部搜索方法提高了算法的整体搜索能力。为验证该方法的有效性,本文以50对序列作为实验,并与常用的ClustalG方法进行了对比。研究结果表明,该方法在多数选取案例中性能优于ClustalG方法。因此,本文提出的多目标进化（EA）算法为个体活动模式相似性评估提供了一种有效方法,也为具有相似活动模式个体中特定群体识别奠定了基础。     

Developing Spatial Weight Matrices for Incorporation into Multiple Linear Regression Models: An Example Using Grizzly Bear Body Size and Environmental Predictor Variables

In this study, we develop spatial autoregressive (SAR) models relating grizzly bear body length to environmental predictor variables in the Alberta Rocky Mountains. We examine the ability of several different spatial neighborhoods to model spatial dependence and compare the estimated parameters and residuals from a standard linear regression model (LRM) with those from three types of SAR models: error, lag, and Durbin. Further, we examine variable selection in the presence of negative dependence by repeating the modeling process using a SAR model. Two findings are that significant negative spatial dependence was present in the residuals of the LRM and that the choice of spatial neighborhood greatly affects the ability to detect spatial dependence. The incorporation of appropriate spatial weights into SAR models improves the fit and increases the significance of the parameter estimates vis‐à‐vis the linear model. The results of this study indicate that negative dependence may not have as severe negative effects on variable selection and parameter estimation as positive dependence. An examination of spatial dependence in regression modeling appears to be an important means of exploring the appropriateness of a sampling framework, predictor variables, and model form. En este estudio desarrollamos modelos espaciales autorregresivos (SAR) que vinculan la longitud del cuerpo de osos grizzli con variables predictivas ambientales en las montañas rocosas de Alberta, Canadá. Examinamos la capacidad de varias vecindades espaciales para modelar la dependencia espacial y la comparación de los parámetros estimados, así como los residuos de un modelo de regresión lineal estándar (LRM) versus tres tipos de modelos SAR: error, retraso (lag) y Durbin. Además, se examina la selección de variables en la presencia de dependencia negativa mediante la repetición del proceso de modelado con un modelo de SAR. El estudio concluye que: 1) existe dependencia espacial negativa significativa en los residuos de la LRM y; 2) la selección de la vecindad espacial afecta en gran medida la capacidad de detectar la dependencia espacial. La incorporación de ponderaciones espaciales correspondientes a los modelos SAR mejora el ajuste y aumenta la importancia de los parámetros estimados versus el modelo lineal. Los resultados de este estudio indican que la dependencia negativa puede no tener los graves efectos negativos en la selección de variables y la estimación de parámetros si se comparan dichos efectos con = la dependencia positiva. Los autores recomiendan un examen de la dependencia espacial en modelos de regresión como medio importante para explorar la conveniencia de un marco de muestreo, de variables de predicción, y de la forma del modelo. 本文构建了阿尔伯达省落基山脉地区的灰熊体态大小与环境预测变量之间的空间自回归模型（SAR）,检验了几种以不同空间邻域矩阵拟合变量的空间相关性,并比较了标准回归模型(LRM)与几种不同类型的SAR模型（空间残差模型、空间滞后模型和空间杜宾模型）的估计参数和残差大小。进而利用一种SAR模型重复模拟过程,进一步测试变量选择对负相关性存在的影响。研究表明,显著的空间负相关存在于LRM的残差中,且空间邻域权重的选择很大程度上影响模型空间相关性的探测能力。将适当的空间权重引入SAR模型中可提高拟合精度,增加相对于线性模型参数估计的显著性。研究结果表明,负相关性在变量选择和参数估计上严重负影响的程度不如正相关性强。回归模型中空间相关性检验似乎是采样结构、预测变量和模型形式适用性分析的一个重要途径。     

Hierarchical Bayesian Spatiotemporal Analysis of Childhood Cancer Trends

In this article, generalized additive mixed models are constructed for the analysis of geographical and temporal variability of cancer ratios. In this class of models, spatially correlated random effects and temporal components are adopted. Spatio‐temporal models that use intrinsic conditionally autoregressive smoothing across the spatial dimension and B‐spline smoothing over the temporal dimension are considered. We study the patterns of incidence ratios over time and identify areas with consistently high ratio estimates as areas for further investigation. A hierarchical Bayesian approach using Markov chain Monte Carlo techniques is employed for the analysis of the childhood cancer diagnoses in the province of Alberta, Canada during 1995–2004. We also evaluate the sensitivity of such analyses to prior assumptions in the Poisson context. En este artículo los autores construyen modelos aditivos generalizados mixtos (generalized additive mixed models) con el fin de analizar la variabilidad geográfica y temporal en las tasas de incidencia de cáncer. Este tipo de modelos emplean efectos aleatorios correlacionados espacialmente y componentes temporales. Los modelos espacio‐temporales emplean un suavizado condicional intrínseco autorregresivo (conditionally autoregressive smoothing) a través de la dimensión espacial y un suavizado de tipo B‐spline sobre la dimensión temporal. Los autores examinan los patrones de las tasas de incidencia a través del tiempo e identifican las áreas con valores consistentemente altos con el fin de sugerir áreas de investigación para el futuro. El estudio utiliza un enfoque jerárquico bayesiano (hierarchical bayesian) que usa una cadena de Markov Monte Carlo para evaluar los diagnósticos de cáncer infantil en la provincia de Alberta, Canadá durante el periodo 1995–2004. Asimismo, también se evalúa la sensibilidad de este tipo de análisis con respecto a los supuestos a‐priori, en el contexto de los modelos tipo Poisson. 本文提出了广义可加和混合模型进行癌症比率的地理和时间变化分析。在这类模型中引入了空间相关的随机效应和时间组分。时空模型在空间维度上采用本征自回归条件平滑,而在时间维度上则使用了B样条平滑。本文研究了疾病发生率模式,并识别出一直具有高比率估计的地区作为进一步调查区。在1995–2004年加拿大亚伯达省儿童癌症的诊断中,采用了基于马尔科夫链‐蒙特卡罗模型的分层贝叶斯方法,并且在泊松先验假设条件下评估了该类分析的敏感性。     

The Spatial Extent of Criminogenic Places: A Changepoint Regression of Violence around Bars

Crime scientists have long known that crime clusters near certain places such as drinking establishments, although the spatial parameters of that clustering are less established. This article proposes a methodology to estimate a distance beyond which there is significantly less evidence of a correlation between locations and concentrations of crime. The technique uses changepoints derived from a segmented regression applied to spatial buffers emanating from around particular crime‐generating land uses. Geographic information system techniques are used to create a series of buffers to determine the density of crime around sites. A changepoint Poisson regression of the buffer midpoints is used to estimate the distance beyond which crime densities do not appear to decline significantly with increasing distance. A case study of violent crime around 1,282 bars in Philadelphia, Pennsylvania, for 2008 reveals that violence is highly clustered within 25.9 m (85 feet) then dissipates rapidly, a pattern that is not replicated using control sites (fire stations). This is an estimate of the spatial extent of violence around bars, and the technique could be used to estimate the extent of other crimes around a variety of crime‐generating locations. Expertos en el estudio del crimen saben desde hace tiempo que los delitos violentos se concentran cerca de algunos lugares tales como establecimientos de bebidas, aunque los parámetros espaciales de dichas aglomeraciones son menos conocidos. Este artículo propone una metodología para estimar la distancia máxima a partir de la cual hay significativamente menos evidencia de una correlación entre puntos de interés y las aglomeraciones de crimen. La técnica empleada utiliza puntos de cambio (changepoints) derivados de una regresión segmentada (segmented regression) aplicada a las zonas de amortiguamiento (buffers) generadas en torno a usos del suelo particulares asociados a delincuencia. Técnicas SIG (Sistema de Información Geográfica) son utilizadas para crear una serie de buffers y determinar la densidad de delitos en torno a la ubicación de cada establecimiento (bar). Una regresión Poisson de tipo changepoint de los puntos medios de los buffers es empleada para estimar la distancia a partir de la cual las densidades del crimen no disminuyen significativamente con la distancia. Un estudio de caso de los delitos violentos en torno a 1.282 bares en Filadelfia, Pennsylvania en 2008 revela que la violencia está muy concentrada dentro de un radio de 25.9 m (85 pies) y luego se disipa rápidamente, un patrón que no se replica cuando el análisis es aplicado a sitios de control (estaciones de bomberos). El resultado es una estimación de la extensión espacial de la violencia alrededor de bares y la conclusión que la técnica podría ser utilizada para estimar la extensión de otros delitos en torno a una gran variedad de lugares asociados con la generación de la delincuencia. 犯罪学家早已明晰犯罪集聚于某些特定区域（如酒吧）的周围,尽管较少地构建这类聚集的空间参数。本文提出了一种方法可估算在一定距离之外,区位与犯罪集聚程度间相关性呈显著减少的证据。将从分段回归中获得的变异点应用于犯罪发生地的空间缓冲区。地理信息系统（GIS）技术用于产生一系列缓冲区以确定地点周围的犯罪密度。缓冲区中点的变异点泊松回归用于估算超出犯罪密度区不呈现随距离增加而显著衰退的距离。本文以宾夕法尼亚州费城1282个酒吧周围暴力犯罪为案例进行研究,揭示出2008年暴力犯罪集聚于25.9m的范围内,并在该距离之外的迅速消失,而当控制点选为消防站时该格局不再出现。实验表明,这是一种估算酒吧暴力犯罪空间范围的方法,并且该技术可用于估算不同类型犯罪产生地的距离范围。     

Exploratory Analysis of Polygons Distributed with Overlap

Spatial analysis often treats a set of polygons distributed with overlap in a region. Daily variation of the market area of a supermarket is represented by a set of overlapping polygons in geographic information systems and analyzed to grasp consumers' purchase behavior. The distribution of ethnic groups is also represented as a set of polygons, and relations among these polygons are discussed in cultural geography. To assist the exploratory analysis of polygons, this article proposes a new method of analyzing the relations among polygons. The method permits us to understand the entire structure of the relations among polygons and to detect local spatial patterns in these relations. The Hasse diagram and its extension are introduced to describe the topological structure of polygons. These tools are effective especially for describing hierarchical relations among polygons. The extension of the Hasse diagram also provides a means of classifying polygons, which is useful for discussing similarities and differences among them. The method is employed to analyze the results of an experiment in environmental psychology. Discussion focuses on technical soundness of and empirical findings obtained with the method. En análisis espacial, es común representar y analizar conjuntos de polígonos con traslapes en una región. La variación diaria de la zona de mercado de un supermercado, por ejemplo, es representada por medio de un conjunto de polígonos traslapados en un Sistema de Información Geográfico (GIS), el cual es analizado para comprender el comportamiento de compra de consumidores. La distribución de los grupos étnicos también es representada como un conjunto de polígonos, mientras que las relaciones entre dichos polígonos es abordado con conceptos de geografía cultural. Para apoyar al análisis exploratorio de los polígonos como éstos, los autores proponen un nuevo método de análisis de las relaciones entre polígonos. El método permite comprender toda la estructura de las relaciones entre polígonos y detectar patrones espaciales locales en estas relaciones. Los autores utilizan el diagrama de Hasse y proponen ciertas ampliaciones al método para describir la estructura topológica de polígonos. Estas herramientas son efectivas especialmente para la descripción de las relaciones jerárquicas entre los polígonos. La ampliación del diagrama de Hasse también proporciona un método para clasificar polígonos que es útil para tratar de similitudes y diferencias entre ellos. El método propuesto es empleado para analizar los resultados de un experimento en psicología ambiental. La discusión se centra en la solidez técnica y en los resultados empíricos obtenidos con el método. 空间分析通常将同个区域的一组多边形分布的图层进行叠加。在GIS系统中对一组超市市场区日常变化的多边形进行叠加分析可掌握消费者的购买行为。种族分布同样可由一组多边形表征,并可在文化地理学框架下进行这些多边形要素间的相关性分析。为支持多边形的探索性分析,本文提出了一种用于分析多边形间关系的新方法。该方法有助于理解多边形间关系的整体结构,并可探测局部空间格局。Hasse图及其扩展用于描述多边形的拓扑结构,对描述图层的层次关系尤为有效。Hasse图的扩展也提供了多边形分类的手段,这对于讨论它们之间的异同性是十分有效。最后,将该方法用于环境心理学试验结果的分析,讨论主要侧重于技术的合理性分析及从该方法得到的经验认识。     

Assessing Quality of Spatial Models Using the Structural Similarity Index and Posterior Predictive Checks

Model assessment is one of the most important aspects of statistical analysis. In geographical analysis, models represent spatial processes, where variability in mapped output results from uncertainty in parameter estimates. Slight spatial misalignments can cause inflated error scores when comparing maps of observed and predicted variables using traditional error metrics at the level of individual spatial units. We conceptualize spatial model assessment as a continuous value map comparison problem and employ methods from image analysis to score model outputs. The structural similarity index, a measure that attempts to replicate the human visual system using a local region approach, is used as an exploratory map comparison statistic. The measure is implemented within a Bayesian spatial modeling framework as a discrepancy measure in a posterior predictive check of model fit. Results are reported for simulation studies representing a variety of spatial processes in a spatial and space–time context. A case study of rainfall mapping in Sri Lanka demonstrates the proposed methodology applied to assessment of Bayesian kriging interpolations. Both simulation studies as well as the case study demonstrate that the approach reveals hidden spatial structure not uncovered by traditional methods. The spatially sensitive assessment methodology provides a diagnostic tool to support spatial modeling and analysis. La evaluación de modelos es uno de los aspectos más importantes de análisis estadístico. En el análisis geográfico, los modelos representan procesos espaciales en los que la variabilidad en los outputs es el resultado de la incertidumbre en los parámetros estimados. Leves desajustes espaciales pueden inflar los valores de error en la comparación entre los mapas de las observaciones y los mapas de las predicciones de las variables si es que se usan medidas tradicionales de medición de error al nivel de unidades espaciales individuales. Los autores conceptualizan la evaluación de modelos espaciales como un problema de comparación mapas de valor continuo y emplea métodos de análisis de imágenes para cuantificar los resultados del modelo. Se utiliza el índice de similitud estructural (SSIM), una medida que intenta replicar el sistema visual humano utilizando un enfoque de región local, como técnica de exploratoria comparación estadística de mapas. El índice es implementado dentro de un marco de modelización espacial bayesiano como medida de discrepancia en la comprobación posterior de predicción del desempeño del modelo. Los resultados se presentan para varios casos de simulación que representan una serie de procesos espaciales en un contexto espacio‐temporal y espacial. Un estudio de caso de mapeo de lluvias en Sri Lanka sirve como demostración de la metodología propuesta y su aplicación a la evaluación de las interpolaciones tipo krigeage (kriging) bayesianas. Tanto los estudios de simulación, así como el estudio de caso demuestran que el enfoque propuesto revela la estructura espacial oculta no evidenciada por métodos tradicionales. La metodología de evaluación espacialmente sensible que se presenta en este artículo proporciona una herramienta de diagnóstico para apoyar la elaboración de modelos y análisis espacial. 模型评估是统计分析中最为重要的内容之一。在地理分析中用模型表达空间过程,参数估计的不确定性会导致地图输出结果的可变性。当采用传统误差指标度量,在个体空间单位水平上进行观测和预测变量的地图比较时,微小的空间错位就可能导致误差的倍增。为此,本文通过将空间模型评估指标概念化为一个连续值图比较问题,并利用图像分析方法来评定模型输出。一种尝试以局域方法仿制人类视觉系统的度量指标——结构相似指数(SSIM),被用作为探索性地图的比较统计量。在贝叶斯空间模型框架下实现其量算,并将其作为一个偏差度量应用于模型拟合的后预测校验。仿真研究的结果显示出空间及时空环境下多类空间过程。以斯里兰卡降雨过程图为案例,展示了上述方法对贝叶斯克里格插值的准确性评估。仿真研究与实证结果均证明本文提出的方法可揭示以往传统方法掩盖的空间结构特征,空间敏感性评价为本研究的空间建模和分析提供了一个诊断工具。     

A Spatiotemporal Analysis of Groundwater Level Changes in Relation to Urban Growth and Groundwater Recharge Potential for Waukesha County,Wisconsin

The main objective of this study was to analyze spatially and temporally groundwater level changes using geographic information systems and spatial analysis with respect to urban development, groundwater withdrawal, and groundwater recharge potential. The study focused on Waukesha County in southeastern Wisconsin, where urban development has been accelerating while groundwater has been declining over the last several decades. We analyzed data about groundwater withdrawal, groundwater level, land use/land cover, and precipitation utilizing correlation analysis, geographically weighted regression, land use change analysis, and map overlay. Our findings include the following: the extent of urban areas expanded faster in regions with high recharge potential than in those with low or moderate recharge potential; the correlation of urban growth with groundwater level is highly variable over time and space; and changes in groundwater level are strongly related to the spatial distribution of groundwater withdrawal. The study pinpointed the need to consider the spatial unevenness of groundwater withdrawal in understanding the changes in groundwater level and groundwater recharge potential for better management of groundwater resources. Un análisis espacio‐temporal de los cambios de nivel de aguas subterráneas en relación con el crecimiento urbano y de las aguas subterráneas de recarga para el Condado de Waukesha, Wisconsin El objetivo principal de este estudio es analizar espacial y temporalmente los cambios del nivel de aguas subterráneas utilizando sistemas de información geográfica (SIG) y análisis espacial con respecto al desarrollo urbano, la extracción de agua subterránea, y el potencial de recarga. El área de estudio es el condado de Waukesha, en el sureste de Wisconsin, donde el desarrollo urbano se ha acelerado mientras que las aguas subterráneas han disminuido en las últimas décadas. Se analizaron los datos sobre la extracción y nivel de aguas subterráneas, la cobertura vegetal/uso del suelo, y la precipitación utilizando un análisis de correlación, una regresión ponderada geográficamente (geographically weigted regression ‐GWR), un análisis de cambio de uso del suelo, y la superposición de mapas (map overlay). Los resultados muestran que: (1) la extensión de las zonas urbanas se incrementó más rápidamente en regiones con alto potencial de recarga que en áreas con potencial de recarga de bajo o moderado; (2) la correlación entre el crecimiento urbano y el nivel de aguas subterráneas es muy variable tanto temporal como espacialmente; y, (3) los cambios en el nivel de aguas subterráneas están fuertemente relacionados con la distribución espacial de las actividades de extracción del recurso. El estudio identificó la necesidad de considerar la irregularidad del patrón espacial de la extracción de aguas subterráneas con el fin de mejorar la comprensión de los cambios en el nivel y el potencial de recarga y así promover una mejor gestión de los recursos hídricos subterráneos. 本研究主要目的在于通过地理信息系统(GISs)和空间分析的时空视角,对城市发展的地下水水位变化,地下水开采和地下水潜在补给进行分析。本文研究区为威斯康星州东南部的沃基肖县,该区在过去近几十年里城市快速增长而地下水急剧减少。通过地下水开采、地下水水位、LUCC（土地利用/土地覆盖LULC）和降水的数据分析,采用方法包括相关分析、地理加权回归模型（GWR）、土地利用变化分析和地图叠加分析。研究结果包括：（1）高地下水补给的城市地域范围的扩大幅度快于中等或低补给地区;（2）城市增长与地下水水位的相关性在时间和空间上是高度变化的;（3）地下水开采的空间分布与地下水位的变化呈强相关。研究精确地指出,为更好管理地下水资源,在理解地下水位变化和地下水潜在补给时需考虑地下水开采的空间不均匀。     

Food Inflation in the European Union: Distribution Analysis and Spatial Effects

In the European Union (EU), homogenous inflation forces are expected to prevail because of increased economic integration, especially after the creation of a single currency area. This expectation is directly related to the issue of inflation convergence, which has gained increasing attention by both academics and policy makers in Europe. Although the examination of core inflation is of great importance for macroeconomic policy, the role of disaggregate inflation indices, and especially food inflation, has also been emphasized in the literature. However, the issue of food inflation convergence has been largely ignored to date in empirical studies. This study explores the evolving distribution of food inflation rates in the EU‐25 member states using distribution dynamics analysis and covering the period from January 1997 to March 2011. This analysis assumes that each country represents an independent observation, providing unique information that can be used to estimate the transition dynamics of inflation. We show that spatial autocorrelation prevails inside the EU‐25, and, therefore, the independency assumption is violated. To ensure spatial independence, the Getis spatial filter is implemented prior to a distribution dynamics analysis. The results of this analysis confirm the existence of convergence trends, which are even clearer after the spatial filtering procedure, indicating, on the one hand, the influence of spatial effects on food inflation and, on the other hand, the effectiveness of the Getis spatial filtering technique. En la Unión Europea (UE), se espera que las fuerzas de inflación homogéneas prevalezcan debido a la mayor integración económica, sobre todo después de la creación de la zona de moneda única. Esta expectativa se relaciona directamente con el tema de la convergencia de la inflación, que ha ganado cada vez más atención por parte de los investigadores académicos y los decisores políticos europeos. Aunque el análisis de la inflación subyacente es de gran importancia para la política macroeconómica, el papel de los índices de inflación a niveles desagregados, sobre todo en el caso de la inflación de alimentos, ha sido un tema enfatizado por la literatura especializada. Sin embargo, la cuestión de la convergencia de la inflación de alimentos carece hasta la fecha de estudios empíricos. El artículo presente estudia la evolución de la distribución de las tasas de inflación de alimentos en los estados miembros de la UE‐25, utilizando el método de análisis de la dinámica de distribución (distribution dynamic analysis) y abarca el período comprendido entre enero de 1997 a marzo de 2011. Este análisis supone que cada país representa una observación independiente que proporciona información única que se puede utilizada para estimar la dinámica de transición inflacionaria. El presente estudio demuestra que la autocorrelación espacial prevalece dentro de los estados UE‐ 25, y por lo tanto la hipótesis de independencia estadística de las observaciones es violada. Para garantizar la independencia espacial, el método de filtrado espacial Getis (Getis Spatial Filter) es implementado antes de proceder con el análisis de la dinámica de distribución. Los resultados del análisis confirman la existencia de las tendencias de convergencia, que son aún más claras después de la aplicación del filtrado espacial. Estos resultados evidencian por un lado, la influencia de los efectos espaciales en la inflación de alimentos, y por otro lado, la eficacia de la técnica de filtrado espacial de Getis. 在欧盟中,由于区域经济一体化进程的推进,同质商品的通货膨胀盛行,特别是单一货币区域建立后,该趋势更为明显。这种演化态势和通货膨胀的收敛问题直接相关,已经引起了欧洲学术界和政策制定者的关注。虽然核心通胀的检测对宏观经济政策十分重要,但是分解通胀指数的作用,尤其是本文所强调的食品通胀也有文献中提及。然而,食品通胀的收敛问题却在实证研究中很大程度上被忽视。本文利用分布动力学方法对欧盟25国1997年1月至2011年3月的食品通胀率的演化分布进行分析。假设每个国家代表一个独立的观察个体并提供唯一的信息,可用于估计通胀的转变动态,研究发现空间自相关在欧盟25国中普遍存在,因此独立假设不成立。为了保证空间独立性,在进行分布动力学分析之前,先使用Getis空间滤波技术进行处理。分析结果证实了收敛趋势的存在,且该趋势在空间滤波程序处理后更为明显。它一方面显示出食品通胀空间效应的影响,另一方面表现出Getis空间滤波技术的有效性。     

Evaluating Polygon Overlay to Support Spatial Optimization Coverage Modeling

Minimizing costs and maximizing coverage are important goals in many planning contexts. These goals often necessitate an abstraction of a continuous demand region, resulting in potential errors when applying traditional coverage models. To reduce coverage errors caused by spatial abstraction, a number of spatial representation schemes have been proposed and applied. A new representation scheme using polygon overlay recently received much attention because potentially it can eliminate representation errors in coverage modeling. However, this overlay‐based approach is computationally challenging in terms of both the generation of demand units and the complexity of the resulting coverage model. This article investigates the operational and computational challenges of polygon overlay for delineating continuous demand in coverage models, an issue that has yet to be fully explored. We present a theoretical evaluation of the computational complexity associated with representation using polygon overlay in coverage modeling. Evaluations of two study regions provide empirical support for the computational complexity analysis. The analysis results provide insight regarding expected problem size and computational requirements if polygon overlay is relied upon to delineate demand unit boundaries in coverage modeling. La minimización de costos y la maximización de la cobertura espacial son objetivos importantes en muchos contextos de planificación. Estas metas a menudo requieren una abstracción de una región continua de demanda , dando lugar a posibles errores en la aplicación de modelos de cobertura tradicionales. Para reducir los errores de cobertura provocadas por la abstracción, la comunidad académica ha propuesto y aplicado una serie de esquemas de representación espacial. Recientemente un nuevo esquema de representación que utiliza la superposición de polígonos ha recibido mucha atención porque potencialmente puede eliminar los errores de representación en el modelado de la cobertura. Sin embargo, este enfoque es computacionalmente difícil, tanto en términos de la generación de unidades de demanda, como en la complejidad del modelo de cobertura resultante. Este artículo investiga los retos operacionales y de cómputo de la superposición de polígonos para delinear la región continua de demanda en los modelos de cobertura, un problema que aún no se ha explorado a fondo. Se presenta una evaluación teórica de la complejidad computacional asociada a la representación mediante superposición de polígonos en el modelado de cobertura espacial. Se presentan evaluaciones de dos regiones de estudio como apoyo empírico para el análisis de la complejidad computacional. Los resultados del análisis proporcionan información sobre el tamaño del problema esperado y los requerimientos computacionales en los casos en que el método de superposición de polígonos es usado para delinear límites de la región de demanda para el modelado de cobertura espacial 最小成本和最大区域覆盖是许多规划情境研究中的重要目标。实现这些目标通常需要对连续需求区域进行抽象,而这又会导致在应用传统覆盖模型时出现潜在误差。为减小由空间抽象引起的覆盖误差,已提出了一系列空间表达方案并得到应用。一种新型的利用多边形覆盖的表达方案,因其或可消除覆盖建模过程的表达误差,近来得到较多的关注。然而,这种基于覆盖的方法在需求单元生成与覆盖模型结果复杂性等方面面临着计算挑战。本文提出了一种在覆盖建模中采用多边形叠加表征的计算复杂度的理论评估方法。两个研究区域的评估为计算复杂度分析提供了经验支撑,如果在覆盖建模中多边形叠加依赖于描述需求单元边界时,该分析结果有助于深入考察预期问题规模大小及计算需求。     

A New Spatial Optimization Perspective on the Delineation of Metropolitan and Micropolitan Statistical Areas

The current method for delineating U.S. metropolitan and micropolitan statistical areas (collectively referred to as core‐based statistical areas—CBSAs) is to cluster counties based on the strength of commuting interactions between outlying areas and urban cores. The social and economic integration of a metropolitan or a micropolitan statistical area is operationalized using what we believe is an outdated, monocentric perspective on functional spatial structure: one that fixates solely on inward and reverse commuting. We propose a new spatial optimization model for delineating CBSAs that can better account for polycentric urban structure by considering all intercounty commuting linkages. Our model seeks to find the boundaries that maximize the containment of entire webs of intra‐CBSA intercounty commuting. We apply the proposed method to delineate alternative sets of metropolitan and micropolitan statistical areas for comparison with the federal government's currently defined (2010 Office of Management and Budget [OMB] standards) official CBSAs. El método actual para delimitar áreas estadísticas micropolitanas y metropolitana en los Estados Unidos (denomidado genéricamente áreas estadísticos nucleares o core‐statistical areas‐ CBSA) consiste en agrupar condados en base a la fuerza de las interacciones de los desplazamientos residencia‐trabajo entre las zonas periféricas y los núcleos urbanos. La integración social y económica de un área estadística metro y micro politana se pone en operación usando un enfoque funcional que los autores consideran obsoleto que se basa en una estructura espacial urbana monocéntrica. Dicho enfoque se limita únicamente a los desplazamientos hacia el interior y de vuelta. Los autores proponen un nuevo modelo de optimización espacial para delinear CBSAs que respresenta mejor una estructura urbana policéntrica que considera todos los desplazamientos entre condados. El modelo busca encontrar los límites que maximizan la inclusión de redes enteras de desplazamientos intra CBSA y entre condados. El método propuesto es aplicado para delinear conjuntos alternativos de áreas estadísticas metropolitanas y micropolitanas y es comparado con los límites oficiales actuales de CBSA proporcionados por gobierno federal (de acuerdo a los lineamientos de la Oficina de Administración y Presupuesto u Office of Management Budget‐OMB) en 2010). 目前用于描述美国大都市区与居住区统计区域（共同地被称为基于核心的统计区域—CBSAs）的方法,主要是基于边远地区和城市核心区间的相互作用强度对郡县进行聚类。大都市区或居住区统计区域社会与经济一体化的运作方式,采用的是一种功能空间结构上落伍的单中心视角,即仅关注内部交换或反向交换。本文提出了一种新的空间优化模型用于描述CBSAs,将所有内部通勤联系考虑在内,能更好地解释多中心城市结构。该模型力图发现能最大化包含内部CBSA通勤联系网的边界。最后,将该模型应用于大都市区和居住区的统计区域描述,并与联邦政府目前定义的官方CBSAs（管理与预算 [OMB]标准2010办公室）进行了对比研究。     

Spatial Patterns of March and September Streamflow Trends in Pacific Northwest Streams, 1958–2008

Summer streamflow is a vital water resource for municipal and domestic water supplies, irrigation, salmonid habitat, recreation, and water‐related ecosystem services in the Pacific Northwest (PNW) in the United States. This study detects significant negative trends in September absolute streamflow in a majority of 68 stream‐gauging stations located on unregulated streams in the PNW from 1958 to 2008. The proportion of March streamflow to annual streamflow increases in most stations over 1,000 m elevation, with a baseflow index of less than 50, while absolute March streamflow does not increase in most stations. The declining trends of September absolute streamflow are strongly associated with seven‐day low flow, January–March maximum temperature trends, and the size of the basin (19–7,260 km²), while the increasing trends of the fraction of March streamflow are associated with elevation, April 1 snow water equivalent, March precipitation, center timing of streamflow, and October–December minimum temperature trends. Compared with ordinary least squares (OLS) estimated regression models, spatial error regression and geographically weighted regression (GWR) models effectively remove spatial autocorrelation in residuals. The GWR model results show spatial gradients of local R ² values with consistently higher local R ² values in the northern Cascades. This finding illustrates that different hydrologic landscape factors, such as geology and seasonal distribution of precipitation, also influence streamflow trends in the PNW. In addition, our spatial analysis model results show that considering various geographic factors help clarify the dynamics of streamflow trends over a large geographical area, supporting a spatial analysis approach over aspatial OLS‐estimated regression models for predicting streamflow trends. Results indicate that transitional rain–snow surface water‐dominated basins are likely to have reduced summer streamflow under warming scenarios. Consequently, a better understanding of the relationships among summer streamflow, precipitation, snowmelt, elevation, and geology can help water managers predict the response of regional summer streamflow to global warming. Patrones espaciales de las tendencias de los caudales de marzo y septiembre en el Pacífico Noroccidental. Los caudales (streamflows) de verano son recursos hídricos vitales para el abastecimiento de agua municipal y domestico así como para el riego agrícola, el hábitat de los salmónidos, la recreación, y para varios servicios de los ecosistemas en el Pacífico Noroccidental (Pacific Nortwest‐PNW) de los Estados Unidos. Este estudio identifica tendencias negativas considerables en los caudales absolutos de septiembre en la mayoría de las 68 estaciones de medición situadas en ríos y arroyos no regulares del PNW entre 1958 y 2008. La proporción del caudal de marzo con respecto al caudal anual aumenta en la mayoría de las estaciones situadas a más de 1000 metros de altitud, que tienen un índice de caudal base (base flow index‐BFI) de menos de 50, pero se mantiene estable en el resto (la mayoría) de las estaciones. Las tendencias decrecientes de los caudales absolutos de septiembre están fuertemente asociadas con el caudal mínimo para siete días (seven‐day low), con las tendencias de temperatura máxima entre enero y marzo, y con el tamaño de la cuenca (19‐7,260 km2). Las tendencias crecientes de la proporción del caudal total correspondiente a marzo están asociadas con la elevación, con un equivalente a agua de la nieve de abril (one April snow wáter equivalent ‐SWE), con la precipitación de marzo, el center timing (TC) de los caudales, y con las tendencias de temperatura mínima entre octubre y diciembre. En comparación con los estimados de modelos de regresión de tipo mínimos cuadrados ordinarios (ordinary least squares‐OLS), los modelos de regresión de error espacial (spatial error regression‐SER) y de regresión ponderada geográficamente (geographically weighted regression‐GWR) eliminan eficazmente la autocorrelación espacial en los residuos. Los resultados del modelo GWR producen mapas con gradientes espaciales donde los valores de los R² locales son consistentemente más altos en las cascadas del norte. Este resultado pone de manifiesto que diferentes factores hidrológicos del paisaje, tales como la geología y la distribución estacional de la precipitación, también influyen en las tendencias de los caudales en el PNW. Adicionalmente, los resultados del modelo de análisis espacial muestran que la inclusión de diversos factores geográficos ayuda a aclarar la dinámica de las tendencias de los caudales en un área geográfica grande, corroborando la mayor utilidad de modelos con enfoque de espacial sobre modelos. Los resultados indican también que en cuencas transicionales (entre lluvia y nieve) donde predominan las aguas superficiales existe una probabilidad más alta de reducccion de caudal de verano en el contexto de escenarios de calentamiento. En consecuencia, una mejor comprensión de las relaciones entre caudal de verano, la precipitación, el derretimiento de nieve, la elevación y la geología puede ayudar a los gestores del agua a predecir la respuesta de los caudales de verano en un escenario de calentamiento global. 夏季径流是美国太平洋西北部地区(PNW)市政与居民水供应、灌溉、鱼类栖息、娱乐及水相关生态系统服务的重要来源。本研究通过1958–2008年PNW地区68个地理位置上未调节的径流测量站分析揭示出九月绝对径流量与该地区水来源呈显著的负相关趋势。三月径流占年际径流量的比例在大多数海拔超过1000米的地区是增加的,然而大多数地区基本径流指数（BFI）少于50,表明三月绝对径流量在多数地区并没有增加。九月绝对径流量的下滑趋势与年均为期7天的低流量,1月–3月最大温度趋势及流域面积（19–7,260 km²）呈强相关性,而三月绝对径流量微量增加的趋势则与海拔高度,四月一日的雪水当量（SWE）,三月降水,径流中心时序（CT）和十月至十二月最小温度趋势相关。与OLS(普通最小二乘法估计)回归模型相比,空间滞后回归（SER）和地理加权回归模型（GWR）能有效剔除空间自相关的残差而更有效。GWR结果显示局部R²值在空间上渐变梯度,并且在北美洲喀斯喀特山脉（Cascade）：北部地区高于其他地区。该发现表明不同水文景观因子,如地质、降水的季节分布,也会对PNW的径流趋势产生影响。另外,我们的空间分析模型结果显示,考虑多种地理因素可解析大面积的地理区域中径流量空间分布趋势的动力机制, 为预测径流趋势提供一种优于OLS空间估计回归模型的空间分析方法。结果表明气候变暖背景下,季节性降水、地表水主导的流域极可能减少夏季径流量。因此,更好地理解夏季径流量与降水、融雪水、海拔和地质的关系,可帮助水资源管理者预测区域夏季径流对全球变暖的响应。     

A Dynamic Spatial Weight Matrix and Localized Space–Time Autoregressive Integrated Moving Average for Network Modeling

Various statistical model specifications for describing spatiotemporal processes have been proposed over the years, including the space–time autoregressive integrated moving average (STARIMA) and its various extensions. These model specifications assume that the correlation in data can be adequately described by parameters that are globally fixed spatially and/or temporally. They are inadequate for cases in which the correlations among data are dynamic and heterogeneous, such as network data. The aim of this article is to describe autocorrelation in network data with a dynamic spatial weight matrix and a localized STARIMA model that captures the autocorrelation locally (heterogeneity) and dynamically (nonstationarity). The specification is tested with traffic data collected for central London. The result shows that the performance of estimation and prediction is improved compared with standard STARIMA models that are widely used for space–time modeling. En los últimos años, se han propuesto diversas especificaciones de modelado estadístico para describir procesos espacio‐temporales. Esto incluye el modelo espacio‐temporal autorregresivo integrado de media móvil (STARIMA) y sus varios derivados. Estas especificaciones de modelo asumen que la correlación de los datos puede ser adecuadamente descrita por parámetros que se fijan a nivel global en el espacio y/o tiempo. Dichos parámetros son inadecuados para los casos en los que las correlaciones entre los datos son dinámicas y heterogéneas, como en el contexto de los datos de la red. El objetivo de este artículo es describir la autocorrelación en los datos de red con una matriz de ponderación espacial dinámica y un modelo STARIMA localizado (LSTARIMA) que captura la autocorrelación local (heterogeneidad) de forma dinámica (no estacionariedad). La especificación del modelo es evaluada con datos de tráfico recolectados en el centro de Londres. Los resultados demuestran que los rendimientos de estimación y predicción mejoran con el método propuesto en comparación con los modelos STARIMA estándar que son ampliamente utilizados para el modelado de espacio‐temporal. 通过设定多种统计模型来描述地理时空过程已提出多年,包括时空自回归移动平均(STARIMA)及其变形。此类模型通过假设数据相关性可由在时间域或者空间域上全局不变的参数加以充分描述。因此,上述模型不适用于具有动态或异质相关性的数据,如网络数据。本文试图采用一个动态空间权重矩阵与局部时空自回归移动平均(LSTARIMA)模型来描述数据的自相关程度,以此捕捉局域自相关(异质性)和动态自相关(非平稳性)。以伦敦市中心的交通数据作为模型的实证案例的测试结果显示,相对于广泛应用于时空过程分析的标准STARIMA模型,本文的模型在参数估计和预测性能上均有提升。     

Incorporating Spatial Dynamics and Temporal Dependency in Land Use Change Models

This article formulates an empirical discrete land use model within a spatially explicit economic structural framework for land use change decisions. The underlying framework goes beyond mechanistically fitting models for the spatial process of land use change to more closely link landowner decision behavior to land use patterns. At the same time, the article explicitly considers spatial spillover effects in the decisions of landowners of proximately located parcels. These spillover or peer influences may be due to strategic or collaborative partnerships between landowners, and can be associated with variables observable to the analyst (such as accessibility to city centers and market places) and variables unobservable to the analyst (such as perhaps soil quality and neighborhood attitudes/politics). In addition to spatial spillover effects, heterogeneity is also likely to exist in the decision‐making process of different landowners because of differential responsiveness to various signals relevant to decision making. This leads to correlation in land uses across time that is stationary for the same spatial unit. The article accommodates these technical considerations by formulating a random coefficients spatial lag discrete choice model using a fine resolution for the spatial unit of analysis. Time‐varying random effects are also considered to capture the effects of time‐varying unobserved factors (for instance, unobserved landowner attitudes regarding specific land uses may shift over time). The model is estimated using Bhat's maximum approximate composite marginal likelihood inference approach. The analysis is undertaken using the City of Austin parcel‐level land use database for multiple years (1995, 2000, 2003, and 2006). The estimation results indicate that proximity to highways and other roadways, distance from floodplains, parcel location in the context of existing development, and distance from schools are important determinants of land use. As importantly, the results provide very strong evidence of temporal dependency and spatial dynamics in land use decisions. There is also a suggestion that major highways may not only physically partition regions, but may also act as social barriers for didactic interactions among individuals. Este artículo presenta un modelo empírico discreto de uso de tierra dentro de un marco económico estructural espacialmente explícito para la toma de decisiones de cambio de uso de suelo. El marco utilizado por los autores va más allá del ajuste mecánico de modelos al proceso espacial de cambio de uso de suelo pues vincula más estrechamente el comportamiento y decisión del terrateniente a los patrones de uso observados. Al mismo tiempo, el estudio considera explícitamente los efectos espaciales de difusión en las decisiones de los propietarios de las parcelas cercanas. Esta difusión (spillover) o influencia de los pares puede deberse a alianzas estratégicas o de colaboración entre los terratenientes. También pueden estar asociadodos a variables observables (como la accesibilidad a los centros de las ciudades y plazas de mercado), así como a las variables no observables (por ejemplo, la calidad del suelo y las actitudes o tendencias políticas de los residentes). Además de los efectos de spillover, también es probable que exista heterogeneidad en el proceso de toma de decisiones de los diversos terratenientes, debido a su diferente capacidad de respuesta a las distintas señales que influencian la toma de decisiones. Esto conduce a que exista correlación en los usos del suelo a través del tiempo que es estacionaria para la misma unidad espacial. En el estudio todas las consideraciones técnicas mencionadas son tomadas en cuenta mediante la formulación de un modelo de elección discreta con rezago espacial con coeficientes aleatorios (random coefficients spatial lag discrete choice model) usando unidades espaciales de alta resolución. Para capturar los efectos de los factores no observados que varían temporalmente los autores utilizan efectos aleatorios (random effects) (por ejemplo, las actitudes de los terratenientes con respecto a usos específicos de tierras pueden cambiar con el tiempo). El modelo propuesto es estimado utilizando el enfoque inferencial de similitudes marginales de aproximaciones máximas compuestas de Bhat (2011) (maximum approximate composite marginal likelihood ‐MACML). El análisis se lleva a cabo usando una base de datos de a nivel de parcela base de uso de la tierra de la ciudad de Austin, Texas para varios años (1995, 2000, 2003 y 2006). Los resultados de la estimación indican que la proximidad a las autopistas y otras carreteras, la distancia de las llanuras de inundación, la ubicación de parcela en el contexto del desarrollo urbano existente, y la distancia a las escuelas, son factores importantes para el uso de suelo. Adicionalmente, los resultados proporcionan evidencia muy clara de la dependencia temporal y la dinámica espacial en las decisiones de uso de suelo. El estudio también sugiere que las carreteras principales dividen las regiones no sólo físicamente, sino que también pueden actuar como barreras sociales para las interacciones entre los individuos. 本文提出了土地利用变化决策下空间经济结构框架的一种经验离散土地利用模型。土地利用变化的空间过程与土地所有者决策行为更为紧密连接的土地利用模式,使得潜在的框架超出已有拟合模型。同时,本文明确认为毗邻地块区位在地主决策中具有空间溢出效应。这种溢出或对等效应可能源于与地主有战略上或合作伙伴的关系,并且这种效应分为分析中可测的变量关联（如城市中心和市场区的可达性）和不可测的变量关联（如可能为土地质量和邻居态度/政策）两种。除空间溢出效应,异质性也可能存在于不同地主决策制定的过程中,因为多种信号的不同响应与决策制定相关。这种异质性导致了土地利用在相同空间单元时间演化过程中的相关性是平稳的。本文在理解这些含义的基础上对此采纳的技术考虑是,通过利用随机高分辨率的空间滞后离散选择模型来分析空间单元。随时间变化的随机效应也被考虑用来捕获不可测的时变因子效应（如不可测的关于特殊土地利用的地主态度可能随时间改变）。本文的模型通过Bhat’s（2011）提出的近似最大联合边缘似然推断方法（MACML）进行估计。分析数据来源于奥斯汀市土地利用关于地块尺度多年的数据库（1995,2000,2003,2006）。估计结果表明,与高速公路和其它道路的近邻性,与洪泛平原的距离,现状开发条件下的地块区位,以及跟学校的距离对于土地利用模式的确定是非常重要的决定因素。同样重要的是,该结果为空间动力学与时间依赖性的土地利用决策提供了非常强的证据。本结果也显示主要高速公路不仅是区域分割的物质表征,而且在个体之间的交流中扮演着社会障碍的角色。     

Structural Determinism in the Interlocking World City Network

Taylor's interlocking specification of the world city network has offered researchers a theoretically informed way to measure the world city network using readily available firm location data. However, the number and size of firms that are viewed as linking cities to one another impose a form of structural determinism on the world city network. Specifically, when a relatively small number of firms is used to define a network among a relatively large set of cities, or when only larger firms are used, this approach is unable to reveal a wide range of structures features that may actually be present in the world city network it is intended to measure. Through a series of examples, I demonstrate how specific features of firm location data predetermine the structure, number and size of cliques, and density of world city networks derived using the interlocking approach. Concluding comments discuss the implications of this structural determinism by focusing on the case of the commonly used Globalization and World Cities data set 11, offer some suggestions concerning how the interlocking approach can be employed while avoiding structural determinism, and identify some alternative approaches to mapping the world city network. El modelo interconexión cerrada (interlocking) de la red mundial de ciudades (world city network) desarrollado por Taylor (2001) ha proporcionado a los investigadores un marco teórico para medir dicha red con datos fácilmente disponibles acerca de la localización de empresas. Sin embargo, el número y tamaño de las empresas que son consideradas como ‘vinculadoras’ (de ciudades) imponen una forma de determinismo estructural sobre la red. En concreto, cuando un número relativamente pequeño de firmas es utilizado para definir la red entre un conjunto relativamente grande de ciudades, o, cuando sólo se utilizan las empresas más grandes, este enfoque (interlocking) no es capaz de identificar o medir una amplia gama de características estructurales presentes potencialmente en la red de ciudades mundiales. A través de una serie de ejemplos, el autor demuestra cómo varias características específicas de los datos de localización de las empresas predeterminan la estructura, número y tamaño de los cliques y la densidad de las redes de ciudades mundiales resultantes de la aplicación de modelos de interlocking. Luego el autor presenta las implicancias del determinismo estructural mencionado, centrándose en el caso la base de datos de uso común GaWC # 11. Finalmente el articulo ofrece algunas sugerencias sobre cómo el enfoque de interlocking puede ser empleado sin caer en un determinismo estructural, y se identifican algunos enfoques alternativos para el mapeo de la red mundial de ciudades. 泰勒（2001）采用相对容易获得的公司地理位置数据提出的世界城市网络连锁度量方法,为研究者提供了一种理论途径。然而,将公司的数量和规模视为城市之间的联接是一种强加的世界城市网络结构决定论形式。特别是当用数量相对较少的公司或仅用大公司的数量来定义相对较大的城市网络集时,该方法不能揭示目前实际情况中大范围世界城市网络需要度量的结构特征。通过一系列的案例分析,本文显示了公司区位数据的特有特征预先决定了集团的结构、数量和规模,以及由泰勒的连锁方法推导出世界城市网络的密度。最后,通过曾被广泛采用的GaWC（全球化和世界城市）数据集的11个城市公司数据讨论了结构决定方法的应用,并就应用该方法如何避免结构决定论提出一些建议,进而识别一些可选用的方法来绘制世界城市网络图。     

A Hierarchical Flow Capturing Location Problem with Demand Attraction Based on Facility Size,and Its Lagrangian Relaxation Solution Method

This article presents a hierarchical flow capturing location problem (HFCLP) and proposes an effective Lagrangian heuristic solution method. The original flow capturing location problem (FCLP) aims to locate a given number of facilities on a network to maximize the total flow that can be serviced at facilities along their preplanned routes, such as daily commute to work. We extend the original model to allow a decision maker to select the size of facilities among m different size alternatives. Larger facilities are assumed to be more attractive and, therefore, can attract more customers, but they cost more to construct than smaller ones. Customers deviate from their preplanned routes to access a facility's service when the size of the facility is sufficiently large. The degree of deviation from the original path is measured by the additional distance customers have to go to access facilities, and the acceptable deviation distance becomes larger as the size of a facility increases. This article presents a new problem in which the number of facilities of each size and their locations are simultaneously determined so as to capture as much flow as possible within the total budget available for locating all facilities. We present an integer programming formulation of the problem and devise a Lagrangian relaxation solution method. The proposed algorithm is tested using road networks with 300 and 500 nodes. The results show that the method produces high‐quality solutions in a fairly short time. Este artículo presenta un problema de localización de captura de flujo jerárquico (hierarchical flow capturing location problem‐HFCLP) y propone un método heurístico eficiente de tipo Lagrange (lagrangian). En su formulación original el HFLCP tiene como objetivo localizar un número determinado de instalaciones en una red con el fin de maximizar el flujo total que puede ser atendido por las instalaciones existentes a lo largo de rutas preestablecidas, como en el caso por ejemplo, de los desplazamientos diarios del lugar de residencia al de trabajo. Los autores amplían el modelo original para permitir que el tomador de decisiones seleccione el tamaño de las instalaciones entre “m” alternativas. Se asume que las instalaciones más grandes son más atractivas que las más pequeñas y, por lo tanto, pueden atraer a más clientes, pero a la vez, son también más costosas de construir. Los clientes se desvían de su ruta preestablecida para acceder al servicio de una instalación cuando el tamaño de la instalación es lo suficientemente grande. El grado de desviación de las rutas se mide por la distancia adicional que los clientes viajan para acceder a las instalaciones. La distancia de desviación aceptable se hace más grande en relación al tamaño de la instalación. En este artículo se presenta un nuevo modelo para el HFLCP en el que el número de las instalaciones de cada tamaño y su ubicación son determinadas simultáneamente con el fin de capturar la mayor cantidad de flujo dentro del presupuesto total disponible para la localización de todas las instalaciones. Los autores presentan una formulación de programación entera (integer programming) del HFCLP e implementan un método que relaja la solución lagrangiana. El algoritmo propuesto es evaluado utilizando redes viales con 300 y 500 nodos. Los resultados muestran que el nuevo método produce soluciones de alta calidad y en tiempos de computación relativamente cortos. 本文介绍了一种分层的截流选址问题 (HFCLP),提出了一个有效的拉格朗日启发式解决方法。最初的截流选址问题(FCLP)目标是在网络上布局给定数量的设施使总流量最大,使按预定路线的行进流可以获得最大的服务,如每日的工作通勤。本文对原始模型进行扩展,让决策者可在不同的设施规模选择方案中进行规模选择。假设更大规模设施具有更大的吸引力,因此也能够吸引更多的客户,但同时也需要更多的建造成本。当设施规模足够大时,消费者会选择偏离预定路径而进入该设施的服务范围。对原始路径的偏离程度可通过用户进入该设施所增加的额外距离度量。可接受的偏差距离随着设施规模增大而增大。本文提出了在总预算确定条件下,同步确定不同规模设施数量及其位置以实现截取最大流量的新问题,并给出了该问题的整数规划方法,设计了拉格朗日松弛解法。通过300和500个节点的网络测试,结果显示该算法可在相当短时间内获得高质量的解决方案。     

Computing the Jacobian in Gaussian Spatial Autoregressive Models: An Illustrated Comparison of Available Methods

When estimating spatial regression models by maximum likelihood using spatial weights matrices to represent spatial processes, computing the Jacobian, ln(|I ? λW|), remains a central problem. In principle, and for smaller data sets, the use of the eigenvalues of the spatial weights matrix provides a very rapid resolution. Analytical eigenvalues are available for large regular grids. For larger problems not on regular grids, including those induced in spatial panel and dyadic (network) problems, solving the eigenproblem may not be feasible, and a number of alternatives have been proposed. This article surveys selected alternatives, and comments on their relative usefulness, covering sparse Cholesky and sparse LU factorizations, and approximations such as Monte Carlo, Chebyshev, and using lower‐order moments with interpolation. The results are presented in terms of component‐wise differences between sets of Jacobians for selected data sets. In conclusion, recommendations are made for a number of analytical settings. Al estimar modelos de regresión espacial con el método del máxima verosimilitud (máximum likelihood) y usando matrices de pesos espaciales para representar procesos espaciales, cálculo del término jacobiano (jabobian)—ln(| I ?λ W |)‐ sigue siendo un problema central. En principio, y para bases de datos más pequeñas, el uso de los valores propios (eigenvalues) de la matriz de pesos espaciales proporciona una solución muy rápida. Los eigenvalues analíticos para retículas o grillas grandes y regulares son ya conocidos. Para problemas más grandes, que no se presentan en mallas regulares ‐incluyendo aquellos que se inducen en problemas de paneles espaciales y en problemas de (redes) diádicas‐, es posible que resolver el eigenproblem no sea posible. Este artículo estudia una selección de alternativas y comenta acerca de su relativa utilidad. Se cubren las facorizaciones de tipo Cholesky disperso (sparse Cholesky) y de tipo LU dispersas (sparse LU), las aproximaciones Monte Carlo, y Chebyshev, así mismo se utiliza momentos de bajo‐orden (lower‐order) con interpolación. Los resultados se presentan en términos de diferencias de componentes entre sets de términos jacobianos para bases de datos seleccionadas. En conclusión, se hacen recomendaciones para una serie de contextos analíticos. 当采用表征空间过程的空间权重矩阵对空间回归模型进行最大似然估计时,雅可比矩阵ln(|I?λW|)的计算仍是核心问题。对于小数据集,原则上可利用空间权重矩阵的特征值提供一种快速的解决方案,对于大型规则格网数据特征值分析同样有效。但对于不规则格网大型问题,包括从空间面板和二元（网络）问题中引伸的问题,利用特征值的解决方案可能不适用,对此学术界提出了多种可选替代方案。本文选取已有的几种替代方案并评论各自的相对有效性,其中包括稀疏Cholesky分解和稀疏LU分解法,Monte Carlo和 Chebyshev近似模拟法以及低阶矩插值法。结果以所选数据集雅可比矩阵间特定组份的差异方式显示。最后,推荐了一些分析设定。