Ich arbeite an der grafischen Darstellung der vorhergesagten Werte aus einem Mehrebenenmodell (unter Verwendung des lme4-Pakets). Ich kann dies erfolgreich mit demEffect() Funktion. Wie nachfolgend dargestellt

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

Aber ich möchte diese Daten als einzelne Zeichnung mit einer Legende darstellen. Mit ggplots kann ich das machen; Aber ich verliere die Fehlerbalken, wie unten gezeigt:

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

Wie kann ich den ersten Plot (mit Fehlerbalken) als einzelnen Plot neu erstellen?