[R-sig-Geo] Back-transforming krigged variables

Edzer Pebesma edzer.pebesma at uni-muenster.de
Mon Mar 1 13:26:23 CET 2010


Ebrahim Jahanshiri wrote:
> I am wondering if there is any function or option to back-transform the
> result of the lognormal kriging in gstat yet.
>   
You can use krigeTg, which does trans-Gaussian kriging for all Box-Cox 
transforms; passing it a lambda parameter with value 0, log-normal 
kriging is used. Pls look at the examples, the back-transformed values 
have funny names like var1TG.pred.

The whole think has been implemented for ordinary kriging only; if 
someone can point me out a source that describes how to do this for 
universal / external drift kriging, I'd be happy to look into how to 
implement it.
> Although for simple kriging it is fairly straight forward, that is, it is
> possible to do back-transform using the vector of predicted values and
> kriging variances, but for ordinary kriging, it might be difficult because
> we dont have lagrange multiplier (according to my web help search) to do the
> back-transformation. I am wondering if there has been any advancement in
> that matter?
>
> thank you
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>   

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster 
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251 
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de 
http://www.52north.org/geostatistics      e.pebesma at wwu.de



More information about the R-sig-Geo mailing list