Discussion :

[draculla]

Bonjour tout le monde :-)

je voudrais sélectionner parmi des modèles logistiques, en utilisant le BIC, le meilleur modèle.

je voudrai savoir s'il existe un package pour ça (le package metomet n'est plus valable), ou bien alors il faut le calculer manuellement pour tous les modèles, et y'en a beaucoup.


Merci les amis.

[chautho]

Bonjour,

Je ne sais pas si par "sélection de modèle" tu entends une méthode types step forward ou backward mais si c'est le cas tu peux faire en sorte que l'algorithle montant ou descendant se fasse en fonction du critère BIC tu peux ajouter à la fonction step() l'option "k=log(n)" pour que le critère soit le critère BIC et non plus AIC

[draculla]

Bonjour, merci de la réponse.

Vous voulez dire que si je fais ça

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

step(modele_initial,scope=list(lower="Y~1", upper=as.formula(modele_additif,trace=TRUE,direction="forward",k=log(nrow(data)))))

Avec l'option k=log(nrow du data), le critére c'est le BIC?

Il ne s'affiche pas par contre, c'est le AIC qui s'affiche !

[chautho]

Tu as comparé entre l'exécution classique et l'exécution avec l'option pour faire le critère BIC ?

Parce qu'il me semble que R continue à t'afficher comme si c'était l'AIC mais en faite c'est bien le BIC qui est calculé...

[draculla]

Oui, je viens de vérifier, et c'est le même, il sélectionne avec le AIC dans les deux cas.

penses-tu que ça marche bien avec le package (BMA)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
glm.out.FT <- bic.glm(X, Y, strict = T, OR = 20,glm.family=binomial(logit), factor.type=TRUE)
 
 
summary(glm.out.FT)
imageplot.bma(glm.out.FT)

[chautho]

Je n'ai jamais utilisé ce package.

Cependant, je ne comprends pas pourquoi cela ne marche sachant que si tu vas dans l'aide de la fonction step() :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

?step

Il y a bien une partie qui précise :

k : the multiple of the number of degrees of freedom used for the penalty. Only k = 2 gives the genuine AIC: k = log(n) is sometimes referred to as BIC or SBC.

Etrange :/