[R-br] Teste de Kruskal Wallis

Prezados, bom dia. Preciso fazer uma comparação entre4 tratamentos via teste de Kruskal Wallis. Utilizei o comando: kruskal.test(Resposta~trat, dados) Kruskal-Wallis rank sum test data: Resposta by trat Kruskal-Wallis chi-squared = 14.602, df = 3, p-value = 0.00219 indicando haver diferenças significativas. Como fazer comparações múltiplas entre os grupos? Desde já agradeço. Att., Ana Paula