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 ae...@duke.edu ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.