经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » R语言 » 查看文章
KEGG富集分析散点图.md
来源:cnblogs  作者:raisok  时间:2019/6/13 8:58:27  对本文有异议

输入数据格式

  1. pathway = read.table("kegg.result",header=T,sep="\t")
  2. pp = ggplot(pathway,aes(richFactor,Pathway)) #Pathwy是ID,richFactor是富集的基因数目除以背景的基因数目
  3. # 改变点的大小
  4. pp + geom_point(aes(size=R0vsR3)) # 以基因的数目表示点大小
  5. pbubble = pp + geom_point(aes(size=R0vsR3,color=-1*log10(Qvalue))) # 显著性表示颜色
  6. # 自定义渐变颜色
  7. pbubble + scale_colour_gradient(low="green",high="red")
  8. # 绘制pathway富集散点图
  9. pr = pbubble + scale_colour_gradient(low="green",high="red") + labs(color=expression(-log[10](Qvalue)),size="Gene number",x="Rich factor",y="Pathway name",title="Top20 of pathway enrichment")
  10. # 改变图片的样式(主题)去除背景色
  11. pr + theme_bw()
  12. #去除网格线
  13. p_remove_grid <- pr +
  14. theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())
  15. #网格线颜色
  16. panel.grid=element_line(color='grey')
  1. 首先需要安装ggplot2
  2. library(ggplot2) #导入ggplot2
  3. x <- read.table("c:/Users/yueyao/Desktop/pathwayenrichment.txt",head = T, sep = "\t") #读入文件,我的文本文件在桌面
  4. pdf(file="c:/Users/yueyao/Desktop/pathway_enrichment.pdf",width=10,height=10)#生成输出文件,双引号里面为路径及文件名,可自行设置
  5. png(file="c:/Users/yueyao/Desktop/pathway_enrichment.png",width=800,height=800)
  6. p <- ggplot(x,aes(x$Rich.Factor,x$Pathway))#作图利用的两列数据
  7. map = p + geom_point(aes(size=x$Genes,colour=x$Qvalue))+theme(axis.text=element_text(color='black'),axis.text.y=element_text(size=14),axis.text.x=element_text(size=14),panel.background=element_rect(fill='transparent'),panel.grid=element_line(color='grey'),panel.border=element_rect(fill='transparent',color='black'),axis.title=element_text(size=16)) +labs(color="Qvalue",size="Gene number",x="Rich factor",y="Pathway name",title="Top20 of pathway enrichment")
  8. map
  9. dev.off()

输出图片

原文链接:http://www.cnblogs.com/raisok/p/11014006.html

 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号