[R-sig-Geo] Log-normal block krige

[Edzer Pebesma]

Javier,

simulation is relatively easy: simulate point support fields on a fine 
grid, take the exponent of everything, compute the block mean of each 
simulation (e.g. using sp::aggregate), and take the statistics of the 
block means for each block. (There is just not yet a simple function 
doing this for you. Or does geoR have this?)

if you want to avoid simulation, you might want to take a look at the 
constrainedKriging package:

http://cran.r-project.org/web/packages/constrainedKriging/index.html


On 02/05/2013 11:30 AM, Javier Elio wrote:
> Dear all,
>
> I am trying to calculate the average value in a study area. I carried out the sampling according to a grid, and the data are lognormal. I would like to do with block krige (ordinary krige), I do the block krige with logarithms but I don't know how to back-transforming.
>
> I read Cressie 2006, but I don't know how to do in gstat (or other package, I aslo try constrainedKriging). I use krigeTg to back-transforming the predictions and for gaussian simulations, but as I know it can't be used for block krige.
>
> Please, can anyone help me?
>
> Best regards,
>
> Javier de Elío Medina
> Universidad Politécnica de Madrid
> Escuela de Ingenieros de Minas
>
> _______________________________________________
> 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
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