<div dir="ltr"><div><font face="arial, helvetica, sans-serif">Helder, boa tarde!</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">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'.</font></div>
<div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">Modifiquei e comento a seguir:</font></div><div><br></div><div><div><font face="courier new, monospace">### <BEGIN></font></div>
<div><font face="courier new, monospace">library(gstat) ### para carregar data(jura), vgm() e krige()</font></div><div><font face="courier new, monospace">library(sp)    ### para usar o coordinates()</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">data(jura)      ### irá adicionar os seguintes objetos</font></div><div><font face="courier new, monospace">#c("jura.grid", "jura.pred", "jura.val", "juragrid.dat","prediction.dat", "transect.dat", "validation.dat")</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">jura.grid <- juragrid.dat         ### 'reinicia' o objeto jura.grid copiando o juragrid.dat</font></div><div><font face="courier new, monospace">coordinates(jura.grid)=~Xloc+Yloc ### transforma em SpatialPixelsDataFrame - SPDF</font></div>
<div><font face="courier new, monospace">gridded(jura.grid) = TRUE         ### informa que o SPDF é um grid</font></div><div><font face="courier new, monospace">class(jura.grid)</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">### os dados pra predição tb. precisam ser da classe SPDF</font></div><div><font face="courier new, monospace">coordinates(jura.pred)=~Xloc+Yloc ### transforma em SpatialPixelsDataFrame - SPDF</font></div>
<div><span style="font-family:'courier new',monospace">class(jura.pred)</span><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">m <- vgm(0.12645599,"Exp", 0.1593586, 0.04550649)</font></div>
<div><font face="courier new, monospace">x <- krige(log(Pb)~1, jura.pred, jura.grid, model=m)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Visualizando...</font></div>
<div><font face="courier new, monospace">class(x)</font></div><div><font face="courier new, monospace">image(x)</font></div><div><font face="courier new, monospace">spplot(x[1], main='K.O. Predição')</font></div><div>
<font face="courier new, monospace">spplot(x[2], main='K.O. Erro')</font></div><div><font face="courier new, monospace">### <END></font></div></div><div><br></div><div class="gmail_extra"><br clear="all"><div>
<div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a href="mailto:omunello.eder@gmail.com" target="_blank">omunello.eder@gmail.com</a>> <br>Dourados, MS - [22 16.5'S, 54 49'W]<br>
</div></div>
<br><br></div></div>