Eu tenho dados de 2 populações. Gostaria de obter o histograma e o gráfico de densidade de ambos no mesmo gráfico. Com uma cor para uma população e outra cor para a outra.

Eu tentei isso (exemplo):

library(ggplot2)
AA <- rnorm(100000, 70,20)
BB <- rnorm(100000,120,20)

valores <- c(AA,BB)
grupo <- c(rep("AA", 100000),c(rep("BB", 100000)))
todo <- data.frame(valores, grupo)


ggplot(todo, aes(x=valores, fill=grupo, color=grupo)) +
  geom_histogram(aes(y=..density..), binwidth=3)+  geom_density(aes(color=grupo))

Mas estou apenas obtendo um gráfico com uma única linha e uma única cor.

Eu gostaria de ter cores diferentes para as duas linhas de densidade. E se possível os histogramas também.

Eu fiz isso com o ggplot2, mas a base R também ficaria bem.

ou não sei o que mudei e agora entendi:

ggplot(todo, aes(x=valores, fill=grupo, color=grupo)) + 
geom_histogram( position="identity", binwidth=3, alpha=0.5)+ 
geom_density(aes(color=grupo))

mas as linhas de densidade não foram plotadas.

ou até coisas estranhas como