我已经生成了一个HCPC对象,并希望返回观察结果的列表以及它们所属的集群。我找不到命令来做这件事,有人知道吗?
FYI,这是我的工作:
file <- read.csv("file", header=TRUE)
library(FactoMineR)
res.mca = MCA(file, graph=FALSE)
# manually cut tree according to inertia gain
res.hcpc = HCPC(res.mca, nb.clust=0) 编辑:一个可复制的例子:
library(FactoMineR)
data(tea)
res.mca <- MCA(tea, quanti.sup = 19, quali.sup = 20:36)
res.hcpc <- HCPC(res.mca, nb.clust = 5)发布于 2014-04-09 15:29:48
输出$data.clust提供带有输入数据的数据和个人所属的集群的列(最后一列)。
library(FactoMineR)
data(tea)
res.mca <- MCA(tea, quanti.sup = 19, quali.sup = 20:36)
res.hcpc <- HCPC(res.mca, nb.clust = 5)
res.hcpc$data.clust请注意,如果希望个体与输入数据集中的顺序相同,则应在HCPC中使用参数order=FALSE:
res.hcpc <- HCPC(res.mca, nb.clust = 5, order=FALSE)
res.hcpc$data.clust发布于 2016-08-05 23:23:34
虽然这篇文章是三年前写的,但我的回答对读者是有帮助的。当我使用HCPC检测异常值时,我只是在寻找同一个问题的解决方案,我的想法是:基于字段“集群”,从hcpc()发布的dataframe $data.clust中提取一个子集。outliers<-subset(DHCPC$data.clust,clust==4)
我在这里找到了答案:how to extract a subset of a data frame based on a condition involving a field?
https://stackoverflow.com/questions/14713100
复制相似问题