[Rd] Error strata in aov (PR#1315)
jzacks@artsci.wustl.edu
jzacks@artsci.wustl.edu
Sun, 17 Feb 2002 21:39:11 +0100 (MET)
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._