{"id":430,"date":"2012-05-04T12:15:56","date_gmt":"2012-05-04T16:15:56","guid":{"rendered":"http:\/\/web.colby.edu\/mgimond\/?p=430"},"modified":"2012-05-04T12:16:50","modified_gmt":"2012-05-04T16:16:50","slug":"430","status":"publish","type":"post","link":"https:\/\/web.colby.edu\/mgimond\/2012\/05\/04\/430\/","title":{"rendered":"Is Moran’s I robust?"},"content":{"rendered":"

Moran’s I (both global and local) is a measure of spatial autocorrelation. \u00a0One form of its equation looks like this:<\/p>\n I = \\frac{N}{\\sum_{i} \\sum_{j} w_{ij}} \\frac{\\sum_{i} \\sum_{j} w_{ij} (X_i - \\bar X)(X_j - \\bar X)}{\\sum_{i} (X_i - \\bar X )^2}\n

Note the residual and spread terms. As such, Moran’s I can be sensitive to outliers. The following example demonstrates this point:<\/p>\n

\r\nlibrary(spdep)\r\nlibrary(maptools)\r\n# Load NC data\r\nNC <- readShapePoly(system.file("etc\/shapes\/sids.shp", package="spdep")[1],\r\n      ID="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))\r\nrn <- sapply(slot(NC, "polygons"), function(x) slot(x, "ID"))\r\nNB <- read.gal(system.file("etc\/weights\/ncCR85.gal", package="spdep")[1], \r\n      region.id=rn)\r\nn  <- length(NB)\r\n# Define some variables\r\nset.seed(4956)\r\nx.norm <- rnorm(n)     # Normal distribution\r\nrho    <- 0.3          # autoregressive parameter\r\nW      <- nb2listw(NB) # Generate spatial weights\r\n# Generate spatially autocorrelated datasets \r\n# (one normally distributed the other skewed)\r\nx.norm.auto <- invIrW(W, rho) %*% x.norm # Generate autocorrelated values\r\nx.skew.auto <- exp(x.norm.auto) # Transform orginal data to create a 'skewed' version\r\n<\/pre>\n
You can map the simulated data as follows:<\/p>\n
\r\n# Append data to NC polygon object\r\nauto_sk <- as.vector(x.skew.auto)\r\nauto_no <- as.vector(x.norm.auto)\r\nrn      <- as.integer(rn)\r\nNCa     <- spCbind(NC,auto_sk)\r\nNCa     <- spCbind(NCa,auto_no)\r\nNCa     <- spCbind(NCa,rn)   # Needed if the GAL file is to be used in GeoDa\r\n# Plot map (skewed data in this example)\r\nbrks <- quantile(auto_sk, seq(0,1,1\/10))\r\ncols <- grey.colors(length(brks)-1, 0.95, 0.55, 2.2)\r\nplot(NCa, col=cols[findInterval(auto_sk, brks, all.inside=TRUE)])\r\n<\/pre>\n
\"\"<\/a><\/p>\n
\r\n# Compute Moran's I for both datasets\r\nI.norm = moran.test(x.norm.auto, listw=W)\r\nI.skew = moran.test(x.skew.auto, listw=W)\r\nI.norm$p.value;I.skew$p.value\r\n# The above approach makes assumptions about the Moran's I when in fact, it may not be.\r\n# To correct for this, we can generate an empirical sampling distribution function\r\n# using a Monte Carlo procedure.\r\nMCI.norm <- moran.mc(x.norm.auto, listw=W, nsim=9999)\r\nMCI.skew <- moran.mc(x.skew.auto, listw=W, nsim=9999)\r\n# Compare the P-values\r\nMCI.norm$p.value;MCI.skew$p.value\r\n<\/pre>\n
The p-values from the Monte Carlo simulations are 0.0133 and 0.1672 for the normal and skewed data respectively. The test results under different data distributions provide us with two distinct conclusions at a 5% significance level!<\/p>\n
\r\noopar <- par(mfrow=c(1,2), mar=c(4,4,3,2)+0.1)\r\nmoran.plot(x.norm.auto,W)\r\nmoran.plot(x.skew.auto,W)\r\npar(oopar)\r\n<\/pre>\n
\"\"<\/a> <\/p>\n
NOTE: To export the spatial data to a shapefile for further analysis in GeoDa:<\/p>\n
\r\n# Write map object to shapefile\r\nwritePolyShape(NCa,"NCa",factor2char = F)\r\n<\/pre>\n
–> If using Geoda, the weight file can be found under C:\\Program Files\\R\\R-2.12.2\\library\\spdep\\etc\\weights<\/p>\n","protected":false},"excerpt":{"rendered":"
Moran’s I (both global and local) is a measure of spatial autocorrelation. \u00a0One form of its equation looks like this: Note the residual and spread terms. As such, Moran’s I can be sensitive to outliers. The following example demonstrates this point: You can map the simulated data as follows: The p-values from the Monte Carlo […]<\/p>\n","protected":false},"author":1199,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"ngg_post_thumbnail":0,"footnotes":""},"categories":[12807,12808],"tags":[],"_links":{"self":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/430"}],"collection":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/users\/1199"}],"replies":[{"embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/comments?post=430"}],"version-history":[{"count":59,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/430\/revisions"}],"predecessor-version":[{"id":520,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/posts\/430\/revisions\/520"}],"wp:attachment":[{"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/media?parent=430"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/categories?post=430"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/web.colby.edu\/mgimond\/wp-json\/wp\/v2\/tags?post=430"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}