calcul de l'AICcorrigé

Questions sur les fonctions statistiques de R

Modérateur : Groupe des modérateurs

verochevalier
Messages : 2
Enregistré le : 15 Avr 2005, 15:40

calcul de l'AICcorrigé

Messagepar verochevalier » 17 Oct 2005, 15:38

Bonjour

Comment fait-on pour calculer un AIC corrigé à partir de l'AIC , pour un modèle de Poisson (GLM) ?

Merci d'avance
Véronique

François Bonnot
Messages : 537
Enregistré le : 10 Nov 2004, 15:19
Contact :

Messagepar François Bonnot » 19 Oct 2005, 06:56

Bonjour,

En principe on peut calculer l'AIC corrigé par

AICc = AIC + 2*k*(k+1)/(n-k-1)
avec
k = nombre de paramètres à estimer
n = nombre d'observations

Je pense qu'il doit y avoir tous les détails dans :
Burnham, K.P. & Anderson, D.R. (2002). Model Selection and Multimodel Inference: a Practical Information - Theoretic Approach. Sringer-Verlag, New-York.
(bien que je ne dispose pas de cet ouvrage)

Donc avec glm j'écrirais :

Code : Tout sélectionner

res <- glm(...)
k <- res$df.null-res$df.residual+1
n <- res$df.null+1
cor <- 2*k*(k+1)/(n-k-1)
aicc <- res$aic+cor

A confirmer quand-même par d'autres utilisateurs...

verochevalier
Messages : 2
Enregistré le : 15 Avr 2005, 15:40

Messagepar verochevalier » 19 Oct 2005, 15:22

Merci , je vais essayer comme ça
Véro

Renaud Lancelot
Messages : 2484
Enregistré le : 16 Déc 2004, 08:01
Contact :

Messagepar Renaud Lancelot » 01 Déc 2005, 16:24

Pour information: le package metomet, disponible sur ce site, contient maintenant une fonction sic() (some information criteria) permettant de calculer l'AIC, l'AICc et le BIC pour des modèles lm et glm.

Amicalement,

Renaud


Retourner vers « Archives : Fonctions statistiques »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité