[R-br] geoestatística - AIC

Paulo Justiniano paulojus em leg.ufpr.br
Sexta Novembro 8 13:41:46 BRST 2013


Invertendo a pergunta... por que deveriam ser iguais??

likfit estimou um modelo gaussiano para os dados
com um termo espacial (efeito aleatório espacial)
No caso um modelo com 3 parametros, um de mdia e 2 de covariancia.

glm dos dados contra os preditos estimou uma regressao linear destes, 
e nao vvejo pq deveriam retornar os mesmo valores
com 2 parametros de média e um de variancia

Nao vejo pq deveriam ter as mesmas verossimilhancas ou AIC's


On Fri, 8 Nov 2013, Hélio Gallo Rocha wrote:

> Caros listeiros.
> Como já haviam citado,
> 
> Prof. Paulo:
> "Pensar se a validacao cruzada é de fato nbecessária. Embora seja 
> usada é bem menos informativa e util que outras medidadas de avaliacao do 
> ajuste baseadas na verossimilhanã e análise dos residuos" 
> 
> Éder citando Davis, B. M., 1987. Uses and abuses of cross-validation in
> geostatistics. Math. Geol., v. 17, p. 563–586.
> "Cross-validation  cannot  confirm  that  a  particular  model  is  or  is  not  the
>  optimum.  It  is  a method  to  better  examine  and  understand  the  phenomenon
>  under  study  using the  available  data"
> 
> 
> 
> 
> situação:
> Caculando o AIC, tive dois resultados diferentes, como segue CRM:
> 
> require(geoR);data(s100)
> vModel <- likfit(s100, ini=c(1,0.5), fix.nugget=T)
> xv.ml<-xvalid(s100,model=vModel)
> 
> AIC(vModel)                                              #173.1391
> reglin_vModel=glm(xv.ml$data~xv.ml$pred)
> AIC(reglin_vModel)                                    #142.4229
> 
> Qual a razão destes valores diferentes?
> 
> 
> 
> Desde já agradeço qualquer elucidação
> 
> 
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
> 
>


Mais detalhes sobre a lista de discussão R-br