GO注释和富集分析

GO注释和富集分析使用TBtools完成,具体步骤可以参考TBtools作者在腾讯课堂开设的一系列视频课程

本文使用的数据是甜樱桃叶绿体蛋白编码基因做GO注释,然后挑部分基因做富集分析,挑选的基因是

rpoC1

rpoB

rpoA

rpoC2

atpI

atpF

atpE

atpH

atpB

atpA

accD

rbcL

rpl22

rpl23

rpl20

rps8

rps7

rps16

rps15

rps14

rps18

做完富集分析得到文件GOenrichmentOutput.txt..GO.Enrichment.final.xls 根据GOplot包的示例数据挑选出其中的5列

Class GO_Name GO_ID GenesOfSelectedSetInGOterm corrected p-value(BH method)

作为数据集1 数据集2包括

ID,logFC,AveExpr,t,P.Value,adj.P.Val,B

数据集2的列变量应该都是转录组数据分析的结果 比如logFC应该是倍数变化Fold change 然后取log AveExpr应该是平均表达量等 然后模仿帮助文档的例子构造数据集

help(package="GOplot")

library(GOplot)

data(EC)

file1

file2

df1

df2

colnames(df1)

df

dim(df)

df气泡图

GOBubble(df)

- 加一些参数的气泡图

GOBubble(df,table.legend = F,table.col = T,ID=T,

display = "multiple")+

scale_x_continuous(limits=c(0,5))+

scale_y_continuous(limits = c(0,10))

- 弦图

pdf("chordpra.pdf",height=15,width = 15)

chord

GOChord(chord)

dev.off()

现在基本可以根据自己的数据来构造GOplot的输入文件,但是作图的具体细节还需要调整 欢迎大家关注我的公众号

大家如果需要文中提到的file1和file2,在我的公众号小明的数据分析笔记本留言就可以了

r语言进行go富集分析_GO富集分析可视化:R语言GOplot包——准备自己的数据相关推荐

  1. r语言进行go富集分析_GO富集柱状图

    分析模块,输入差异基因GO富集分析结果,由分析模块"GO Enrichment Analysis"生成,输出GO富集结果柱状图. 输入: 差异基因GO富集分析结果文件,由分析模块& ...

  2. r语言进行go富集分析_GO富集分析可视化:R语言GOplot包

    找R语言做弦图的教程的时候发现了这个包:GOplot.其主要功能是可视化GO富集分析的结果.自己应该会用得到. 第一步是学习其帮助文档中的例子,然后学习如何准备自己的数据,并利用这个包中的函数来绘图 ...

  3. R语言可视化散点图、ggrepel包geom_text_repel函数避免数据点之间的标签互相重叠(设置nudge_x参数、hjust参数和direction参数配置数据点标签在可视化图像左右两侧边缘

    R语言ggplot2可视化:可视化散点图并为散点图数据点添加文本标签.使用ggrepel包的geom_text_repel函数避免数据点之间的标签互相重叠(设置nudge_x参数.hjust参数和di ...

  4. 使用R语言进行Metroplis-in-Gibbs采样和MCMC运行分析

    全文链接:http://tecdat.cn/?p=12200 对于许多模型,例如逻辑模型,没有共轭先验分布.因此,吉布斯采样不适用(点击文末"阅读原文"获取完整代码数据). 这篇文 ...

  5. 使用R语言绘制富集条形图,轻松分析基因表达数据

    一.引言 富集分析(enrichment analysis)是一种生物信息学方法,它可以帮助我们识别基因或其他的生物实体在某个特定的类别中过度表示的趋势.通俗来说,富集分析通过将基因分类到特定的集合中 ...

  6. R软件--GEO数据分析教程:差异性分析、富集分析(GO\KEGG\DO)

    title: "GSE93798" output: word_document knitr::opts_chunk$set(echo = TRUE) R Markdown 关于GS ...

  7. 富集分析和基因表达花样可视化

    GOplot包介绍 GOplot包用于生物数据的可视化.更确切地说,该包将表达数据与功能分析的结果整合并进行可视化.但是要注意该包不能用于执行这些分析,只能把分析结果进行可视化.在所有科学领域,由于空 ...

  8. 转录组GO富集与微生物相关性分析

    转录组GO富集与微生物相关性分析 原始数据格式 使用TCC包进行差异基因分析 使用topGO包进行GO富集分析 计算转录组与微生物组相关性 根据GO号提取相关性矩阵并绘图 提取相关性矩阵 绘制热图 输 ...

  9. R语言——GO分析和KEGG分析

    目录 介绍 计算原理 实际操作过程 加载相关包 读取文件 格式转化 GO富集分析 KEGG分析 保存相关结果 介绍 在DEseq2的基础上,进行GO分析和KEGG分析 # 设置国内镜像,安装时运行一次 ...

最新文章

  1. VS013的单元测试去哪里了
  2. python【蓝桥杯vip练习题库】ALGO-50 数组查找及替换
  3. 【AI】Win10-Tensorflow
  4. python语言介绍-00-python语言介绍
  5. android串口通讯奇偶校验,串口通讯奇偶数校验及CRC校验如何使用详解
  6. 线上Tomcat支持Eclipse远程调试的方法
  7. 电子设计基础——电源的各项指标及测量方法
  8. Bubble Sort Aizu - ALDS1_2_A
  9. 不作35岁的程序员?
  10. AlexNet--CNN经典网络模型详解(pytorch实现)
  11. macbook装双系统多分区其实很简单,你只要把macbook当作一台普通pc就可以了!
  12. WEEX|初始化工程
  13. javascript函数传参方式
  14. 关于程序员的面试于自我介绍模板
  15. 网易音乐社交产品 MUS 上线
  16. python格式化输出函数之format
  17. 逆向易语言程序执行流程
  18. 网易2018校园招聘:相反数 [python]
  19. 计算电费(厦大PAT)
  20. 认沽期权是什么意思:认沽期权怎么交易

热门文章

  1. 在Vs Code中搭建JSP开发环境
  2. 电子书包教师端界面简介
  3. java openjdk_OpenJDK的安装与配置
  4. openwrt 最简设置ipv6
  5. windows使用scoop安装gcc报404错误的解决方法
  6. FDK算法的实现过程
  7. VMware Tools 其实很重要
  8. 求矩阵的1,和2范数
  9. 深度学习网络和特征图可视化的工具介绍
  10. 目标管理的培训要点是什么?