É importante voce fazer algums analises explortorios para verificar o compotamento dos dados anted de variograms ou likfit. em qualquer caso abaixo a assimetria é nitida e gtransformação é recomendada plot(jura_dados.Cu, low=T) plot(jura_dados.Cu, trend=~Landuse, low=T) plot(jura_dados.Cu, trend=~Rock, low=T) plot(jura_dados.Cu, trend=~Rock+Landuse, low=T) require(MASS) boxcox(jura_dados.Cu, trend=~Rock+Landuse) plot(jura_dados.Cu, trend=~Rock+Landuse, low=T, lambda=0) Siga as analises considerando esta (ou alguma outra TRANSFORMACAO) likfit() retorna logLik, AIC, BIC para os odelos ajustados e portanto é só estratir e comparar, em especial apra avaliar os efeitos das potencias covariaveis (landuse e rock) por algum teste Lembre que na pagina da geoR tem uma sessao de tutotias e varios scripts com exemplo, alem do vignette do pacote On Sun, 12 Jun 2011, nani.design@bol.com.br wrote:
Boa noite,
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. Estou utilizando o conjunto de dados "Jura", incluido no pacote "gstat", 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: Landuse: 1 = floresta; 2 = pastagem; 3 = campo; 4 = lavoura
Rock: 1 = Argovian; 2 = Kimmeridgian; 3 = Sequanian; 4 = Portlandian; 5 = Quaternary
Utilizando o pacote "geoR", gerei os variogramas de cada metal através da função "variog". Por exemplo, o variograma do metal Cobre através do R:
library(gstat) library(geoR)
data(jura) prediction.dat # conjun to de dados Jura para o estudo
jura_dados.Cu <- as.geodata(prediction.dat, data.col = 8, covar.col = 3:4) varioJura.Cu <- variog(jura_dados.Cu, max.dist=4.5) # variograma do Cobre plot(varioJura.Cu, main="Variograma - Cobre", type="b")
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.
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?
Não sei se consegui ser bem clara na minha dúvida, mas se alguém souber me ajudar ficarei imensamente gra ta. Obrigada,
Ariane