ComplexHeatmap
我正在尝试使用库创建带有外部树状图的热图。
df <- data.frame(genes=c("G1","G2","G3","G4","G5","G6","G7","G8","G9","G10",
"G11","G1","G12","G4","G15","G6","G17","G8","G1","G2"),
rel=c("A111SD","G422ER","A112SA","B457EE","B33","N124A","F124A",
"G900GG","I332LP","I332LO",
"M332LP","A322TR","C14SA","B467ET","Z653","R124T","F334A",
"G901GZ","R330TP","L982LP"))
df_dist <- stringdist::stringdistmatrix(df$rel,useNames = T)
df_hclust <- hclust(df_dist)
plot(df_hclust)
df$Id <- seq.int(nrow(df))
df <- spread(df,genes,Id)
rownames(df) <- df$rel
df$rel<- NULL
df[!is.na(df)] <- 1
df[is.na(df)] <- 0
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("ComplexHeatmap")
library(ComplexHeatmap)
Heatmap(as.matrix(df),cluster_rows = df_hclust,cluster_columns = F)
创建的热图与使用plot()
.