Re: [R-sig-Geo] st_sample with raster

2019-08-10 Thread Edzer Pebesma
sf::st_sample has a size argument, which can be specified as a vector to operate per feature. If you give that a value equal to st_area(obj) * obj$required_density and maybe set exact = FALSE, you should get the densities you want. On 7/30/19 12:25 AM, Tyler Frazier via R-sig-Geo wrote: > Hi All

Re: [R-sig-Geo] regression-kriging and co-kriging

2019-08-10 Thread Edzer Pebesma
Hard to tell from your script. Maybe give a reproducible example? On 8/6/19 1:07 PM, Emanuele Barca wrote: > Dear r-sig-geo friends, > > I produced two maps garnered in the following way: > > # for regression-kriging > Piezo.map <-autoKrige(LivStat ~ Z, input_data = mydata.sp, new_data > = cov

[R-sig-Geo] Error when running a Spatial Durbin Error Model without intercept

2019-08-10 Thread Rolando Valdez
Dear list, I am trying to run a Spatial Durbin Error Model without intercept, however, I do get an error: In spatialreg::errorsarlm(formula = formula, data = data, listw = listw, : model configuration issue: no total impacts This is what I am doing: library(spatialreg) library(spdep) data(o

Re: [R-sig-Geo] Error when running a Spatial Durbin Error Model without intercept

2019-08-10 Thread Roger Bivand
On Fri, 9 Aug 2019, Rolando Valdez wrote: Dear list, I am trying to run a Spatial Durbin Error Model without intercept, however, I do get an error: In spatialreg::errorsarlm(formula = formula, data = data, listw = listw, : model configuration issue: no total impacts I can reproduce the pro

Re: [R-sig-Geo] Error when running a Spatial Durbin Error Model without intercept

2019-08-10 Thread Roger Bivand
And for completeness (https://github.com/r-spatial/spatialreg/issues/7): > dem <- lmSLX(eq, data = COL.OLD, listw = W) Error in lmSLX(eq, data = COL.OLD, listw = W) : object 'dirImps' not found In addition: Warning message: In lmSLX(eq, data = COL.OLD, listw = W) : model configuration issue: