[R-br] Krigagem lognormal

[Paulo Justiniano]

me parece correta, mas talvez fazendopasso a passo coisas que a 
krige.conv() já faz 
internamente

O autor fez isto para rodar a xvalid() mas se for para predicao a ida e 
volta dele nao é necessária

Por exemplo:

library(geoR)
data(s100)
s100$data <- exp(s100$data)
reml <- likfit(s100, ini = c(0.5, 0.5), fix.nug = TRUE, lik.met = "REML",
     lambda = 0.1, fix.lambda = F)
grid.pred <- pred_grid(c(0,1), c(0,1), by=0.1)

preds <- krige.conv(s100, loc=grid.pred, krige=krige.control(obj.m=reml))

já vai retornar predicoes na escala original dos dados


On Wed, 26 Nov 2014, Felinto COSTA wrote:

> Professor:
> Essa rotina (não sei o autor) lhe parece apropriada?  
> http://rstudio-pubs-static.s3.amazonaws.com/743_cfad3b800ac140ba97740f4c678b6772.html
> 
> FCosta
> 
> ps. não há menção de seu autor
> 
> __________________________________________________________________________________________________________________________________________________________________
>
> 
> On 26/11/2014 13:21, Paulo Justiniano wrote:
>       Sim e sim
>
>       Existem expressoes analiticas para transformar
>       as predicoes e erros de predição para escala original
>       (procure log-normal kriging e volte a escrever se nao encontrar rapidamente)
>
>       Na geoR se voce usa
>       krige.conv()
>       com os dados na escala original e define lambda = 0
>       (ou passa um objeto com este ajuste)
>       a krigagem é feita na escala transformada e
>       retornada automaticamente para original
>
>       Em outras palavras, a geoR vai sempre te dar os resultados da krigagem na escala na qual voce passa os dados
> 
> 
>
>       On Wed, 26 Nov 2014, Helder Gramacho wrote:
>
>             Bom dia,
>
>             Quando numa variável foi aplicada a transformação logarítmica como é possível realizar a krigagem
>             lognormal no geo-r? e após a obtenção da estimativa como deve ser feita a transformação para a escala
>             original? 
>             Alguém possui um exemplo de como isso pode ser realizado?
>
>             Desde já muito obrigado,
>
>             Hélder Gramacho 
>             Recife-PE / agrohelder em gmail.com
> 
> 
> 
> 
> 
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
> 
> 
> 
>