matlab的火山图,ggplot绘制火山图
火山(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绘制火山图相关推荐
- ggplot绘制小提琴图
我们既往已经讲过使用ggplot绘制箱型图,今天我们来聊聊ggplot绘制小提琴图,小提琴图在SCI论文中非常常见,特别是基因类的文章. 小提琴图其实可以算是箱型图的变种(等于:箱型图+核密度图),箱 ...
- 利用思维导图软件绘制鱼骨图怎样做
思维导图的类型有哪些? 圆圈图 主要用于把一个主题展开来,联想或描述细节.小圈圈是主题,而外面的大圈圈里放的是和这个主题有关的细节或特征 气泡图 由很多泡泡组成,中间一个主题泡泡描述核心主题,周围的属 ...
- R语言ggplot绘制鸡冠花图(南丁格尔玫瑰图)
南丁格尔玫瑰图也就是鸡冠花图,一种极坐标的柱状图.图形的大小跟变量分组的数值呈正比.从极坐标延伸出来的每一环可以当做标尺来使用,用来表示分段的大小并代表数值的高低. 首先绘制柱状图 bar<- ...
- 圆形和多边形雷达图python-Matplotlib绘制雷达图和三维图的示例代码
1.雷达图 程序示例 '''1.空白极坐标图''' import matplotlib.pyplot as plt plt.polar() plt.show() '''2.绘制一个极坐标点''' im ...
- 棒棒糖也能画成图—— ggplot2绘制棒棒糖图/哑铃图
欢迎关注公众号查看原文 Part 1 :棒棒糖图 棒棒糖图因其形状和棒棒糖相似而得名,具体来看实际上是一个散点和一条线段的组合.棒棒糖图是散点图的一种变体,又与柱状图非常相似,但其在清晰展示数据的同时 ...
- python话雷达图-PYTHON绘制雷达图代码实例
这篇文章主要介绍了PYTHON绘制雷达图代码实例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下 1.雷达图 import matplotlib.py ...
- 利用python绘制雪景图_python绘制雪景图
本文实例为大家分享了python绘制雪景图的具体代码,供大家参考,具体内容如下 绘制雪景图,应用到turtle和random. from turtle import * from random imp ...
- python彩虹图_python绘制彩虹图
本文实例为大家分享了python绘制彩虹图的具体代码,供大家参考,具体内容如下 from turtle import * #控制彩虹路径 def path(pen, r, g, b): pen.pen ...
- python分析数据走势图_python绘制趋势图的示例
import matplotlib.pyplot as plt #plt用于显示图片 import matplotlib.image as mping #mping用于读取图片 import date ...
- python热图_Python绘制热图
1.matplitlib绘制热图 import random from matplotlib import pyplot as plt from matplotlib import cm from m ...
最新文章
- 动态加载flex皮肤.
- DbHelperSQL.cs
- ECMAScript6 规范
- 1047. Student List for Course (25)
- 使用Visual Studio Code调试运行在SAP云平台上处于运行状态的nodejs应用
- 柠檬汁制成的电池可以开动超100千克的车子吗?
- 天池 在线编程 放小球(动态规划)
- 关于字节对齐以及内存占用
- Linux系统如何安装VMware(虚拟机)
- js base64图片太大_JS实现base64图片下载 简易方法
- 图片像素对比OpenCV实现,实现人工分割跟算法分割图像结果的对比
- 计算机excel素材,excel操作练习题素材.doc
- 155款安卓开源项目源码整理+20个Android必备第三方框架
- 将子龙山人的box2d碰撞小球修改到了cocos2d-x下,提供源码
- 每日新闻早报简报十条和1条微语分享
- 用计算机完成下表的视距测量计算公式,2012测量学计算题库及参考答案
- 「硬核讲解」通达信跨周期引用均线指标公式
- ARFoundation入门到精通 - 1.8 远程调试
- mysql到底是什么_MySQL数据库之MySQL到底是什么
- 大家都看得见学计算机的决定,卢须振二丫-图书检索系统.PDF
热门文章
- python读取lmdb文件_python中读写LMDB数据库
- 【FICO】SAP中的银行
- Guice依赖注入(接口多实现)
- VMware Workstation 12 Pro的安装
- 数字化转型背景下的“新IT职业教育” | 职业教育系列报告(四)
- 安徽省2019c语言二级答案,2019-C语言二级考试题库(含答案).docx
- Java实现 LeetCode 289 生命游戏
- 苹果系统摩尔庄园是什么服务器,摩尔庄园手游官服和渠道服有什么区别_可以一起玩吗_官服和渠道服详细介绍...
- VMware12虚拟机中安装win7(安装好的解压版Win7_64-cl1.vmdk)
- 微信公众号开发之微信支付开发