<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=windows-1252">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="-1">Prezados.<br>
<br>
Por razões didáticas gostaria de decompor a variável regionalizada
"v" no exemplo abaixo em duas partes:<br>
(i) uma parte atribuída a uma tendência associada às suas
coordenadas (x,y)<br>
(ii) outra que seriam seus resíduos suavizados pela mediana
(median polish)<br>
<br>
</font><font size="-1">library(geoR)<br>
dados.mat=matrix(nrow=8, ncol=3,byrow=F,<br>
c(486259,483200,487021,478100,482646,481735,480789,484181,<br>
7428046,7428217,7424717,7420978,7426862,7427070,7427206,7427091,<br>
159,350,220,265,450,110,224,530))<br>
colnames(dados.mat)=c("x","y","v")<br>
dados.mp=medpolish(dados.mat)<br>
dados.res=as.vector(dados.mp$residuals[,3])<br>
dados.mat.mp[,3]=dados.res<br>
dados.mat.mp<br>
geo.dados.mat.mp=as.geodata(dados.mat.mp,coords.col =
1:2,data.col=3)<br>
<br>
Minhas dúvidas são:<br>
<br>
(i) a "extração" dos resíduos está correta?<br>
(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<br>
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?<br>
(ou teoricamente não há razão para isso?)<br>
<br>
<br>
<br>
Antecipadamente agradeço a colaboração<br>
e desejo a todos um ótimo 2015.<br>
<br>
Felinto <br>
<br>
<br>
<br>
<br>
</font>
<pre class="moz-signature" cols="72"> </pre>
</body>
</html>