Barplot con diferencias significativas e interacciones?

Me gustaría visualizar mis datos y estadísticas ANOVA. Es común hacer esto utilizando una gráfica de barras con líneas agregadas que indican diferencias e interacciones significativas. ¿Cómo haces una trama como esta usando R?

Esto es lo que me gustaría:

Diferencias significativas:

Interacciones significativas:

Fondo

Actualmente estoy usandobarplot2{ggplots} para trazar barras e intervalos de confianza, pero estoy dispuesto a usar cualquier paquete / procedimiento para hacer el trabajo. Para obtener las estadísticas que estoy usando actualmenteTukeyHSD{stats} opairwise.t.test{stats} para las diferencias y una de las funciones anova (aov, ezANOVA{ez}, gls{nlme}) para interacciones.

Solo para darte una idea, esta es mi trama actual:

Respuestas a la pregunta(3)

Su respuesta a la pregunta