[R] Post Hoc Groupings

[Jarrett Byrnes]

Indeed, the following works as well
On Oct 26, 2005, at 5:23 PM, P Ehlers wrote:

> fm1 <- aov(breaks ~ wool*tension, data = warpbreaks)
> TukeyHSD(fm1, c("wool","tension", "wool:tension"))

However, when working with my own dataset,  I get the following errors. 
  I have some inkling this may be due to a slightly unbalanced sample 
size, but am uncertain of this.

 > rich.aov <- aov(X.open ~ Dock*Slip, data=dock_2004_data)
 > TukeyHSD(rich.aov, c("Dock", "Slip"))

Error in rep.int(n, length(means)) : unimplemented type 'NULL' in 'rep'
In addition: Warning messages:
1: non-factors ignored: Slip in: replications(paste("~", xx), data = mf)
2: non-factors ignored: Dock, Slip in: replications(paste("~", xx), 
data = mf)

I am pleased to know that these errors are not quite the