Re: [R-sig-Geo] gstat::variogram - distance calculation is incorrect; kriging on sphere

2008-10-31 Thread Greg Lee
No sign of sp 0.9-27 at this end (Melbourne, Australia) yet. Interesting that it should take so long. I temporarily reset my repository to access the updates. Croatia already has sp 0.9-28, along with gstat 0.9-53. (Does Tomislav receive special treatment?) Your example now works perfectly. Thank

Re: [R-sig-Geo] gstat::variogram - distance calculation is incorrect; kriging on sphere

2008-10-31 Thread Edzer Pebesma
That's right; my fault that I sent the email and updated the sp package on CRAN almost at the same time. The new sp (0.9-27 would suffice) should now have propagated to your mirror, simply update your sp package again and it should work. Let me know if it doesn't. -- Edzer Greg Lee wrote: Hel

Re: [R-sig-Geo] gstat::variogram - distance calculation is incorrect; kriging on sphere

2008-10-30 Thread Greg Lee
Hello Edzer, I was curious to run the example you provided (using the latest CRAN versions of all packages), but as written the line > idw.spdf = as(idw.out, "SpatialPolygonsDataFrame") produces Error in as(idw.out, "SpatialPolygonsDataFrame") : no method or default for coercing "SpatialPixels

Re: [R-sig-Geo] gstat::variogram - distance calculation is incorrect; kriging on sphere

2008-10-30 Thread Edzer Pebesma
Dear all, In continuation of this thread, I've spent some time looking at kriging on the sphere, corrected some bugs, and need further help. First of all, distances for covariances on the sphere were computed incorrectly as well in gstat, so I hope not too many people have been relying on th

Re: [R-sig-Geo] gstat::variogram - distance calculation is incorrect

2008-09-14 Thread Edzer Pebesma
Thanks for the bug report; indeed, distance computation for longlat data variograms went wrong so far. I just uploaded a new gstat version to /incoming on CRAN that seems to repair this bug. Next thing I'm not 100% sure of is the direction selections for directional variograms of longlat data.

[R-sig-Geo] gstat::variogram - distance calculation is incorrect

2008-09-14 Thread Timothy W. Hilton
Hello, I am trying to use 'variogram' from the gstat package to compute an empirical semivariogram for a dataset that spans the North American continent. I have been unable to get gstat to calculate the great circle distances correctly. For example, in the example below, the first two lines o