如何在R中的距离矩阵旁边绘制树状图？

With regards to doing this in ggplot2; I wrote a function at some point that helps with this, though it is not without flaws. It takes an hclust object and uses that to plot a dendrogram as the axis guide. First we'll grab the dendrogram from the heatmap you had before.

library(gplots)
#> Warning: package 'gplots' was built under R version 4.0.2
#> 
#> Attaching package: 'gplots'
#> The following object is masked from 'package:stats':
#> 
#>     lowess
library(ggplot2)
library(ggh4x) #devtools::install_github("teunbrand/ggh4x")

set.seed(2)
x <- matrix(rnorm(100), nrow = 5)
dist.fn <- function(x) as.dist(1-cor(t(x)))
hclust.com <- function(x) hclust(x, method="complete")
h.ori <- heatmap.2(x, trace="none", distfun=dist.fn, hclustfun=hclust.com,dendrogram = "row",main = "Fig1")
h.ori$rowInd
#> [1] 1 3 5 4 2

然后hclust，将其格式化为一个对象，然后将其放入秤。比例尺（理论上）应根据聚类自动对变量进行排序。

我只是在图的每一侧添加树状图，因此您可以选择真正想要的树状图。

# Plot prep: making the distance and hclust objects
clust <- as.hclust(h.ori$rowDendrogram)
df <- reshape2::melt(cor(t(x)))

ggplot(df, aes(Var1, Var2, fill = value)) +
  geom_raster() +
  scale_fill_gradient2(low = "red", mid = "white", high = "yellow")+
  scale_x_dendrogram(hclust = clust) +
  scale_y_dendrogram(hclust = clust) +
  guides(
    x.sec = guide_dendro(dendro = ggdendro::dendro_data(clust), position = "top"),
    y.sec = guide_dendro(dendro = ggdendro::dendro_data(clust), position = "right")
  ) +
  coord_equal()

需要注意的是，对标签的控制还没有很好。如果您在使用该功能时遇到任何麻烦，请告诉我，以便我改善它。

祝好运！