Juste pour une précision, je suis dans une situation ou le calcul de l'AICc est impossible car n-k-1 = 0 --> division par zero impossible / résultat AICc infini.
Or, le rapport n/k est inférieur à 40, je devrais donc préférentiellement utiliser AICc à AIC.
Une idée ? utiliser l'AIC en spécifiant dans ma démarche que le calcul de l'AICc est non faisable ? Je pose cette question, car les résultats des variations d'AIC et AICc ne me donnent pas du tout les même choses...
Exemple:
Code : Tout sélectionner
> DATA <- fcbe
>
> # modèle nul
> m00 = ASRF3~ta*co*po-ta*co*po
> aov00 <- aov(m00, data=DATA)
> monAIC(aov00)
n k nsk AIC AICc
aov00 10 2 5 17.15587 18.87016
>
> # modèle complet
> m0 <- ASRF3~ta*co*po
> aov0 <- aov(m0, data=DATA)
> monAIC(aov0)
n k nsk AIC AICc
aov0 10 9 1.111111 15.7834 Inf
>
> # modèle sans l'interaction double
> m1 <- ASRF3~ta+co+po+ta:co+ta:po+co:po
> aov1 <- aov(m1, data=DATA)
> monAIC(aov1)
n k nsk AIC AICc
aov1 10 9 1.111111 15.7834 Inf
>
> # modèles sans l'interaction triple, avec retrait d'une intéraction simple à chaque fois
> m1a <- ASRF3~ta+co+po+ta:co+ta:po #+co:po
> aov1a <- aov(m1a, data=DATA)
>
> m1b <- ASRF3~ta+co+po+ta:co+co:po #+ta:po
> aov1b <- aov(m1b, data=DATA)
>
> m1c <- ASRF3~ta+co+po+ta:po+co:po #+ta:co
> aov1c <- aov(m1c, data=DATA)
>
> monAIC(aov1, aov1a, aov1b, aov1c)
n k nsk AIC AICc
aov1 10 9 1.111111 15.7834 Inf
aov1a 10 9 1.111111 15.7834 Inf
aov1b 10 9 1.111111 15.7834 Inf
aov1c 10 9 1.111111 15.7834 Inf
> m2 <- ASRF3~ta+co+po
> aov2 <- aov(m2, data=DATA)
> monAIC(aov2)
n k nsk AIC AICc
aov2 10 8 1.25 14.06430 158.0643
>
> m2a <- ASRF3~ta+co #+po
> aov2a <- aov(m2a, data=DATA)
>
> m2b <- ASRF3~ta+po #+co
> aov2b <- aov(m2b, data=DATA)
>
> m2c <- ASRF3~co+po #+ta
> aov2c <- aov(m2c, data=DATA)
>
> monAIC(aov2, aov2a, aov2b, aov2c)
n k nsk AIC AICc
aov2 10 8 1.250000 14.06430 158.06430
aov2a 10 7 1.428571 16.02761 72.02761
aov2b 10 5 2.000000 17.28014 32.28014
aov2c 10 7 1.428571 13.00565 69.00565
--> selon les variations de l'AIC, je devrais conserver dans mon modèle les effets principaux po et co
--> selon les variations de l'AICc, je devrais partir sur un modèle nul
Auriez vous un conseil à me donner ?
Merci d'avance ;-)