[R] TukeyHSD troubles

[Amy Freitag]

I've tried to run a Tukey post-hoc but keep getting this weird error,
whether the aov was significant or not. treat_code is a dummy
variable, but that shouldn't matter.  Any suggestions?
Thanks
Amy

> summary(aov(EtoH~treat_code, mydata))
             Df Sum Sq Mean Sq F value   Pr(>F)
treat_code    1  16.44   16.44  11.027 0.001014 **
Residuals   287 427.91    1.49
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> TukeyHSD(aov(EtoH~treat_code, mydata))
Error in TukeyHSD.aov(aov(EtoH ~ treat_code, mydata)) :
  no factors in the fitted model
In addition: Warning message:
In replications(paste("~", xx), data = mf) :
  non-factors ignored: treat_code


-- 
Amy Freitag
Ph.D. student, Marine Science and Conservation
Nicholas School for the Environment
Duke University
aef10 at duke.edu