<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Corrigindo os limites x-y:<br>
<br>
n = 30;<br>
coo = cbind(runif(n)*2, runif(n))<br>
y = rnorm(n)<br>
<br>
x0 = seq(0, 2, length=201)<br>
y0 = seq(0, 1, length=101)<br>
grid = as.matrix(expand.grid(x0, y0))<br>
<br>
p = 5<br>
w = 1/(outer(grid[,1], coo[,1], '-')^2 + <br>
outer(grid[,2], coo[,2], '-')^2)^(p/2)<br>
w = w/rowSums(w)<br>
<br>
result = matrix(w%*%y, length(x0))<br>
<br>
require(fields)<br>
image.plot(list(x=x0, y=y0, z=result), asp=1)<br>
points(coo, cex=0.1+(y-min(y))/diff(range(y)))<br>
<br>
<div class="moz-cite-prefix">On 29/05/14 08:21, Elias T. Krainski
wrote:<br>
</div>
<blockquote cite="mid:5386D1F8.7070609@yahoo.com.br" type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Hélio, <br>
<br>
Se nao ha' dependencia espacial nao faz sentido interpolar. O IDW
pode ser resumido numa conta matricial. Exemplo: <br>
<br>
n = 30;<br>
coo = cbind(runif(n)*2, runif(n))<br>
y = rnorm(n)<br>
<br>
x0 = seq(0, 1, length=101)<br>
y0 = seq(0, 2, length=201)<br>
grid = as.matrix(expand.grid(x0, y0))<br>
<br>
p = 5<br>
w = 1/(outer(grid[,1], coo[,1], '-')^2 + <br>
outer(grid[,2], coo[,2], '-')^2)^(p/2)<br>
w = w/rowSums(w)<br>
<br>
result = matrix(w%*%y, length(x0))<br>
<br>
require(fields)<br>
image.plot(list(x=x0, y=y0, z=result), asp=1)<br>
points(coo, cex=0.1+(y-min(y))/diff(range(y)))<br>
<br>
Elias.<br>
<br>
<div class="moz-cite-prefix">On 28/05/14 21:21, Hélio Gallo Rocha
wrote:<br>
</div>
<blockquote
cite="mid:CA+ppPW59H8BTHBNWJfMeo_twbMqiToqn2tfauh0_-Ksumenqtw@mail.gmail.com"
type="cite">
<div dir="ltr">Caros Listeiros.
<div><br>
</div>
<div>Estou utilizando o pacote geoR para análise
geoestatística. </div>
<div>Numa das áreas em estudo não apresentou dependência
espacial.</div>
<div>Sei que no pacote geosat contempla o método IDW.</div>
<div>Alguém teria um CRM com este método?</div>
<div><br>
</div>
<div>Agradeço desde já.</div>
<div> <br clear="all">
<div><br>
</div>
-- <br>
Hélio Gallo Rocha<br>
IFSULDEMINAS - Câmpus Muzambinho<br>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true" 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>
<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>