將組信息添加到 pvclust 樹狀圖而不是樣本名稱

Question

這就是我正在做的

mat <-  read.table("Model_pvclust/Model18_FAB_M5_vs_M0_MAP_TF.txt",sep = "\t",strip.white = FALSE,check.names = FALSE,header=TRUE,row.names=1)
drop <- c("gene","baseMean","log2FoldChange","lfcSE","stat","pvalue","padj","UP_DOWN")

d1 <- select(mat, -one_of(drop))

####### 讀取元數據

metadata <-read.table("Model_hmap_meta/FAB_table.txt",sep = "\t",strip.white = FALSE,check.names = FALSE,header=TRUE,row.names=1)
head(metadata)



head(metadata)
             prior_malignancy FAB    Risk_Cyto
TCGA-AB-2856               no  M4 Intermediate
TCGA-AB-2849               no  M0         Poor
TCGA-AB-2971               no  M4 Intermediate
TCGA-AB-2930               no  M2 Intermediate
TCGA-AB-2891               no  M1         Poor
TCGA-AB-2872               no  M3         Good

繪制集群

pvc <- pvclust( data = d1 , method.dist = "correlation", method.hclust = "complete",parallel = T)
plot(pvc,las=2,hang = -0.5)
pvrect(pvc, alpha = 0.9)

我得到的圖像是我的樣本名稱被標記的地方。 而不是那些樣本名稱，我想根據與樣本名稱順序匹配的 FAB 列來標記它們。

我的數據文件墊和我的元數據

更新的答案

Answer 1

您可以使用dendextend::labels函數替換當前標簽。

library("dendextend")
labels(pvc$hclust) <- metadata$FAB
plot(pvc,las=2,hang = -0.5)
pvrect(pvc, alpha = 0.9)