Comment puis-je tester les effets dans une ANOVA Split-Plot en utilisant des comparaisons de modèles appropriées à utiliser avec les arguments X
et M
de anova.mlm()
dans R? Je connais bien ?anova.mlm
Dalgaard (2007) [1]. Malheureusement, il ne brosse que les conceptions Split-Plot. Faire cela dans une conception entièrement randomisée avec deux facteurs intra-sujets:
N <- 20 # 20 subjects total
P <- 3 # levels within-factor 1
Q <- 3 # levels within-factor 2
DV <- matrix(rnorm(N* P*Q), ncol=P*Q) # random data in wide format
id <- expand.grid(IVw1=gl(P, 1), IVw2=gl(Q, 1)) # intra-subjects layout of data matrix
library(car) # for Anova()
fitA <- lm(DV ~ 1) # between-subjects design: here no between factor
resA <- Anova(fitA, idata=id, idesign=~IVw1*IVw2)
summary(resA, multivariate=FALSE, univariate=TRUE) # all tests ...
Les comparaisons de modèles suivantes conduisent aux mêmes résultats. Le modèle restreint n'inclut pas l'effet en question mais tous les autres effets du même ordre ou inférieur, le modèle complet ajoute l'effet en question.
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw2, test="Spherical") # IVw1
anova(fitA, idata=id, M=~IVw1 + IVw2, X=~IVw1, test="Spherical") # IVw2
anova(fitA, idata=id, M=~IVw1 + IVw2 + IVw1:IVw2,
X=~IVw1 + IVw2, test="Spherical") # IVw1:IVw2
Une conception Split-Splot avec un facteur intra et un entre sujets:
idB <- subset(id, IVw2==1, select="IVw1") # use only first within factor
IVb <- gl(2, 10, labels=c("A", "B")) # between-subjects factor
fitB <- lm(DV[ , 1:P] ~ IVb) # between-subjects design
resB <- Anova(fitB, idata=idB, idesign=~IVw1)
summary(resB, multivariate=FALSE, univariate=TRUE) # all tests ...
Ce sont les anova()
commandes pour répliquer les tests, mais je ne sais pas pourquoi ils fonctionnent. Pourquoi les tests des comparaisons de modèles suivants conduisent-ils aux mêmes résultats?
anova(fitB, idata=idB, X=~1, test="Spherical") # IVw1, IVw1:IVb
anova(fitB, idata=idB, M=~1, test="Spherical") # IVb
Deux facteurs intra-sujets et un facteur inter-sujets:
fitC <- lm(DV ~ IVb) # between-subjects design
resC <- Anova(fitC, idata=id, idesign=~IVw1*IVw2)
summary(resC, multivariate=FALSE, univariate=TRUE) # all tests ...
Comment reproduire les résultats ci-dessus avec les comparaisons de modèles correspondantes à utiliser avec les arguments X
et M
de anova.mlm()
? Quelle est la logique derrière ces comparaisons de modèles?
EDIT: suncoolsu a souligné qu'à toutes fins pratiques, les données de ces conceptions devraient être analysées à l'aide de modèles mixtes. Cependant, j'aimerais toujours comprendre comment reproduire les résultats de summary(Anova())
avec anova.mlm(..., X=?, M=?)
.
[1]: Dalgaard, P. 2007. Nouvelles fonctions pour l'analyse multivariée. R News, 7 (2), 2-7.
lme4
package pour s'adapter au modèle ET NONlm
. Mais cela peut être une vue très spécifique basée sur un livre. Je laisserai les commentaires des autres à ce sujet. Je peux donner un exemple basé sur la façon dont je l'interprète qui est différent du vôtre.