Wie füllt man ein Histogramm mit Farbverlauf?

Question

Ich habe ein einfaches Problem. Wie man Histogramm mit ggplot2 mit örtlich festgelegtem binwidth plant und mit Regenbogenfarben (oder irgendeiner anderen Palette) füllte?

Können sagen, ich habe eine Daten wie folgt aus:

myData <- abs(rnorm(1000)) 

Ich möchte Histogramm darstellen, zum Beispiel unter Verwendung von binwidth=.1. Das wird jedoch unterschiedliche Anzahl von Behältern verursachen, abhängig von Daten:

ggplot() + geom_histogram(aes(x = myData), binwidth=.1) 

Wenn ich Anzahl der Bins wusste (zB n=15) ich so etwas wie verwenden würde:

ggplot() + geom_histogram(aes(x = myData), binwidth=.1, fill=rainbow(n)) 

Aber mit wechselnder Anzahl von Behältern bin ich irgendwie auf diesem einfachen Problem fest.

Answer 1

Wenn Sie wirklich die Anzahl der Bins flexibel wollen, hier ist meine kleine Abhilfe:

library(ggplot2) 

gg_b <- ggplot_build(
    ggplot() + geom_histogram(aes(x = myData), binwidth=.1) 
) 

nu_bins <- dim(gg_b$data[[1]])[1] 

ggplot() + geom_histogram(aes(x = myData), binwidth=.1, fill = rainbow(nu_bins)) 

Answer 2

Bei der binwidth fixiert ist, ist hier eine alternative Lösung, die die interne Funktion verwendet ggplot2:::bin_breaks_width(), um die Anzahl der Bins vor Erstellen der Grafik zu erhalten. Es ist immer noch eine Abhilfe, sondern vermeidet geom_histogram() doppelt so in the other solution zu nennen:

# create sample data 
set.seed(1L) 
myData <- abs(rnorm(1000)) 
binwidth <- 0.1 

# create plot  
library(ggplot2) # CRAN version 2.2.1 used 
n_bins <- length(ggplot2:::bin_breaks_width(range(myData), width = binwidth)$breaks) - 1L 
ggplot() + geom_histogram(aes(x = myData), binwidth = binwidth, fill = rainbow(n_bins)) 

---

Als dritte Alternative kann die Aggregation außerhalb von ggplot2 erfolgen. Dann werden geom_col() Nocken statt geom_histogram() verwendet:

# start binning on multiple of binwidth 
start_bin <- binwidth * floor(min(myData)/binwidth) 
# compute breaks and bin the data 
breaks <- seq(start_bin, max(myData) + binwidth, by = binwidth) 
myData2 <- cut(sort(myData), breaks = breaks, by = binwidth) 

ggplot() + geom_col(aes(x = head(breaks, -1L), 
         y = as.integer(table(myData2)), 
         fill = levels(myData2))) + 
    ylab("count") + xlab("myData") 

Beachten Sie, dass breaks auf der x-Achse aufgetragen wird statt levels(myData2) die x-Achse kontinuierlich zu halten. Andernfalls würde jedes Faktoretikett geplottet werden, das die x-Achse durcheinander bringen würde. Beachten Sie auch, dass anstelle der rainbow() die integrierte Farbpalette ggplot2 verwendet wird.