<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Éder, essa e' a ultima versao... <br>
<br>
<div class="moz-cite-prefix">On 29/05/14 17:16, Éder Comunello
wrote:<br>
</div>
<blockquote
cite="mid:CABmC8gnR7Z_qU1Ccv_xbY2TtWDb14WRtjwTULVqkbGZfbf=V5g@mail.gmail.com"
type="cite">
<div dir="ltr">Elias, bom dia!
<div><br>
</div>
<div>Reconheço o mérito das suas colocações, só complementei! </div>
<div><br>
</div>
<div>Aproveitando a deixa... A exemplo do material do Jackson
Aquino, mantenho também para consulta um material seu,
"Introdução à análise de dados espacialmente referenciados".
Tenho aqui uma versão de 2012. Tens alguma edição mais atual
ou intenção de transformar esse material em livro?</div>
<div><br>
</div>
<div>Grato,</div>
<div><br>
<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">Em 29 de maio de 2014 07:18, Elias
T. Krainski <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:eliaskrainski@yahoo.com.br"
target="_blank">eliaskrainski@yahoo.com.br</a>></span>
escreveu:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px
0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> humilhou... :)
<div>
<div class="h5"><br>
<br>
<div>On 29/05/14 12:36, Éder Comunello wrote:<br>
</div>
</div>
</div>
<blockquote type="cite">
<div>
<div class="h5">
<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" 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
moz-do-not-send="true"
href="http://grid.pt" 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
moz-do-not-send="true"
href="http://idw.pt" target="_blank">idw.pt</a>
<- idw(chuva~1, pr, <a
moz-do-not-send="true"
href="http://grid.pt" 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
moz-do-not-send="true"
href="http://idw.pt" 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 moz-do-not-send="true"
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></fieldset>
<br>
</div>
</div>
<div class="">
<pre>_______________________________________________
R-br mailing list
<a moz-do-not-send="true" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a moz-do-not-send="true" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</div>
</blockquote>
<br>
</div>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br"
target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a moz-do-not-send="true"
href="http://www.leg.ufpr.br/r-br-guia"
target="_blank">http://www.leg.ufpr.br/r-br-guia</a>)
e forneça código mínimo reproduzível.<br>
</blockquote>
</div>
<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>