Helder, boa tarde! Acredito que faltou um procedimento no script, que foi definir o objeto jura (origem dos dados para a K.O.). Avaliando o conjunto de dados carregado pelo data(), acredito que o correto seria usar o objeto jura.pred ou copiá-lo como jura. É necessário que os dados utilizem classes da {sp}, nesse caso 'SpatialPixelsDataFrame'. Modifiquei e comento a seguir: ### <BEGIN> library(gstat) ### para carregar data(jura), vgm() e krige() library(sp) ### para usar o coordinates() data(jura) ### irá adicionar os seguintes objetos #c("jura.grid", "jura.pred", "jura.val", "juragrid.dat","prediction.dat", "transect.dat", "validation.dat") jura.grid <- juragrid.dat ### 'reinicia' o objeto jura.grid copiando o juragrid.dat coordinates(jura.grid)=~Xloc+Yloc ### transforma em SpatialPixelsDataFrame - SPDF gridded(jura.grid) = TRUE ### informa que o SPDF é um grid class(jura.grid) ### os dados pra predição tb. precisam ser da classe SPDF coordinates(jura.pred)=~Xloc+Yloc ### transforma em SpatialPixelsDataFrame - SPDF class(jura.pred) m <- vgm(0.12645599,"Exp", 0.1593586, 0.04550649) x <- krige(log(Pb)~1, jura.pred, jura.grid, model=m) ### Visualizando... class(x) image(x) spplot(x[1], main='K.O. Predição') spplot(x[2], main='K.O. Erro') ### <END> Éder Comunello <c <comunello.eder@gmail.com>omunello.eder@gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]