如何在R中获取热图的每个簇中的基因列表

install.packages("heatmaply")
install.packages("pheatmap")

library(heatmaply)
library(pheatmap)

# 从文件中读取基因表达矩阵数据
expression_matrix <- read.table("expression_matrix.txt", header = TRUE, row.names = 1)

# 对基因表达矩阵进行层次聚类，并生成热图
heatmap_result <- pheatmap(expression_matrix, clustering_distance_rows = "euclidean", clustering_distance_cols = "euclidean")

# 获取行聚类结果
row_clusters <- heatmap_result$tree_row$order

# 将聚类结果划分为簇
cluster_labels <- cutree(heatmap_result$tree_row, k = n_clusters)

# 获取每个簇中的基因列表
gene_lists <- lapply(unique(cluster_labels), function(cluster) {
  rownames(expression_matrix)[cluster_labels == cluster]
})