[R-br] geoestatística - ajuste de modelo

[Paulo Justiniano]

É 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 em 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
> 
> 
>