[R-br] likfit (geoR), normalidade da variável e do resíduo do modelo Trend

[Éder Comunello]

Pessoal, bom dia!

Até onde sei, a normalidade não é uma prerrogativa para os modelos
geoestatísticos. Mas haveria sim uma preocupação com a simetria da
distribuição dos dados, uma vez que a presença de "caudas" poderia
interferir no ajustamento.  Acredito que a ideia do Box-Cox seja melhorar a
distribuição dos dados, sem necessariamente normalizá-los.

Por outro lado é necessário ter estacionariedade. Quando você adiciona o
efeito de tendência no modelo (trend) a ideia é tentar considerar (ou
"retirar") um componente não-estacionário da série, obtendo um resíduo
estacionário (retratando adequadamente  a dependência espacial), que seria
apto para o ajustamento dos modelos geoestatísticos. Assim como o dado
original, acredito que não deve haver pressuposto da normalidade,
devendo-se atentar para a simetria da distribuição.

No caso do pH, lembrando que é uma medida de natureza logarítmica, talvez
seja prudente testar lambda=0 que equivaleria à transformação logarítmica.
O comando likfit() também recebe o argumento lambda, então dá pra testar a
sugestão do Wagner com diferentes lambdas. Os comandos demonstram
preocupação com anisotropia. É provável que o efeito seja minimizado ao
considerar efeitos de tendência.

Pelo que lembro, o krige.conv retorna os valores na escala de entrada. Se
você indicou o trend e lambda no modelo, eles serão considerados no ajuste,
mas o resultado final já sai "transformado de volta".


​
================================================
Éder Comunello
Agronomist (UEM), MSc in Environ. Sciences (UEM)
DSc in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil |<O>|
================================================
GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00




Em 2 de abril de 2016 11:07, Wagner Wolff <wwolff em usp.br> escreveu:

> Olá pessoal a dúvida é a seguinte.
> Preciso interpolar uma variável na qual não consegui normalizá-la pelo
> BoxCox(); entretanto, a análise exploratória de tendência espacial indicou
> um modelo para remoção da mesma; assim, passo a interpolar os resíduos
> deste modelo, correto?
> Sendo assim, posso desconsiderar a não normalidade da variável inicial e
> apenas considerar a normalidade dos resíduos do modelo Trend?
> Outra questão é se realmente preciso normalizar a Variável, é plausível
> aceitar a normalidade dela por este gráfico abaixo? Acabei aceitando esta
> condição e não tive problemas na obtenção dos meus parâmetros.
>
>
>
> ​Att.
> Wagner Wolff
>
> _______________________________________________
> 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/20160404/1ec6ab1d/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: 2016-04-02-115136_1366x768_scrot.png
Tipo: image/png
Tamanho: 35494 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160404/1ec6ab1d/attachment.png>