Facteurs et contrasts régréssion

Bonjour,

Ma régression ne passe pas alors que mes données sont bien des facteurs.
Je vous montre :

Code:

---

1 2 3

data_modelisation$ec <- as.factor(data_modelisation$ec) data_modelisation$tm <- as.factor(data_modelisation$tm) modele=glm(RA ~ C(ec,sum) + C(tm,sum) ,family=binomial(logit), data=data_modelisation)

---

RA étant défini plus haut dans le code.

J'obtiens :

Code:

---

1 2	Erreur dans `contrasts<-`(`*tmp*`, value = contr) : contrasts can be applied only to factors with 2 or more levels

---

Alors que :

Code:

---

1 2 3 4

> levels(data_modelisation$ec) [1] "1" "2" "3" "4" "5" "6" > levels(data_modelisation$tm) [1] "1" "2" "3" "4"

---

Merci à ceux qui pourront m'aider

Citation:

---

Envoyé par Stepsbysteps

Bonjour,

Ma régression ne passe pas alors que mes données sont bien des facteurs.
Je vous montre :

RA étant défini plus haut dans le code.

---

Bonjour,

Quel est l'intérêt d'utiliser un contraste avec sum?
Il se passe quoi quand tu fais attributes(C(data_modelisation$ec,sum)) ?


Peut on avoir le summary ou le table de RA?

Bonjour,

Merci.

Les RA sont des 0 ou des 1.
Je veux avoir la somme des effets égales à 0 et non le première effet pris en référence.

Et j'obtiens

Code:

---

1 2 3 4 5 6 7 8 9

> attributes(C(data_modelisation$ec,sum)) $levels [1] "1" "2" "3" "4" "5" "6"   $class [1] "factor"   $contrasts [1] "contr.sum"

---

Ok bon recapitulons.

Il y a un problème de contrast.
Or ce n'est pas la fonction C() qui fait bugger.

Si tu lances ta regression sans le C() je pense que ca renverra exactement la même erreur (dis moi si je me trompe).

Diagnostique?
Une de tes variables explicatives n'a qu'une modalité.

"Oui mais nan car:"

Code:

---

1 2	> levels(data_modelisation$tm) [1] "1" "2" "3" "4"

---

Sauf que :

Code:

---

1 2 3 4 5 6

facteur<-as.factor(c(1,2,1,1,1)) facteur<-facteur[-2] > facteur [1] 1 1 1 1 Levels: 1 2

---

:aie:

haha sympa le message il m'as fait sourire!

Nan sans le C ça marche bien!

Et bizarre pour le facteur alors faut que je regarde de plus près

Citation:

---

Diagnostique?
Une de tes variables explicatives n'a qu'une modalité.

---

Bizarre car :

Arf dommage, j'y ais cru car c'est exactement la même erreur que quand tu fais une regression avec un facteur d'une modalité, et comme tu arrives à utiliser C() en dehors de la regression...

Heuresement que je suis pas docteur :mouarf:

Je sais pas ce qui fait planter du coup :ccool:

Arf dommage merci quand même

Ca marche ainsi,

Code:

---

modele=glm(RA ~ ec + tm ,family=binomial(logit), data=data_modelisation,contrasts=list(ec="contr.sum",tm="contr.sum"))

---

:ccool: