Dear list users,
I am trying to learn Repeated measures ANOVA using the aov() interface, but
I'm struggling to understand its output.
According to tutorials on the web, formula for a repeated measures design
is:
aov(Y ~ IV+ Error(SUBJECT/IV) )
This formula does work but it returns three strata (Error:SUBJECT, Error:
SUBJECT:IV, Error: Within), when I would expect two strata (Within and
Between subjects). I've seems some tutorials show the exactly same setup,
but returning only the two first strata.
Is it possible to have two or three strata depending on the data?
If there is always three strata, how this would fit the interpretation of
between vs within effects?
Below a reproducible example that gives three strata:
data(beavers)
data=data.frame(id =
rep(c("beaver1","beaver2"),c(nrow(beaver1),nrow(beaver2))),rbind(beaver1,beaver2))
data$activ=factor(data$activ)
#balance dataset to have 6 samples for every combination of beaver and
activity.
balanced = split(data,interaction(data$id,data$activ))
sizes = sapply(balanced,nrow)
selected = lapply(sizes,sample.int,6)
balanced = mapply(function(x,y) {x[y,]}, balanced,selected,SIMPLIFY=F)
balanced = do.call(rbind,balanced)
aov(temp~activ+Error(id/activ),data=balanced)
Thanks,
Jorge
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
