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
documentation of "predict" in gstat it is mentioned as debug level 64
(order relation violations (indicator kriging values before and after
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
any pointers in the right direction
Thank you very much in advance
Best wishes
Dafne
[[alternative HTML version deleted]]