[R-br] Dúvida Tutorial Geoestatística

[Éder Comunello]

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 em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131104/c67a55ca/attachment.html>