转录组——热图（pheatmap）

最新推荐文章于 2024-08-07 16:07:24 发布

wt12138

最新推荐文章于 2024-08-07 16:07:24 发布

阅读量3.7k

点赞数

文章标签：聚类 csv openssh bind yum

本文链接：https://blog.csdn.net/wt141643/article/details/110424853

版权

本文详细介绍了如何利用pheatmap包在R中对转录组数据进行热图绘制，包括数据预处理、csv导入、聚类分析等步骤，帮助读者理解基因表达模式。

摘要由CSDN通过智能技术生成

###########################
# 2020/11/18 heatmap
###########################
setwd(dir = "../muscle")
library(tidyverse)
# 1. 导入并筛选，差异基因
gene_info <- read.csv(file = "zd_gene_info.csv")
names(gene_info) <- c("gene_id","Swissprot_ID","Gene_symbol","Function")
data <- read.csv(file = "../muscle/res",
                 header = T,
                 sep = "\t" )
de_result <- left_join(data, gene_info, by = c("id" = "gene_id")) # 合并表格
all_gene_TMM <- read.csv(file = "gene_all.TMM.EXPR.matrix",sep = "\t")
names(all_gene_TMM) <- c("id","H1","H2","H3","L1","L2","L3")
plot_result <- left_join(de_result,all_gene_TMM, by = c("id"="id"))
deg <- select(plot_result , id , log2Fo