Hélio,

Se nao ha' dependencia espacial nao faz sentido interpolar. O IDW pode ser resumido numa conta matricial. Exemplo:

n = 30;
coo = cbind(runif(n)*2, runif(n))
y = rnorm(n)

x0 = seq(0, 1, length=101)
y0 = seq(0, 2, length=201)
grid = as.matrix(expand.grid(x0, y0))

p = 5
w = 1/(outer(grid[,1], coo[,1], '-')^2 +
outer(grid[,2], coo[,2], '-')^2)^(p/2)
w = w/rowSums(w)

result = matrix(w%*%y, length(x0))

require(fields)
image.plot(list(x=x0, y=y0, z=result), asp=1)
points(coo, cex=0.1+(y-min(y))/diff(range(y)))

Elias.

On 28/05/14 21:21, Hélio Gallo Rocha wrote:

Caros Listeiros.

Estou utilizando o pacote geoR para análise geoestatística.

Numa das áreas em estudo não apresentou dependência espacial.

Sei que no pacote geosat contempla o método IDW.

Alguém teria um CRM com este método?

Agradeço desde já.

--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.