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