Hello, I think that to apply the Neuman-Keuls test under R it's possible with the TukeyHSD function and I know that to do a Neuman-Keuls test I have to sort the means but I can't use the TukeyHSD function I don't understand how to do to apply this function to my data. In the R help they use the TukeyHSD function to a aov() do you think that I have to do teh same? thanks. ___________________________________________________________________________ [[alternative HTML version deleted]]