Wenn ich einen Datenrahmen wie diesen habe:

obs<-rnorm(20)
d<-data.frame(year=2000:2019,obs=obs,pred=obs+rnorm(20,.1))
d$pup<-d$pred+.5
d$plow<-d$pred-.5
d$obs[20]<-NA
d

Und ich möchte, dass die Fehlerbalken für Beobachtung und Modellvorhersage ungefähr so aussehen:

(p1<-ggplot(data=d)+aes(x=year)
 +geom_point(aes(y=obs),color='red',shape=19)
 +geom_point(aes(y=pred),color='blue',shape=3)
 +geom_errorbar(aes(ymin=plow,ymax=pup))
 )

Wie füge ich eine Legende / Skala / einen Schlüssel hinzu, die / der die roten Punkte als Beobachtungen und die blauen Punkte mit Fehlerbalken als Punktvorhersagen mit Bereichen identifiziert?