ANOVA à 2 facteurs et à mesures appariées : test de sphéricité

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Matisse Poupard
Messages : 1
Enregistré le : 04 Avr 2020, 15:20

ANOVA à 2 facteurs et à mesures appariées : test de sphéricité

Messagepar Matisse Poupard » 05 Avr 2020, 11:12

Bonjour tout le monde,

N'ayant que quelques bases sur r je me heurte à un problème dans une analyse de données. J'aimerai effectuer une analyse de variance sur le jeu de données* suivant (5 premières lignes):

sujet;score;heure;sexe
sujet 1;17;Matin;Homme
sujet 2;18;Matin;Homme
sujet 3;17;Matin;Homme
sujet 4;18;Matin;Homme
sujet 5;18;Matin;Homme

Pour cela voici mon code : (j'importe les données depuis un fichier csv, je ne trouve pas comment vous le transmettre)

Code : Tout sélectionner

library(multcompView)
library(car)
library(rstatix)

data_ex5 <- read.csv("exo5.csv",header = TRUE, sep=";")

score<-data_ex5$score
heure<-data_ex5$heure
sexe<-data_ex5$sexe
sujet<-data_ex5$sujet

by(score, list(heure,sexe), shapiro.test)


anova2x2<-aov(score~heure+sexe+heure*sexe+Error(factor(sujet)),data = data_ex5)
summary(anova2x2)


J'aimerai maintenant vérifier l'hypothèse de sphéricité (test de Mauchly). Pour cela j'utilise le package "rstatix", et la fonction anova_test().

Elle fait automatiquement le test de Mauchly et la correction si nécessaire. Cela marche sur un jeu de données à un facteur (par exemple score en fonction de l'heure), pourtant dans ce cas précis (score en fonction de l'heure et du sexe) la fonction ne marche pas :

Code : Tout sélectionner

anova_test(score~heure*sexe+Error(sujet/(heure*sexe)),data=data_ex5)


Erreur :
Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
0 (non-NA) cases


Je ne sais donc pas comment faire pour tester cette sphéricité autrement, ni pourquoi cette méthode ne marche pas...

Merci !

* je ne trouve pas comment vous transmettre mon jeu de données

Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité