[R-sig-Geo] Indicator kriging map creation

Matevž Pavlič matevz.pavlic at gi-zrmk.si
Mon Jul 25 22:32:06 CEST 2011

Hi Edzer, 

sorry for the long time to  reply.

>It might be annoying, and after all, it is your sole responsibility to make sure the right labels end up in the right place. R can only prevent you from making mistakes that are too obvious.

What do you think would be the best thing to do about that? Predict the classes on 3D grid? I mean, i would like to create a code that would generate maps of several classes (variing numbers of them) each depth interval. I think that prediction on 3D grid would be the right thing?

>Not so likely; I think it comes with simple kriging -- rare classes often show weak autocorrelation, resulting in low estimates everywhere.
>Without having seen your data, I suspect your classes 9 and 10 are like that.

Attached are the co-variograms of ten classes. As i can see there are few classes (devet=level 10 , osem=level 9...) that have nugget (pure nugget effect). This is probably because there are only a few point with that class in the area :

devet   dva   ena   nic  osem   pet sedem   set stiri   tri 
    5    50    85   209     3     7    83    71    49    89

Is there a way to fix this without discarding this classes?


-------------- next part --------------
A non-text attachment was scrubbed...
Name: variogram.png
Type: image/png
Size: 26318 bytes
Desc: variogram.png
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110725/487b5c82/attachment.png>

More information about the R-sig-Geo mailing list