Prezados.
Preciso estimar os parâmetros do componente espacial (e do erro) da
função de correlação usando a função <likfit> com a remoção de tendência
da média.
Para isso estou inserindo uma função para esta tendência
<trend=~coords[,1]+coords[,2]+
I(coords[,1]^2)+I(coords[,2]^2)+coords[,1]*coords[,2]>, mas
que parece ser incorreta para esta situação pelo aviso de erro que retorna.
Seria o caso de estimar os resíduos partindo da regressão da matriz de
tendência (via <trend.spatial>) pelos valores observados e daí sim
usá-los como o novo <geodata> a função <likfit>?
Uma vez estimados os parâmetros da função usando os resíduos, a
aplicação dela para realizar estimativas
em pontos não amostrados deverá ser sobre o <geodata> original?
O o <geodata> dos resíduos, somando-se a esta estimativa a parte
estimada pelo polinômio de 2ª ordem?
Antecipadamente agradeço a atenção de todos.
_______________________________________________
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.
