Les posteurs les plus actifs de la semaine
Aucun utilisateur |
Sujets les plus vus
Regression logistique avec que des zéro pour un groupe
3 participants
Page 1 sur 1
Regression logistique avec que des zéro pour un groupe
Bonjour !
Je souhaite comparer la réaction d'individus entre 2 types de plantes.tout en contrôlant pour l'effet groupe d'age (stage_ad) et le sexe.
Or pour une plante, aucun mâle ne consomme, j'ai donc que des zéro pour tous ce groupe (plante1-male) et donc aucune variance. Cependant, au lieu d'avoir quelque chose très significatif, j'ai un "estimate" énorme et un p proche de 1?
Y a-t-il un moyen pour pallier à ce genre de biais?
Ci dessous mes lignes de codes et résultats sur R:
model1<-glmer(Peat~(as.factor(Ses) + Type + Sexe + as.factor(Stage_ad) + Type:Sexe +(1|Idi) ), family=binomial,na.action=na.exclude,data1)
> summary(model1)
Generalized linear mixed model fit by the Laplace approximation
Formula: Psmel ~ (as.factor(Ses) + Type + Sexe + as.factor(Stage_ad) + Type:Sexe + (1 | Idi))
Data: data1
AIC BIC logLik deviance
169.6 199.1 -75.82 151.6
Random effects:
Groups Name Variance Std.Dev.
Idi (Intercept) 0.61361 0.78333
Number of obs: 195, groups: Idi, 13
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.3665 0.9456 -2.5027 0.0123 *
as.factor(Ses)2 -1.0035 0.5128 -1.9570 0.0503 .
as.factor(Ses)3 -0.8535 0.5005 -1.7055 0.0881 .
TypePlante 0.9230 0.8906 1.0364 0.3000
SexeM -16.8616 2349.1866 -0.0072 0.9943
as.factor(age)2 1.8705 0.7494 2.4960 0.0126 *
as.factor(age)3 -0.3371 1.1848 -0.2846 0.7760
TypePlante:SexeM 16.6156 2349.1865 0.0071 0.9944
Merci d'avance
dao
Je souhaite comparer la réaction d'individus entre 2 types de plantes.tout en contrôlant pour l'effet groupe d'age (stage_ad) et le sexe.
Or pour une plante, aucun mâle ne consomme, j'ai donc que des zéro pour tous ce groupe (plante1-male) et donc aucune variance. Cependant, au lieu d'avoir quelque chose très significatif, j'ai un "estimate" énorme et un p proche de 1?
Y a-t-il un moyen pour pallier à ce genre de biais?
Ci dessous mes lignes de codes et résultats sur R:
model1<-glmer(Peat~(as.factor(Ses) + Type + Sexe + as.factor(Stage_ad) + Type:Sexe +(1|Idi) ), family=binomial,na.action=na.exclude,data1)
> summary(model1)
Generalized linear mixed model fit by the Laplace approximation
Formula: Psmel ~ (as.factor(Ses) + Type + Sexe + as.factor(Stage_ad) + Type:Sexe + (1 | Idi))
Data: data1
AIC BIC logLik deviance
169.6 199.1 -75.82 151.6
Random effects:
Groups Name Variance Std.Dev.
Idi (Intercept) 0.61361 0.78333
Number of obs: 195, groups: Idi, 13
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -2.3665 0.9456 -2.5027 0.0123 *
as.factor(Ses)2 -1.0035 0.5128 -1.9570 0.0503 .
as.factor(Ses)3 -0.8535 0.5005 -1.7055 0.0881 .
TypePlante 0.9230 0.8906 1.0364 0.3000
SexeM -16.8616 2349.1866 -0.0072 0.9943
as.factor(age)2 1.8705 0.7494 2.4960 0.0126 *
as.factor(age)3 -0.3371 1.1848 -0.2846 0.7760
TypePlante:SexeM 16.6156 2349.1865 0.0071 0.9944
Merci d'avance
dao
dao06- Nombre de messages : 3
Date d'inscription : 14/08/2009
Re: Regression logistique avec que des zéro pour un groupe
Il est normal que le coefficient soit très grand et qu'il ne soit pas significatif. Comme tu intègres une interaction pour entre le type de plante et le sexe et que pour l'une des combinaison tu as très très peu d'individus, alors l'incertitude que tu as quand tu estimes ton coefficient est très grande, d'où un écart-type très élevé et une valeur non significative. De plus, suivant si pour cette combinaison plante-sexe tu as un succès ou un échec alors tu vas avoir un coefficient très fort (en valeur absolue) qui représente le rapport des chances de cette combinaison.
Après il faut savoir si c'est un biais expérimental qui fait que tu n'as pas de mâle sur cette plante ou si c'est une réalité biologique.
Après il faut savoir si c'est un biais expérimental qui fait que tu n'as pas de mâle sur cette plante ou si c'est une réalité biologique.
droopy- Nombre de messages : 1156
Date d'inscription : 04/09/2009
Re: Regression logistique avec que des zéro pour un groupe
Merci beaucoup pour une aussi rapide réponse !
Il ne s'agit pas d'un biais d'échantillonnage à priori. En fait les donnée sont codée 0 (n'a pas consommé) et 1 (a consommé).
Du coup le fait qu'aucun male n'a consommé la plante1 a une signification biologique en relation avec mon hypothèse de départ. Mais toutes leurs valeurs étant codée 0, il n'ya aucune variance et le test le met comme non significatif. Si je truque mes donnée en rajoutant un 1 pour un individu, le test sort très significatif,
mais y a-t-il moyen d'y arriver sans truquer ces données?
Merci d'avance
Il ne s'agit pas d'un biais d'échantillonnage à priori. En fait les donnée sont codée 0 (n'a pas consommé) et 1 (a consommé).
Du coup le fait qu'aucun male n'a consommé la plante1 a une signification biologique en relation avec mon hypothèse de départ. Mais toutes leurs valeurs étant codée 0, il n'ya aucune variance et le test le met comme non significatif. Si je truque mes donnée en rajoutant un 1 pour un individu, le test sort très significatif,
mais y a-t-il moyen d'y arriver sans truquer ces données?
Merci d'avance
dao06- Nombre de messages : 3
Date d'inscription : 14/08/2009
Re: Regression logistique avec que des zéro pour un groupe
Hello,
Je m'étais déjà posé une question dans le genre je pense et droopy m'avais fournis une réponse très claire ici:
https://statistiques.forumpro.fr/questions-generales-f9/regression-logistique-avec-variable-independante0-t1128.htm
Bon courage
Je m'étais déjà posé une question dans le genre je pense et droopy m'avais fournis une réponse très claire ici:
https://statistiques.forumpro.fr/questions-generales-f9/regression-logistique-avec-variable-independante0-t1128.htm
Bon courage
jigouen- Nombre de messages : 54
Date d'inscription : 04/09/2009
Re: Regression logistique avec que des zéro pour un groupe
Merci beaucoup !
Je pense que c'est la même situation en effet !
Je pense que c'est la même situation en effet !
dao06- Nombre de messages : 3
Date d'inscription : 14/08/2009
Sujets similaires
» régression logistique binaire avec R
» Validation code regression logistique Binaire avec R
» regression logistique avec variable independante=0%
» Régression logistique généralisée, var avec 10 modalités
» Modéliser un taux avec une régression logistique
» Validation code regression logistique Binaire avec R
» regression logistique avec variable independante=0%
» Régression logistique généralisée, var avec 10 modalités
» Modéliser un taux avec une régression logistique
Page 1 sur 1
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum