[R-br] geoestatística - AIC

Hélio Gallo Rocha heliogallorocha em gmail.com
Sexta Novembro 8 22:16:50 BRST 2013


Entendi...

Muitíssimo grato prof. Paulo,por mais essa!

Hélio


Em 8 de novembro de 2013 18:17, Paulo Justiniano [via R-br] <
ml-node+s2285057n4660900h68 em n4.nabble.com> escreveu:

> 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
> >
> >
> _______________________________________________
> R-br mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=4660900&i=0>
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
>
> http://r-br.2285057.n4.nabble.com/R-br-geoestatistica-AIC-tp4660891p4660900.html
>  To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
> .
> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>



-- 
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131108/f01b92e6/attachment.html>


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