主页 / coding / 问题 / 78941527
Accepted
Bogaso
Asked: 2024-09-03 01:25:50 +0800 CST

在 ggplot 中的直方图中添加正态密度图

  • 772

这是我上一篇文章《在单个图的顶部和底部绘制 2 个直方图》中的后续问题

我画了以下一对直方图

library(ggplot2)
set.seed(1)
dat = rbind(data.frame('val' = rnorm(100), 'met' = 'Metric1'), data.frame('val' = rt(100, 2), 'met' = 'Metric12'))

ggplot(dat, aes(x = val, fill = met)) + 
  geom_histogram(data    = dat[dat$met == 'Metric1',], 
                 breaks  = seq(-10, 10, 0.5),
                 mapping = aes(y = after_stat(density)),
                 colour  = "black",
                 alpha   = 0.3) +
  stat_function(fun = dnorm, args = list(mean = 0, sd = 1)) +
  geom_histogram(data    = dat[dat$met == 'Metric12', ],
                 breaks  = seq(-10, 10, 0.5),
                 mapping = aes(y = -after_stat(density)),
                 colour  = "black",
                 alpha   = 0.3)+
  coord_cartesian(xlim = c(-10, 10))

用这个我可以绘制normal density plot第一个直方图，但无法叠加第二个直方图的正常密度曲线。

关于如何绘制normal density plot第二个直方图的任何建议都很棒。

感谢您的时间。

1 个回答

  1. Best Answer
    ggplot() + 
  geom_histogram(data    = dat[dat$met == 'Metric1',], 
                 breaks  = seq(-10, 10, 0.5),
                 mapping = aes(x = val, y = after_stat(density)),
                 colour  = "black",
                 alpha   = 0.3) +
  geom_histogram(data    = dat[dat$met == 'Metric12', ],
                 breaks  = seq(-10, 10, 0.5),
                 mapping = aes(x = val, y = -after_stat(density)),
                 colour  = "black",
                 alpha   = 0.3)+
  coord_cartesian(xlim = c(-10, 10))+ 
  geom_line(
    aes(
      x = dat[dat$met == 'Metric1',]$val,
      y = dnorm(dat[dat$met == 'Metric1',]$val),
      color = as.factor(1))
  )+
  geom_line(
    aes(
      x = dat[dat$met == 'Metric12',]$val,
      y = -dnorm(dat[dat$met == 'Metric12',]$val),
      color = as.factor(2))
  )
    • 2

相关问题