Estoy trabajando en graficar los valores predichos de un modelo multinivel (usando el paquete lme4). Puedo hacer esto con éxito usando elEffect() función. Como se muestra abajo:

library(lme4)
library(effects)
m1=lmer(price~depth*cut+(1|cut),diamonds)
plot(Effect(c("cut","depth"),m1))

Pero, quiero presentar estos mismos datos como una sola trama con una leyenda. Usando ggplots, puedo hacer esto; pero pierdo las barras de error, como se muestra a continuación:

ggplot(data.frame(Effect(c("cut","depth"),m1)),
       aes(x=depth,y=fit,color=cut,group=cut))+
  geom_line()

¿Cómo puedo recrear la primera trama (con barras de error) como una sola trama?