<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
humilhou... :)<br>
<br>
<div class="moz-cite-prefix">On 29/05/14 12:36, Éder Comunello
wrote:<br>
</div>
<blockquote
cite="mid:CABmC8gmK3g5YZkEc5FNrTSGYuz101bu+O6gvRuLy8HE2eYYCdA@mail.gmail.com"
type="cite">
<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
moz-do-not-send="true" href="http://grid.pt">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
moz-do-not-send="true" href="http://grid.pt">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
moz-do-not-send="true" href="http://idw.pt">idw.pt</a>
<- idw(chuva~1, pr, <a moz-do-not-send="true"
href="http://grid.pt">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
moz-do-not-send="true" href="http://idw.pt">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="cid:part6.02060300.09060607@yahoo.com.br"
style="margin-right: 0px;"><br>
</div>
<br>
</div>
<div class="gmail_extra"><br clear="all">
<div>
<div dir="ltr">Éder Comunello <<a moz-do-not-send="true"
href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a
moz-do-not-send="true"
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 class="gmail_quote"><br>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
</body>
</html>