<br />Boa noite,<br /><br /> Estive tentando estudar um pouco de Geoestatística, conhecer sobre o assunto e tenho algumas dúvidas com relação ao R na hora de fazer as devidas modelagens.<br /> Estou utilizando o conjunto de dados "Jura", incluido no pacote <em>"gstat"</em>, para um pequeno estudo. O conjunto de dados contém as coordenadas espaciais X e Y (Xloc e Yloc), duas variáveis categóricas Landuse e Rock, e a leitura de 7 metais em ppm. Onde:<br /> <em> Landuse:</em> 1 = floresta; <br /> 2 = pastagem;<br /> 3 = campo; <br /> 4 = lavoura<br /><br /> <em> Rock: </em>1 = Argovian;<br /> 2 = Kimmeridgian;<br /> 3 = Sequanian;<br /> 4 = Portlandian; <br /> 5 = Quaternary<br /><br /> Utilizando o pacote<em> "geoR"</em>, gerei os variogramas de cada metal através da função "variog". Por exemplo, o variograma do metal Cobre através do R:<br /><span style="color: #333399;"><br />library(gstat)<br />library(geoR)<br /><br />data(jura)<br />prediction.dat    # conjun
 to de dados Jura para o estudo<br /><br />jura_dados.Cu <- as.geodata(prediction.dat, data.col = 8, covar.col = 3:4)<br />varioJura.Cu <- variog(jura_dados.Cu, max.dist=4.5)   # variograma do Cobre<br />plot(varioJura.Cu, main="Variograma - Cobre", type="b")<br /><br /><span style="color: #000000;"> Gostaria agora de ajustar um modelo para este variograma através da função "likfit", mas não consigo. Não sei se há a necessidade de fazer algumas transformações antes, ou definir de alguma forma o modelo a ser ajustado, etc.<br /><br /> Além disto, gostaria de fazer testes através da razão de máxima verossimilhança de modelos reduzidos (por exemplo, testando o efeito da rocha "Argovian" vamos supor, entre outros). Como faço isto na função "likfit"? Existe outras funções melhores neste caso para fazer os devidos ajustes?<br /> <br /></span></span> Não sei se consegui ser bem clara na minha dúvida, mas se alguém souber me ajudar ficarei imensamente gra
 ta.<br /> Obrigada,<br /><br /><br />Ariane<br /><br />