[R-sig-Geo] using gwr for interpolation

[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

On Friday 29 May 2009 11:36:08 am 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$SDF["gwr.e"])
>
> But it doesn't give what I'd expected - somehow nicely interpolated
> cadmium values. I probably misused fit.points, but I couldn't see any
> predict method. Is this possible at all?