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 for your contribution! Regards, Tom Paul Hiemstra wrote: > 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, he is my supervisor). If people are interested in this > package, it can be downloaded from > http://intamap.geo.uu.nl/~paul/Downloads.html. If you have any comments, > complaints or suggestions please do not hesitate to e-mail me. > > Some sample code: > library(automap) > data(meuse) > coordinates(meuse) = ~x+y > data(meuse.grid) > gridded(meuse.grid) = ~x+y > kr.ok = autoKrige(log(zinc)~1, meuse, meuse.grid) > plot(kr.ok) > kr.uk = autoKrige(log(zinc)~sqrt(dist), meuse, meuse.grid) > plot(kr.uk) > > cheers, > Paul > > -- Thomas E Adams National Weather Service Ohio River Forecast Center 1901 South State Route 134 Wilmington, OH 45177 EMAIL: [EMAIL PROTECTED] VOICE: 937-383-0528 FAX: 937-383-0033 _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo