在 ggplot2 Rstudio 中获取第二个绘图

Question

我希望分布 f 出现在我的情节中。

library(ggplot2)
f <- dnorm(x=2)
x <- c(v6$Zwischenankunftszeit)
counts <- c(v6$ZugangFolgeVrg)
df <- data.frame(x=x, counts=counts)

plt <- ggplot(df + f) + geom_bar(aes(x=x, y=counts), stat="identity")
print(plt)

Answer 1

如果你想在顶部绘制一个正态分布的直方图，这是一种更简单的方法：

library(ggplot2)

df <- data.frame(
  x = rnorm(100, mean = 50, sd = 5)
)

ggplot(df, aes(x)) +
  geom_histogram(aes(y = stat(density))) +
  stat_function(fun = dnorm, args = list(mean = mean(df$x), sd = sd(df$x)))
#> `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

^{由reprex package (v0.3.0) 于 2020 年 5 月 26 日创建}

如果要将 y 轴保持为计数单位，则必须包装 dnorm 函数并按 binwidth 和观察次数对其进行缩放：

wrap_dnorm <- function(x, mean = 0, sd = 1, binwidth = 1, count = 1) {
  dnorm(x, mean = mean, sd = sd) * binwidth * count
}

ggplot(df, aes(x)) +
  geom_histogram(binwidth = 1) +
  stat_function(fun = wrap_dnorm, 
                args = list(mean = mean(df$x), sd = sd(df$x), 
                            binwidth = 1, count = nrow(df)))