Prof Brian Ripley <[EMAIL PROTECTED]> writes: > > My questions is: > > which formula should I use if I want to use a random effects model? > > I think I should use "lme", but I don't know how. > > aov() is the function for anova modelling. so you want > > aov(A ~ B, data = DATA) # fixed effects > aov(A ~ Error(B), data = DATA) # random effects > > but a 1-way ANOVA with fixed or random effects is the same analysis, and > only the interpretation differs.
Not quite. Only if the design is balanced and you're not looking to compute things like the variance of the estimated overall mean. In the unbalanced case, I don't believe aov() is doing the right thing (it might not be very wrong, though, if the imbalance is slight) and lme/lmer is closer on the mark. -- O__ ---- Peter Dalgaard Øster Farimagsgade 5, Entr.B c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
