Prezados.
Tenho dois modelos teóricos (exponencial e esférico ) que foram ajustados aos dados por máxima verossimilhança.
As coordenadas numéricas (UTM) dos dados foram reduzidas pela subtração da coordenada média segundo cada eixo.
Em ambos considero a existência de uma tendência polinomial de 2ª ordem na média. 
Eles foram ajustados sem problemas, assim como a validação cruzada de ambos pode ser feita.
O problema está ao querer estimar algum ponto(s) usando a função krige.conv.  Apenas o modelo
exponencial é aceito. Ao usar o esférico a função retorna a seguinte mensagem:

krige.conv: model with mean given by a 2nd order polynomial on the coordinates
krige.conv: anisotropy correction performed
krige.conv: performing the Box-Cox data transformation
Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) :  sistema é computacionalmente singular: condição recíproca número = 1.76896e-16

Alguém saberia como contornar esse problema?

Antecipadamente agradeço a colaboração
Felinto COSTA