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_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.