Scale geom_density zu geom_bar mit Prozent auf y

Question

Da ich über die Mathematik last time I tried asking this verwirrt war, hier ist ein weiterer Versuch. Ich möchte ein Histogramm mit einer geglätteten Verteilungsanpassung kombinieren. Und ich möchte, dass die y-Achse in Prozent ist.

Ich kann keine gute Weg finden, um dieses Ergebnis zu erhalten. Letztes Mal habe ich einen Weg gefunden, die geom_bar auf den gleichen Maßstab wie geom_density zu skalieren, aber das ist das Gegenteil von dem, was ich wollte.

Meine aktuellen Code erzeugt diese Ausgabe:

ggplot2::ggplot(iris, aes(Sepal.Length)) + 
    geom_bar(stat="bin", aes(y=..density..)) + 
    geom_density() 

Die Dichte und die Bar y-Werte übereinstimmen, aber die Skalierung ist unsinnig. Ich will Prozent auf den y-Achsen, nicht gut, die Dichte.

Einige neue Versuche. Wir beginnen mit einem Balkendiagramm geändert Prozentsätze zu zeigen, statt zählt:

gg = ggplot2::ggplot(iris, aes(Sepal.Length)) + 
    geom_bar(aes(y = ..count../sum(..count..))) + 
    scale_y_continuous(name = "%", labels=scales::percent) 

Dann versuchen wir, einen geom_density hinzufügen, dass es irgendwie richtig zu skalieren erhalten:

gg + geom_density() 

gg + geom_density(aes(y=..count..)) 

gg + geom_density(aes(y=..scaled..)) 

gg + geom_density(aes(y=..density..)) 

Gleiche wie die erste.

gg + geom_density(aes(y = ..count../sum(..count..))) 

gg + geom_density(aes(y = ..count../n)) 

scheint off-Faktor von etwa 10 zu sein ...

gg + geom_density(aes(y = ..count../n/10)) 

gleiche wie:

gg + geom_density(aes(y = ..density../10)) 

Aber Ad-hoc-Zahlen eingefügt scheint wie eine schlechte Idee.

Ein nützlicher Trick besteht darin, die berechneten Werte des Diagramms zu überprüfen. Diese werden normalerweise nicht im Objekt gespeichert, wenn eines gespeichert wird. Allerdings kann man verwenden:

gg_data = ggplot_build(gg + geom_density()) 
gg_data$data[[2]] %>% View 

Da wir die Dichte wissen passen um x = 6 sollte etwa 0,04 (4%), können wir um für ggplot2 berechneten Werte suchen, die uns dorthin zu gelangen, und die einzige Was ich sehe, ist Dichte/10.

Wie bekomme ich geom_density passen auf die gleiche y-Achse skalieren wie die modifizierte geom_bar?

Bonus Frage: Warum sind die Gruppierungen der Bars unterschiedlich? Die aktuelle Funktion hat keine Leerzeichen zwischen Balken.

Answer 1

ist eine einfache Lösung:

library(scales) # ! important 
library(ggplot2) 
ggplot(iris, aes(Sepal.Length)) + 
stat_bin(aes(y=..density..), breaks = seq(min(iris$Sepal.Length), max(iris$Sepal.Length), by = .1), color="white") + 
geom_line(stat="density", size = 1) + 
scale_y_continuous(labels = percent, name = "percent") + 
theme_classic() 

Plot is here

Answer 2

Versuchen Sie, diese

ggplot2::ggplot(iris, aes(x=Sepal.Length)) + 
geom_histogram(stat="bin", binwidth = .1, aes(y=..density..)) + 
geom_density()+ 
scale_y_continuous(breaks = c(0, .1, .2,.3,.4,.5,.6), 
     labels =c ("0", "1%", "2%", "3%", "4%", "5%", "6%")) + 
ylab("Percent of Irises") + 
xlab("Sepal Length in Bins of .1 cm") 

Ich denke, Ihr erstes Beispiel ist das, was Sie wollen, Sie wollen einfach nur die Etiketten ändern, um es scheint, wie es Prozent, also genau das tun und nicht herumspielen. Hier