Review on Spatial Data
|
|
- Raymond Heath
- 6 years ago
- Views:
Transcription
1 Week 12 Lecture: Spatial Autocorrelation and Spatial Regression Introduction to Programming and Geoprocessing Using R GEO GEO Point data Review on Spatial Data Area/lattice data May be regular or irregular Recall the key components of spatial data: Spatial information (coordinates, CRS, etc.) Attributes Some of this lecture based on a workshop presented by Elisabeth Root, Department of Geography, University of Colorado, Boulder 1
2 Tobler s First Law of Geography "Everything is related to everything else, but near things are more related than distant things." Spatial autocorrelation is the formal property that measures the degree to which near and distant things are related The key: A statistical ttiti ltest tof match thbt between locational similarity and attribute similarity Positive, negative or zero relationship Spatial Regression Steps in determining the extent of spatial autocorrelation in your data and running a spatial regression: 1. Choose a neighborhood criterion Which areas are linked? 2. Assign weights to the areas that are linked Create a spatial weights matrix 3. Run statistical test to examine spatial autocorrelation 4. Run an OLS regression 5. Run a spatial regression(s) By applying weights matrices 2
3 Assessing Spatial Autocorrelation Is a spatial model worth the trouble? Spatial Weights Matrices Neighborhoods can be defined in multiple ways Contiguity it (common boundary) But what is a shared boundary? Distance (distance band, K nearest neighbors) How many neighbors to include, what distance do we use? General weights (social distance, exponential decay) 3
4 Importing shapefiles into R and constructing neighborhood sets STEP 1: CHOOSE A NEIGHBORHOOD CRITERION R Libraries for Spatial Analyses Can install our libraries individually using install.packages Or, we can install the entire Spatial set of libraries, called a view. > install.packages( ctv ) lib ( t ) > library( ctv ) > install.views( Spatial ) 4
5 Spatial View Maintained by Bivand, and described here: project.org/web/views/spatial.html Other views: project.org/web/views/ 5
6 R Libraries for Spatial Analyses If we install the entire Spatial view, then all our various spatial libraries i will then be available: > install.packages( ctv ) > library( ctv ) > install.views( Spatial ) > library(maptools) > library(rgdal) > library(spdep) 6
7 Loading an ESRI Shapefile > library(maptools) > getinfo.shape("c:/tmp/data/sids.shp") Shapefile type: Polygon, (5), # of Shapes: 100 > sids <- readshapepoly("c:/tmp/data/sids.shp") > class(sids) [1] "SpatialPolygonsDataFrame" attr(,"package") Loading a Shapefile With Projection Information > library(rgdal) > sids <- readogr( dsn="c:/tmp/data/", layer="sids ) OGR data source with driver: ESRI Shapefile Source: "C:/tmp/data/", layer: "sids" with 100 features and 18 fields Feature type: wkbpolygon with 2 dimensions > class(sids) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" 7
8 Projection of the Shapefile If the shapefile has no.prj file associated with it, you need to assign a coordinate system: > proj4string(sids)<-crs("+proj=longlat ellps=wgs84") And then we can transform the map into any projection: > sids_nad <- sptransform(sids, CRS("+init=epsg:3358")) > sids_sp <- sptransform(sids, CRS("+init=ESRI:102719")) For a list of applicable CRS codes: Easiest with the EPSG and ESRI pre defined codes 8
9 Contiguity Based Neighbors Areas sharing any boundary point ( queen s ) are taken as neighbors, using the poly2nb function, which acceptsa a SpatialPolygonsDataFrame > library(spdep) > sids_nbq<-poly2nb(sids) If contiguity is defined as areas sharing more than one boundary point ( rook s ), the queen= argument is set to FALSE > sids_nbr<-poly2nb(sids, queen=false) > coords<-coordinates(sids) > plot(sids) > plot(sids_nbq, coords, add=t) Queen s Rook s 9
10 Distance Based Neighbors Using K Nearest Neighbors (KNN) > coords <- coordinates(sids_sp) > IDs <- row.names(as(sids_sp, "data.frame")) > sids_kn1<-knn2nb(knearneigh(coords, k=1), row.names=ids) > sids_kn2<-knn2nb(knearneigh(coords, k=2), row.names=ids) > sids_kn4<-knn2nb(knearneigh(coords, k=4), row.names=ids) > plot(sids_sp) > plot(sids_kn2, coords, add=t) k=2 k=3 k=1 k=1 k=2 k=4 10
11 Distance Based Neighbors We might also assign neighbors based on a specified distance: > dist <- unlist(nbdists(sids_kn1, coords)) > summary(dist) Min. 1st Qu. Median Mean 3rd Qu. Max > max_k1 <- max(dist) > sids_kd1<-dnearneigh(coords, d1=0, d2 = 0.75*max_k1, row.names=ids) > sids_kd2<-dnearneigh(coords, d1=0, d2 = 1*max_k1, row.names=ids) > sids_kd3<-dnearneigh(coords, dnearneigh(coords d1=0, d2 = 1.5*max_k1, row.names=ids) dist = 0.75*max_k1 dist = 1*max_k1 = dist = 1.5 * max_k1 11
12 Creating spatial weights matrices using neighborhood lists STEP 2: ASSIGN WEIGHTS TO THE AREAS THAT ARE LINKED Spatial weights matrices Once we define which observations are neighbors, we assign weights iht to each link Can be binary or variable Even when the values are binary (0/1), the issue ofwhat to do with no neighborneighbor observations arises 12
13 Row standardized Weights Matrix > sids_nbq_w<- nb2listw(sids_nbq, style="w") > sids_nbq_w Characteristics of weights list: Neighbour list object: Number of regions: 100 Number of nonzero links: 490 Percentage nonzero weights: 4.9 Average number of links: 4.9 Weights style: W Weights constants summary: n nn S0 S1 S2 W Row standardization is used to create proportional weights when features have unequal numbers of neighbors Divide each neighbor weight for a feature by the sum of all neighbor weights Obs i has 3 neighbors, each has a weight of 1/3 Obs j has 2 neighbors, each has a weight of 1/2 Must be used if you want comparable spatial parameters across different data sets with different connectivity structures > sids_nbq_wb <- nb2listw(sids_nbq, style="b") > sids_nbq_wb Characteristics ti of weights list: Neighbour list object: Number of regions: 100 Number of nonzero links: 490 Percentage nonzero weights: 4.9 Average number of links: 4.9 Binary weights Row standardised weights increase the influence of links from observations with few neighbours Weights style: B Weights constants summary: n nn S0 S1 S2 B Binary weights vary the influence ceof observations o s Those with many neighbours are upweighted compared to those with few 13
14 Binary vs. Row Standardized A binary weights matrix looks like: A1 A2 A3 A4 A A A A A row standardized matrix it looks like: A1 A2 A3 A4 A A A A Regions With No Neighbors Error in nb2listw(filename): Empty neighbor sets found You have some regions that have no neighbors Could be: A problem in the GIS topology (digitizing errors, slivers, etc) Induced by hard coded distance thresholds, and are true islands (e.g., Hawaii) May want to use k nearest neighbors Or add zero.policy=t to the nb2listw call > sids_nbq_w<-nb2listw(sids_nbq, zero.policy=t) 14
15 How to use the spatial weights, assess degree of spatial autocorrelation STEP 3: EXAMINE SPATIAL AUTOCORRELATION Spatial Autocorrelation Test for the presence of spatial autocorrelation ti Global Moran s I Geary s C Local (LISA Local Indicators of Spatial Autocorrelation) Local Moran s I and Getis G i * 15
16 Moran s I in R > moran.test(sids_nad$sidr79, listw=sids_nbq_w, alternative= two.sided ) two.sided H A : I I 0 greater H A : I > I 0 Moran's I test under randomisation data: sids_nad$sidr79 weights: sids_nbq_w Moran I statistic standard deviate = , p-value = alternative hypothesis: greater sample estimates: Moran I statistic Expectation Variance Moran s I in R > moran.test(sids_nad$sidr79, listw=sids_nbq_wb) Moran's I test under randomisation data: sids_nad$sidr79 weights: sids_nbq_wb Moran I statistic standard deviate = , p-value = alternative hypothesis: greater sample estimates: Moran I statistic Expectation Variance
17 Spatial Regression And how do we do it in R? Spatial Autocorrelation in Residuals: Spatial Error Model Incorporates spatial effects through error term Where: ε is λ is the spatial error coefficient y = xβ + ε ε = λwε + ξ the vector of error terms, spatially weighted using the weights matrix (W) ξ is a vector of uncorrelated error terms If there is no spatial correlation between the errors, then λ = 0 17
18 Spatial Autocorrelation in Obs. Values: Spatial Lag Model Incorporates spatial effects by including a spatially lagged dependent variable as an additional predictor Where: Wy is ε is a vector of error terms y = ρ Wy + xβ + ε the spatially lagged DVs for weights matrix W x is a matrix of observations on the explanatory variables ρ is the spatial coefficient If there is no spatial dependence, and y does no depend on neighboring y values, ρ = 0 Good books Bivand, R., et al. (2007) Applied Spatial Data Analysis with R. New York: Springer. Ward, M.D. and K.S. Gleditsch (2008) Spatial Regression Models. Thousand Oaks, CA: Sage. 18
19 CRIM Spatial Regression in R Example: Housing Prices in Boston per capita crime rate by town ZN proportion of residential land zoned for lots over 25,000 ft 2 INDUS proportion of non retail business acres per town CHAS Charles River dummy variable (=1 if tract bounds river; 0 otherwise) NOX Nitrogen oxide concentration (parts per 10 million) RM average number of rooms per dwelling AGE proportion of owner occupied units built prior to 1940 DIS weighted distances to five Boston employment centres RAD index of accessibility to radial highways TAX full value property tax rate per $10,000 PTRATIO pupil teacher ratio by town B 1000(Bk 0.63) 2 where Bk is the proportion of blacks by town LSTAT % lower status of the population MEDV Median value of owner occupied homes in $1000's Spatial Regression in R 1. Read in a shapefile (boston.shp) 2. Define neighbors (k nearest w/point data) 3. Create weights matrix 4. Moran s test of observed, Moran scatterplot 5. Run OLS regression 6. Check residuals for spatial dependence 7. Determine which SR model to use w/lm tests 8. Run spatial regression model 19
20 Define Neighbors and Create Weights Matrix > boston<-readogr(dsn= C:/tmp/data/",layer="boston") > class(boston) > boston$logmedv<-log(boston$cmedv) > coords<-coordinates(boston) > IDs<-row.names(as(boston, "data.frame")) > bost_ kd1<-dnearneigh(coords, d e ds, d1=0, d2=3.973, row.names=ids) s) > plot(boston) > plot(bost_kd1, coords, add=t) > bost_kd1_w<- nb2listw(bost_kd1) Moran s I on the Observed Median House Values > moran.test(boston$logmedv, listw=bost_kd1_w) Moran's I test under randomisation data: boston$logmedv weights: bost_kd1_w Moran I statistic standard deviate = , p-value < 2.2e-16 alternative hypothesis: greater sample estimates: Moran I statistic Expectation Variance
21 Moran Plot for the DV > moran.plot(boston$logmedv, bost_kd1_w, labels=as.character(boston$id)) OLS Regression bostlm<-lm(logmedv~rm + LSTAT + CRIM + ZN + CHAS + DIS, data=boston) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr(> t ) (Intercept) < 2e-16 *** RM e-11 *** LSTAT < 2e-16 *** CRIM < 2e-16 *** ZN *** CHAS *** DIS e-06 *** --- Residual standard error: on 499 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 6 and 499 DF, p-value: < 2.2e-16 21
22 Checking residuals for spatial autocorrelation > boston$lmresid<-residuals(bostlm) > lm.morantest(bostlm, bost_kd1_w) Global Moran's I for regression residuals Moran I statistic standard deviate = , p-value = 2.396e-09 alternative hypothesis: greater sample estimates: Observed Moran's I Expectation Variance Determining the type of dependence > lm.lmtests(bostlm, bost_kd1_w, test="all") Lagrange multiplier diagnostics for spatial dependence LMerr = , df = 1, p-value = 3.201e-07 LMlag = , df = 1, p-value = 8.175e-12 RLMerr = , df = 1, p-value = RLMlag = , df = 1, p-value = 4.096e-07 SARMA = , df = 2, p-value = 5.723e-12 Robust tests used to find a proper alternative Only use robust forms when BOTH LMErr and LMLag are significant 22
23 One more diagnostic > install.packages( lmtest ) > library(lmtest) > bptest(bostlm) studentized Breusch-Pagan test data: bostlm BP = , df = 6, p-value = 2.651e-13 Indicates errors are heteroskedastic Not surprising since we have spatial dependence Running a spatial lag model > bostlag<-lagsarlm(logmedv~rm + LSTAT + CRIM + ZN + CHAS + DIS, data=boston, bost_kd1_w) Type: lag Coefficients: (asymptotic standard errors) Estimate Std. Error z value Pr(> z ) ) (Intercept) < 2.2e-16 RM e-10 LSTAT < 2.2e-16 CRIM < 2.2e-16 ZN CHAS DIS e-11 Rho: , LR test value:37.426, p-value:9.4936e-10 Asymptotic standard error: z-value: , p-value: e-11 Wald statistic: , p-value: e-11 Log likelihood: for lag model ML residual variance (sigma squared): , (sigma: ) AIC: , (AIC for lm: ) 23
24 A few more diagnostics LM test for residual autocorrelation test value: , p-value: > bptest.sarlm(bostlag) studentized Breusch-Pagan test data: BP = , df = 6, p-value = 4.451e-11 LM test suggests there is no more spatial autocorrelation in the data BP test indicates remaining heteroskedasticity in the residuals Most likely due to misspecification Running a spatial error model > bosterr<-errorsarlm(logmedv~rm + LSTAT + CRIM + ZN + CHAS + DIS, data=boston, listw=bost_kd1_w) Type: error Coefficients: (asymptotic standard errors) Estimate Std. Error z value Pr(> z ) ) (Intercept) < 2.2e-16 RM e-09 LSTAT < 2.2e-16 CRIM < 2.2e-16 ZN CHAS DIS Lambda: , LR test value: , p-value: e-07 Asymptotic standard error: z-value: , p-value: e-12 Wald statistic: 46.35, p-value: e-12 Log likelihood: for error model ML residual variance (sigma squared): , (sigma: ) AIC: , (AIC for lm: ) 24
GRAD6/8104; INES 8090 Spatial Statistic Spring 2017
Lab #5 Spatial Regression (Due Date: 04/29/2017) PURPOSES 1. Learn to conduct alternative linear regression modeling on spatial data 2. Learn to diagnose and take into account spatial autocorrelation in
More informationMultiple Regression Part I STAT315, 19-20/3/2014
Multiple Regression Part I STAT315, 19-20/3/2014 Regression problem Predictors/independent variables/features Or: Error which can never be eliminated. Our task is to estimate the regression function f.
More informationData Mining Techniques. Lecture 2: Regression
Data Mining Techniques CS 6220 - Section 3 - Fall 2016 Lecture 2: Regression Jan-Willem van de Meent (credit: Yijun Zhao, Marc Toussaint, Bishop) Administrativa Instructor Jan-Willem van de Meent Email:
More informationLecture 6: Hypothesis Testing
Lecture 6: Hypothesis Testing Mauricio Sarrias Universidad Católica del Norte November 6, 2017 1 Moran s I Statistic Mandatory Reading Moran s I based on Cliff and Ord (1972) Kelijan and Prucha (2001)
More information<br /> D. Thiebaut <br />August """Example of DNNRegressor for Housing dataset.""" In [94]:
sklearn Tutorial: Linear Regression on Boston Data This is following the https://github.com/tensorflow/tensorflow/blob/maste
More informationHW1 Roshena MacPherson Feb 1, 2017
HW1 Roshena MacPherson Feb 1, 2017 This is an R Markdown Notebook. When you execute code within the notebook, the results appear beneath the code. Question 1: In this question we will consider some real
More informationIntroduction to PyTorch
Introduction to PyTorch Benjamin Roth Centrum für Informations- und Sprachverarbeitung Ludwig-Maximilian-Universität München beroth@cis.uni-muenchen.de Benjamin Roth (CIS) Introduction to PyTorch 1 / 16
More informationStat588 Homework 1 (Due in class on Oct 04) Fall 2011
Stat588 Homework 1 (Due in class on Oct 04) Fall 2011 Notes. There are three sections of the homework. Section 1 and Section 2 are required for all students. While Section 3 is only required for Ph.D.
More informationBayesian Model Averaging (BMA) with uncertain Spatial Effects A Tutorial. Martin Feldkircher
Bayesian Model Averaging (BMA) with uncertain Spatial Effects A Tutorial Martin Feldkircher This version: October 2010 This file illustrates the computer code to use spatial filtering in the context of
More informationSupervised Learning. Regression Example: Boston Housing. Regression Example: Boston Housing
Supervised Learning Unsupervised learning: To extract structure and postulate hypotheses about data generating process from observations x 1,...,x n. Visualize, summarize and compress data. We have seen
More informationOn the Harrison and Rubinfeld Data
On the Harrison and Rubinfeld Data By Otis W. Gilley Department of Economics and Finance College of Administration and Business Louisiana Tech University Ruston, Louisiana 71272 (318)-257-3468 and R. Kelley
More informationSpatial Regression. 10. Specification Tests (2) Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 10. Specification Tests (2) Luc Anselin http://spatial.uchicago.edu 1 robust LM tests higher order tests 2SLS residuals specification search 2 Robust LM Tests 3 Recap and Notation LM-Error
More informationIntroduction to Spatial Statistics and Modeling for Regional Analysis
Introduction to Spatial Statistics and Modeling for Regional Analysis Dr. Xinyue Ye, Assistant Professor Center for Regional Development (Department of Commerce EDA University Center) & School of Earth,
More informationLecture 1: Introduction to Spatial Econometric
Lecture 1: Introduction to Spatial Econometric Professor: Mauricio Sarrias Universidad Católica del Norte September 7, 2017 1 Introduction to Spatial Econometric Mandatory Reading Why do We Need Spatial
More informationCOLUMN. Spatial Analysis in R: Part 2 Performing spatial regression modeling in R with ACS data
Spatial Demography 2013 1(2): 219-226 http://spatialdemography.org OPEN ACCESS via Creative Commons 3.0 ISSN 2164-7070 (online) COLUMN Spatial Analysis in R: Part 2 Performing spatial regression modeling
More informationCreating and Managing a W Matrix
Creating and Managing a W Matrix Carlos Hurtado Department of Economics University of Illinois at Urbana-Champaign hrtdmrt2@illinois.edu Junel 22th, 2016 C. Hurtado (UIUC - Economics) Spatial Econometrics
More informationON CONCURVITY IN NONLINEAR AND NONPARAMETRIC REGRESSION MODELS
STATISTICA, anno LXXIV, n. 1, 2014 ON CONCURVITY IN NONLINEAR AND NONPARAMETRIC REGRESSION MODELS Sonia Amodio Department of Economics and Statistics, University of Naples Federico II, Via Cinthia 21,
More informationSPACE Workshop NSF NCGIA CSISS UCGIS SDSU. Aldstadt, Getis, Jankowski, Rey, Weeks SDSU F. Goodchild, M. Goodchild, Janelle, Rebich UCSB
SPACE Workshop NSF NCGIA CSISS UCGIS SDSU Aldstadt, Getis, Jankowski, Rey, Weeks SDSU F. Goodchild, M. Goodchild, Janelle, Rebich UCSB August 2-8, 2004 San Diego State University Some Examples of Spatial
More informationStat 401B Final Exam Fall 2016
Stat 40B Final Exam Fall 0 I have neither given nor received unauthorized assistance on this exam. Name Signed Date Name Printed ATTENTION! Incorrect numerical answers unaccompanied by supporting reasoning
More informationMeasuring The Benefits of Air Quality Improvement: A Spatial Hedonic Approach. Chong Won Kim, Tim Phipps, and Luc Anselin
Measuring The Benefits of Air Quality Improvement: A Spatial Hedonic Approach Chong Won Kim, Tim Phipps, and Luc Anselin Paper prepared for presentation at the AAEA annual meetings, Salt Lake City, August,
More informationAnalyzing spatial autoregressive models using Stata
Analyzing spatial autoregressive models using Stata David M. Drukker StataCorp Summer North American Stata Users Group meeting July 24-25, 2008 Part of joint work with Ingmar Prucha and Harry Kelejian
More informationSTK 2100 Oblig 1. Zhou Siyu. February 15, 2017
STK 200 Oblig Zhou Siyu February 5, 207 Question a) Make a scatter box plot for the data set. Answer:Here is the code I used to plot the scatter box in R. library ( MASS ) 2 pairs ( Boston ) Figure : Scatter
More informationDISCRIMINANT ANALYSIS: LDA AND QDA
Stat 427/627 Statistical Machine Learning (Baron) HOMEWORK 6, Solutions DISCRIMINANT ANALYSIS: LDA AND QDA. Chap 4, exercise 5. (a) On a training set, LDA and QDA are both expected to perform well. LDA
More informationSpatial Investigation of Mineral Transportation Characteristics in the State of Washington
Khachatryan, Jessup 1 Spatial Investigation of Mineral Transportation Characteristics in the State of Washington Hayk Khachatryan Graduate Student Email: hkhachatryan@wsu.edu Eric L. Jessup Assistant Professor
More informationSKLearn Tutorial: DNN on Boston Data
SKLearn Tutorial: DNN on Boston Data This tutorial follows very closely two other good tutorials and merges elements from both: https://github.com/tensorflow/tensorflow/blob/master/tensorflow/examples/skflow/boston.py
More informationSpatial Econometrics with R - Spatial Data Analysis of the 5-Region Script Example
1 Prof. Dr. Reinhold Kosfeld October 2017 Table of contents 1. Introduction Spatial Econometrics with R - Spatial Data Analysis of the 5-Region Script Example 2. Creating spatial weights matrices, variable
More informationSpatial Regression. 9. Specification Tests (1) Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 9. Specification Tests (1) Luc Anselin http://spatial.uchicago.edu 1 basic concepts types of tests Moran s I classic ML-based tests LM tests 2 Basic Concepts 3 The Logic of Specification
More informationLecture 3: Exploratory Spatial Data Analysis (ESDA) Prof. Eduardo A. Haddad
Lecture 3: Exploratory Spatial Data Analysis (ESDA) Prof. Eduardo A. Haddad Key message Spatial dependence First Law of Geography (Waldo Tobler): Everything is related to everything else, but near things
More informationST430 Exam 2 Solutions
ST430 Exam 2 Solutions Date: November 9, 2015 Name: Guideline: You may use one-page (front and back of a standard A4 paper) of notes. No laptop or textbook are permitted but you may use a calculator. Giving
More informationA Robust Approach of Regression-Based Statistical Matching for Continuous Data
The Korean Journal of Applied Statistics (2012) 25(2), 331 339 DOI: http://dx.doi.org/10.5351/kjas.2012.25.2.331 A Robust Approach of Regression-Based Statistical Matching for Continuous Data Sooncheol
More informationSpatial Regression. 1. Introduction and Review. Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 1. Introduction and Review Luc Anselin http://spatial.uchicago.edu matrix algebra basics spatial econometrics - definitions pitfalls of spatial analysis spatial autocorrelation spatial
More information22s:152 Applied Linear Regression. In matrix notation, we can write this model: Generalized Least Squares. Y = Xβ + ɛ with ɛ N n (0, Σ)
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y Ordinary Least Squares Estimation The classical models we have fit so far with a continuous response
More informationSpatial Stochastic frontier models: Instructions for use
Spatial Stochastic frontier models: Instructions for use Elisa Fusco & Francesco Vidoli June 9, 2015 In the last decade stochastic frontiers traditional models (see Kumbhakar and Lovell, 2000 for a detailed
More informationSpatial Autocorrelation
Spatial Autocorrelation Luc Anselin http://spatial.uchicago.edu spatial randomness positive and negative spatial autocorrelation spatial autocorrelation statistics spatial weights Spatial Randomness The
More informationSpatial Autocorrelation (2) Spatial Weights
Spatial Autocorrelation (2) Spatial Weights Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline
More informationLecture 3: Exploratory Spatial Data Analysis (ESDA) Prof. Eduardo A. Haddad
Lecture 3: Exploratory Spatial Data Analysis (ESDA) Prof. Eduardo A. Haddad Key message Spatial dependence First Law of Geography (Waldo Tobler): Everything is related to everything else, but near things
More informationHeteroskedasticity. Part VII. Heteroskedasticity
Part VII Heteroskedasticity As of Oct 15, 2015 1 Heteroskedasticity Consequences Heteroskedasticity-robust inference Testing for Heteroskedasticity Weighted Least Squares (WLS) Feasible generalized Least
More informationSpatial Analysis 2. Spatial Autocorrelation
Spatial Analysis 2 Spatial Autocorrelation Spatial Autocorrelation a relationship between nearby spatial units of the same variable If, for every pair of subareas i and j in the study region, the drawings
More information22s:152 Applied Linear Regression. Returning to a continuous response variable Y...
22s:152 Applied Linear Regression Generalized Least Squares Returning to a continuous response variable Y... Ordinary Least Squares Estimation The classical models we have fit so far with a continuous
More informationExploratory Spatial Data Analysis (ESDA)
Exploratory Spatial Data Analysis (ESDA) VANGHR s method of ESDA follows a typical geospatial framework of selecting variables, exploring spatial patterns, and regression analysis. The primary software
More informationNon-independence due to Time Correlation (Chapter 14)
Non-independence due to Time Correlation (Chapter 14) When we model the mean structure with ordinary least squares, the mean structure explains the general trends in the data with respect to our dependent
More informationSparse polynomial chaos expansions as a machine learning regression technique
Research Collection Other Conference Item Sparse polynomial chaos expansions as a machine learning regression technique Author(s): Sudret, Bruno; Marelli, Stefano; Lataniotis, Christos Publication Date:
More informationModeling the Ecology of Urban Inequality in Space and Time
Modeling the Ecology of Urban Inequality in Space and Time Corina Graif PhD Candidate, Department Of Sociology Harvard University Presentation for the Workshop on Spatial and Temporal Modeling, Center
More informationThe Use of Spatial Weights Matrices and the Effect of Geometry and Geographical Scale
The Use of Spatial Weights Matrices and the Effect of Geometry and Geographical Scale António Manuel RODRIGUES 1, José António TENEDÓRIO 2 1 Research fellow, e-geo Centre for Geography and Regional Planning,
More informationSpatial Econometrics
Spatial Econometrics Lecture 5: Single-source model of spatial regression. Combining GIS and regional analysis (5) Spatial Econometrics 1 / 47 Outline 1 Linear model vs SAR/SLM (Spatial Lag) Linear model
More informationEksamen på Økonomistudiet 2006-II Econometrics 2 June 9, 2006
Eksamen på Økonomistudiet 2006-II Econometrics 2 June 9, 2006 This is a four hours closed-book exam (uden hjælpemidler). Please answer all questions. As a guiding principle the questions 1 to 4 have equal
More informationGeoDa and Spatial Regression Modeling
GeoDa and Spatial Regression Modeling June 9, 2006 Stephen A. Matthews Associate Professor of Sociology & Anthropology, Geography and Demography Director of the Geographic Information Analysis Core Population
More informationSPATIAL ECONOMETRICS: METHODS AND MODELS
SPATIAL ECONOMETRICS: METHODS AND MODELS STUDIES IN OPERATIONAL REGIONAL SCIENCE Folmer, H., Regional Economic Policy. 1986. ISBN 90-247-3308-1. Brouwer, F., Integrated Environmental Modelling: Design
More informationSummary of OLS Results - Model Variables
Summary of OLS Results - Model Variables Variable Coefficient [a] StdError t-statistic Probability [b] Robust_SE Robust_t Robust_Pr [b] VIF [c] Intercept 12.722048 1.710679 7.436839 0.000000* 2.159436
More informationLinear Probability Model
Linear Probability Model Note on required packages: The following code requires the packages sandwich and lmtest to estimate regression error variance that may change with the explanatory variables. If
More informationESE 502: Assignments 6 & 7
ESE 502: Assignments 6 & 7 Claire McLeod April 28, 2015 Many types of data are aggregated and recorded over a sub- region. Such areal data is particularly common in health and census data, where metrics
More informationFinding Hot Spots in ArcGIS Online: Minimizing the Subjectivity of Visual Analysis. Nicholas M. Giner Esri Parrish S.
Finding Hot Spots in ArcGIS Online: Minimizing the Subjectivity of Visual Analysis Nicholas M. Giner Esri Parrish S. Henderson FBI Agenda The subjectivity of maps What is Hot Spot Analysis? Why do Hot
More informationLuc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign
GIS and Spatial Analysis Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline GIS and Spatial Analysis
More informationRoger S. Bivand Edzer J. Pebesma Virgilio Gömez-Rubio. Applied Spatial Data Analysis with R. 4:1 Springer
Roger S. Bivand Edzer J. Pebesma Virgilio Gömez-Rubio Applied Spatial Data Analysis with R 4:1 Springer Contents Preface VII 1 Hello World: Introducing Spatial Data 1 1.1 Applied Spatial Data Analysis
More informationsphet: Spatial Models with Heteroskedastic Innovations in R
sphet: Spatial Models with Heteroskedastic Innovations in R Gianfranco Piras Cornell University Abstract This introduction to the R package sphet is a (slightly) modified version of Piras (2010), published
More informationBayesian Classification Methods
Bayesian Classification Methods Suchit Mehrotra North Carolina State University smehrot@ncsu.edu October 24, 2014 Suchit Mehrotra (NCSU) Bayesian Classification October 24, 2014 1 / 33 How do you define
More informationSingle Index Quantile Regression for Heteroscedastic Data
Single Index Quantile Regression for Heteroscedastic Data E. Christou M. G. Akritas Department of Statistics The Pennsylvania State University SMAC, November 6, 2015 E. Christou, M. G. Akritas (PSU) SIQR
More informationMultiple Regression Analysis: Heteroskedasticity
Multiple Regression Analysis: Heteroskedasticity y = β 0 + β 1 x 1 + β x +... β k x k + u Read chapter 8. EE45 -Chaiyuth Punyasavatsut 1 topics 8.1 Heteroskedasticity and OLS 8. Robust estimation 8.3 Testing
More informationProblem #1 #2 #3 #4 #5 #6 Total Points /6 /8 /14 /10 /8 /10 /56
STAT 391 - Spring Quarter 2017 - Midterm 1 - April 27, 2017 Name: Student ID Number: Problem #1 #2 #3 #4 #5 #6 Total Points /6 /8 /14 /10 /8 /10 /56 Directions. Read directions carefully and show all your
More informationA SPATIAL ANALYSIS OF A RURAL LAND MARKET USING ALTERNATIVE SPATIAL WEIGHT MATRICES
A Spatial Analysis of a Rural Land Market Using Alternative Spatial Weight Matrices A SPATIAL ANALYSIS OF A RURAL LAND MARKET USING ALTERNATIVE SPATIAL WEIGHT MATRICES Patricia Soto, Louisiana State University
More informationGMM - Generalized method of moments
GMM - Generalized method of moments GMM Intuition: Matching moments You want to estimate properties of a data set {x t } T t=1. You assume that x t has a constant mean and variance. x t (µ 0, σ 2 ) Consider
More informationTime-Series Regression and Generalized Least Squares in R*
Time-Series Regression and Generalized Least Squares in R* An Appendix to An R Companion to Applied Regression, third edition John Fox & Sanford Weisberg last revision: 2018-09-26 Abstract Generalized
More informationStatistical Inference. Part IV. Statistical Inference
Part IV Statistical Inference As of Oct 5, 2017 Sampling Distributions of the OLS Estimator 1 Statistical Inference Sampling Distributions of the OLS Estimator Testing Against One-Sided Alternatives Two-Sided
More informationChapter 4 Dimension Reduction
Chapter 4 Dimension Reduction Data Mining for Business Intelligence Shmueli, Patel & Bruce Galit Shmueli and Peter Bruce 2010 Exploring the data Statistical summary of data: common metrics Average Median
More informationCase of single exogenous (iv) variable (with single or multiple mediators) iv à med à dv. = β 0. iv i. med i + α 1
Mediation Analysis: OLS vs. SUR vs. ISUR vs. 3SLS vs. SEM Note by Hubert Gatignon July 7, 2013, updated November 15, 2013, April 11, 2014, May 21, 2016 and August 10, 2016 In Chap. 11 of Statistical Analysis
More informationRepresentation of Geographic Data
GIS 5210 Week 2 The Nature of Spatial Variation Three principles of the nature of spatial variation: proximity effects are key to understanding spatial variation issues of geographic scale and level of
More informationISSN Article
Econometrics 2014, 2, 217-249; doi:10.3390/econometrics2040217 OPEN ACCESS econometrics ISSN 2225-1146 www.mdpi.com/journal/econometrics Article The Biggest Myth in Spatial Econometrics James P. LeSage
More informationGraduate Econometrics Lecture 4: Heteroskedasticity
Graduate Econometrics Lecture 4: Heteroskedasticity Department of Economics University of Gothenburg November 30, 2014 1/43 and Autocorrelation Consequences for OLS Estimator Begin from the linear model
More informationGeoDa-GWR Results: GeoDa-GWR Output (portion only): Program began at 4/8/2016 4:40:38 PM
New Mexico Health Insurance Coverage, 2009-2013 Exploratory, Ordinary Least Squares, and Geographically Weighted Regression Using GeoDa-GWR, R, and QGIS Larry Spear 4/13/2016 (Draft) A dataset consisting
More informationDr Arulsivanathan Naidoo Statistics South Africa 18 October 2017
ESRI User Conference 2017 Space Time Pattern Mining Analysis of Matric Pass Rates in Cape Town Schools Dr Arulsivanathan Naidoo Statistics South Africa 18 October 2017 Choose one of the following Leadership
More informationWorking Paper No Introduction to Spatial Econometric Modelling. William Mitchell 1. April 2013
Working Paper No. 01-13 Introduction to Spatial Econometric Modelling William Mitchell 1 April 2013 Centre of Full Employment and Equity The University of Newcastle, Callaghan NSW 2308, Australia Home
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 15: Examples of hypothesis tests (v5) Ramesh Johari ramesh.johari@stanford.edu 1 / 32 The recipe 2 / 32 The hypothesis testing recipe In this lecture we repeatedly apply the
More informationSpatial analysis. Spatial descriptive analysis. Spatial inferential analysis:
Spatial analysis Spatial descriptive analysis Point pattern analysis (minimum bounding box, mean center, weighted mean center, standard distance, nearest neighbor analysis) Spatial clustering analysis
More informationBasics of Geographic Analysis in R
Basics of Geographic Analysis in R Spatial Autocorrelation and Spatial Weights Yuri M. Zhukov GOV 2525: Political Geography February 25, 2013 Outline 1. Introduction 2. Spatial Data and Basic Visualization
More informationUsing AMOEBA to Create a Spatial Weights Matrix and Identify Spatial Clusters, and a Comparison to Other Clustering Algorithms
Using AMOEBA to Create a Spatial Weights Matrix and Identify Spatial Clusters, and a Comparison to Other Clustering Algorithms Arthur Getis* and Jared Aldstadt** *San Diego State University **SDSU/UCSB
More informationAttribute Data. ArcGIS reads DBF extensions. Data in any statistical software format can be
This hands on application is intended to introduce you to the foundational methods of spatial data analysis available in GeoDa. We will undertake an exploratory spatial data analysis, of 1,387 southern
More informationIntroduction. Part I: Quick run through of ESDA checklist on our data
CSDE GIS Workshop Series Spatial Regression Chris Fowler csfowler@uw.edu Introduction The goal of these exercises is to give you a chance to put the concepts we have just discussed into practice. Keep
More informationFinding Outliers in Models of Spatial Data
Finding Outliers in Models of Spatial Data David W. Scott Dept of Statistics, MS-138 Rice University Houston, TX 77005-1892 scottdw@stat.rice.edu www.stat.rice.edu/ scottdw J. Blair Christian Dept of Statistics,
More informationSpatial autocorrelation: robustness of measures and tests
Spatial autocorrelation: robustness of measures and tests Marie Ernst and Gentiane Haesbroeck University of Liege London, December 14, 2015 Spatial Data Spatial data : geographical positions non spatial
More informationEconometrics. 9) Heteroscedasticity and autocorrelation
30C00200 Econometrics 9) Heteroscedasticity and autocorrelation Timo Kuosmanen Professor, Ph.D. http://nomepre.net/index.php/timokuosmanen Today s topics Heteroscedasticity Possible causes Testing for
More informationSpatial Econometric Analysis of the Hungarian Border Crossings
Spatial Econometric Analysis of the Hungarian Border Crossings Zsombor Szabó 1,*, Tibor Sipos 1, and Árpád Török 1 1 Budapest University of Technology an Economics, Department of Transport Technology and
More informationWeek 11 Heteroskedasticity and Autocorrelation
Week 11 Heteroskedasticity and Autocorrelation İnsan TUNALI Econ 511 Econometrics I Koç University 27 November 2018 Lecture outline 1. OLS and assumptions on V(ε) 2. Violations of V(ε) σ 2 I: 1. Heteroskedasticity
More informationEconometrics. Week 8. Fall Institute of Economic Studies Faculty of Social Sciences Charles University in Prague
Econometrics Week 8 Institute of Economic Studies Faculty of Social Sciences Charles University in Prague Fall 2012 1 / 25 Recommended Reading For the today Instrumental Variables Estimation and Two Stage
More informationCHAPTER 3 APPLICATION OF MULTIVARIATE TECHNIQUE SPATIAL ANALYSIS ON RURAL POPULATION UNDER POVERTYLINE FROM OFFICIAL STATISTICS, THE WORLD BANK
CHAPTER 3 APPLICATION OF MULTIVARIATE TECHNIQUE SPATIAL ANALYSIS ON RURAL POPULATION UNDER POVERTYLINE FROM OFFICIAL STATISTICS, THE WORLD BANK 3.1 INTRODUCTION: In regional science, space is a central
More informationLecture 7 Autoregressive Processes in Space
Lecture 7 Autoregressive Processes in Space Dennis Sun Stanford University Stats 253 July 8, 2015 1 Last Time 2 Autoregressive Processes in Space 3 Estimating Parameters 4 Testing for Spatial Autocorrelation
More informationIntroduction to Econometrics. Heteroskedasticity
Introduction to Econometrics Introduction Heteroskedasticity When the variance of the errors changes across segments of the population, where the segments are determined by different values for the explanatory
More informationUsing Spatial Statistics Social Service Applications Public Safety and Public Health
Using Spatial Statistics Social Service Applications Public Safety and Public Health Lauren Rosenshein 1 Regression analysis Regression analysis allows you to model, examine, and explore spatial relationships,
More informationHypothesis Testing hypothesis testing approach
Hypothesis Testing In this case, we d be trying to form an inference about that neighborhood: Do people there shop more often those people who are members of the larger population To ascertain this, we
More informationUnivariate analysis. Simple and Multiple Regression. Univariate analysis. Simple Regression How best to summarise the data?
Univariate analysis Example - linear regression equation: y = ax + c Least squares criteria ( yobs ycalc ) = yobs ( ax + c) = minimum Simple and + = xa xc xy xa + nc = y Solve for a and c Univariate analysis
More informationOutline. Overview of Issues. Spatial Regression. Luc Anselin
Spatial Regression Luc Anselin University of Illinois, Urbana-Champaign http://www.spacestat.com Outline Overview of Issues Spatial Regression Specifications Space-Time Models Spatial Latent Variable Models
More informationBig Value from Big Data: SAS/ETS Methods for Spatial Econometric Modeling in the Era of Big Data
Paper SAS535-2017 Big Value from Big Data: SAS/ETS Methods for Spatial Econometric Modeling in the Era of Big Data Guohui Wu and Jan Chvosta, SAS Institute Inc. ABSTRACT Data that are gathered in modern
More informationAn Application of Spatial Econometrics in Relation to Hedonic House Price Modelling. Liv Osland 1 Stord/Haugesund University College
An Application of Spatial Econometrics in Relation to Hedonic House Price Modelling Liv Osland 1 Stord/Haugesund University College 1 Bjørnsonsgt. 45, 5528 Haugesund, Norway (e-mail: liv.osland@hsh.no,
More information2008 ESRI Business GIS Summit Spatial Analysis for Business 2008 Program
A GIS Framework F k to t Forecast F t Residential Home Prices By Mak Kaboudan and Avijit Sarkar University of Redlands School of Business 2008 ESRI Business GIS Summit Spatial Analysis for Business 2008
More informationMacroeconomic Drivers of Global Food Security; A Spatial Econometrics Analysis
Macroeconomic Drivers of Global Food Security; A Spatial Econometrics Analysis Working Paper on Spatial effects on Global Food security. Presented by Philip Munyua 2008 commodity booms reignited the global
More informationCircle the single best answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 6, 2017 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 32 multiple choice
More informationAirport Noise in Atlanta: Economic Consequences and Determinants. Jeffrey P. Cohen, Ph.D. Associate Professor of Economics. University of Hartford
Airport Noise in Atlanta: Economic Consequences and Determinants Jeffrey P. Cohen, Ph.D. Associate Professor of Economics University of Hartford Cletus C. Coughlin, Ph.D. Vice President and Deputy Director
More informationEconometrics Homework 1
Econometrics Homework Due Date: March, 24. by This problem set includes questions for Lecture -4 covered before midterm exam. Question Let z be a random column vector of size 3 : z = @ (a) Write out z
More informationDealing with Heteroskedasticity
Dealing with Heteroskedasticity James H. Steiger Department of Psychology and Human Development Vanderbilt University James H. Steiger (Vanderbilt University) Dealing with Heteroskedasticity 1 / 27 Dealing
More informationModeling Spatial Relationships Using Regression Analysis
Esri International User Conference San Diego, California Technical Workshops July 24, 2012 Modeling Spatial Relationships Using Regression Analysis Lauren M. Scott, PhD Lauren Rosenshein Bennett, MS Answering
More informationAnalysis of multi-step algorithms for cognitive maps learning
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 62, No. 4, 2014 DOI: 10.2478/bpasts-2014-0079 INFORMATICS Analysis of multi-step algorithms for cognitive maps learning A. JASTRIEBOW
More information