Use the Running Reality desktop app or website to create your own videos, and customize how the map looks.

This is a book for Economists who regularly use spatial datasets for their econometric work.

The 'TopologyException: Input geom 1 is invalid' self-intersection error which arises from invalid polygon geometries has been widely discussed. However, I haven't found a convenient solution on the web that solely relies on R functionality. For instance, I have managed to create a 'SpatialPolygons' object from the output of map("state", ...) following Josh O'Brien's nice answer here. library(maps) library(maptools) map_states = map("state", fill = TRUE, plot = FALSE) IDs = sapply(strsplit(map_states$names, ":"), "[[", 1) spydf_states = map2SpatialPolygons(map_states, IDs = IDs, proj4string = CRS("+init=epsg:4326")) plot(spydf_states) The problem with this widely applied dataset is now that self-intersection occurs at the point given below. rgeos::gIsValid(spydf_states) [1] FALSE Warning message: In RGEOSUnaryPredFunc(spgeom, byid, "rgeos_isvalid") : Self-intersection at or near point -122.22023214285259 38.060546477866055 Unfortunately, this problem prevents any further use of 'spydf_states', e.g. when calling rgeos::gIntersection. How can I solve this issue from within R?

I am working on running an OLS on a shapefile and running tests on the residuals of that OLS where I am getting an error. library(maptools) library(sp) library(spdep) library(lmtest) library(spatialreg) library(rgdal) library(spstat) library(spatstat) library(spgwr) library(RColorBrewer) library(sf) library(dplyr) #read in my shapefile fires_OGR = readOGR(dsn = "C:/Users/westi/Desktop/unm_grad_21/580/project", layer="grid_1acre_fires") class(fires_OGR) class(fires_OGR$FIRE) = "Numeric" class(fires_OGR$CA_GAP) = "Numeric" class(fires_OGR$MEAN_PPT) = "Numeric" class(fires_OGR$MEAN_SLOPE) = "Numeric" class(fires_OGR$MEAN_ELEV) = "Numeric" class(fires_OGR$MEAN_VPD_M) = "Numeric" class(fires_OGR$MEAN_VPD_M.1) = "Numeric" class(fires_OGR$MEAN_ASPEC) = "Numeric" #row standardized "listw" object (default style = "W") fires_queen_nb = poly2nb(fires_OGR, queen=TRUE) fires_w = nb2listw(fires_queen_nb, zero.policy = TRUE) #ols for fire count per grid fires_ogr_OLS = lm(FIRE~CA_GAP + MEAN_PPT + MEAN_ELEV + MEAN_SLOPE + MEAN_VPD_M + MEAN_VPD_M.1 + MEAN_ASPEC, data = fires_OGR) summary(fires_ogr_OLS) #saving residuals and fitted values fires_OGR$olsresid = residuals(fires_ogr_OLS) fires_OGR$ols_fitted = fitted(fires_ogr_OLS) #moran test for residuals lm.morantest(fires_ogr_OLS, fires_w) #lagrange multiplier tests for residuals #fits a linear regression model for the squared residuals and examines #whether the fitted model is significant lm.LMtests(fires_ogr_OLS, fires_w, test = "all") #breusch pagan test for heteroskedasticity bptest(fires_ogr_OLS) Here is the error in the R terminal: I have tried making all the attributes I am using in my regression numeric since the error suggests that "y is not numeric" and when I run the "class(fires_OGR)" So I thought changing all the attributes used to numeric would fix the problem, but it did not. I guess the problem is I don't know what "y" is in this situation. Here is the summary from fires_OGR: Object of class SpatialPolygonsDataFrame Coordinates: min max x -120.69009 -115.78261 y 32.49987 35.02487 Is projected: FALSE proj4string : [+proj=longlat +datum=WGS84 +no_defs] Data attributes: CA_GAP MEAN_PPT MEAN_VPD_M MEAN_VPD_M.1 MAX_SLOPE MEAN_SLOPE MEAN_ASPEC MEAN_ELEV Length:2346 Min. : 61.84 Min. :0.697 Min. : 4.269 Min. : 0.00 Min. : 0.000 Min. : 0.0 Min. : -66.5 Class :character 1st Qu.: 285.70 1st Qu.:2.427 1st Qu.:18.219 1st Qu.:19.45 1st Qu.: 3.708 1st Qu.:160.9 1st Qu.: 220.7 Mode :character Median : 389.03 Median :3.544 Median :22.199 Median :29.23 Median : 7.937 Median :182.0 Median : 552.6 Mean : 404.13 Mean :3.923 Mean :22.441 Mean :26.65 Mean : 8.741 Mean :178.1 Mean : 687.7 3rd Qu.: 497.21 3rd Qu.:4.997 3rd Qu.:25.526 3rd Qu.:34.57 3rd Qu.:12.912 3rd Qu.:200.3 3rd Qu.:1077.2 Max. :1237.33 Max. :9.588 Max. :40.049 Max. :57.39 Max. :27.985 Max. :279.1 Max. :2904.3 MAJORITY_E FIRE mn_ppt mn_vpd_max mx_slope64 mn_slope mn_aspect Length:2346 Length:2346 Length:2346 Length:2346 Length:2346 Length:2346 Length:2346 Class :character Class :character Class :character Class :character Class :character Class :character Class :character Mode :character Mode :character Mode :character Mode :character Mode :character Mode :character Mode :character mn_elev maj_elev mn_vpd_min Length:2346 Length:2346 Length:2346 Class :character Class :character Class :character Mode :character Mode :character Mode :character

A port of the 'spReapportion' package, using Simple Features in order to lose the dependencies to the retired 'maptools' and 'rgeos' packages.

Most cloud-based mapping tools ask you to choose. You can have something simple that your whole team can use, or you can have something powerful that requires specialized training. GIS analysts know this tradeoff well. They’ve spent years learning desktop software that does exactly what they need, and they’ve watched colleagues struggle with those same […] The post Is Maptive Powerful Enough for GIS Experts & Analysts? first appeared on Flowster.

> > Couldn't reproduce. > Do you have any mods installed which replace the player's hand? (besides default from MTG) > There was a bugfix which can break wrongly defined hand replacements. https://github.com/luanti-org/luanti/pull/15743 > Nodecore also had a problem, I addressed this about two months ago, I hope they fixed it. > > If I understand right what mods do you mean, yes. There are: > terraform > treecapitator > tunnelmaker > creative (in mtg) > replacer > worldedit > maptools > markers > tape_measure maptools removed, all other mods tested in singleplayer, without problems... this will be hard swimming among 250 mods... Statistics: Posted by MatyasP — Sun May 25, 2025 14:56 * * *

> Couldn't reproduce. > Do you have any mods installed which replace the player's hand? (besides default from MTG) > There was a bugfix which can break wrongly defined hand replacements. https://github.com/luanti-org/luanti/pull/15743 > Nodecore also had a problem, I addressed this about two months ago, I hope they fixed it. If I understand right what mods do you mean, yes. There are: terraform treecapitator tunnelmaker creative (in mtg) replacer worldedit maptools markers tape_measure Statistics: Posted by MatyasP — Sun May 25, 2025 12:05 * * *

A port of the 'spReapportion' package, using Simple Features in order to lose the dependencies to the retired 'maptools' and 'rgeos' packages.

(4) Curso sobre el software estadístico R: La librería ma...