Prezados.

Por razões didáticas gostaria de decompor a variável regionalizada "v" no exemplo abaixo em duas partes:
(i) uma parte atribuída a uma tendência associada às suas coordenadas (x,y)
(ii) outra que seriam seus resíduos suavizados pela mediana (median polish)

library(geoR)
dados.mat=matrix(nrow=8, ncol=3,byrow=F,
c(486259,483200,487021,478100,482646,481735,480789,484181,
7428046,7428217,7424717,7420978,7426862,7427070,7427206,7427091,
159,350,220,265,450,110,224,530))
colnames(dados.mat)=c("x","y","v")
dados.mp=medpolish(dados.mat)
dados.res=as.vector(dados.mp$residuals[,3])
dados.mat.mp[,3]=dados.res
dados.mat.mp
geo.dados.mat.mp=as.geodata(dados.mat.mp,coords.col = 1:2,data.col=3)

Minhas dúvidas são:

(i) a "extração" dos resíduos está correta?
(ii) numa análise geoestatística, ao se modelar os resíduos (semivariograma experimental > ajuste dos parâmetros de uma função teórica) e usar esta função para
estimar os valores numa grade de pontos (krigagem ordinária) ainda deve ser "adicionada" a parte atribuída à tendência removida dos dados experimentais no começo?
(ou teoricamente não há razão para isso?)

 

Antecipadamente agradeço a colaboração
e desejo a todos um ótimo 2015.

Felinto