Boa noite, sou iniciante na Geoestatística e no R, estava praticando numa espécie de tutorial situado no seguinte link: http://jurachumbo.blogspot.com.br/p/resultados_06.html Estava executando os comandos até que apareceu um erro que não estou conseguindo entender, seguem os comandos e o erro: data(jura) jura.grid=juragrid.dat coordinates(jura.grid)=~Xloc+Yloc gridded(jura.grid) = TRUE m<-vgm(0.12645599,"Exp", 0.1593586, 0.04550649) x<-krige(log(Pb)~1,jura, jura.grid, model=m) Error in krige(log(Pb) ~ 1, jura, jura.grid, model = m) : erro na avaliação do argumento 'locations' na seleção do método para a função 'krige': 'Error: object 'jura' not found ' Alguém saberia me explicar qual é o problema? Desde já agradeço, Hélder Gramacho dos Santos Engenheiro Agrônomo Mestrando Ciênc. Geod. Tec. Geoinformação-UFPE *agrohelder@gmail.com <agrohelder@hotmail.com>*
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]
Parece que o posto no blog estṕa incompleto, sugiro contatar o autor do mesmo TALVEZ o desejado necessite incluir jura <- prediction.dat coordinates(jura) = ~Xloc+Yloc x <- krige(log(Pb)~1, jura, jura.grid, model = m) On Sat, 2 Nov 2013, Helder Gramacho wrote:
Boa noite,
sou iniciante na Geoestatística e no R, estava praticando numa espécie de tutorial situado no seguinte link: http://jurachumbo.blogspot.com.br/p/resultados_06.html
Estava executando os comandos até que apareceu um erro que não estou conseguindo entender, seguem os comandos e o erro:
data(jura) jura.grid=juragrid.dat coordinates(jura.grid)=~Xloc+Yloc gridded(jura.grid) = TRUE m<-vgm(0.12645599,"Exp", 0.1593586, 0.04550649) x<-krige(log(Pb)~1,jura, jura.grid, model=m) Error in krige(log(Pb) ~ 1, jura, jura.grid, model = m) : erro na avaliação do argumento 'locations' na seleção do método para a função 'krige': 'Error: object 'jura' not found ' Alguém saberia me explicar qual é o problema?
Desde já agradeço,
Hélder Gramacho dos Santos Engenheiro Agrônomo Mestrando Ciênc. Geod. Tec. Geoinformação-UFPE agrohelder@gmail.com
participantes (3)
-
Helder Gramacho -
Paulo Justiniano -
Éder Comunello