This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Add annotation color bar to ggplot objects

there is a picture below, I want to use my data to create a similar plot like it, but I don't know how to add annotation color bar for it! color bar in the red circle

r ggplot
• 22,843 views • link
updated by 董泽 • 董泽 • written by Yun Yun

Did you manage to plot other bits already? Share the code/errors/warnings.

Please provide example data, link to the paper, so we know what this plot is about.

• 0 views • link
written by zx8754 zx8754

Sorry for my blurred question, the picture is from paper https://www.ncbi.nlm.nih.gov/pubmed/29420467

there is a similar question here Add annotation color bar to ggplot or ggvis barplot, but the answer is for a heatmap object, I want a way to add a annotation color bar in a ggplot object

I have write some code like this

test <- data.frame(sample = paste0("sample", 1:100), tumor = rep(c("breast", "lung", "bladder", "colorectal", "prostate"), 20), allele = rep(c("exon1", "exon2", "exon3", "exon4", "exon5", "exon6", "exon7", "exon8", "exon9", "exon10"), 10), values = rnorm(100))
ggplot(test)+geom_bar(mapping = aes(x = sample, y = values, fill = allele), stat = "identity")

and now I want to use tumor variable in data.frame test to create a annotation color bar in the top of the picture

• 0 views • link
updated by Kevin Blighe Kevin Blighe • written by Yun Yun

Looking at that figure, my first guess was that it was generated via ComplexHeatmap, or it could be just a bunch of separate figures pieced together in a third-party software. That's probably not what you wanted to hear!

• 0 views • link
written by Kevin Blighe Kevin Blighe

Thanks for your reply, I'm read the reference manual for Complexheatmap to see if it can help me to make such figure. I have heard that I can use cowplot package to make it, Can it be a workable solution?

• 0 views • link
written by Yun Yun

I've been on this journey - you will ultimately end up at ComplexHeatmap, because whatever cowplot can do here, ComplexHeatmap can do it in a more straightforward manner.

• 0 views • link
written by Ram Ram

It definitely is ComplexHeatmap (looks like a bunch of Heatmaps with rowAnnotations added in) or a bunch of regular plots combined using third party software.

• 0 views • link
written by Ram Ram

Thanks you very much. I'll study the package ComplexHeatmap

• 0 views • link
written by Yun Yun

why don't you ask the author? my personal impression is that it is a combined plot.

• 0 views • link
written by H.Hasani H.Hasani

I don't have the contact of the author, I'll try the ComplexHeatmap, Thanks!!!

• 0 views • link
written by Yun Yun

You do - it's public information. All publications have a "corresponding author", who is the generally the person to reach out to with questions about the publication. In this case, it is David Hyman (hymand[at]mskcc[dot]org). See the picture below:

• 0 views • link
written by Ram Ram

Oh, Thanks, I'll try the package ComplexHeatmap firstly. Thanks for your help

• 0 views • link
written by Yun Yun

Although this plot has been depicted, I find the plots created by above code have some color which insinuate two different values. It is difficult for us to set complete different color aesthethic in two different plots ,especially for discrete variables.

is there any palette for this circumstance?

• 0 views • link
written by Yun Yun

It's probably best to re-phrase this as a new question. Include the code for the plot and use the image itself to point out which colors you'd like to display differently. You can, of course, also link to this post to keep the context.

• 0 views • link
written by Friederike Friederike

1 answer

Thanks for the people helping me in completing this plot. I have created this plot use two different ways, the one uses ggplot2 and cowplot, the other uses ComplexHeatmap

here is my code

for the data code

 test <- data.frame(sample = paste0("sample", 1:100), 
               tumor = rep(c("breast", "lung", "bladder", "colorectal", "prostate"), 20), 
               allele = rep(c("exon1", "exon2", "exon3", "exon4", "exon5", "exon6", "exon7", "exon8", "exon9", "exon10"),10), 
               values = rnorm(100))
 test <- test[order(test$tumor),]
 test$sample <- factor(test$sample, levels = test$sample)

for the ggplot2 and cowplot code:

h1 <- ggplot(test)+
       geom_bar(mapping = aes(x = sample, y = values, fill = allele), 
       stat = "identity", 
       position = "dodge")+
       guides(fill = guide_legend(ncol = 2))+
       theme(axis.text.x.bottom = element_blank(), 
    axis.ticks = element_blank(),
    panel.spacing.x = unit(1, "mm"),
    axis.title.x = element_blank(),
    strip.background.x = element_blank(),
    strip.text.x = element_blank())+
    facet_grid(.~tumor, scales = "free_x")
    h2 <- ggplot(test)+
  geom_bar(mapping = aes(x = sample, y = 1, fill = tumor), 
       stat = "identity", 
       width = 1)+
      theme_void()+
      theme(panel.spacing.x = unit(1, "mm"))+
     facet_grid(.~tumor, scales = "free_x")
     legend <- plot_grid(get_legend(h2), get_legend(h1), ncol = 1)
    h1 <- h1 + theme(legend.position = "none")
    h2 <- h2 + theme(legend.position = "none")
   plot <- plot_grid(h2, h1, align = "v", ncol = 1, axis = "tb", rel_heights = c(0.5, 15))
   plot_grid(plot, legend, nrow = 1, rel_widths = c(10, 1.5))

for the ComplexHeatmap code:

allele_col = structure(1:10, names = paste0("exon", 1:10))
    ha_list <- HeatmapAnnotation(bar = anno_barplot(test$values,
gp = gpar(fill = allele_col[test$allele]),
height = unit(6, "cm"))) %v%
   HeatmapAnnotation(tumor = anno_simple(test$tumor))
  draw(ha_list, heatmap_legend_list = Legend(title = "Allele", labels = names(allele_col), legend_gp = gpar(fill = allele_col)))
• 0 views • link
written by Yun Yun

This method still can't tackle the problem. The length of the "tumor type" label is not equal

• 0 views • link
written by 董泽 • 董泽

Log in to answer this question.

More posts like this