Les posteurs les plus actifs de la semaine
louiv
 
gg
 
Eric Wajnberg
 


Regression logistique avec que des zéro pour un groupe

Aller en bas

Regression logistique avec que des zéro pour un groupe

Message par dao06 le Mar 23 Fév 2010 - 11:07

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

dao06

Nombre de messages : 3
Date d'inscription : 14/08/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Regression logistique avec que des zéro pour un groupe

Message par droopy le Mar 23 Fév 2010 - 12:55

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.
avatar
droopy

Nombre de messages : 1097
Date d'inscription : 04/09/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Regression logistique avec que des zéro pour un groupe

Message par dao06 le Mar 23 Fév 2010 - 14:47

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

dao06

Nombre de messages : 3
Date d'inscription : 14/08/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Regression logistique avec que des zéro pour un groupe

Message par jigouen le Mar 23 Fév 2010 - 19:01

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:
http://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

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Regression logistique avec que des zéro pour un groupe

Message par dao06 le Mer 24 Fév 2010 - 9:41

Merci beaucoup !
Je pense que c'est la même situation en effet !

dao06

Nombre de messages : 3
Date d'inscription : 14/08/2009

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Re: Regression logistique avec que des zéro pour un groupe

Message par Contenu sponsorisé


Contenu sponsorisé


Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum