BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图
Volcano plot | 别再问我这为什么是火山图 一文解释了火山图如何解读。不太难看懂,而一旦看懂了,图也就知道怎么绘制了。
假设我们已经有了一个差异基因鉴定后的表格文件 590e7b6b-c279-40da-b1d2-1017464cea02.untrt.vs.trt.results.txt (看到这一串无规律的符号做文件名就知道这是我们平台输出的差异分析结果了,之前介绍过可以翻一翻),那怎么绘制火山图呢?
上传文件
首先注册个账户。这里不是强制大家注册,只是火山图文件一般比较大,直接粘贴进文本框会导致浏览器卡顿,为了更好的用户体验,建议注册个账户,在个人中心上传 (支持断点续传,大文件也不怕)。如果想直接粘贴,也没问题。
绘制火山图
进入到火山图绘制页面http://www.ehbio.com/Cloud_Platform/front/#/analysis?page=b%27Nw%3D%3D%27。
Input way
选择Select uploaded file
;选择上传的文件,如果没找到,则刷新下页面再重复此操作。
选择的文件会在文本域中显示预览(不可修改,但可以复制),点击Check Data
,核对数据格式没有问题,会激活下方的必选选项。
选择对应的列进行绘制,绘图就是把数据的结构用几何形状表示出来,并用颜色、大小等代表特定的属性展示。
从下拉框选择
log2FoldChange
指定为Fold change column
;为什么选这一列,因为这个参数值跟参数名字太像了;
从下拉框选择
padj
指定为Statistical significance column
;为什么选这一列,因为
pval
和padj
都是统计显著性的指示指标,为了解决总被审稿人提起的多重假设检验校正是什么?,我们选padj
列;padj
列数值没有进行过转换,这里选择Log10 transform significance value
。Padj
越小转换后的值越大,越在图的上方;Gene expression change status variable
这是一个可选参数,是说文件中是否已经根据某个阈值做了差异基因标记,哪些上调了,哪些下调了,就是这里面的level
列,我们选择上;如果之前没做过筛选,没有
level
列也没关系,DE genes filtering threshold
参数可以设置筛选阈值,默认为0.05,1
,第一个数字0.05
表示统计pvalue<0.05
或padj<0.05
(取决于Statistical significance column
的选择);第二个数字1
表示变化倍数,变化倍数取完log2
后的绝对值大于1
。
点击提交,出来一个结果; 颜色不对,需要调整下。
先设置下level
变量的顺序,然后按顺序设置颜色
Status variable order
依次选择下拉内容为trt up, untrt UP, NoDiff
;Customized point colors
用颜色选择器选择3个颜色
提交后获得结果,正确了,可以下载PDF格式和相关的R代码
上图中的两条垂直虚线和一条水平虚线是参数DE genes filtering threshold
控制的,如果你筛选差异基因的标准(生成level
列中哪些上调、哪些下调的标准)不是默认标志,则需要修改这个值为你设置的阈值,从而调整线的位置。
也可以设置参数Coordinate flip
旋转火山图。
用的不多,但有用户提过这个需求,就加上了。
另外一个常见需求是在火山图上标记基因,现在的策略是要在数据表中加一列,不太方便,回头更新一版新的处理方式再介绍。
如果你想了解哪个图的绘制,请留言告诉我们,我们优先推出教程。
测试数据获取:https://gitee.com/ct5869/bic
点击阅读原文访问网站
往期精品(点击图片直达文字对应教程)
机器学习
后台回复“生信宝典福利第一波”或点击阅读原文获取教程合集
BIC无代码绘制差异基因火山图相关推荐
- ImageGP/BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图 Volcano plot | 别再问我这为什么是火山图 一文解释了火山图如何解读.不太难看懂,而一旦看懂了,图也就知道怎么绘制了. 假设我们已经有了一个差异基因鉴定后的表格文 ...
- ggplot2绘制差异表达基因火山图
一.前置环境 1.1 R 语言 下载对应系统的R软件 R: The R Project for Statistical Computing (r-project.org) 以win11为演示 http ...
- 一行代码绘制高分SCI火山图
一.概述 在近半年中,我读了很多的高分SCI文章,很多文章中都有多种不同的火山图,包括「普通的火山图.渐变火山图.以及包含GO通路信息的火山图」! 经过一段时间的文献阅读和资料查询,终于找到了一个好用 ...
- 无代码绘制基因表达箱线图
给定一个基因表达矩阵和样本分组信息,如何绘制样品整体表达箱线图.单个或多个基因表达箱线图. 获取示例数据 我们从中截取前面4行作为演示例子. 利用工具Wide to Long把表达矩阵转换为长表格 把 ...
- 转录组-差异基因热图
top_de_exp<-dplyr::slice(de_result2,1:20)%>%#挑取差异最大的select(-c(2:8))%>%#去掉2-8列column_to_rown ...
- 差异表达基因-火山图和聚类图解释
想研究某现象的分子机制,老板豪气的来一句,先测个转录组吧,看下差异表达基因. 是否在心里窃喜,制个样就完事了,太easy有木有.等大堆数据回来的时候,是不是傻眼了? 从何下手挑选差异表达基因呢? 今天 ...
- 【Python】pandas一行代码绘制26种美图
本文目录 1.单组折线图 2.多组折线图 3.单组条形图 4.多组条形图 5.堆积条形图 6.水平堆积条形图 7.直方图 8.分面直方图 9.箱图 10.面积图 11.堆积面积图 12.散点图 13. ...
- 差异表达基因火山图(ggplot函数)
1. 读入数据 差异表达基因来自limma分析结果. # read the file data <- read.csv("diff_expr_genes.csv",row.n ...
- pandas一行代码绘制26种美图
↑↑↑关注后"星标"简说Python 人人都可以简单入门Python.爬虫.数据分析 简说Python推荐 来源:pythonic生物人作者:pythonic生物人 本文目录 1. ...
最新文章
- 使用Hyperopt实现机器学习自动调参
- 【LeetCode笔记】347. 前K个高频元素(Java、优先队列、小顶堆、HashMap)
- 安徽关节式焊接机器人_机器人自动焊接技术的优势及应用介绍
- rc脚本(类的定义与脚本的结构、start方法、stop和status方法、以daemon方式启动)...
- elasticsarch6.5.4安装插件 searchguard和elasticsearh-head插件安全性问题
- 子类重写方法aop切不到_SpringBoot源码之旅——AOP
- Docker学习总结(21)——Docke网络bridge详解
- python如何输入多行数据合并_Python如何将多行数据合并成一行|python如何实现excle数据合并...
- CAD手机精准看图隐私政策
- ant design-Table中的分页实现(包括点击跳页,输入跳页,以及指定每页项目数)
- case when 失效,看了这篇文章就明白了
- Latex 中文Beamer模板
- 手机话费充值页面HTMLcss3+html5模板
- echarts 环形图 不同区域背景色自定义
- 【5GC】什么是5G切片?5G切片如何工作?
- js 数字转为千分符格式
- 大前端求人不如求己系列工具:如何PSD文件自动标注切图和真机预览
- 运行时绘制Gizmo——关于unity3D的GL图像库的使用(非常详细)
- 计算机网络中怎么共享文件,图文教你如何设置局域网里网络邻居中的共享文件...
- 电脑:软件(win)