AIC(lm(Fertility ~ ., data=swiss))
[1] 326.0716
ok, ya que el AIC se calcula como
-2*logLik(lm(Fertility ~ ., data=swiss)) + 2*7
¿Por qué el stepAIC devuelve un AIC menor?
stepAIC(lm(Fertility ~ ., data=swiss))
Start: AIC=190.69
AIC(lm(Fertility ~ ., data=swiss))
[1] 326.0716
ok, ya que el AIC se calcula como
-2*logLik(lm(Fertility ~ ., data=swiss)) + 2*7
¿Por qué el stepAIC devuelve un AIC menor?
stepAIC(lm(Fertility ~ ., data=swiss))
Start: AIC=190.69
Leer ?AIC
. En la sección de detalles tenemos:
The log-likelihood and hence the AIC is only defined up to an
additive constant. Different constants have conventionally be
used for different purposes and so ‘extractAIC’ and ‘AIC’ may give
different values (and do for models of class ‘"lm"’: see the help
for ‘extractAIC’).
Así que está viendo el efecto de diferentes constantes aditivas, como stepAIC()
(en el paquete MASS) está utilizando extractAIC()
para calcular el AIC de los modelos.
I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.