• 新疆铁路货运单日装车突破7000车大关 2019-04-09
  • 端午小长假杭州接待游客454.8万人次 9成民宿预订一空 2019-04-02
  • 溃疡恶变也会腹痛 这些肿瘤最容易被误诊 2019-03-27
  • 为藏药材可持续发展垦出“良田沃土” 2019-03-27
  • 吴尊送儿子上学晒自拍 父子颜值超高Max笑容灿烂 2019-03-24
  • 月薪3000男会计豪掷930万打赏 女主播该退钱吗? 2019-03-23
  • 我国成功组织航天员沙漠野外生存训练 2019-03-23
  • 网红健康险——红在哪 险在哪 2019-03-18
  • 升级是硬道理 MQB平台全新一代宝来驾乘体验更舒适 2019-03-18
  • 马特乌斯评德国队首战失利:“几乎什么都欠缺” 2019-03-02
  • 世界自然遗产旅行地中国国家地理网 2019-03-02
  • 计划经济走向完美,看搞好公有制计划经济存在不好说来看看?是统得太死还是要这样做才科学? 2019-02-17
  • 马思纯《我是大侦探》送礼物 吴磊:我们大老爷们好糙啊 2019-02-17
  • 回复@海之宁:不给你分配生产资料你咋尽所能? 2019-01-28
  • Conférence de presse du Premier ministre chinois 2019-01-28
  • 画图专题 | 差异表达基因结果展示系列之柱状图

    如何将好几张功能富集结果中的柱状图的横坐标的范围全部调整为一样的?

    一般画这个柱状图都是用Y叔的clusterprofiler包中的barplot函数对使用这个包的功能富集结果进行一键绘图,超级简单方便。但是当我去查找这个函数的调整坐标的参数时:

    Usage
    ##?S3?method?for?class?‘enrichResult’
    barplot(height,?x?=?“Count”,?colorBy?=?“pvalue”,
    ?showCategory?=?5,?font.size?=?12,?title?=?“”,?…)
    Arguments
    height?????????????enrichResult?object
    x?????????????????one?of?‘Count’?and?‘GeneRatio’
    colorBy???????????one?of?‘pvalue’,?‘p.adjust’,?‘qvalue’
    showCategory?????number?of?categories?to?show
    font.size??????????font?size
    title??????????????plot?title
    …???????????????other?parameter,?ignored

    other?parameter,?ignored,other?parameter,?ignored,other?parameter,?ignored……

    也就是说这个函数里面没有调整坐标范围的参数,不知道是不是因为我没找到还是什么,反正我找了好久没有找到,我有点方。

    后面我觉得还是自己用ggplot2画吧,反正这个包里面一键绘图也是根据ggplot2来的。

    这个图与一般的函数barplot画出来的不一样的地方在于它的颜色,这张图里面的颜色反应的是fdr的大小,是一个连续值,ggplot2可以将连续值映射到到颜色上,横坐标是通路中感兴趣基因的个数。

    用来画图的数据示例:

    #加载ggplot2包
    library(ggplot2)?

    #一键清除
    rm(list=ls())

    #设置文件路径
    setwd(“D:/zhangj/DZH/enrichment”)?

    #读取功能富集结
    enrich?<-?read.table(“S01_S03_S05_vs_S02_S04_S06_Biological_Process_enrich.list”,?header=T,sep=”\t”,stringsAsFactors=F,comment.char=””,quote=”\””)

    #对富集结果按照qvalue进行从小到大排序,保证最显著的通路在前
    enrich1?<-?enrich[order(enrich$qvalue),]?

    #这里画图只展示top10的通路
    enrich2?<-?enrich1[1:10,]

    #提取每条通路里面差异表达的基因数
    count?<-?as.numeric(unlist(strsplit(enrich2$GeneRatio,”/5498″,fixed=T)))?
    enrich3?<-?data.frame(enrich2[,2],count,enrich2[,7])
    colnames(enrich3)?<-?c(“ID”,”count”,”qvalue”)

    #fill=qvalue?fill颜色填充,使用连续值qvalue
    p?<-?ggplot(data=enrich3,aes(x=ID,y=count,fill=qvalue))

    #coord_flip()颠倒坐标轴
    p1?<-?p?+?geom_bar(stat=”identity”)?+?coord_flip()
    p2?<-?p1?+?theme(panel.background=element_rect(fill=’transparent’,color=’gray’),
    ???????????axis.text.y=element_text(color=”black”,size=12))

    #ylim(0,20)?更改横坐标的范围这里坐标轴颠倒了,虽然看起来是x轴,但其实是y轴
    p3?<-?p2?+?ylim(0,20)?+?scale_fill_gradient(low=”red”,high=”blue”)

    p4?<-?p3?+?scale_x_discrete(limits=rev(enrich3[,1]))?+labs(x=””,y=””,title=”Biological_Process”)

    #输出为png格式的图片
    png(“S01_S03_S05_vs_S02_S04_S06_Biological_Process_enrich.png”,width=680,height=480)
    print(p4)
    dev.off()

    #输出为pdf的文件
    pdf(“S01_S03_S05_vs_S02_S04_S06_Biological_Process_enrich.pdf”,width=9)
    print(p4)
    dev.off()

    最后结果图:

     

     

    最近文章
    qq游戏赢三张牌下载
  • 新疆铁路货运单日装车突破7000车大关 2019-04-09
  • 端午小长假杭州接待游客454.8万人次 9成民宿预订一空 2019-04-02
  • 溃疡恶变也会腹痛 这些肿瘤最容易被误诊 2019-03-27
  • 为藏药材可持续发展垦出“良田沃土” 2019-03-27
  • 吴尊送儿子上学晒自拍 父子颜值超高Max笑容灿烂 2019-03-24
  • 月薪3000男会计豪掷930万打赏 女主播该退钱吗? 2019-03-23
  • 我国成功组织航天员沙漠野外生存训练 2019-03-23
  • 网红健康险——红在哪 险在哪 2019-03-18
  • 升级是硬道理 MQB平台全新一代宝来驾乘体验更舒适 2019-03-18
  • 马特乌斯评德国队首战失利:“几乎什么都欠缺” 2019-03-02
  • 世界自然遗产旅行地中国国家地理网 2019-03-02
  • 计划经济走向完美,看搞好公有制计划经济存在不好说来看看?是统得太死还是要这样做才科学? 2019-02-17
  • 马思纯《我是大侦探》送礼物 吴磊:我们大老爷们好糙啊 2019-02-17
  • 回复@海之宁:不给你分配生产资料你咋尽所能? 2019-01-28
  • Conférence de presse du Premier ministre chinois 2019-01-28
  • 中国体彩超级大乐透 三分彩计划网站 混合过关中有一场输平手 双色球规则更改 时时彩开彩结果 北京赛车pk10技巧5码公式 新时时彩二星和值遗漏 15选5 大奖复式 手机投注新浪彩票 pk10大小单双在线计划 北京时时彩规律大全 pk10必胜单双技巧 pc蛋蛋兑奖有限制吗 中大奖技巧 广东快乐十分官网