我已经创建了一个矩阵,现在我想使用 pheatmap 来绘制热图,同时保留矩阵行的顺序。我想关闭集群。目前,当我运行以下脚本时,pheatamp 正在对行进行聚类:
tissuedata<-read.table('OrderedByGeneCounts.csv', header=TRUE, sep=',')
tissuedata2<-data.matrix(tissuedata)
tissuedata3<-log10(tissuedata2)
is.na(tissuedata3) <- do.call(cbind,lapply(tissuedata3, is.infinite))
tissuedata3[is.na(tissuedata3)] <- 0
tissue_exons<-scan(file="OrderedExonsList.txt", what="character", strip.white=TRUE)
rownames(tissuedata3)<-tissue_exons
pdf(file="pheatmapENSG00000026025.pdf", paper="a4", height=15)
pheatmap(tissuedata3,
color = colorRampPalette(rev(c("#D73027", "#FC8D59", "#FEE090", "#FFFFBF", "#E0F3F8", "#91BFDB", "#4575B4")))(100),
cellwidth = 25, cellheight = 3.8,
treeheight_row = 0,
treeheight_column=0,
kmeans_k = NA,
show_rownames = T, show_colnames = T,
fontsize=4,
scale="none",
legend=FALSE,
)
dev.off()