Forum de Statistiques
Vous souhaitez réagir à ce message ? Créez un compte en quelques clics ou connectez-vous pour continuer.
Les posteurs les plus actifs de la semaine
Aucun utilisateur

Le Deal du moment : -17%
Casque de réalité virtuelle Meta Quest 2 ...
Voir le deal
249.99 €

Régression logistique et anova sur R

2 participants

Aller en bas

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

Message par Amorphis 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

Revenir en haut Aller en bas

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

Message par droopy 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 : 1156
Date d'inscription : 04/09/2009

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