如何将高斯曲线添加到使用 qplot 创建的直方图？

Question

我的问题可能类似于 将密度曲线拟合到直方图R。使用 qplot，我使用此命令创建了 7 个直方图：

 (qplot(V1, data=data, binwidth=10, facets=V2~.)   

对于每个切片，我想添加一条拟合高斯曲线。当我尝试使用 lines() 方法时，出现错误：

Error in plot.xy(xy.coords(x, y), type = type, ...) : 
plot.new has not been called yet

正确执行此操作的命令是什么？

Answer 1

您尝试过stat_function吗？

+ stat_function(fun = dnorm)

您可能需要使用 aes(y = ..密度..) 绘制直方图，以便绘制密度值而不是计数。

许多有用的信息可以在这个中找到 问题，包括在不同方面绘制不同正态曲线的一些建议。

以下是一些示例：

dat <- data.frame(x = c(rnorm(100),rnorm(100,2,0.5)), 
                  a = rep(letters[1:2],each = 100))

在每个面上覆盖单个法线密度：

ggplot(data = dat,aes(x = x)) + 
  facet_wrap(~a) + 
    geom_histogram(aes(y = ..density..)) + 
    stat_function(fun = dnorm, colour = "red")

从我链接到的问题中，创建具有不同正态曲线的单独数据框：

grid <- with(dat, seq(min(x), max(x), length = 100))
normaldens <- ddply(dat, "a", function(df) {
  data.frame( 
    predicted = grid,
    density = dnorm(grid, mean(df$x), sd(df$x))
  )
})

并使用 geom_line 分别绘制它们：

ggplot(data = dat,aes(x = x)) + 
    facet_wrap(~a) + 
    geom_histogram(aes(y = ..density..)) + 
    geom_line(data = normaldens, aes(x = predicted, y = density), colour = "red")

Answer 2

ggplot2 使用与基础图形不同的图形范例。 （虽然您可以使用 grid 图形，但最好的方法是向绘图添加新的 stat_function 图层。ggplot2 代码是 请注意

，我无法使用 qplot 来实现此功能，但到 ggplot 的转换相当简单，最重要的区别是您的数据必须。 采用 data.frame 格式。

另请注意 y 美学 aes=aes(y=..密度..)) 的显式映射 - 这有点不寻常，但会获取 stat_function 结果并将其映射到数据：

library(ggplot2)
data <- data.frame(V1 <- rnorm(700), V2=sample(LETTERS[1:7], 700, replace=TRUE))
ggplot(data, aes(x=V1)) + 
  stat_bin(aes(y=..density..)) + 
  stat_function(fun=dnorm) + 
  facet_grid(V2~.)