<div dir="ltr">Caros Elias e Éder.<div><br></div><div>Grato pelas opiniões e exemplos.</div><div><br></div><div>Vou testar com meus dados.</div><div><br></div><div>Mais uma vez, agradeço imensamente!!!</div><div><br></div>
<div>Hélio</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 29 de maio de 2014 07:39, Eder Comunello [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4662232h60@n4.nabble.com" target="_blank">ml-node+s2285057n4662232h60@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

        <div dir="ltr"><div>Hélio, bom dia!</div><div><br></div><div>Faz sentido interpolar? Acredito que sim, partindo do princípio que a dependência espacial pode existir, mas seu método não permitiu detectá-la. Acredito que a causa mais comum seja o esquema de amostragem utilizado, sobretudo no que se refere à distância entre amostras. No caso particular, você conseguiu modelar em algumas áreas e não em outras, dando margem à essa interpretação.</div>

<div><br></div><div>O colega Elias já postou uma solução, mas acrescento o código abaixo, caso ainda tenha interesse em utilizar o {gstat}.</div><div><br></div><div>Em termos gerais, a ideia de operação no {gstat} é similar a do {geoR}, no sentido em que você vai precisar criar um grid pra receber o resultado da interpolação. O espaçamento da grade será a resolução da interpolação. O ponto principal é que tem que trabalhar com objetos da classe 'sp'.</div>

<div><br></div><div>Verifique os parâmetros número de vizinhos (nmax) e peso/potência da  distância (idp). No {gstat} por default utilizam-se todos os pontos e o idp=2.</div><div><br></div><div>Atte.,</div><br><div><div>
<font face="courier new, monospace">### <code r></font></div>
<div><font face="courier new, monospace">sapply(c("gstat", "sp", "geoR", "RColorBrewer"), require, character.only=T)</font></div><div><font face="courier new, monospace"><br></font></div>

<div><font face="courier new, monospace">data(parana); names(parana)</font></div><div><font face="courier new, monospace">hist(parana$data, col=3)</font></div><div><font face="courier new, monospace">points(parana, pt.divide='quart')</font></div>

<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Criar objeto 'sp'</font></div><div><font face="courier new, monospace">pr <- data.frame(x=parana$coords[,1], y=parana$coords[,2], chuva=parana$data)</font></div>

<div><font face="courier new, monospace">names(pr); coordinates(pr) <- ~x+y</font></div><div><font face="courier new, monospace">class(pr)</font></div><div><font face="courier new, monospace">plot(pr, asp=1, axes=T, pch=20); polygon(parana$borders, border=2)</font></div>

<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Criar 'grid'</font></div><div><font face="courier new, monospace">bbox(pr)</font></div><div><font face="courier new, monospace">lim   <- round(bbox(pr)+c(-1,-1, 1, 1)*50); lim ### amplia a área do bbox</font></div>

<div><font face="courier new, monospace">grid  <- expand.grid(x=seq(lim[1,1],lim[1,2], by=10), y=seq(lim[2,1],lim[2,2], by=10))</font></div><div><font face="courier new, monospace">grid  -> <a href="http://grid.pt" rel="nofollow" link="external" target="_blank">grid.pt</a></font></div>

<div><font face="courier new, monospace">gridded(grid) = ~x+y        ### SpatialPixels</font></div><div><font face="courier new, monospace">coordinates(<a href="http://grid.pt" rel="nofollow" link="external" target="_blank">grid.pt</a>) = ~x+y ### SpatialPoints</font></div>

<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### IDW Default</font></div><div><font face="courier new, monospace">idw    <- idw(pr$chuva~1, pr, grid)</font></div><div>

<font face="courier new, monospace"><a href="http://idw.pt" rel="nofollow" link="external" target="_blank">idw.pt</a> <- idw(chuva~1, pr, <a href="http://grid.pt" rel="nofollow" link="external" target="_blank">grid.pt</a>)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Visualização</font></div>
<div><font face="courier new, monospace">spplot(idw,    "var1.pred", main = "IDW Default")</font></div><div><font face="courier new, monospace">spplot(<a href="http://idw.pt" rel="nofollow" link="external" target="_blank">idw.pt</a>, "var1.pred", main = "IDW default - Mapa Pontuado")</font></div>

<div><font face="courier new, monospace">image(idw); polygon(parana$borders); points(pr)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Variando parâmetros</font></div>

<div><font face="courier new, monospace">idw.data <- data.frame(</font></div><div><font face="courier new, monospace">          default = idw(chuva ~ 1, pr, grid)$var1.pred, ### nmax=todos & idp=2</font></div><div>

<font face="courier new, monospace">          idw6    = idw(chuva ~ 1, pr, grid, nmax=6)$var1.pred,</font></div><div><font face="courier new, monospace">          idw9    = idw(chuva ~ 1, pr, grid, nmax=9)$var1.pred,</font></div>

<div><font face="courier new, monospace">          idp1    = idw(chuva ~ 1, pr, grid, idp = 1)$var1.pred,</font></div><div><font face="courier new, monospace">          idp4    = idw(chuva ~ 1, pr, grid, idp = 4)$var1.pred,</font></div>

<div><font face="courier new, monospace">          idp8    = idw(chuva ~ 1, pr, grid, idp = 8)$var1.pred)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">grid.data <- SpatialPixelsDataFrame(grid, idw.data)</font></div>

<div><font face="courier new, monospace">dput(names(grid.data))</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">spplot(grid.data, c("default", "idw6", "idw9", "idp1", "idp4", "idp8"), main = "IDW", col.regions=rainbow(16)</font></div>

<div><font face="courier new, monospace"># idp:  numeric; specify the inverse distance weighting power</font></div><div><font face="courier new, monospace"># nmax: the number of nearest observations that should be used </font></div>

<div><font face="courier new, monospace">### </code></font></div></div><div class="gmail_extra"><br></div><div class="gmail_extra"><div><img src="http://r-br.2285057.n4.nabble.com/attachment/4662232/0/Rplot_idw.png" style="margin-right:0px"><br>
</div>
<br></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr">Éder Comunello <<a href="http://user/SendEmail.jtp?type=node&node=4662232&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a><a href="http://user/SendEmail.jtp?type=node&node=4662232&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> <br>

Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div>
<br><br><div class="gmail_quote"><br></div></div></div><div class="">
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4662232&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.

        
        
        
        <br>
        <br>
        <hr noshade size="1" color="#cccccc">
        </div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div class="">
                <div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
                </div><a href="http://r-br.2285057.n4.nabble.com/R-br-IDW-tp4662226p4662232.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-IDW-tp4662226p4662232.html</a>
        </div><div class="HOEnZb"><div class="h5">
        <div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
                
                To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>

                <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
        </div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>