生信分析矫正P值_好看的富集分析图GOplot
###########安装#install.packages('GOplot')#install_github('wencke/wencke.github.io')
###########加载library(GOplot)
数据准备:
data(EC)head(EC$david)###DAVID富集分析结果head(EC$genelist)###差异分析结果###大家也可以根据上面的格式将自己的文件读入circ head(circ)
其中,category代表GO term categoryID代表GO term ID term代表GO term namecount代表该GO term中的基因数量genes代表GO term中相关的基因名logFC代表相关基因差异表达分析中的logFC值adj_pval代表GO term的矫正后p值
其中up,down,count,分别代表在相关GO term中上调,下调,总计的基因数准备好上面的数据,就可以开始绘制各种GO富集分析图啦~figure-1 GOBar
GOBar(circ, display = 'multiple', title = 'Z-score coloured barplot', ##设置标题 zsc.col = c('firebrick3', 'white', 'royalblue3')##设置颜色 )
figure-2 GOBubble
OBubble(circ, ,##设置标题
GOBubble(circ, title = 'GOBubble plot', labels = 3,
figure-3 GOCircle
GOCircle(circ, rad1=2, rad2=3,##内径、外径设置 zsc.col=c('firebrick3', 'white', 'royalblue3'),#z-score颜色设置 lfc.col=c('firebrick3', 'royalblue3'),##上调下调颜色设置 label.size=5,label.fontface='bold',##字体大小格式设置 table.legend=TRUE##右侧表格设置 )
GOCircle(circ, nsub = 5)###选择GO term 进行展示(即选取前n行展示)
IDs GOCircle(circ, nsub = IDs)##用户自定义选择GO term 进行展示
figure-4 GOChord
head(EC$genes)
figure-5 GOHeat
GOHeat(chord[,-ncol(chord)], ##去除最后一列logFC nlfc = 0,#定义logFC列的数量(默认值为0) )
GOHeat(chord, nlfc = 1, fill.col = c('firebrick3', 'white','royalblue3')#颜色设置 )
figure-6 GOCluster
GOCluster(data=circ, process=EC$process,#选择GO term metric='euclidean',#选择距离度量方法 clust='average', #选择聚类方法 clust.by = 'term',#指定是否应该对基因表达模式或功能类别进行聚类。term(default) or logFC term.width = 2,term.col=brewer.pal(length(EC$process), "Set3"),#GO term 宽度,颜色设置 nlfc=FALSE, lfc.col=c('firebrick3', 'white','royalblue3')#是否包含多个logFC列,颜色设置 )
figure-7 GOVenn(max 3sets)
l1 l2 l3 GOVenn(l1,l2,l3, lfc.col=c('firebrick3', 'white','royalblue3'),#logFC颜色设置 circle.col=brewer.pal(3, "Dark2"),#Venn颜色设置 label = c('heart development', 'plasma membrane', 'tissue morphogenesis')#标签设置 )
可共参选的样式还是很多的,大家可以尝试起来喽~参考网址:https://wencke.github.io/
欢迎关注生信人
TCGA | 小工具 | 数据库 |组装| 注释 | 基因家族 | Pvalue
基因预测 |bestorf | sci | NAR | 在线工具 | 生存分析 | 热图
生信不死 | 初学者 | circRNA | 一箭画心| 十二生肖 | circos
舞台|基因组 | 黄金测序 | 套路 | 杂谈组装 | 进化 | 测序简史
生信分析矫正P值_好看的富集分析图GOplot相关推荐
- 生信分析矫正P值_生物信息分析:从入门到精(fang)通(qi) 第1期FASTQ! BAM! VCF! 傻傻分不清楚...
生信小白:肉哥,上次听完你的介绍,我满脑子跟这张图片一样...凌乱?!我们为什么不一次性把整本书读取了,非要把这本书撕碎呢? 西克孚肉:这主要受限于技术,测序仪一次只能读取几百.几千.几万个碱基,这与 ...
- 生信分析矫正P值_TCGA数据分析系列(二):数据库之GEPIA2
所谓工欲善其事,必先利其器,从今天开始,我们来介绍TCGA数据库的使用.今天我们来介绍一款非常容易上手的数据库:GEPIA2 链接点这里 话不多说,直接进入正题 GEP1A2数据来源 如图所示,GEP ...
- docker卸载命令_使用docker完成生信分析环境搭建
生信开发人员最头疼的问题,可能就是平台搭建和软件安装了.部署和迁移上要费很大力气.本文讲述使用docker制作一个镜像,后续通过导入自己定制的镜像,复制文件完成分析流程的部署和迁移. 如何使用dock ...
- 属于服务器端运行的程序_生信分析云平台产品开发 - 5 生信分析pipeline服务器端运行...
在上文 [生信分析云平台产品开发 - 4 生信分析pipeline的图形化] 讨论了生信分析pipeline的图形化,如何用图形的方式显示生信pipeline,但是pipeline脚本按照变量的形式保 ...
- mirna富集分析_经验之谈丨生信分析文章套路原来这么简单!
近两年,不做实验或者仅需要少量实验的生物信息学分析文章,发表量越来越多.如果利用数据库检索,高效的发出一篇文章.是科研工作者关注的话题,今天我们就用一篇生信分析的文章作为切入点,来谈谈生信分析文章的套 ...
- r语言python0基础_万丈高楼平地起,从零开始学R语(生信分析向,R语言基础)...
一.R语言的安装 我所做的R语言学习主要是为应对工作上的生信分析,如果你是往这个方向走,那么这篇文章很大可能会帮助你.如果你也是有python基础的话,那我写的可能会更対你口味. 不过现在看到这篇文章 ...
- cluego使用说明_生信分析绘图神器,你值得拥有!
GO和KEGG分析是最常用的生信分析方法,在SCI论文中也经常见到,那么你能想到的GO和KEGG分析结果的展示方法有哪些? 条形图? 条形图2? 饼状图? 表格? 相比于上面这些,这样的网络图展示起来 ...
- perl语言入门第七版 电子版_百迈客带您走近生信分析【入门篇】
年末促销倒计时:59天 百迈客推出年末活动促销,发文有礼.推广有礼以及多种产品钜惠来袭,百迈客为您倾情打造科研福利,您还在等什么?快快行动起来领取您的超级奖励吧!(详情请见"决战2020!品 ...
- 2020.09.30【RNA-seq流程】丨转录组生信分析全流程
RNA-Seq生信分析全流程 摘要 第一部分 step.1 下载数据 step.2 数据质控 第二部分 step.3序列比对 step.4 计算基因表达量 step.5 插入片段长度检验 step.6 ...
最新文章
- ckeditor 在dwz里面使用
- 史上最全SpaceX火箭数据开源,核心、组员舱、起落架、发射信息全都有!
- python通过代理发送邮件_Python实现SMTP发送邮件详细教程
- GNU make manual 翻译( 一百零九)
- GDCM:读取gdcm::Curve的测试程序
- Shell(6)——array的删改unset
- werkzeug create_environ
- Go 自动构建工具 dogo 代码已托管到 Git@OSC
- 嵌入式linux led驱动有几种写法,嵌入式Linux字符驱动LED灯设计
- Python数据类型(数字)
- ExecutorService的shutDown和shutDownNow方法的区别
- how to set up github blog
- python零基础编程入门先学什么_零基础怎么学Python编程,新手常犯哪些错误
- 读《亿级流量网站架构核心技术》
- 无线电发射器 Android,简单无线电app
- 两个表格数据合并怎么完成?
- 大三下了,一点感概,仅作共勉
- 头歌c语言实训作业题解
- 在 Linux 上安装 chm 文件阅读器
- apply方法和unapply方法