<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="cid:ii_hv70i1ad0_145fc2c91b43f43b" 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="cid:ii_hv70m0h71_145fc2f677d3739e" width="455" height="234.13541666666666"><br>O segundo Gaussiano por ML:</div><div><img src="cid:ii_hv70nzv42_145fc30cf0341f87" 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="mailto:agrohelder@hotmail.com" target="_blank">agrohelder@gmail.com</a></i></div></i></div><div style="text-align:center"><br></div></div></div>
</div>