Les posteurs les plus actifs de la semaine
Lalie17
Régression logistique et anova sur R Vote_lcapRégression logistique et anova sur R Voting_barRégression logistique et anova sur R Vote_rcap 


Régression logistique et anova sur R

Aller en bas

Régression logistique et anova sur R Empty Régression logistique et anova sur R

Message par Amorphis le Dim 15 Aoû 2010 - 10:17

Bonjour,
N'étant pas douée du tout en statistique, je suis un peu paumée dans mes calculs...
En fait, j'essaie d'expliquer une variable de type binaire (présence=1 ou absence=0) en fonction de 3 variables explicatives :
X1 : de 1 à 7
X2 : de 1 à 3
X3 : de 1 à 5

exemple de mon tableau :

Y X1 X2 X3
1 6 2 4
1 6 2 3
0 6 2 2
1 7 3 4
1 7 3 2
0 7 3 1


j'ai utilisé la formule glm mais lorsque je veux rajouter la fonction "family=binomial" ou "family=quasibinomial", R me met :
Erreur dans eval(expr, envir, enclos) :
les valeurs de y doivent être 0 <= y <= 1

Or, mon Y varie bien entre 0 et 1

Du coup je l'ai pas marqué dans la formule.. mais j'imagine que c'est important de le marquer...

Puis j'ai vu sur le net qu'il fallait faire une anova pour analyser les données mais ça ne marche pas.
J'ai tapé :
glm1=glm(Y~X1+X2+X3) (bon je sais qu'il manque la fonction family)
anova(glm1, test = "Chisq")
Et là, R me sort :
Erreur dans anova(glm1, test = "Chisq") :
argument(s) inutilisé(s) (test = "Chisq")


Du coup, qu'est-ce qu'il aurait fallu faire?
Merci d'avance de votre aide.
Sophie

Amorphis

Nombre de messages : 1
Date d'inscription : 15/08/2010

Voir le profil de l'utilisateur

Revenir en haut Aller en bas

Régression logistique et anova sur R Empty Re: Régression logistique et anova sur R

Message par droopy le Lun 16 Aoû 2010 - 7:32

bonjour, c'est surement un problème de virgule floattante. C'est à dire que ce que tu penses être un 1 est en fait 1.0000000000000000000235 quelque chose dans ce gout la. Du coup quand tu démarres le glm en family binomial il contrôle tes données et tu ne passes pas ce stade. Une façon d'en être sûre :
Code:
sprintf("%.30f", y)
# ou
all(y<=1 & y>=0)
Et tu verras si tu as des valeurs qui sortent de l'intervalle 0-1.
droopy
droopy

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

Voir le profil de l'utilisateur

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