indicator kriging with R
Yes, this is what indicator kriging does. What you want is probably going back from the [0..1] estimated probabilities to the original Z values. Gstat doesn't help here, but it shouldn't be too hard to write a function for this. You need to make assumptions about the distribution between different cutoff values, in particular about both extreme (open) classes. I think the GSLIB and Goovaerts books would be the first ones to look into.
On 22/10/15 10:19, Dafne Berg wrote:
Dear list, I am trying to do a indicator kriging of a continuous variable. I am especially interested in performing the "order relation correction". In the documentation of "predict" in gstat it is mentioned as debug level 64 (order relation violations (indicator kriging values before and after order relation correction). However, when running predict I can only get the values in the range 0 to 1, as expected for indicator kriging. I am probably missing something very basic here, but I will be grateful for any pointers in the right direction Thank you very much in advance Best wishes Dafne [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster, Heisenbergstra?e 2, 48149 M?nster, Germany; +49 251 83 33081 Journal of Statistical Software: http://www.jstatsoft.org/ Computers & Geosciences: http://elsevier.com/locate/cageo/ Spatial Statistics Society http://www.spatialstatistics.info -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20151022/bb0c5089/attachment.bin>