Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-04-09 Thread RICHARD M. HEIBERGER
On Wed, Apr 7, 2010 at 9:25 PM, Eric Scott wrote: > Thank you for your reply. The WoodEnergy example helped a lot. I > understand now that it is inappropriate to make all pairwise comparisons > with an interaction present and better to make comparisons between levels of > one factor within a co

Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-04-07 Thread Eric Scott
Thank you for your reply. The WoodEnergy example helped a lot. I understand now that it is inappropriate to make all pairwise comparisons with an interaction present and better to make comparisons between levels of one factor within a constant level of the second factor. As I understand it, the s

Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-03-15 Thread RICHARD M. HEIBERGER
In addtition to the example I mentioned previously, demo("MMC.WoodEnergy-aov", "HH") Please also see demo("MMC.WoodEnergy", "HH") In this example, since anova(energy.aov.4), shows that the Wood factor and Stove:Wood interaction are significant, all possible pairwise comparisons of the 12 Stove:Wo

Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-03-15 Thread Eric Scott
Thanks for the example, but I'm still not sure from this example how to see the pairwise comparisons for the interaction. For example, if I have two factors, X and Y; and X has 2 levels, A and B, and Y has 3 levels, 1, 2, and 3, a TukeyHSD would give the following comparisons with p-values for eac

Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-03-15 Thread RICHARD M. HEIBERGER
Please see the maiz example in ?MMC in the HH package. maiz is the last example in the help file. Keep going all the way to the end of the help file. See also the demo("MMC.WoodEnergy-aov", "HH") These examples show how to use glht in the presence of interactions and covariates. Rich

Re: [R] Multiple comparisons for a two-factor ANCOVA

2010-03-15 Thread Dieter Menne
Eric Scott-3 wrote: > > I'm trying to do an ANCOVA with two factors (clipping treatment with two > levels, and plot with 4 levels) and a covariate (stem diameter). The > response variable is fruit number. The minimal adequate model looks like > this: > > model3<-lm(fruit~clip + plot + st.dia

[R] Multiple comparisons for a two-factor ANCOVA

2010-03-15 Thread Eric Scott
I'm trying to do an ANCOVA with two factors (clipping treatment with two levels, and plot with 4 levels) and a covariate (stem diameter). The response variable is fruit number. The minimal adequate model looks like this: model3<-lm(fruit~clip + plot + st.dia + clip:plot) I'd like to get some mu