3

「pheatmap」パッケージを使用してヒートマップを生成するための次のコードがあります。ヒートマップの作成には成功しましたが、列に注釈を付けるのに非常に苦労しています。9 つの細胞型があり、各細胞型には 3 つの複製があります。

これがデータです。サンプルには 01_01,01_02,01_03 から 09_01, 09_02,09_03 までのラベルを付けました。

列の注釈を「A」、「B」、「C」、「D」、「E」、「F」、「G」、「H」、「I」にしたいと思います。行と列をクラスタリングする pheatmap

クラスター列に注釈を付けるために最善を尽くしましたが、成功しませんでした。誰かが列の注釈をカスタマイズする方法を説明できれば、本当に感謝しています。

データ: 私のデータファイル

library(readxl)
library(pheatmap)
library(tibble)
library(dplyr)

ファイルを読む

test<-read_excel("~test.xlsx",1)

test%>% 
   column_to_rownames("Metabolite") %>%
   as.data.frame()

test<-as.data.frame(test)
rownames(test)<-test$Metabolite
test$Metabolite<-NULL


    pheatmap(test,cluster_rows=TRUE,cluster_cols=TRUE,breaks=NA,scale="none",legend=TRUE,color=colorRampPalette(c("navy","white","firebrick3"))(50),margins=c(5,10),fontsize_row=4,fontsize_col=4,cellwidth=5,cellheight=5)
4

0 に答える 0