### <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>