Question

以下代码适用于使用基于模型的方法对分类数据进行聚类

dat <- read_xlsx("C:/CLU/mydata.xlsx")

dat$studies <- as.factor(dat$studies)
dat$Age <- as.factor(dat$Age)
dat$CAT <- as.factor(dat$CAT)
dat$DIV <- as.factor(dat$DIV)
dat$FLEX <- as.factor(dat$FLEX)
dat$INH <- as.factor(dat$INH)
dat$MEN <- as.factor(dat$MEN)
dat$PLA <- as.factor(dat$PLA)

dat <- data.frame(dat$id, dat$studies, dat$Age, dat$CAT, dat$DIV, dat$FLEX,
                  dat$INH, dat$MEN, dat$PLA)

#Model-based 
install.packages("Rmixmod")
library(Rmixmod)
out_mix <- mixmodCluster(dat, nbCluster=2)
summary(out_mix)

这是mydata http://www.mediafire.com/file/03w2e335biz826h/mydata.xlsx/file

现在我正在尝试绘制轮廓

我发现了类似的问题 Silhouette plot in R 但是对于基于模型的我无法找出如何使用摘要（out_mix）中的结果进行工作

之所以可能是因为它们在非常相似的情况下在此处显示了剪影图 http://amsdottorato.unibo.it/4302/1/Dissertation_Anderlucci.pdf

但无法弄清楚如何。

Answer 1

我在这里找到了建议的解决方案 https://support.bioconductor.org/p/106627/

基于模型的聚类的轮廓图

1 个答案: