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 <comunello.eder@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]