lenda para estatísticas resumidas em ggplot2

Aqui está o código para o enredo

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)

Eu quero ter uma legenda para este enredo que irá identificar os valores de dados e valores médios alguma coisa assim

Black point = Data
Red point   = Mean.

Qualquer ponteiro para obter o resultado desejado será altamente apreciado. obrigado

questionAnswers(2)

yourAnswerToTheQuestion