R help - Oct 2009 - Multiple comparison test split plot analysis in R

Hi all!


I have a split plot analysis in which the model parameters are:



Efficacy = y

Rep = block

Product and rateind = main plot

Variety = subplot

I can do the split plot analysis with:

summary(model1<-aov(efficacy~rep+product*rateind*variety+Error(rep:product:rate),
data=ht28))


How do I multiple comparison tests of the main plot factor, subplot factor and
interaction? I cannot see how to do this  using TukeyHSD and use the correct
error term.

Thanks

Regards

Jo

Joanne Stringer
Senior Biometrician
BSES Limited
Postal Address: PO Box 86
Street Address: 50 Meiers Road
Indooroopilly QLD 4068 Australia
Ph: +61 (0) 7 33313337
Fax: +61 (0) 7 38710383
Mobile: +61 (0) 408011205
Email: jstringer@bses.org.au<blocked::mailto:jstringer@bses.org.au>
Web: www.bses.org.au

BSES Limited Disclaimer

--------------------------------------------------------------------------------

This email and any files transmitted with it are confide...{{dropped:15}}

Please see the maiz example in ?HH:::glht.mmc

You may need to install HH first with

install.packages("HH")