Início / coding / Perguntas / 78941527
Accepted
Bogaso
Asked: 2024-09-03 01:25:50 +0800 CST

Adicionando gráfico de densidade normal ao histograma no ggplot

  • 772

Esta é uma pergunta complementar da minha postagem anterior em Plotando 2 histogramas nas partes superior e inferior de um único gráfico

Desenhei abaixo o par de histogramas

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))

Com isso, consegui desenhar normal density ploto primeiro histograma, mas não consegui sobrepor a curva de densidade normal para o segundo histograma.

Qualquer sugestão de como posso desenhar o segundo normal density plotpara o segundo histograma seria ótima.

Obrigado pelo seu tempo.

1 respostas

  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

relate perguntas