[R-sig-Geo] Gaussian transformation

Carlo Innocenti carlo.innocenti at isprambiente.it
Thu Apr 30 14:16:57 CEST 2015


Thank you for your answers.

In the meantime, I found the package RGeostats developed by Geosciences
Research Center of MINES ParisTech on http://rgeostats.free.fr and at the
moment I'm studying it.
It seems promising with 3D geostatistics capability and function for
gaussian anamorphosis transformation and back transformation.

Carlo




2015-04-29 15:11 GMT+02:00 Paulo Justiniano Ribeiro Jr <
paulojus at c3sl.ufpr.br>:

> Perhaps you can also look at the Box-Cox family of transformation for
>> which some cases such as log is a particular case
>>
>> You may look at MASS::boxcox()
>> and/or boxcox related functions in geoR
>>
>> If log works well you benefit from the results of log-normal kriging.
>>
>>
>
> Paulo Justiniano Ribeiro Jr
> LEG (Laboratorio de Estatistica e Geoinformacao)
> Universidade Federal do Parana
> Caixa Postal 19.081
> CEP 81.531-990
> Curitiba, PR  -  Brasil
> Tel: (+55) 41 3361 3573
> VOIP: (+55) (41) (3361 3600) 1053 1066
> Fax: (+55) 41 3361 3141
> e-mail: paulojus AT  ufpr  br
> http://www.leg.ufpr.br/~paulojus
>
>
> On Wed, 29 Apr 2015, Carlo Innocenti wrote:
>
>  Thank you Edzer.
>>
>> It works, but I need the back transformation as well.
>> What about the gaussian anamorphosis with hermite polynomials? I know that
>> is a method used to normalize the data before to apply the conditional
>> gaussian simulation.
>> There is a package in R to transform and back transform with this method?
>> Another method that is used in these cases should be the normal scores
>> transformation: is it implemented in some package?
>>
>> Carlo
>>
>>
>>
>>
>> 2015-04-28 22:32 GMT+02:00 Edzer Pebesma <edzer.pebesma at uni-muenster.de>:
>>
>>  one way is the rank transform:
>>>
>>> x = rgamma(100, 0.1)
>>> hist(x)
>>> y = qnorm((order(x)/(length(x) + 0.5)))
>>> qqnorm(y)
>>>
>>> but transforming your conditionally simulated values back will be
>>> tricky. And I'm not sure about the 0.5.
>>>
>>> On 04/28/2015 06:09 PM, Carlo Innocenti wrote:
>>> > I would like to perform a gaussian conditional simulation on a variable
>>> > that is highly skewed.
>>> > Could someone suggest me how to transform my variable in a gaussian
>>> > distribution?
>>> >
>>> > Thanks,
>>> > Carlo
>>> >
>>> >
>>>
>>> --
>>> 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
>>>
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> R-sig-Geo at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>>
>>>
>>
>> --
>>
>>
>> Carlo Innocenti, PhD
>>
>> ISPRA
>> Istituto Superiore per la Protezione e la Ricerca Ambientale
>> High Institute for Environmental Protection and Research
>>
>> Via Vitaliano Brancati, 60
>> 00144 Roma
>> Italy
>>
>> phone: +39 06 5007 4639
>> email: carlo.innocenti at isprambiente.it
>>
>>         [[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
>
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>


-- 


Carlo Innocenti, PhD

ISPRA
Istituto Superiore per la Protezione e la Ricerca Ambientale
High Institute for Environmental Protection and Research

Via Vitaliano Brancati, 60
00144 Roma
Italy

phone: +39 06 5007 4639
email: carlo.innocenti at isprambiente.it

	[[alternative HTML version deleted]]



More information about the R-sig-Geo mailing list