[R-br] Estimação dos parâmetros <likfit>

[Felinto COSTA]

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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141113/d8ce5757/attachment.html>