I have a between-subjects factor, A, with three levels, and a
within-subjects factor with 4 levels. There are 6 subjects at each A
level. I get the correct F tests using the aov function with the stack
form of the data frame but I want the epsilon adjustment
(Greenhouse-Geisser, Huynh-Feldt). for nonspheriicty as well. I got the
SPSS results using the standard broad data frame (a row for each subject
with two within-subject variables) in a simple repeated measures design
but get an error message when one factor is between subjects.. I tried
what I thought was a simple extension of this approach with one between-
and one within-subject factor. I tried
B = factor(c("B1","B2","B3","B4"))
A=factor(rep(c("A1","A2","A3"),each=6))
Bdf=data.frame(B)
attach(Bdf)
model.lm = lm(cbind(B1,B2,B3,B4) ~A,data=mix.data)
mix.model=Anova(model.lm, idata=Bdf, idesign=
B,Type="III",data=mix.data)
summary(mix.model,multivariate=FALSE)
This gives some output but it yields an error message:
"Error in x$terms : $ operator is invalid for atomic vectors
In addition: Warning message:
In Anova.mlm(model.lm, idata = Bdf, idesign = B, Type = "III") :
the model contains only an intercept: equivalent Type III test
substituted"
I've tried variations of the lm parameters (e.g., ~1 inserted) but it
still yields an error message. I'd appreciate suggestions.
Jerry Myers
--
Jerome L. Myers
Department of Psychology, Tobin Hall
University of Massachusetts
Box 37710
Amherst, MA 01003-7710
Phone: (413) 545-2331
Fax: (413) 545-0996
http://www.psych.umass.edu/people/jeromemyers/home/
______________________________________________
R-help@r-project.org mailing list
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.
