Hello ! Thanks for your advice. I tried it, but the output is the same: > HSD.test(anova.typabunmit, "typ", group=TRUE) Name: typ ds.typabunmit$typ I don't get the values...!?!? -- View this message in context: http://r.789695.n4.nabble.com/Anova-and-tukey-grouping-tp4644485p4644513.html Sent from the R help mailing list archive at Nabble.com.