Discussion :

[mrdeep]

Bonjour à tous,

Je suis nouveau avec R et dans les stats aussi ...

Je souhaiterais effectuer une ANOVA avec 2 facteurs, un facteur indépendant (Groupe) x un facteur à mesure répété (Bras).

J'ai trouvé plein de renseignements sur comment faire, j'ai essayé plusieurs méthodes différentes, notamment avec le package afex.

Voici ce que j'ai essayé :

-Avec afex :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

aov_ez(id="Groupe", dv="myVD", data=myDATA, between = "Groupe", within = "Bras")

Et j'obtiens :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
Univariate Type III Repeated-Measures ANOVA Assuming Sphericity
 
                 SS num Df  Error SS den Df      F Pr(>F)
(Intercept) 0.75747      1 0.0216294      1 35.021 0.1066
Bras          0.01060      1 0.0014469      1  7.326 0.2253

Est-ce que cela vous parait normal ? On dirait qu'il manque des données d'analyse ...

-Avec une méthode trouvé sur internet :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
anova_2f = aov(data=myDATA, myVD~Bras*Groupe + Error(Groupe/(Bras*Groupe)) )
summary(anova_2f)

Je ne suis pas sur de la méthode ... Quelqu'un peut-il m'éclairer ou m'indiquer comment procéder à mon ANOVA ?

Merci par avance pour votre aide.

Cordialement,

[faubry]

Il existe de nombreuses manières pour traiter les Anova à mesures répétées. Les modèles sous-jacents sont différents ; de plus es-tu vraiment sur de ton modèle d'erreur dans aov (relis la doc) :
- liste générale recensant les approches :
1) http://www.r-statistics.com/2010/04/...h-r-tutorials/
2) https://gribblelab.wordpress.com/200...anova-using-r/
- utilisation de aov
1) https://cran.r-project.org/doc/contrib/Baron-rpsych.pdf
- modèle mixte http://www.uni-kiel.de/psychologie/r...novaMixed.html

[Alain.Blanchet]

Le package AFEX est excellent pour ce type d'analyse, mais votre formule n'est pas correcte :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
x<-aov_ez("sujets", "variable mesurée", between = "groupe", within = "bras")
summary(x)
anova(x)

Comme ceci, ça devrait marcher.

Cordialement