zoukankan      html  css  js  c++  java
  • KEGG富集分析散点图.md

    输入数据格式

    pathway = read.table("kegg.result",header=T,sep="	")
     
    pp = ggplot(pathway,aes(richFactor,Pathway)) #Pathwy是ID,richFactor是富集的基因数目除以背景的基因数目
    # 改变点的大小
    pp + geom_point(aes(size=R0vsR3)) # 以基因的数目表示点大小
     
    pbubble = pp + geom_point(aes(size=R0vsR3,color=-1*log10(Qvalue))) # 显著性表示颜色
    # 自定义渐变颜色
    pbubble + scale_colour_gradient(low="green",high="red")
     
    # 绘制pathway富集散点图
    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")
    # 改变图片的样式(主题)去除背景色
    pr + theme_bw()
    #去除网格线
    p_remove_grid <- pr +
    theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank())
    #网格线颜色
    panel.grid=element_line(color='grey')
    
    首先需要安装ggplot2
    library(ggplot2) #导入ggplot2
    x <- read.table("c:/Users/yueyao/Desktop/pathwayenrichment.txt",head = T, sep = "	") #读入文件,我的文本文件在桌面
    pdf(file="c:/Users/yueyao/Desktop/pathway_enrichment.pdf",width=10,height=10)#生成输出文件,双引号里面为路径及文件名,可自行设置
    png(file="c:/Users/yueyao/Desktop/pathway_enrichment.png",width=800,height=800)
    p <- ggplot(x,aes(x$Rich.Factor,x$Pathway))#作图利用的两列数据
    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")
    map
    dev.off()
    

    输出图片

  • 相关阅读:
    HTML常用标签(自用,可能不严谨,勿怪)
    Nginx负载均衡和反向代理设置
    Django的列表反序
    Python装饰器通用样式
    WCF、Web API、WCF REST、Web Service的区别
    C++11 标准新特性: 右值引用与转移语义
    在windows下vs使用pthread
    部分浏览器记住密码后可能会带来的问题
    SQL Server、 My SQL、PG Sql、Oracle、 Access 不同数据库sql差异
    sql中select语句的逻辑执行顺序
  • 原文地址:https://www.cnblogs.com/raisok/p/11014006.html
Copyright © 2011-2022 走看看