<div dir="ltr">Caro Eder,<div><br></div><div>Já havia feito a retirada de tendencia com o atributo altura da planta, que está ligado com a produção.</div><div><br></div><div>O variograma gerado ficou diferente se usar polinomial de 1st, mas não sei como validar isso depois...</div>
<div><br></div><div>Abraço</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 31 de outubro de 2013 13:15, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Boa tarde, Helio!<div><br></div><div>Vi seu email, mas vou precisar reservar um tempo pra tratar dele. Você chegou a usar as glebas como covariável na modelagem? Para isso cada ponto deveria ter uma var a mais identificando a gleba origem.</div>
<div><br></div><div>Atte.,</div><div><br></div><div> <br><div class="gmail_extra"><br clear="all"><div><div dir="ltr">Éder Comunello <<a href="mailto:comunello.eder@gmail.com" target="_blank">c</a><a 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 outubro de 2013 22:25, Hélio Gallo Rocha <span dir="ltr"><<a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Caro Eder,<div><br></div><div>Só hoje vi seu post com CRM anexo.<br><div><br></div><div>Muito obrigado pela revisão do CRM.</div><div><br></div><div>Como também está a procura de uma maior entendimento da geoestatística, vamos seguindo.<br>
</div><div><br></div><div><br></div><div>Vou numerar para separar as observações:</div><div><br></div><div>1. Pra entender melhor a área, realmente é café, o recorte esta meio estranho, mas na realidade são três talhões de 0.5 ha cada.</div>
<div>Cada um tem uma característica bem diferenciada com relação a produção:alta, baixa e média.</div><div><br></div><div>Segundo li, numa situação em que a média é flutuante, como parece ser o caso, teríamos de usar krigagem ordinária.</div>
<div><br></div><div><br></div><div><div>A maioria das observações tinha notado, a linha:</div><div> <span style="font-size:13px;font-family:'courier new',monospace">plot(dados, lowess=T, trend='1st')</span></div>
<div><span style="font-family:'courier new',monospace;font-size:13px">achei interessante, mostrando o resultado da retirada de tendencia.</span><br></div><div><br></div><div>2. no final :</div><div><br></div><div>
<div style="font-family:arial,sans-serif;font-size:13px"><font face="courier new, monospace">### Lembra do resíduo que separei anteriormente?</font></div><div style="font-family:arial,sans-serif;font-size:13px"><font face="courier new, monospace">### Se fosse usar...</font></div>
<div style="font-family:arial,sans-serif;font-size:13px"><font face="courier new, monospace">va.r = variog(dados.res, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical")</font></div><div style="font-family:arial,sans-serif;font-size:13px">
<font face="courier new, monospace">va.r.env <- variog.mc.env(dados.res, obj.var = va.r)</font></div><div style="font-family:arial,sans-serif;font-size:13px"><font face="courier new, monospace">plot(va.r, xlab="Dist. (m)", ylim=c(0,summary(va.r$v)[6]*1.1), env=va.r.env, main="Residual")</font></div>
</div><div style="font-family:arial,sans-serif;font-size:13px"><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace"><font face="arial, sans-serif">Foi possível gerar o semivariograma sem tendencia, com os </font>resíduos<font face="arial, sans-serif">, fazer o ajuste assentimento usando: </font></font><span style="font-family:arial,sans-serif">eyefit(va.r)</span></div>
<div><font face="arial, sans-serif"><br></font></div><div><font face="courier new, monospace"><font face="arial, sans-serif">Dai fazer a validação, fazendo a validação ou não e gerar o mapa sem tendencia.</font></font></div>
<div><font face="courier new, monospace"><font face="arial, sans-serif"><br></font></font></div><div><font face="courier new, monospace"><font face="arial, sans-serif">3. Quanto ao envelopamento, não sei se estou certo, mas ele simula valores máximos e minimos para dependencia, mas como é uma simulação, a cada vez que vc. roda, pode ter pontos com dependencia, nem que seje um, e as vezes não dá nenhum, mas de qualquer forma indica que o valor do alcance é pequeno em relação ao todo e baixa dependencia espacial, como nos meus dados.</font></font></div>
<div><font face="courier new, monospace"><font face="arial, sans-serif"><br></font></font></div><div><font face="courier new, monospace"><font face="arial, sans-serif">Esse item do envelopamento </font></font><span style="font-family:arial,sans-serif">vale</span><span style="font-family:arial,sans-serif"> </span><span style="font-family:arial,sans-serif">discussão, principalmente em situações que a variancia sobe além do limite superior e desce, isso significa algo importante, mas o que?</span><br>
</div><div><span style="font-family:arial,sans-serif"><br></span></div><div><span style="font-family:arial,sans-serif">4. Quanto aos links do Landim já tinha visto, usando o surfer, também vale atenção.</span></div><div>
<span style="font-family:arial,sans-serif">Fiz o recomendado por ele, mas dobrou os valores do mapa., vc. tentou?</span><br>
</div></div></div><div><span style="font-family:arial,sans-serif"><br></span></div><div><span style="font-family:arial,sans-serif">5. Quanto ao ajuste por verosimilhança colocada pelo Prof. Paulo:</span></div><div><span style="font-family:arial,sans-serif">Na primeira vez que havia rodado os dados achei muito estranho, pois o semivariograma subia indefinidamente e o semivariograma teórico por verossimilhança, </span><span style="font-family:arial,sans-serif">usando likfit, baixou muito o patamar.</span></div>
<div><span style="font-family:arial,sans-serif">Interessante se usar likfit com trend="1st", o semivariograma teórico baixa mais ainda, das dai não consigo krigar</span></div><div><span style="font-family:arial,sans-serif"><br>
</span></div><div><span style="font-family:arial,sans-serif">6. Não sei se posso fazer o ajuste assentimento usando o semivariograma com tendencia mas usar os parametros sem tendencia.</span></div><div><span style="font-family:arial,sans-serif"><br>
</span></div><div><span style="font-family:arial,sans-serif">7. Tenho de retirar a tendencia sempre?</span></div><div><span style="font-family:arial,sans-serif">Tem casos que há uma patamar claro, mas se rodar o variog com trend, fica efeito pepita puro...e ai?</span></div>
<div><br></div><div>8. Outra coisa que não sei se estou entendendo:</div><div>após a retirada da tendencia, fazer o mapa com o residuo, mas os valores não tem nada a ver com os reais, onde estou entendendo errado?</div><div>
<br>
</div><div>Abraço,</div><div><br></div><div>Hélio</div><div><span style="font-family:arial,sans-serif"><br></span></div><div><span style="font-family:arial,sans-serif"><br></span></div><div><span style="font-family:arial,sans-serif"><br>
</span></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 26 de outubro de 2013 11:05, Eder Comunello [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4660724h12@n4.nabble.com" target="_blank">ml-node+s2285057n4660724h12@n4.nabble.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Senhores, bom dia!<div><br></div><div>Agradeço novamente ao Prof. Paulo pela atenção dispensada.</div><div><br></div><div>Hélio, quanto ao script que eu tinha ficado de olhar, vou colocar minhas considerações abaixo.</div>
<div><br></div><div>Procurei não mudar muito, mas acabei alterando um pouco. Os comentários estão no corpo do script.</div><div><br></div><div>Lembrando sempre que também estou aprendendo sobre o assunto, então pode ter muita coisa errada ou mal interpretada aí. Mas acho que pode te ajudar.</div>
<div><br></div><div>Outra coisa é que não tenho muito ideia do que os dados são, como foram obtidos e as condições da área de onde vieram, o que pra mim é fundamental pra análise. Pelos indícios nos dados (localização, estrutura), vou arriscar que são dados da produção de café ou citros e que o terreno está numa pendente da direita pra esquerda, com maior produção na parte baixa. O que chamou a atenção é o formato da área que é bem particular e não parece muito 'natural'.</div>
<div><br></div><div>Mas como o foco principal é o uso do R/geoR, vou tentar me deter nisso.</div><div><br></div><div>Use por sua própria conta e risco! :D</div><div><br></div><div><br></div><div><div><font face="courier new, monospace">### <BEGIN></font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">require(geoR) #;require(tcltk2)</font></div><div><font face="courier new, monospace">setwd("C:/LAB/RGIS/trend"); dir()</font></div>
<div><font face="courier new, monospace">par.ori <- par(no.readonly = TRUE)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">url1 <- '<a href="https://dl.dropboxusercontent.com/u/117618178/helio/dados.txt" rel="nofollow" link="external" target="_blank">https://dl.dropboxusercontent.com/u/117618178/helio/dados.txt</a>'</font></div>
<div><font face="courier new, monospace">url2 <- '<a href="https://dl.dropboxusercontent.com/u/117618178/helio/recorte.txt" rel="nofollow" link="external" target="_blank">https://dl.dropboxusercontent.com/u/117618178/helio/recorte.txt</a>'</font></div>
<div>
<font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Como criar uma pasta pública no Dropbox...</font></div><div><font face="courier new, monospace">browseURL('<a href="https://www.dropbox.com/help/16/pt_BR" rel="nofollow" link="external" target="_blank">https://www.dropbox.com/help/16/pt_BR</a>')</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### 1. Carregando dados ------------------------------------------------------</font></div><div><font face="courier new, monospace">dados <- read.geodata(url1, coords.col=1:2, data.col=3, sep="", header=T)</font></div>
<div><font face="courier new, monospace">dados$borders <- read.table(url2, header=T)</font></div><div><font face="courier new, monospace">summary(dados); length(dados$data)</font></div><div><font face="courier new, monospace">ls(dados);ls.str(dados)</font></div>
<div><font face="courier new, monospace"> </font></div><div><font face="courier new, monospace">### 2. Visualizando dados ----------------------------------------------------</font></div><div><font face="courier new, monospace"><br>
</font></div><div><font face="courier new, monospace">### Visão geral</font></div><div><font face="courier new, monospace">par(mfrow=c(1,1)); par()$mfrow</font></div><div><font face="courier new, monospace">points(dados)</font></div>
<div><font face="courier new, monospace">points(dados, pt.div="quart")</font></div><div><font face="courier new, monospace">points(dados, pt.div="quint")</font></div><div><font face="courier new, monospace">### tem pelos três sub-áreas bem marcadas...</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Verificando tendências com o plot.geodata --------------------------------</font></div><div><font face="courier new, monospace">### Antes de partir pros variogramas, use o plot.geodata para pré-avaliar tendências</font></div>
<div><font face="courier new, monospace">### Bom para esclarecer opções de trend, pois notei certa confusão (cte,1st,2nd)</font></div><div><font face="courier new, monospace">### Recomendo ler o item 'trend' na ajuda da função plot.geodata</font></div>
<div><font face="courier new, monospace">?plot.geodata </font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(dados, lowess=T)</font></div><div><font face="courier new, monospace">plot(dados, lowess=T, trend='cte') ### mesmo que anterior</font></div>
<div><font face="courier new, monospace">### linhas de comando anteriores são equivalentes: trend='cte' é 'default' da geoR!</font></div><div><font face="courier new, monospace">### nos gráf. 2 e 3, plota-se as coord vs. dados. Observar as distorções na linha 'lowess'</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(dados, lowess=T, trend=~coords)</font></div><div><font face="courier new, monospace">plot(dados, lowess=T, trend='1st') </font></div>
<div><font face="courier new, monospace">plot(dados, lowess=T, trend=~dados$coords[,1]+dados$coords[,2])</font></div><div><font face="courier new, monospace">### as três linhas anteriores são equivalentes</font></div><div>
<font face="courier new, monospace">### o uso da polinomial de primeira ordem melhora o comportamento do histograma</font></div><div><font face="courier new, monospace">### mas a linha de lowess continua com comportamento anômalo</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(dados, lowess=T, trend='2nd') </font></div><div><font face="courier new, monospace">plot(dados, lowess=T, trend=~I(coords)+I(coords^2)+I(coords[,1]*coords[,2]))</font></div>
<div><font face="courier new, monospace">### as duas linhas anteriores são equivalentes (polinomal de segunda ordem)!</font></div><div><font face="courier new, monospace">### não aparenta agregar nenhuma melhoria</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Mais sobre modelos de tendência com trend.spatial ------------------------</font></div><div><font face="courier new, monospace">### Pra melhorar o entendimento dos modelos, dá pra dar uma olhada nas matrizes</font></div>
<div><font face="courier new, monospace">?trend.spatial</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial('cte', dados))) ### constante</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~coords, dados))) ### polinomal 1a ordem, usando coords</font></div>
<div><font face="courier new, monospace">head(unclass(trend.spatial('1st', dados))) ### mesmo que anterior</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial('2nd', dados))) ### polinomal 2a ordem, usando coords</font></div>
<div><font face="courier new, monospace">head(unclass(trend.spatial(~coords[,1], dados))) ### em função de X</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~coords[,2], dados))) ### em função de Y</font></div>
<div><font face="courier new, monospace">head(unclass(trend.spatial(~coords+I(coords[,1]^2), dados))) ### coords + X**2</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~I(coords)+I(coords^2), dados)))</font></div>
<div><font face="courier new, monospace">head(unclass(trend.spatial(~I(coords)+I(coords^2)+I(coords[,1]*coords[,2]), dados))) ### mesmo que 2nd</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~I(coords)+I(coords[,1]^2), dados)))</font></div>
<div><font face="courier new, monospace">head(unclass(trend.spatial(~I(coords^2), dados)))</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~coords[,1]*coords[,2], dados)))</font></div><div>
<font face="courier new, monospace">head(unclass(trend.spatial(~I(coords[,1]*coords[,2]), dados)))</font></div><div><font face="courier new, monospace">head(unclass(trend.spatial(~coords[,1]+I(coords[,1]^2)+coords[,2], dados)))</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Separando o drift com base na matriz trend.spatial -----------------------</font></div><div><font face="courier new, monospace">trend.matrix <- unclass(trend.spatial(trend='1st', geodata=dados)); nrow(trend.matrix)</font></div>
<div><font face="courier new, monospace">trend.lm <- lm(dados$data ~ trend.matrix)</font></div><div><font face="courier new, monospace">trend.res <- lm(dados$data ~ trend.matrix)$residuals</font></div><div><font face="courier new, monospace">trend.fit <- lm(dados$data ~ trend.matrix)$fitted.values</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### 'Empacotando' os resíduos em um geodata-----------------------------------</font></div><div><font face="courier new, monospace">dados.res <- dados</font></div>
<div><font face="courier new, monospace">names(trend.res) <- NULL; trend.res</font></div><div><font face="courier new, monospace">dados.res$data <- trend.res; dados.res$data</font></div><div><font face="courier new, monospace">str(dados.res)</font></div>
<div><font face="courier new, monospace">### Pode usar dados.res pra fazer a análise variográfica no lugar de dados, </font></div><div><font face="courier new, monospace">### caso tenha problemas com krige.conv()!!!</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Embora o uso da polinomial de 1a ordem tenha representado alguma melhoria,</font></div><div><font face="courier new, monospace">### eu não fiquei muito entusiamado por conta do 'lowess' apresentado</font></div>
<div><font face="courier new, monospace">### Talvez fosse o caso de testar covariáveis. </font></div><div><font face="courier new, monospace">### Uma fácil de obter e que pode ter relação é a altitude...</font></div><div>
<font face="courier new, monospace">### Mas vamos investigar o efeito nos variogramas</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### 3. Distancia máxima e parâmetros para semivariograma ---------------------</font></div>
<div><font face="courier new, monospace">hmax <- summary(dados)[[3]][[2]]*.8 ; hmax ### 80% da máxima distância</font></div><div><font face="courier new, monospace">uvec <- 12; pairs <- 20 ### pairs.min</font></div>
<div><font face="courier new, monospace"> </font></div><div><font face="courier new, monospace">### 4. Variogramas -----------------------------------------------------------</font></div><div><font face="courier new, monospace">### Não precisaria rodar todos, afinal como visto, tá repetindo modelos!!!</font></div>
<div><font face="courier new, monospace">### v0:média constante; v1:1st; v2:2nd; v3:cte; v4:coords</font></div><div><font face="courier new, monospace">v0 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical")</font></div>
<div><font face="courier new, monospace">v1 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical", trend="1st")</font></div><div><font face="courier new, monospace">v2 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical", trend="2nd")</font></div>
<div><font face="courier new, monospace">v3 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical", trend="cte")</font></div><div><font face="courier new, monospace">v4 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical", trend=~coords)</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### envelopes</font></div><div><font face="courier new, monospace">v0.env <- variog.mc.env(dados, obj.var = v0)</font></div>
<div><font face="courier new, monospace">v1.env <- variog.mc.env(dados, obj.var = v1)</font></div><div><font face="courier new, monospace">v2.env <- variog.mc.env(dados, obj.var = v2)</font></div><div><font face="courier new, monospace">v3.env <- variog.mc.env(dados, obj.var = v3)</font></div>
<div><font face="courier new, monospace">v4.env <- variog.mc.env(dados, obj.var = v4)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### 5. Visualizando variogramas ----------------------------------------------</font></div>
<div><font face="courier new, monospace">par(par.ori)</font></div><div><font face="courier new, monospace">par(mfrow=c(2,3))</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">points(dados, pt.div="quartiles",xlab="Leste", ylab="Norte")</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(v0, xlab="Dist. (m)", ylim=c(0,summary(v0$v)[6]*1.1), env=v0.env, main="Isotrópico com tendência")</font></div>
<div><font face="courier new, monospace">text(v0$u,v0$v,round(v0$n,1),pos=1);text(v0$u,v0$v,round(v0$u,1),pos=3);lines(v0$u,v0$v);abline(h=mean(v0$v))</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">plot(v1, xlab="Dist. (m)", ylim=c(0,summary(v1$v)[6]*1.1), env=v1.env, main="Isotrópico trend 1st")</font></div><div><font face="courier new, monospace">text(v1$u,v1$v,round(v1$n,1),pos=1);text(v1$u,v1$v,round(v1$u,1),pos=3);lines(v1$u,v1$v);abline(h=mean(v1$v))</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(v2, xlab="Dist. (m)", ylim=c(0,summary(v2$v)[6]*1.1), env=v2.env, main="Isotrópico trend 2nd")</font></div>
<div><font face="courier new, monospace">text(v2$u,v2$v,round(v2$n,1),pos=1);text(v2$u,v2$v,round(v2$u,1),pos=3);lines(v2$u,v2$v);abline(h=mean(v2$v))</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">plot(v3, xlab="Dist. (m)", ylim=c(0,summary(v3$v)[6]*1.1), env=v3.env, main="Isotrópico trend cte")</font></div><div><font face="courier new, monospace">text(v3$u,v3$v,round(v3$n,1),pos=1);text(v3$u,v3$v,round(v3$u,1),pos=3);lines(v3$u,v3$v);abline(h=mean(v3$v))</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">plot(v4, xlab="Dist. (m)", ylim=c(0,summary(v4$v)[6]*1.1), env=v4.env, main="Isotrópico trend coords")</font></div>
<div><font face="courier new, monospace">text(v4$u,v4$v,round(v4$n,1),pos=1);text(v4$u,v4$v,round(v4$u,1),pos=3);lines(v4$u,v4$v);abline(h=mean(v4$v))</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">### Analisando os variogramas ------------------------------------------------</font></div><div><font face="courier new, monospace">### Realmente o comportamento do histograma de trend='cte' parece conter uma tendência,</font></div>
<div><font face="courier new, monospace">### mas os modelos de 'retirada' usados deixam os dados sem uma estrutura de </font></div><div><font face="courier new, monospace">### variação. Dificilmente vai conseguir ajustar um bom modelo nesses casos, pois </font></div>
<div><font face="courier new, monospace">### o comportamento é muito errático. É possível que O erro que você estava obtendo </font></div><div><font face="courier new, monospace">### seja decorrente da fragilidade do modelo ajustado.</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Realmente não ficou bom! Mas e aí? O que dá pra fazer? Tenta ssim mesmo?</font></div><div><font face="courier new, monospace">### Ainda tem que testar a anisotropia, que não tem no seu script original.</font></div>
<div><font face="courier new, monospace">?variog4</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">va <- variog4(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical")</font></div>
<div><font face="courier new, monospace">par(par.ori)</font></div><div><font face="courier new, monospace">plot(va)</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Pronto! Apareceu anisotropia...</font></div>
<div><font face="courier new, monospace">### Tem duas sugestões de leitura abaixo pra dar uma ajuda...</font></div><div><font face="courier new, monospace">browseURL('<a href="http://www.rc.unesp.br/igce/aplicada/DIDATICOS/LANDIM/tkrigagem.pdf" rel="nofollow" link="external" target="_blank">http://www.rc.unesp.br/igce/aplicada/DIDATICOS/LANDIM/tkrigagem.pdf</a>')</font></div>
<div><font face="courier new, monospace">browseURL('<a href="https://www.soils.org/publications/sssaj/abstracts/61/1/SS0610010298" rel="nofollow" link="external" target="_blank">https://www.soils.org/publications/sssaj/abstracts/61/1/SS0610010298</a>')</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Você vai ter que estudar adequadamente a questão da anisotropia em seus dados</font></div><div><font face="courier new, monospace">### Isso não te isenta do efeito combinado anisotropia + tendência</font></div>
<div><font face="courier new, monospace">### Só pra testar as possibilidades, vamos fazer o variograma na direção 135 graus</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Vai ter que entrar a direções, sendo mais fácil usar PI/n</font></div>
<div><font face="courier new, monospace">### Tabela das direções (PI/n) para variogramas</font></div><div><font face="courier new, monospace">denom <- c(1:4,6,8,10,12,24,36,45,60); k <- (pi/180) ### Denominadores</font></div>
<div><font face="courier new, monospace">data.frame(rad=paste0('pi/', denom), rad2=round(pi/denom,4), deg=pi/denom/k)</font></div><div><font face="courier new, monospace">### 135 é 3*45 ou 3*pi/4</font></div><div>
<font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">va.135 <- variog(dados, uvec=uvec, pairs.min=pairs, dir=3*pi/4, max.dist=hmax, estimator="classical")</font></div>
<div><font face="courier new, monospace">plot(va.135)</font></div><div><font face="courier new, monospace">### Aqui parece ser mais fácil ajustar. Talvez tenha que ajustar o 'lag'</font></div><div><font face="courier new, monospace">### Avaliar a questão da distância entre amostras...</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### Lembra do resíduo que separei anteriormente?</font></div><div><font face="courier new, monospace">### Se fosse usar...</font></div>
<div><font face="courier new, monospace">va.r = variog(dados.res, uvec=uvec, pairs.min=pairs, max.dist=hmax, estimator="classical")</font></div><div><font face="courier new, monospace">va.r.env <- variog.mc.env(dados.res, obj.var = va.r)</font></div>
<div><font face="courier new, monospace">plot(va.r, xlab="Dist. (m)", ylim=c(0,summary(va.r$v)[6]*1.1), env=va.r.env, main="Residual")</font></div><div><font face="courier new, monospace"><br></font></div>
<div><font face="courier new, monospace">### Vou parar por aqui, porque antes de prosseguir vai precisar ver </font></div><div><font face="courier new, monospace">### adequadamente essa questão da anisotropia.</font></div>
<div><font face="courier new, monospace">### No caso do efeito ser só da anisotropia, vai simplificar a krigagem e você</font></div><div><font face="courier new, monospace">### não deve mais ter o problema que vinha tendo.</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">### <END></font></div></div><div><br></div><div class="gmail_extra"><div><div dir="ltr">Éder Comunello <<a href="http://user/SendEmail.jtp?type=node&node=4660724&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a><a href="http://user/SendEmail.jtp?type=node&node=4660724&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> <br>
<div>
Dourados, MS - [22 16.5'S, 54 49'W]<br></div></div></div>
<br><br><div class="gmail_quote"><br></div></div></div><div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4660724&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</div><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><div>
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</div><a href="http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pacote-geoR-tp4660619p4660724.html" target="_blank">http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pacote-geoR-tp4660619p4660724.html</a>
</div><div><div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a><span><font color="#888888">
</font></span></div></div></div></blockquote></div><span><font color="#888888"><br><br clear="all"><span class="HOEnZb"><font color="#888888"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</font></span></font></span></div><span class="HOEnZb"><font color="#888888">
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a 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 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></font></span></blockquote></div><br></div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>