Caros da lista,

Com certeza tem algo nos dados, pois além do erro na krigagem apareceu este ao estimar por verossimillhança, com comando normal, sempre usei:

LIKFIT.reml <- likfit(gd, ini = ini.grid ,lik.method ="REML",nug =grid.pepita )

onde:
ini.grid é uma matriz com combinações do paramar e alcance
grid.pepita são várias "chutes" de nugge

gera este
 erro:

Erro em optim(par = c(1000, 0.04), fn = function (pars, fp, ip, temp.list)  : 
  valor não finito fornecido por optim

SE alguém puder ter uma ideia...


--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho