using gwr for interpolation
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?