Oto kod fabuły

library(ggplot2)
df <- data.frame(gp = factor(rep(letters[1:3], each = 10)), y = rnorm(30))
library(plyr)
ds <- ddply(df, .(gp), summarise, mean = mean(y), sd = sd(y))
ggplot(df, aes(x = gp, y = y)) +
   geom_point() +
   geom_point(data = ds, aes(y = mean), colour = 'red', size = 3)

Chcę mieć legendę dla tego wykresu, która zidentyfikuje wartości danych i oznacza wartości takie jak ta

Black point = Data
Red point   = Mean.

Każdy wskaźnik, aby uzyskać pożądany wynik, zostanie bardzo doceniony. Dzięki