Criar um histograma de densidade em ggplot2?

Question

Eu quero criar o seguinte histograma de densidade de terreno com ggplot2.No "normal" forma (da base de dados de pacotes) é muito fácil:

set.seed(46)
vector <- rnorm(500)  
breaks <- quantile(vector,seq(0,1,by=0.1))
labels = 1:(length(breaks)-1)
den = density(vector)
hist(df$vector,
     breaks=breaks,
     col=rainbow(length(breaks)),
     probability=TRUE)
lines(den)

Com ggplot eu ter chegado a este tão longe:

seg <- cut(vector,breaks,
           labels=labels,
           include.lowest = TRUE, right = TRUE)
df = data.frame(vector=vector,seg=seg)

ggplot(df) + 
     geom_histogram(breaks=breaks,
                    aes(x=vector,
                        y=..density..,
                        fill=seg)) + 
     geom_density(aes(x=vector,
                      y=..density..))

Mas o "y" escala tem o errado dimensão.Tenho observado que a próxima execução fica com o "y" escala de direito.

 ggplot(df) + 
     geom_histogram(breaks=breaks,
                    aes(x=vector,
                    y=..density..,
                    fill=seg)) + 
     geom_density(aes(x=vector,
                      y=..density..))

Eu só não entendo isso. y=..density.. é lá, que deve ser a altura.Então, por que na terra minha escala é modificado quando eu tento preenchê-lo?

Eu preciso das cores.Eu só quero um histograma onde as quebras e as cores de cada bloco são unilateralmente definidas de acordo com o padrão ggplot de preenchimento de cores.

Answer 1

Manualmente, eu adicionei cores para o seu percentil bares.Veja se isso funciona para você.

library(ggplot2)

ggplot(df, aes(x=vector)) +   
   geom_histogram(breaks=breaks,aes(y=..density..),colour="black",fill=c("red","orange","yellow","lightgreen","green","darkgreen","blue","darkblue","purple","pink")) + 
   geom_density(aes(y=..density..)) +
   scale_x_continuous(breaks=c(-3,-2,-1,0,1,2,3)) +
   ylab("Density") + xlab("df$vector") + ggtitle("Histogram of df$vector") +
   theme_bw() + theme(plot.title=element_text(size=20),
                      axis.title.y=element_text(size = 16, vjust=+0.2),
                      axis.title.x=element_text(size = 16, vjust=-0.2),
                      axis.text.y=element_text(size = 14),
                      axis.text.x=element_text(size = 14),
                      panel.grid.major = element_blank(),
                      panel.grid.minor = element_blank())

Answer 2

fill=seg resultados no agrupamento.Na verdade, você está recebendo uma diferente histograma para cada valor de seg.Se você não precisa das cores, você poderia usar:

ggplot(df) + 
  geom_histogram(breaks=breaks,aes(x=vector,y=..density..), position="identity") + 
  geom_density(aes(x=vector,y=..density..))

Se você precisa de cores, pode ser mais fácil para calcular os valores de densidade de fora da ggplot2.

Answer 3

A confusão sobre a interpretação y-axis pode ser devido a densidade é plotado em vez de contar.Assim, os valores em y-axis são proporções do total da amostra, onde o soma das barras é igual a 1.