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@gmail.com





_______________________________________________
R-br mailing list
R-br@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.

_______________________________________________
R-br mailing list
R-br@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.