data(MaskedPriming)
str(MaskedPriming)
plot(MaskedPriming[sample(6381,100), 7:10])
library(lme4)
cmat <- matrix(c(-1, 1, 0,
-1, -1, 2), 3, 2,
dimnames=list(c("BM1", "BM2", "SK"),
c(".BM1-2", ".BM-SK")))
m0 <- lmer(rrt ~ p*f*e + (1 | subjects) + (0 + p | subjects) +
(0 + f | subjects) + (1 | words), contrasts=list(e=cmat),
data=d)
m1p <- lmer(rrt ~ p*f*e + (p | subjects) + (0+f | subjects) + (1 | words),
contrasts=list(e=cmat)
m2 <- lmer(rrt ~ p*f*e + (p + f | subjects) + (1 | words),
contrasts=list(e=cmat), data=d)
anova(m0, m1p, m2)
Run the code above in your browser using DataLab