How to fill histogram with color gradient?

Question

I have a simple problem. How to plot histogram with ggplot2 with fixed binwidth and filled with rainbow colors (or any other palette)?

Lets say I have a data like that:

myData <- abs(rnorm(1000))

I want to plot histogram, using e.g. binwidth=.1. That however will cause different number of bins, depending on data:

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

If I knew number of bins (e.g. n=15) I'd use something like:

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

But with changing number of bins I'm kind of stuck on this simple problem.

So if I understand you correctly, you want each bin of the histogram to be coloured differently according to a rainbow gradient? — sebastian-c, Oct 27 '16 at 12:12
@user20650 I saw this answer, but over there number of bins is fixed so its not solving my problem. — Art, Oct 27 '16 at 12:15
Why not just force the number of bins through the `bins` argument? — Roman Luštrik, Oct 27 '16 at 12:22
@RomanLuštrik: Because I need to have fixed bin width not fixed number of bins. If I could do that I wouldn't ask :) — Art, Oct 27 '16 at 12:24
This seems to work... `n <- round(((max(myData)-min(myData))/.1)+1)` — Rupert, Oct 27 '16 at 12:26

score 8 · Accepted Answer · answered Oct 27 '16 at 12:24

8

If you really want the number of bins flexible, here is my little workaround:

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

answered Oct 27 '16 at 12:24

J_F

9,956
2
31
55

Yep, that is exactly what I need. I was just thinking (and hoping) that there is some more straightforward method to do that :) Thanks! – Art Oct 27 '16 at 12:29

Uwe · Answer 2 · 2017-06-16T22:16:15.237

In case the binwidth is fixed, here is an alternative solution which is using the internal function ggplot2:::bin_breaks_width() to get the number of bins before creating the graph. It's still a workaround but avoids to call geom_histogram() twice as in the other solution:

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

As a third alternative, the aggregation can be done outside of ggplot2. Then, geom_col() cam be used instead of geom_histogram():

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

Note that breaks is plotted on the x-axis instead of levels(myData2) to keep the x-axis continuous. Otherwise each factor label would be plotted which would clutter the x-axis. Also note that the built-in ggplot2 color palette is used instead of rainbow().

How to fill histogram with color gradient?

2 Answers2

Linked