<div dir="ltr"><div><div>Olá André te aconselho além de mexer na opção max.dist, mexer também nas opções uvec e trend, esse último é interessante caso tenha observado uma tendência. Segue um exemplo.<br><br>v.pluviANO <-  variog(Pluvi.geo,max.dist=200,uvec=seq(0, 200, by=10),trend="1st").<br><br></div>Outra questão é usar a função likfit() que ajusta os parâmetros dos modelos diretamente aos dados, sem a necessidade prévia de fazer um variograma empírico. Assim, você utilizaria o AIC como critério de escolha do melhor modelo. Segue um exemplo.<br><br>Pluvilf<-list()<br>Pluvilf$exp <- likfit(Pluvi.geo, ini=c(2.80e-6,24.32), nug=3.5e-7,trend="1st")<br>Pluvilf$gau <- likfit(Pluvi.geo, cov.model = "gau", ini=c(2.80e-6,36.49), nug=3.5e-7,trend="1st")  <br>Pluvilf$sph <- likfit(Pluvi.geo, cov.model = "sph", ini=c(2.80e-6,72.97), nug=3.5e-7,trend="1st") <br>Pluvilf$cir <- likfit(Pluvi.geo, cov.model = "cir", ini=c(2.80e-6,72.97), nug=3.5e-7,trend="1st")<br>Pluvilf$kappa1.5 <- likfit(Pluvi.geo, cov.model = "mat", kappa= 1.5, ini=c(2.80e-6,12.16), nug=3.5e-7,trend="1st")  <br>Pluvilf$kappa2.5 <- likfit(Pluvi.geo, cov.model = "mat", kappa= 2.5, ini=c(2.80e-6,12.16), nug=3.5e-7,trend="1st") <br><br></div>## Escolha pelo AIC o menor é o melhor!!<br><div>sort(sapply(Pluvilf4, AIC))<br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 3 de abril de 2016 18:12, Andre Oliveira <span dir="ltr"><<a href="mailto:andreolsouza@yahoo.com.br" target="_blank">andreolsouza@yahoo.com.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:bookman old style,new york,times,serif;font-size:16px"><div dir="ltr">Pessoal aproveitando os tópico e sou novato no tema também.  Estou fazendo tentativas para melhor simivariograma empírico e já que o pH é uma va que não é normal Fiz o seguinte.</div><div dir="ltr"><br></div><div dir="ltr">require(geoR)<br></div><div dir="ltr">boxcox(dados$data~1)    #  Achei lambda = 0.7 para o pH <br></div><div dir="ltr"><br></div><div dir="ltr">vario.emp.1 <- variog(dados,max.dist=30 direction="omnidirectional",lam=0.7)</div><div dir="ltr">vario.emp.2 <- variog(dados,max.dist=70, direction="omnidirectional",lam=0.7)</div><div dir="ltr">vario.emp.3 <- variog(dados,max.dist=100, direction="omnidirectional",lam=0.7)</div><div dir="ltr">vario.emp.4 <- variog(dados,max.dist=150 , direction="omnidirectional",lam=0.7)<br></div><div dir="ltr"><br></div><div dir="ltr"><br></div><div dir="ltr">Não consegui em nenhum dos casos  um variograma que se estabilizasse,  fiz outras variações para  o parâmetro max.dist! O que eu deveria fazer para resolver meu problema além de suar lambda =0.7? Uma outra dúvida quando uso lambda = 0.7 a função krige.conv faz a interpolação em escala original?</div><div dir="ltr"><br></div><div dir="ltr"> <br></div><div dir="ltr"><br></div><div dir="ltr"><br> </div><div><span></span></div><div> </div><div>André Oliveira Souza. 
 Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo.  IFES</div> <div><br><br></div><div style="display:block"> <div style="font-family:bookman old style,new york,times,serif;font-size:16px"> <div style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"><div><div class="h5"> <div dir="ltr"><font face="Arial" size="2"> Em Domingo, 3 de Abril de 2016 10:26, Wagner Wolff <<a href="mailto:wwolff@usp.br" target="_blank">wwolff@usp.br</a>> escreveu:<br></font></div>  <br><br> </div></div><div><div><div class="h5"><div><div><div dir="ltr">Olá César </div>
<div dir="ltr">A questão aqui é de interpolação espacial, envolvendo questões de modelagem geoestatistica, krigagem e tudo mais. Realmente não especifique isto na pergunta é ficou confuso, mas mesmo assim muito obrigado.</div>
<div dir="ltr">Abraço</div>
<div>Em 2 de abr de 2016 8:11 PM, "Cesar Rabak" <<a rel="nofollow" shape="rect" href="mailto:cesar.rabak@gmail.com" target="_blank">cesar.rabak@gmail.com</a>> escreveu:<br clear="none"><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Wagner,<div><br clear="none"></div><div>Está um pouco confusa a pergunta colocada por você.</div><div><br clear="none"></div><div>Para interpolar uma variável você precisa de uma relação entre ela e uma outra (independente).</div><div><br clear="none"></div><div>A normalidade ou não dessas duas variáveis não tem a menor importância para interpolação, o que importa, isso sim, é que a regressão feita (para encontrar uma lei que as "correlacionem") tenha resíduos normais.</div><div><br clear="none"></div><div>Mas não me lembro de nenhum texto (e nem vejo sentido físico) em interpolar os resíduos, que por definição indicariam ruído na sua regressão (seria algo como tentar interpolar folhas de chá após apreciar a bebida!).</div><div><br clear="none"></div><div>HTH</div></div><div><br clear="none"><div>2016-04-02 12:07 GMT-03:00 Wagner Wolff <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:wwolff@usp.br" target="_blank">wwolff@usp.br</a>></span>:<br clear="none"><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Olá pessoal a dúvida é a seguinte. <br clear="none">Preciso interpolar uma variável na qual não consegui normalizá-la pelo BoxCox(); entretanto, a análise exploratória de tendência espacial indicou um modelo para remoção da mesma; assim, passo a interpolar os resíduos deste modelo, correto? <br clear="none">Sendo assim, posso desconsiderar a não normalidade da variável inicial e apenas considerar a normalidade dos resíduos do modelo Trend?  <br clear="none"></div>Outra questão é se realmente preciso normalizar a Variável, é plausível aceitar a normalidade dela por este gráfico abaixo? Acabei aceitando esta condição e não tive problemas na obtenção dos meus parâmetros.<br clear="none"><div><br clear="none"><br clear="none"><img src="cid:fEz9D4ZjRj4PcJ14rW31" height="224" width="442"><br clear="none"></div><div>​Att. <br clear="none"><span><font color="#888888">Wagner Wolff<br clear="none"></font></span></div></div>
<br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" 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 clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" 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.<div><br clear="none"></div></blockquote></div><div><br clear="none"></div></div><div>
<br clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" 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 clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" 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 clear="none"></div></blockquote></div></div></div><br></div></div><div><div><div class="h5">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" 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 clear="none"></div></div>Leia o guia de postagem (<a shape="rect" href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne� c�igo m�imo reproduz�el.</div><br><br></div>  </div> </div>  </div></div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" 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="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br></div>