Re: [R] Contrast SS

2008-11-28 Thread Ista Zahn
Thank you! > m.1 <- aov(y ~ x) > summary(m.1, split = list(x = list("groups 1 and 2 vs group 3" = 1, "group 1 > vs group 2" = 2))) does what I want. On Fri, Nov 28, 2008 at 12:02 PM, Richard M. Heiberger <[EMAIL PROTECTED]> wrote: > I believe you are looking for the summary.lm() function and i

Re: [R] Contrast SS

2008-11-28 Thread Richard M. Heiberger
I believe you are looking for the summary.lm() function and its split= argument. Look at example(summary.aov) for several examples. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:/

Re: [R] Contrast SS

2008-11-28 Thread Gabor Grothendieck
Is your objective to break down the sum of squares for x into one for the group 1 vs group 3 contrast and another for the group 2 vs group 3 contrast? If that's it try this. We create the model matrix, make the two columns generated from x into separate explanatory variables, x1 and x2, and then

[R] Contrast SS

2008-11-28 Thread Ista Zahn
Hi all, Hopefully this is an easy question and I've overlooked something simple. I want to present specific contrasts in an ANOVA table. I can of course get the significance tests with summary.lm, but my adviser wants the results in ANOVA format. Here is a minimal example: > x <- factor(c(rep("con