转录组DEGs聚类热图和功能富集分析 |
您所在的位置:网站首页 › commen名词 › 转录组DEGs聚类热图和功能富集分析 |
写在前面:经常做转录组分析,就是把差异基因搞个火山图和Venn图看各组差异基因的共有和特有情况。看见有个比较好的选择,能直观比较各种处理带来的影响,如下: image.png 来自Nature plants的一篇文章 Ref: https://github.com/YulongNiu/MPIPZ_microbe-host_homeostasis https://www.nature.com/articles/s41477-021-00920-2 这个图很牛逼啊,表示的信息量很全,值得学习。去扒作者的代码,复现出了大部分所需文件: 总的基因丰度表,即各个基因在每个样品中的丰度 image.png各个样品的基因差异表达情况 举一个例子,这是Deseq2算出来的: image.png 1. 导入数据并做一些处理 setwd("C:/Users/yjk/Desktop/cluster_DEGs") library("dplyr") library("ComplexHeatmap") library("tibble") library("RColorBrewer") library("ggplot2") my_theme() all_genes 1, padj < 0.05 HF12N_Rs % write.table('kmeans10_KEGG.txt', quote = FALSE, sep = "\t")列一个例子,GO富集的生物学过程 image.png到此结束 |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |