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

Utilizando un ejemplo normal estándar simple, aquí hay algunos datos:

x <- rnorm(1000)

Puedo hacer

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

pero esto proporciona la escala del histograma en frecuencias y no en densidades. con..density.. Puedo obtener la escala adecuada 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, comocurve(), pero entonces no como capa?