Re: [R-sig-Geo] kernel density estimation to google earth

2008-01-27 Thread Tomislav Hengl
# Modified by T. Hengl (http://spatial-analyst.net) # 27.01.2008 library(rgdal) library(maptools) library(splancs) # Import the points and study area: data.shp - readOGR(C:/, layer=events) str(data.shp) poly.shp - readOGR(C:/, layer=hull) str(poly.shp) poly -

Re: [R-sig-Geo] kernel density estimation to google earth

2008-01-27 Thread Roger Bivand
On Sun, 27 Jan 2008, Tomislav Hengl wrote: # Modified by T. Hengl (http://spatial-analyst.net) # 27.01.2008 Thanks, very clear - the kernel values at the grid of projected points are not an a grid in geographical coordinates, and need to be interpolated (warped) to a grid in geographical

Re: [R-sig-Geo] Regression kriging

2008-01-27 Thread Edzer Pebesma
Jose, is your model linear, or are you using a generalized linear model? The questions is not so much: model parameters before or after kriging residuals, but rather: model parameters under the assumption of independent observations (the usual regression approach), or model parameters under the

Re: [R-sig-Geo] kernel density estimation to google earth

2008-01-27 Thread G. Allegri
Thanks, a really interesting post. May I ask a stupid question? Where are the following function documented? getPolygonCoordsSlot(getPolygonsPolygonsSlot(getSpPpolygonsSlot(poly.shp )[[1]])[[1]]) They are in the sp package but I can't find references but in some posts on this ML. Thanks,