[R-br] krige.conv
Felinto COSTA
incorpld em onda.com.br
Segunda Fevereiro 23 19:26:26 BRT 2015
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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150223/3e0b6bec/attachment.html>
Mais detalhes sobre a lista de discussão R-br