Full_Name: Jeff Zacks Version: 1.4.0 OS: Mac OS X Submission from: (NULL) (209.16.201.87) Hi, I believe there may be a problem with how Error strata are handled in aov, for designs with multiple error strata involving multiple interactions. The following performs the expected repeated-measures ANOVA when called in S-Plus, gives an error when called in R. In this example, the design is a simple 2x2x2 repeated-measures design with 2 subjects. Hope this is helpful. Best, Jeff Zacks
dat <- read.table("foo", header=T)
dat
Subj within1 within2 within3 depvar 1 A 0 0 1 1 2 A 0 1 0 2 3 A 1 0 0 3 4 A 0 1 1 4 5 A 1 0 1 5 6 A 1 1 0 6 7 B 0 0 1 2 8 B 0 1 0 3 9 B 1 0 0 4 10 B 0 1 1 5 11 B 1 0 1 6 12 B 1 1 0 7
summary(aov(depvar ~ within1*within2*withn3 + Error(Subj:within1 +
Subj:within2 + Subj:within3 + Subj:within1:within2 + Subj:within1+within3 + Subj:within2:within3), data=dat)) Error in model.frame.default(formula = depvar ~ Subj:within1 + Subj:within2 + : object is not a matrix
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._