R ggplot с двумя сериями: точки и ошибки с легендами
Если у меня есть такой фрейм данных:
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
И я хочу, чтобы столбцы ошибок наблюдения и прогнозирования модели выглядели примерно так:
(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))
)
Как добавить легенду / шкалу / ключ, определяющий красные точки в качестве наблюдений и синие плюсы с барами ошибок в качестве точечных прогнозов с диапазонами?