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] using gwr

2008-11-20 Thread Roger Bivand
On Thu, 20 Nov 2008, Virginia wrote: Hi list, I would like to know if R does kernel smooth with borders restriction. Please use the list as described in the posting guide for all R lists: http://www.r-project.org/posting-guide.html You have asked a new question, but have not started a new t

Re: [R-sig-Geo] using gwr

2008-11-20 Thread Virginia
Hi list, I would like to know if R does kernel smooth with borders restriction. Virginia ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Re: [R-sig-Geo] using gwr

2008-11-11 Thread Danlin Yu
I can't read this Brazilian (Portugese :)), but it seems that the distance calculation is in error (possibly some coordinates were missing causing distance calculation yield NA/NaN/Inf value), which leads to distance not being able to calculated. Like Edzer mentioned, if a sample set of your d

Re: [R-sig-Geo] using gwr

2008-11-11 Thread Roger Bivand
On Tue, 11 Nov 2008, Virginia wrote: I have some questions: I got an error. I dont't know why. Somebody to help me? Your code template is not helpful. A more recent version below, which will probably show that your input data are in error: Thanks, Virginia Ragoni require(spgwr) require

Re: [R-sig-Geo] using gwr

2008-11-11 Thread Edzer Pebesma
It is no problem for us to read all this Brazilian (eh, sorry, Portugese), but it is hard to reproduce your error without having access to your, or similar, data. -- Edzer Raphael Saldanha wrote: Olá Virgínia! Legal ver brazileiros nessa lista! Vou dar uma olhada no seu problema e lhe escrevo

Re: [R-sig-Geo] using gwr

2008-11-11 Thread Raphael Saldanha
Olá Virgínia! Legal ver brazileiros nessa lista! Vou dar uma olhada no seu problema e lhe escrevo. On Tue, Nov 11, 2008 at 10:06 AM, Virginia <[EMAIL PROTECTED]> wrote: > I have some questions: > > I got an error. I dont't know why. > > Somebody to help me? > > Thanks, > > Virginia Ragoni > > r

[R-sig-Geo] using gwr

2008-11-11 Thread Virginia
I have some questions: I got an error. I dont't know why. Somebody to help me? Thanks, Virginia Ragoni require(spgwr) require(foreign) x <- read.shape("cidade.shp") plot(x) mappolys <- Map2poly1(x,as.character(x$att.data$ID9100),raw=FALSE) mapbbs=Map2bbs(x) mapcents=get.Pcent(x) tab= x$att