ggplot2热图扩展包（ggalign）的细节

文章预览

前面我们在 人工智能大模型不会告诉你的热图绘制技巧 演示了如何使用ggplot2热图扩展包（ggalign），可以快速替代之前的 pheatmap： 比如我们可以先去geo数据库里面下载 GSE104171_NormalizedMatrix.txt.gz 这个文件， 然后提取里面的指定的基因的表达量矩阵： rm(list = ls()) #清空当前的工作环境 options(stringsAsFactors = F) #不以因子变量读取 options(scipen = 20) #不以科学计数法显示 library(data.table) library(tinyarray)  data < -data.table::fread( "GSE104171_NormalizedMatrix.txt.gz" ,                         data.table = F)  data=data[!duplicated(data $V1 ),] mat < -data[,c( 2:ncol(data))] rownames(mat)=data[,1] mat[1:4,1:4] keep_feature  < - rowSums (mat > 1) > 1 ;table(keep_feature) ensembl_matrix  < - mat[keep_feature, ]   rownames(ensembl_matrix)=rownames(mat)[keep_feature] ensembl_matrix[1:4,1:4]  symbol_matrix=ensembl_matrix symbol_matr ………………………………