Re: [R-sig-Geo] spgrass6 and R 2.9

2009-05-29 Thread Torleif Markussen Lunde
Hi First try updating glibc, and install spgrass6 again. If this does not work you could try export MALLOC_CHECK_=0 before you start R. The error might be related to "somewhere something" is trying to free memory that has already been unallocated. I don't know the details here, but it could

[R-sig-Geo] spgrass6 and R 2.9

2009-05-29 Thread Dylan Beaudette
Hi, Just upgraded to R 2.9.0 this morning, and tried to update all of my packages. Everyhing seemed to work, except spgrass6. Manually removing the old packing and installing fresh resulted in: install.packages('spgrass6', dep=TRUE) trying URL 'http://cran.cnr.Berkeley.edu/src/contrib/spgrass6_

[R-sig-Geo] Spatial regression

2009-05-29 Thread Adrian Toti
Hi everyone, In R, the function lm.LMtests is used to select which spatial model fits better the data, an error or lag one. It looks to me that there are some “limitations” using this function like being constraint to use only row-standardized spatial weights and it does not support weighted lm ob

Re: [R-sig-Geo] using gwr for interpolation

2009-05-29 Thread Roger Bivand
On Fri, 29 May 2009, Torleif Markussen Lunde wrote: Hi I think the error is in .GWR_int line 42: df[i, (m + 3)] <- ei[i] ei will have length equal to dim(data)[1], while i is equal to nrow(fit.points). So as i > dim(data)[1] ei[i] will be NA. This could explain why only the first 155 points a

Re: [R-sig-Geo] using gwr for interpolation

2009-05-29 Thread Roger Bivand
On Fri, 29 May 2009, Edzer Pebesma wrote: Can I use gwr for interpolation? I would for example use something like library(spgwr) data(meuse) coordinates(meuse) = ~x+y data(meuse.grid) gridded(meuse.grid) = ~x+y x = gwr(cadmium ~ dist, meuse, bandwidth = 228, fit.points = meuse.grid) spplot(x$S

Re: [R-sig-Geo] using gwr for interpolation

2009-05-29 Thread Torleif Markussen Lunde
Hi I think the error is in .GWR_int line 42: df[i, (m + 3)] <- ei[i] ei will have length equal to dim(data)[1], while i is equal to nrow(fit.points). So as i > dim(data)[1] ei[i] will be NA. This could explain why only the first 155 points are evaluated (in your example). Best wishes Torleif

[R-sig-Geo] using gwr for interpolation

2009-05-29 Thread Edzer Pebesma
Can I use gwr for interpolation? I would for example use something like library(spgwr) data(meuse) coordinates(meuse) = ~x+y data(meuse.grid) gridded(meuse.grid) = ~x+y x = gwr(cadmium ~ dist, meuse, bandwidth = 228, fit.points = meuse.grid) spplot(x$SDF["gwr.e"]) But it doesn't give what I'd ex

Re: [R-sig-Geo] Problems loading hydrosanity package

2009-05-29 Thread Roger Bivand
On Fri, 29 May 2009, Felix Andrews wrote: 2009/5/29 Roger Bivand : On Thu, 28 May 2009, brwin...@aol.com wrote: Good Morning: I need to construct a set of Thiessen Polygons centered around approximately 4500 longlat positions in the US. After searching the R site, I was going to attempt to