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? _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
