[R-sig-Geo] indicator kriging with R

[Edzer Pebesma]

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>