Ben Meijering
2008-Dec-20 15:01 UTC
[R] How test contrasts/coefficients of Repeated-Measures ANOVA?
Hi all, I'm doing a Repeated-Measures ANOVA, but I don't know how to test its contrasts or where to find the p-values of its coefficients. I know how to find the coefficient estimates of a contrast, but not how to test these estimates. First I do something like: y.aov <- aov(y ~ fac1 * fac2 + Error(subj/(fac1 * fac2)), data=data) Then, with coef(y.aov) I get the coefficients (which are, by default, treatment contrasts?), but p-values are not provided. How do I obtain these p-values? ~ Ben
Dieter Menne
2008-Dec-20 19:29 UTC
[R] How test contrasts/coefficients of Repeated-Measures ANOVA?
Ben Meijering <b.meijering <at> student.rug.nl> writes:> I'm doing a Repeated-Measures ANOVA, but I don't know how to test its > contrasts or where to find the p-values of its coefficients. I know > how to find the coefficient estimates of a contrast, but not how to > test these estimates. > > First I do something like: > > y.aov <- aov(y ~ fac1 * fac2 + Error(subj/(fac1 * fac2)), data=data) > > Then, with > > coef(y.aov) > > I get the coefficients (which are, by default, treatment contrasts?), > but p-values are not provided.See model.table and se.contrast; and consider using lme in package nlme as an alternative which make live much easier, if it can handle your (over?)complex nesting. Dieter