火山(Volcano Plot)图在一张图中显示了两个重要的指标(Fold change/pvalue),可以非常直观且合理地筛选出在两样本间发生差异表达的基因。检验分析出两样本间显著差异表达的基因后,以log2(fold change)为横坐标,以T检验显著性检验P值的负对数-log10(pvalue)为纵坐标,即可得火山图(Volcano Plot)。

> install.packages("ggplot2")

> library(ggplot2)

> R0_vs_R3_isoforms_fitter

> data

> r03

> r03 + geom_point()

运行过程数据及图片:

改变点的颜色:

> r03 + geom_point(color="red")

> r03 + geom_point(aes(color="red"))

> r03 + geom_point(aes(color=significant))

设置标题,坐标轴范围

>r03+geom_point(aes(color=significant))+xlim(-4,4)+ylim(0,30)

>r03xy=r03+geom_point(aes(color=significant))+xlim(-4,4)+ylim(0,30)

> r03xy + labs(title="Volcano plot",x="log2(FC)")

> r03xy + labs(title="Volcanoplot",x=expression(log[2](FC),y=expression(-log[10](FDR))))

自定义颜色

> r03xy + labs(title="Volcanoplot",x=expression(log[2](FC),y=expression(-log[10](FDR))))

> r03xyc = r03xy+labs(title="Volcanoplot",x=expression(log[2](FC)), y=expression(-log[10](FDR)))

> r03xyc + scale_color_manual(values =c("green","black", "red"))

>r03xyc + scale_color_manual(values =c("green","blue", "red"))

> volcano = r03xyc +scale_color_manual(values = c("green","blue","red"))

>volcano+geom_hline(yintercept=1.3)+geom_vline(xintercept=c(-1,1))

>volcano+geom_hline(yintercept=1.3,linetype=4)+geom_vline(xintercept=c(-1,1),linetype=4)

保存图片:

>ggsave("volcano.png")

>ggsave("volcano8.png",volcano,width=8,height=8)

好吧,学习使人疯狂,脑袋疼学起画图一样起劲的忘休息。

matlab的火山图,ggplot绘制火山图相关推荐

  1. ggplot绘制小提琴图

    我们既往已经讲过使用ggplot绘制箱型图,今天我们来聊聊ggplot绘制小提琴图,小提琴图在SCI论文中非常常见,特别是基因类的文章. 小提琴图其实可以算是箱型图的变种(等于:箱型图+核密度图),箱 ...

  2. 利用思维导图软件绘制鱼骨图怎样做

    思维导图的类型有哪些? 圆圈图 主要用于把一个主题展开来,联想或描述细节.小圈圈是主题,而外面的大圈圈里放的是和这个主题有关的细节或特征 气泡图 由很多泡泡组成,中间一个主题泡泡描述核心主题,周围的属 ...

  3. R语言ggplot绘制鸡冠花图(南丁格尔玫瑰图)

    南丁格尔玫瑰图也就是鸡冠花图,一种极坐标的​柱状图.图形的大小跟变量分组的数值呈正比.从极坐标延伸出来的每一环可以当做标尺来使用,用来表示分段的大小并代表数值的高低. 首先绘制柱状图 bar<- ...

  4. 圆形和多边形雷达图python-Matplotlib绘制雷达图和三维图的示例代码

    1.雷达图 程序示例 '''1.空白极坐标图''' import matplotlib.pyplot as plt plt.polar() plt.show() '''2.绘制一个极坐标点''' im ...

  5. 棒棒糖也能画成图—— ggplot2绘制棒棒糖图/哑铃图

    欢迎关注公众号查看原文 Part 1 :棒棒糖图 棒棒糖图因其形状和棒棒糖相似而得名,具体来看实际上是一个散点和一条线段的组合.棒棒糖图是散点图的一种变体,又与柱状图非常相似,但其在清晰展示数据的同时 ...

  6. python话雷达图-PYTHON绘制雷达图代码实例

    这篇文章主要介绍了PYTHON绘制雷达图代码实例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下 1.雷达图 import matplotlib.py ...

  7. 利用python绘制雪景图_python绘制雪景图

    本文实例为大家分享了python绘制雪景图的具体代码,供大家参考,具体内容如下 绘制雪景图,应用到turtle和random. from turtle import * from random imp ...

  8. python彩虹图_python绘制彩虹图

    本文实例为大家分享了python绘制彩虹图的具体代码,供大家参考,具体内容如下 from turtle import * #控制彩虹路径 def path(pen, r, g, b): pen.pen ...

  9. python分析数据走势图_python绘制趋势图的示例

    import matplotlib.pyplot as plt #plt用于显示图片 import matplotlib.image as mping #mping用于读取图片 import date ...

  10. python热图_Python绘制热图

    1.matplitlib绘制热图 import random from matplotlib import pyplot as plt from matplotlib import cm from m ...

最新文章

  1. 动态加载flex皮肤.
  2. DbHelperSQL.cs
  3. ECMAScript6 规范
  4. 1047. Student List for Course (25)
  5. 使用Visual Studio Code调试运行在SAP云平台上处于运行状态的nodejs应用
  6. 柠檬汁制成的电池可以开动超100千克的车子吗?
  7. 天池 在线编程 放小球(动态规划)
  8. 关于字节对齐以及内存占用
  9. Linux系统如何安装VMware(虚拟机)
  10. js base64图片太大_JS实现base64图片下载 简易方法
  11. 图片像素对比OpenCV实现,实现人工分割跟算法分割图像结果的对比
  12. 计算机excel素材,excel操作练习题素材.doc
  13. 155款安卓开源项目源码整理+20个Android必备第三方框架
  14. 将子龙山人的box2d碰撞小球修改到了cocos2d-x下,提供源码
  15. 每日新闻早报简报十条和1条微语分享
  16. 用计算机完成下表的视距测量计算公式,2012测量学计算题库及参考答案
  17. 「硬核讲解」通达信跨周期引用均线指标公式
  18. ARFoundation入门到精通 - 1.8 远程调试
  19. mysql到底是什么_MySQL数据库之MySQL到底是什么
  20. 大家都看得见学计算机的决定,卢须振二丫-图书检索系统.PDF

热门文章

  1. python读取lmdb文件_python中读写LMDB数据库
  2. 【FICO】SAP中的银行
  3. Guice依赖注入(接口多实现)
  4. VMware Workstation 12 Pro的安装
  5. 数字化转型背景下的“新IT职业教育” | 职业教育系列报告(四)
  6. 安徽省2019c语言二级答案,2019-C语言二级考试题库(含答案).docx
  7. Java实现 LeetCode 289 生命游戏
  8. 苹果系统摩尔庄园是什么服务器,摩尔庄园手游官服和渠道服有什么区别_可以一起玩吗_官服和渠道服详细介绍...
  9. VMware12虚拟机中安装win7(安装好的解压版Win7_64-cl1.vmdk)
  10. 微信公众号开发之微信支付开发