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