Professor: Essa rotina (não sei o autor) lhe parece apropriada? http://rstudio-pubs-static.s3.amazonaws.com/743_cfad3b800ac140ba97740f4c678b... 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.

Olá, Muito obrigado pela ajuda Prof. Paulo e Felinto. *Hélder Gramacho * Recife-PE / *agrohelder@gmail.com <agrohelder@hotmail.com>* Em 26 de novembro de 2014 15:46, Paulo Justiniano <paulojus@leg.ufpr.br> escreveu:

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.