[R-br] Krigagem lognormal

[Helder Gramacho]

Olá,

Muito obrigado pela ajuda Prof. Paulo e Felinto.

*Hélder Gramacho *
Recife-PE /
*agrohelder em gmail.com <agrohelder em hotmail.com>*


Em 26 de novembro de 2014 15:46, Paulo Justiniano <paulojus em 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 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.
>>
>>
>>
>>
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141127/4ce6d303/attachment.html>