MUITO obrigada pelas dicas Paulo! <br />Vou olhar melhor essa página do geoR, não sabia que existia esta pagina com scripts, etc<br /><br />Obrigada mesmo. <br />Abraço,<br /> <br />Ariane<br /><br /><br /><br />
<hr style="border-top: 1px solid #ccc;" />
Em 13/06/2011 08:57, <strong>Paulo Justiniano < paulojus@leg.ufpr.br ></strong> escreveu:<br />É importante voce fazer algums analises explortorios para verificar o <br />compotamento dos dados anted de variograms ou likfit.<br /><br />em qualquer caso abaixo a assimetria é nitida e gtransformação é <br />recomendada<br /><br /> plot(jura_dados.Cu, low=T)<br /> plot(jura_dados.Cu, trend=~Landuse, low=T)<br /> plot(jura_dados.Cu, trend=~Rock, low=T)<br /> plot(jura_dados.Cu, trend=~Rock+Landuse, low=T)<br /><br /> require(MASS)<br /> boxcox(jura_dados.Cu, trend=~Rock+Landuse)<br /><br /> plot(jura_dados.Cu, trend=~Rock+Landuse, low=T, lambda=0)<br /><br /> Siga as analises considerando esta (ou alguma outra TRANSFORMACAO)<br /><br /> likfit() retorna logLik, AIC, BIC para os odelos ajustados e portanto<br />é só estratir e comparar, em especial apra avaliar os efeitos das <br />potencias covariaveis (landuse e rock) por algum teste<br /><br />Lembre que na pagina 
 da geoR tem uma sessao de tutotias e varios scripts <br />com exemplo, alem do vignette do pacote<br /><br /><br /><br />On Sun, 12 Jun 2011, nani.design@bol.com.br wrote:<br /><br />> <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<br />> R na hora de fazer as devidas modelagens.<br />> Estou utilizando o conjunto de dados "Jura", incluido no pacote "gstat", para um pequeno estudo. O conjunto de dados<br />> contém as coordenadas espaciais X e Y (Xloc e Yloc), duas variáveis categóricas Landuse e Rock, e a leitura de 7<br />> metais em ppm. Onde:<br />> Landuse: 1 = floresta;<br />> 2 = pastagem;<br />> 3 = campo;<br />> 4 = lavoura<br />> <br />> Rock: 1 = Argovian;<br />> 2 = Kimmeridgian;<br />> 3 = Sequanian;<br />> 4 = Portlandian;<br />> 5 = Quaternary<br />> <br />> Utilizando o pacote "geoR", gerei o
 s variogramas de cada metal através da função "variog". Por exemplo, o variograma<br />> do metal Cobre através do R:<br />> <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 />> Gostaria agora de ajustar um modelo para este variograma através da função "likfit", mas não consigo. Não sei se há<br />> 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,<br />> testando o efeito da rocha "A
 rgovian" vamos supor, entre outros). Como faço isto na função "likfit"? Existe outras<br />> funções melhores neste caso para fazer os devidos ajustes?<br />> <br />> 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 />> <br />><br /><br />