Para fazer comparações múltiplas no R, há várias funções disponíveis. No sistema básico, há a função TukeyHSD(), que faz as comparações múltiplas de Tukey. Para fazer comparações por Bonferroni e Scheffé, conforme descrito em Neter et al (1994), veja Multiple comparisons JIMRC.
Um pacote...