与密度曲线叠加直方图

Question

我正在尝试制作密度值的直方图，并将其与密度函数的曲线（而不是密度估计）叠加。

使用一个简单的标准正态示例，这里有一些数据：

x <- rnorm(1000)

我可以这样做：

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

但这给出了频率而不是密度的直方图的比例。使用 ..密度.. 我可以在直方图上获得正确的比例：

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

但是现在这给出了一个错误：

q + stat_function( fun = dnorm )

有什么我没有看到的吗？

另一个问题，有没有办法绘制函数的曲线，例如curve()，但不是作为图层？

Answer 1

干得好！

# 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

，传递观察到的平均值和标准差，并使用 ggplot 代替 qplot：

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

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

Ramnath 答案的一个更简单的替代方案 sstatic.net/uIE0B.png" rel="noreferrer">

Answer 3

使用ggplot2中的geom_密度()怎么样？像这样：

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

ggplot(df, aes(x)) +
  geom_histogram(aes(y=..density..)) +  # scale histogram y
  geom_density(col = "red")

这也适用于多模式分发，例如：

df <- data.frame(x = c(rnorm(1000, 2, 2), rnorm(1000, 12, 2), rnorm(500, -8, 2)))

ggplot(df, aes(x)) +
  geom_histogram(aes(y=..density..)) +  # scale histogram y
  geom_density(col = "red")

Answer 4

我正在尝试获取虹膜数据集。您应该能够在这些简单的代码中看到您需要的图表：

ker_graph <- ggplot(iris, aes(x = Sepal.Length)) + 
geom_histogram(aes(y = ..density..),
colour = 1, fill = "white") +
geom_density(lwd = 1.2,
linetype = 2,
colour = 2)