在转录组学分析中,识别差异基因后,根据差异基因进一步注释到其富集的通路对深入挖掘差异基因的生物学信息具有重要的意义。差异基因比对不同的数据库富集到不同类型通路,其中重要的是KEGG与Gene_Ontology两个数据库,在上一篇中我们分享可制作kegg气泡图的教程。今天和大家分享的是GO条形图的R语言实现
#设置工作目录,注意R语言识别的路径与传统意义上的电脑储存路径写法不一致,直接“属性”复制的时候需要注意转换“\”成“/”。
setwd("E:/....../GO_class/GO_class")
#加载R包
library(ggplot2)
#加载并查看数据
GO<-read.delim(choose.files(),header= TRUE,check.names=F)
head(GO)
#调整顺序
GO$type_order=factor(rev(as.integer(rownames(GO))),labels=rev(GO$GO_classify2))
#绘图(纵向)
ggplot(data=GO,aes(x=GO$type_order,y=
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)