Pathway富集分析氣泡圖


  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") > pr=pr + theme_bw()

     


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM