用R画GO注释二级分类统计图

原文
来自生信菜鸟团，又添加了三个配色！

## load data
load(url("https://raw.githubusercontent.com/x2yline/Rdata/master/data%20visualization%20R/go_enrich_df.Rdata"))
head(go_enrich_df)
## numbers as data on x axis
go_enrich_df$number <- factor(rev(1:nrow(go_enrich_df)))
## shorten the names of GO terms
shorten_names <- function(x, n_word=4) {
  if (length(strsplit(x, " ")[[1]]) > n_word)
  {
    return(paste(paste(strsplit(x, " ")[[1]][1:n_word],
                       collapse=" "), "...", sep=""))
  } 
  else
  {
    return(x)
  }
}

labels=(sapply(
  levels(go_enrich_df$Description)[as.numeric(go_enrich_df$Description)],
  shorten_names))
names(labels) = rev(1:nrow(go_enrich_df))

## colors for bar
CPCOLS <- c("#66C3A5", "#8DA1CB", "#FD8D62")
#CPCOLS<-c("#CC6666", "#9999CC", "#66CC99")
#CPCOLS<-c("#999999", "#E69F00", "#56B4E9")
#CPCOLS<-alpha(c("blue", "red", "yellow"), .8)
library(ggplot2)
ggplot(data=go_enrich_df, aes(x=number, y=GeneNumber, fill=type)) +
  geom_bar(stat="identity", width=0.8) + coord_flip() + 
  scale_fill_manual(values = CPCOLS) + theme_bw() + 
  scale_x_discrete(labels=labels) +
  xlab("GO term") + 
  theme(axis.text=element_text(face = "bold", color="gray50")) +
  labs(title = "The Most Enriched GO Terms")

image.1
image.2 image.3 image.4

备注：

#CPCOLS<-c("#CC6666", "#9999CC", "#66CC99")
#CPCOLS<-c("#999999", "#E69F00", "#56B4E9")
#CPCOLS<-alpha(c("blue", "red", "yellow"), .8)
## 为新添加的配色， alpha选择透明度。

参考
http://www.biotrainee.com/thread-242-1-1.html
http://www.cookbook-r.com/Graphs/Colors_(ggplot2)/
https://zhuanlan.zhihu.com/p/27093478
https://zhuanlan.zhihu.com/p/25173606