Histograma de superposición con curva de densidad

Question

Estoy tratando de hacer un histograma de valores de densidad y superponerlo con la curva de una función de densidad (no la densidad estimada).

El uso de un ejemplo sencillo normal estándar, aquí es algunos datos:

x <- rnorm(1000) 

que puedo hacer:

q <- qplot(x, geom="histogram") 
q + stat_function(fun = dnorm) 

pero esto da la escala del histograma en frecuencias y no densidades. con ..density.. puedo conseguir la escala apropiada en el histograma:

q <- qplot(x,..density.., geom="histogram") 
q 

Pero ahora esto da un error:

q + stat_function(fun = dnorm) 

¿Hay algo que no estoy viendo?

Otra pregunta, ¿hay alguna manera de trazar la curva de una función, como curve(), pero no como capa?

Answer 1

¡Aquí tienes!

# create some data to work with 
x = rnorm(1000); 

# overlay histogram, empirical density and normal density 
p0 = qplot(x, geom = 'blank') + 
    geom_line(aes(y = ..density.., colour = 'Empirical'), stat = 'density') + 
    stat_function(fun = dnorm, aes(colour = 'Normal')) +      
    geom_histogram(aes(y = ..density..), alpha = 0.4) +       
    scale_colour_manual(name = 'Density', values = c('red', 'blue')) + 
    theme(legend.position = c(0.85, 0.85)) 

print(p0) 

Answer 2

A más escueto alternativo a la respuesta de Ramnath, que pasa a la media observada y la desviación estándar, y el uso de ggplot vez de qplot:

df <- data.frame(x = rnorm(1000, 2, 2)) 

# overlay histogram and normal density 
ggplot(df, aes(x)) + 
    geom_histogram(aes(y = ..density..)) + 
    stat_function(fun = dnorm, 
       args = list(mean = mean(df$x), sd = sd(df$x)), 
       lwd = 2, 
       col = 'red')