Aquí está el código para la trama.

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)

Quiero tener una leyenda para este gráfico que identifique los valores de los datos y los valores medios algo como esto

Black point = Data
Red point   = Mean.

Cualquier puntero para obtener el resultado deseado será altamente apreciado. Gracias