<div dir="ltr">Caros Hélder e Éder.<div><br></div><div>Vou dar minha opinião, mas como o Éder, é para tentar ajudar.</div><div><br></div><div>Como o Éder comentou no item (6), seis pontos parece ser uma quantidade pequena e como a distância minima é de 10 m, então teria de aumentar a distância máxima para uns 80%, teria de testar.</div>
<div><br></div><div>Hélder, vc. selecionou o modelo gaussiano por ML como o melhor, mas foi o que apresentou o mais alto valor de AIC. Se for seguir este indicativo, este seria o pior modelo.</div><div><br></div><div>Pelo que andei lendo, a seleção do modelo seria pela observação do mapa da variância da krigagem. O que apresenta a menor variância é o selecionado.</div>
<div><br></div><div>Observando a escala dos dois mapas da variância, o ajustado por OLS apresenta valores inferiores indicando menor variância, o que seria, melhor ajuste, mas como disse antes, o modelo gaussiano ajustado por ML, analisando pelo AIC, é o pior.</div>
<div><br></div><div>Abraço</div><div><br></div><div>Hélio</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 14 de maio de 2014 16:39, Hélder [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4662125h43@n4.nabble.com" target="_blank">ml-node+s2285057n4662125h43@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"><div>Gostaria de tirar algumas dúvidas, desde já peço desculpas se cometi algum deslize, tenho tentado evoluir na utilização do R me baseando apenas em tutoriais, e nas orientações que recebo aqui do grupo.<br>
</div><div><br></div><div>A validação cruzada no pacote Geo R nos fornece como resultado: o Erro Médio (EM), Erro Médio Reduzido (ER), Desvio padrão do erro médio (SEM), Desvio padrão do erro reduzido (SER) e obtenho por cálculo o Erro Absoluto (EA), também consigo o AIC e BIC nos modelos por Máxima Verossimilhança, quando do ajuste do modelo.<br>
</div><div><br></div><div>1) Todos estes parâmetros possuem igual nível de importância na análise, ou algum deles tem peso maior?</div><div>2) Estou analisando cada um dos parâmetros e o modelo que "vence" na maior parte deles admito que é o melhor modelo, esta forma de analisar está correta?</div>
<div>3) Analisando-se o melhor modelo (escolhido como na pergunta 2) com o Gráfico do Semivariograma Teórico ajustado, nem sempre este parece que é o melhor modelo ajustado, nestes casos, deve prevalecer a análise da validação cruzada ou a visual?</div>
<div>4) Comparando os resultados da krigagem do modelo escolhido pela Validação Cruzada com o obtido pela análise visual, a krigagem utilizando como modelo o semivariograma "melhor" visualmente parece representar adequadamente o variação do fenômeno na área, nestes casos o que vocês sugerem?</div>
<div>5) Em alguns casos a Validação cruzada retornou como resultado NAN e INF, trata-se de algum problema com os dados?</div><div><br></div><div>Caso concreto:</div><div><br></div><div>A tabela abaixo traz os resultados da validação cruzada da variável condutividade hidráulica Ks, por dois métodos (Mínimos Quadrados e Máxima Verossimilhança), pela interpretação dos resultados, o modelo Gaussiano por Máxima Verossimilhança seria o melhor modelo a meu ver.</div>
<div><table frame="VOID" cellspacing="0" cols="10" rules="NONE" border="0">
<colgroup><col width="86"><col width="86"><col width="86"><col width="257"><col width="86"><col width="86"><col width="86"><col width="86"><col width="86"><col width="86"></colgroup>
<tbody>
<tr>
<td width="86" height="18" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="257" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
<td width="86" align="LEFT"><br></td>
</tr>
<tr>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" height="18" align="CENTER"><b>Variável</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>Método</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>Modelo</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>EM</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>ER</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>Sem</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>Ser</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>EA</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>AIC</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0);border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>BIC</b></td>
</tr>
<tr>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" height="18" align="LEFT"><b><br></b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">OLS</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">Esférico</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">-1,27E-16</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">-6,42E-16</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">0,19729</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER"><b>1,01442</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">6,51815</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER"><br></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER"><br></td>
</tr>
<tr>
<td height="18" align="LEFT"><br></td>
<td align="CENTER">OLS</td>
<td align="CENTER">Exponencial</td>
<td align="CENTER"><b>-0,00132</b></td>
<td align="CENTER"><b>-0,00346</b></td>
<td align="CENTER"><b>0,17911</b></td>
<td align="CENTER">0,95090</td>
<td align="CENTER">5,99872</td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
</tr>
<tr>
<td height="18" align="LEFT"><br></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">OLS</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">Gaussiano</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">-0,00097</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">-0,00259</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">0,18444</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER">0,98982</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER"><b>5,87635</b></td>
<td align="CENTER"><br></td>
<td align="CENTER"><br></td>
</tr>
<tr>
<td height="18" align="CENTER"><b>Ks</b></td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">ML</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">Esférico</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">-5,18E-16</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">-2,69E-015</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">0,19729</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">1,02439</td>
<td style="border-top-width:1px;border-top-style:solid;border-top-color:rgb(0,0,0)" align="CENTER">6,51815</td>
<td align="CENTER">-12,18</td>
<td align="CENTER">-5,233</td>
</tr>
<tr>
<td height="18" align="LEFT"><b><br></b></td>
<td align="CENTER">ML</td>
<td align="CENTER">Exponencial</td>
<td align="CENTER">-0,00118</td>
<td align="CENTER">-0,00327</td>
<td align="CENTER">0,18002</td>
<td align="CENTER">1,01035</td>
<td align="CENTER">5,91153</td>
<td align="CENTER">-15,24</td>
<td align="CENTER">-8,29</td>
</tr>
<tr>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" height="18" align="LEFT"><br></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00">ML</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00">Gaussiano</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00"><b>-0,00134</b></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00"><b>-0,00372</b></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00">0,17934</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00"><b>1,00992</b></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00">5,89874</td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00"><b>-15,35</b></td>
<td style="border-bottom-width:1px;border-bottom-style:solid;border-bottom-color:rgb(0,0,0)" align="CENTER" bgcolor="#FFFF00"><b>-8,4</b></td>
</tr>
<tr>
<td height="18" align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
<td align="LEFT"><br></td>
</tr>
</tbody>
</table></div><div>Seguem abaixo os Semivariogramas Teóricos Ajustados pelos dois métodos, o da esquerda é por OLS e o da direita é por ML, o modelo Gaussiano está em verde, analisando visualmente achei que o Modelo Gaussiano por OLS está melhor ajustado:<br>
</div>
<div><img src="http://r-br.2285057.n4.nabble.com/attachment/4662125/0/comp_ols_ml_ks.png" width="455" height="234.13541666666666"></div><div>Nas figuras abaixo estão os resultados da krigagem comparando os dois modelos, o primeiro Gaussiano por OLS:<br>
</div>
<div><img src="http://r-br.2285057.n4.nabble.com/attachment/4662125/2/comp_krig_ks_ols_gaus.png" width="455" height="234.13541666666666"><br>O segundo Gaussiano por ML:</div><div><img src="http://r-br.2285057.n4.nabble.com/attachment/4662125/1/comp_krig_ks_ml_gaus.png" width="455" height="234.13541666666666"><br>
Desde já agradeço a toda ajuda,</div><div><br></div><div><div dir="ltr"><div style="text-align:center"><font size="1"><i><b>Hélder Gramacho </b></i></font></div><div style="text-align:center"><font size="1">Recife-PE / </font><i style="font-size:x-small;color:rgb(51,51,255)"><div style="display:inline!important">
<i><a href="http://user/SendEmail.jtp?type=node&node=4662125&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a></i></div></i></div><div style="text-align:center"><br></div></div></div>
</div>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4662125&i=1" 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 style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-Validacao-Cruzada-e-Krigagem-tp4662125.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Validacao-Cruzada-e-Krigagem-tp4662125.html</a>
</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>
</div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>