[R-es] Error seleccionando modelos por el AIC

[Jorge I Velez]

Hola Javier,

La ayuda de ?aictabe establece que la implementacion del AIC no funciona
para objetos de clase "lmer" (a la que pertenecen los modelos ajustados).
De ahi el error.

La idea es simplemente comparar los modelos basados en el AIC?  O tienes
pensado algo mas?  Si es lo primero, los objetos lmer contienen el AIC --
es solo cuestion de usar str(objetolmer), ubicar la entrada correspondiente
y  extraerla.   Si es lo segundo, ahi tendrias tu que ayudarnos con un poco
mas de informacion.

Saludos cordiales,
Jorge.-


On Fri, Jan 30, 2015 at 6:37 AM, javier bueno enciso <
jbuenoenciso en hotmail.com> wrote:

> Hola, os escribo por un error que me sale al tratar de hacer una tabla de
> selección de modelos según AIC, con la función aictab() del paquete
> AICcmodavg.
> Os copio aquí el script que utilizo:
>
> m1<- lmer(ldate~A+(1|zona/area), data=Data, subset= (Data$brood== 0))
> m2<-lmer(ldate~B+(1|zona/area),  data=Data, subset= (Data$brood== 0))
> m3<-lmer(ldate~A*B+(1|zona/area),data=Data, subset= (Data$brood== 0))
> m4<-lmer(ldate~A+B+(1|zona/area), data=Data, subset= (Data$brood== 0))
> m5<-lmer(ldate~A+A*B+(1|zona/area), data=Data, subset= (Data$brood== 0))
> m6<-lmer(ldate~B+A*B+(1|zona/area), data=Data, subset= (Data$brood== 0))
> m7<-lmer(ldate~A+B+A*B+(1|zona/area), data=Data, subset= (Data$brood== 0))
>
> ## library("AICcmodavg", lib.loc="~/R/win-library/3.1")
> modList <- list(glm1, glm2, glm3, glm4, glm5, glm6, glm7)
> Modnames<- c("1", "2", "3", "4", "5", "6", "7")
>
> aictab(cand.set= modList, modnames= Modnames)
>
> El error es el siguiente:
> Error in aictab.default(cand.set = modList, modnames = ModNames) :
> Function not yet defined for this object classHe buscado por internet pero
> no consigo solucionarlo, cualquier ayuda sería de gran utilidad.
>
> Muchas gracias.
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

	[[alternative HTML version deleted]]