Discussion :

[lmoulin]

Bonjour à tous,

J'essaie d'estimer un modèle à l'aide d'une régression logistique en ayant deux variables explicatives liées. Voici mon modèle :

glm(Y ~ X1 + X2 + X3 * X4 + X3 + X4), family = binomial,data = XXX)

X3 et X4 sont des variables binaires et ce sont elles qui sont linéairement dépendantes. Ainsi X3 représente deux années, disons 2000 et 2001 et X4 une politique économique binaire, par exemple verser une allocation chômage, qui n'existe pas en 2000 mais qui existe en 2001 mais ne touche pas tous les individus.

Je pensais qu'il fallait mettre un signe * entre les variables linéairement dépendantes pour tenir compte de l'indépendance mais apparemment ce n'est pas la bonne solution puisque R m'indique "Coefficients: (1 not defined because of singularities)" tout en me sortant des résultats.

Une idée ?

Merci d'avance

[Sengar]

Hello,

L'erreur renvoyée signifie qu'une variable est combinaison linéaire d'autres variables du modèle.

A ce que j'ai compris:

X3=2000 (0 en binaire) => X4=Pas d'allocation (0 en binaire)
X3=2001 (1 en binaire)=> on ne peut rien conclure sur X4

Du coup en binaire:
X4 = X3:X4
d'où l'erreur puisque dans ton modèle tu met à la fois X4 et X3:X4 qui sont deux variables identiques!

Conclusion: enlève X4 (ou X3:X4 au choix):
glm(Y ~ X1 + X2 + X3 :X4 + X3 ), family = binomial,data = XXX)

NB: Attention à l'operateur * dans une formule.
X3*X4 ca correspond à X3+X4+X3:X4.
Si tu veux juste le croisement, tu utilises ":" .

[lmoulin]

Merci beaucoup pour la réponse je n'avais pas reussi à la trouver sur internet

Juste une dernière question pour être sur d'obtenir ce que je cherche. J'aimerais avoir l'effet de la politique mise en place (X4) en contrôlant par l'année (X3) afin d'être sur que l'effet observé ne provienne pas d'un changement de l'échantillon entre les années.
J'obtiens les résultats suivant :

X3 -1.8e-01 1.3e-01 -1.4 0.16
X3:X4 modalité 1 1.4e-04 6.9e-05 2.0 0.05

La ligne qui m’intéresse (voir si la politique mise en place entre les deux années a un effet) est donc bien X3:X4 ?

En fait je me permets encore deux petites questions :
-Aurais tu une source sur l'utilisation de ":" dans les régressions ?
-Peut on l'utiliser de la même manière dans les régressions logistiques multinomiales ?

Merci d'avance

[Sengar]

Merci beaucoup pour la réponse je n'avais pas reussi à la trouver sur internet

J'obtiens les résultats suivant :

X3 -1.8e-01 1.3e-01 -1.4 0.16
X3:X4 modalité 1 1.4e-04 6.9e-05 2.0 0.05

La ligne qui m’intéresse (voir si la politique mise en place entre les deux années a un effet) est donc bien X3:X4 ?

Oui.
Et effectivement, la politique a un effet significatif puisque la p value est de 5%.
Ce n'est clairement pas le cas de l'année (p value de 16%).

En fait je me permets encore deux petites questions :
-Aurais tu une source sur l'utilisation de ":" dans les régressions ?

Oui: Moi
Tu peux me citer ca me dérange pas.

Aha nan serieusement si on demandait une source à chaque fois que quelqu'un répond à un problème de syntaxe R

-Peut on l'utiliser de la même manière dans les régressions logistiques multinomiales ?
Merci d'avance

Le ":" n'est aucunement une syntaxe propre à la fonction glm.
Il est propre aux objets de type formule.
En clair ca marche avec toute sorte de regression si tu veux tant que tu utilises un objet de type formula.

Au passage la "source" que tu voullais c'est l'aide de R sur formula si ca te tient vraiment à coeur ^^

An expression of the form y ~ model is interpreted as a specification that the response y is modelled by a linear predictor specified symbolically by model. Such a model consists of a series of terms separated by + operators. The terms themselves consist of variable and factor names separated by : operators. Such a term is interpreted as the interaction of all the variables and factors appearing in the term.

In addition to + and :, a number of other operators are useful in model formulae. The * operator denotes factor crossing: a*b interpreted as a+b+a:b.

[lmoulin]

De nouveau merci pour tes réponses qui vont m'être d'une grande utilité

Sujet résolu