AICc non calculable

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

E.H. [compte supprimé]

AICc non calculable

Messagepar E.H. [compte supprimé] » 21 Mai 2007, 12:36

Bonjour tout le monde,

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 ;-)

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité