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