Barplot com diferenças e interações significativas?

Gostaria de visualizar meus dados e estatísticas ANOVA. É comum fazer isso usando um gráfico de barras com linhas adicionadas indicando diferenças e interações significativas. Como você faz um enredo como esse usando R?

Isto é o que eu gostaria:

Diferenças significantes:

Interações significativas:

fundo

Eu estou usando atualmentebarplot2{ggplots} para traçar barras e intervalos de confiança, mas estou disposto a usar qualquer pacote / procedimento para fazer o trabalho. Para obter as estatísticas que estou usando no momentoTukeyHSD{stats} oupairwise.t.test{stats} diferenças e uma das funções anova (aov, ezANOVA{ez}, gls{nlme}) para interações.

Só para você ter uma ideia, esta é minha trama atual: