Re: [R-sig-Geo] Package for automatic interpolation

2008-03-12 Thread Edzer Pebesma
I agree with Thomas that seeing a warning is a worry. If Paul's automap takes care of singular fits by discarding them, it should try to hide the warning for the user, because otherwise it will worry them. I see that the warning is issued by the C code in gstat, and will rethink about how a wra

Re: [R-sig-Geo] Package for automatic interpolation

2008-03-12 Thread Paul Hiemstra
Hi Thomas, The automatic fitting procedure checks a range of possible variogram models, including Spherical, Exponential, Gaussian and Matern models. Each of these models could lead to a singular model fit. The fitting function in gstat (fit.variogram) says the following on singular model fits:

Re: [R-sig-Geo] Package for automatic interpolation

2008-03-12 Thread Thomas Adams
Paul, This is outstanding; thank you very much! BTW, I noticed that when I run your example that I get: Warning: singular model in variogram fit …actually, 2 such warnings with the OK example, and 1 with the UK example. Are you seeing these too? Otherwise, the results look fine… Thanks again

[R-sig-Geo] Package for automatic interpolation

2008-03-11 Thread Paul Hiemstra
Dear r-sig-geo, In the course of my PhD I have written a R package that performs automatic interpolation (kriging) and might be interesting for other people. I work on real-time automatic interpolation of radioactivity levels. It depends heavily on the gstat pacakge written by Edzer (Pebesma,