我正试图为我用98个主题创建的结构主题模型估计R中的聚类树状图。
我首先运行了以下操作,效果良好:
res.hc <- eclust(scale(out_corr$cor), "hclust", nboot = 500)
然后,我尝试使用以下语法可视化树状图:
fviz_dend(res.hc, rect = TRUE)
在这里,我收到以下错误:rect_dendrogram中的错误(dend,k=k,palette=rect_border,rect_fill=rect_fill,:k必须介于2和97 之间
这是因为我的模型中的主题数量是98吗?如果是这样的话,有没有一种方法可以在不将我的主题减少到97的情况下仍然可视化树状图?
谢谢!
以下步骤有助于解决问题:
- 估计聚类树状图
res.hc <- eclust(scale(out_corr$cor), "hclust", nboot = 500)
- 安装denextend
install.packages("dendextend")
library(dendextend)
- 安装dplyr
install.packages("dplyr")
library(dplyr)
- 将聚类估计保存为树状图
dend<-as.dendrogram(res.hc)
- 集群级别中的颜色
par(mar=c(1,1,1,7))
dend %>%
set("labels_col", value = c("skyblue", "red", "grey", "blue"), k=4) %>%
set("branches_k_color", value = c("skyblue", "red", "grey", "blue"), k = 4) %>%
plot(horiz=FALSE, axes=FALSE)
abline(v = 350, lty = 2)