I ran a split-plot anova using aov() and would like to do
pairwise comparisons for an interaction term using TukeyHSD().  However,
when I try this, I get an error.  Can anyone help?  Here is the code:

> attach(soilcom.data)
> soilcom.aov <- aov(Depth ~ Stand*Distance + Error(Site/Stand))
> summary(soilcom.aov)

Error: Site
          Df Sum Sq Mean Sq F value Pr(>F)
Residuals 28  29937  1069.2

Error: Site:Stand
          Df Sum Sq Mean Sq F value   Pr(>F)
Stand      1  32654   32654  132.53 3.94e-12 ***
Residuals 28   6899     246
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Error: Within
                Df  Sum Sq Mean Sq F value   Pr(>F)
Distance         4   531.8  132.95  1.2363 0.296330
Stand:Distance   4  1858.2  464.54  4.3196 0.002193 **
Residuals      224 24089.6  107.54
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

> TukeyHSD(soilcom.aov)

Error in UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied
to an object of class "c('aovlist', 'listof')"



Chris Habeck
Research Associate
Department of Biology
Washington University at St. Louis

	[[alternative HTML version deleted]]

