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.