Thank you all for your **very good** answers: Using aovp(..., perm="Exact") seems to be the way to go for small datasets, and also I should definitely try ?kruskal.test. Juan [[alternative HTML version deleted]]