通路富集分析气泡图(R实现)

  1. data.tsv 
  2. > pathway = read.table("data.tsv",header = T, sep="\t")
    > library(ggplot2) 
    > p = ggplot(pathway,aes(Pvalue,Pathway))
    > p=p + geom_point()  
    > p=p + geom_point(aes(size=Count))
    > pbubble = p+ geom_point(aes(size=Count,color=-1*log10(Pvalue)))
    > pr = pbubble+scale_color_gradient(low="green",high = "red")
    > pr = pr+labs(color=expression(-log[10](Pvalue)),size="Gene number",x="Pvalue",y="Pathway name",title="Pathway enrichment")+theme(plot.title=element_text(hjust=0.5))
    > pr=pr + theme_bw()

猜你喜欢

转载自blog.csdn.net/qq_42458954/article/details/87627765