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 / [email protected]
_______________________________________________
R-br mailing list
[email protected]
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
[email protected]
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.
