Regression logistique

[Chérif Diallo]

Bonjour,

Je réalise une régression logistique pour expliquer une variable Trafic (qui a 3 modalités 0,1,2)

soit donc: regression<-glm(Trafic~.,data=bdd,family=binomial)

Dans mes variables explicatives, j'ai une variable " commune" qui possède les modalités suivantes:
Basse commune, Villecentre , boloc, zairo.
Je voulais prendre zairo en référence dans ma regression logistique pour cette variable "commune", quel est la fonction R que je dois utiliser?

[Alexandre Dangléant]

Bonjour,

Par défaut, le premier niveau du facteur est pris comme référence. Donc il suffit de réordonner les niveaux du facteur, par exemple :

Code : Tout sélectionner

fac <- factor(fac, levels = c( "zairo", "Basse commune", "Villecentre", "boloc"))

Il existe aussi la fonction relevels il me semble, mais j'ai toujours préféré recréer le facteur...

Cela étant dit, attention : vous voulez faire une régression multinomiale puisque votre votre variable réponse à 3 modalités, et non binomiale. De mémoire, glm(..., family = "binomial") n'est pas adaptée car elle fera implicitement la comparaison entre la référence et les autres modalités (mélangées!), cela se voit dans le summary du modèle car il n'y a pas une ligne par modalité de la variable réponse. Il existe des packages spécifiques pour cela, "multinom." de mémoire, mais j'ai un doute...

Bonne journée.

[Logez Maxime]

Bonjour,

dans le même ordre d'idée il existe la fonction relevel :

Code : Tout sélectionner

fac <- c("Basse commune", "Villecentre", "boloc", "zairo")
fac<- factor(fac)
fac
1] Basse commune Villecentre   boloc         zairo       
Levels: Basse commune boloc Villecentre zairo

relevel(fac, 4)
[1] Basse commune Villecentre   boloc         zairo       
Levels: zairo Basse commune boloc Villecentre


Cordialement,
Maxime