我在尝试使用肘部法来寻找名为“data.clustering”的数据中的聚类数量。它有两个特征：年龄和性别。以下是它的头部信息：head(data.clustering);

> head(data.clustering)age gender2 2 13 6 24 2 15 2 16 6 27 6 1

我在“data.clustering”数据框上寻找k-聚类值的代码如下：

# 包含库require(stats)library(GMD)library(ggplot2)library(parallel)# 包含函数source('~/Workspaces/Projects/RProject/MovielensCluster/readData.R');###elbow.k <- function(mydata){## 使用肘部法确定一个“好”的k值dist.obj <- dist(mydata);hclust.obj <- hclust(dist.obj);css.obj <- css.hclust(dist.obj,hclust.obj);elbow.obj <- elbow.batch(css.obj);# 打印elbow.objk <- elbow.obj$kreturn(k)}# 包含文件filePath <- "dataset/u.user";data.original <- readtext.tocsv(filePath);data.convert <- readtext.convert(filePath);data.clustering <- data.convert[,c(-1,-4)];# 寻找k值no_cores <- detectCores();cl<-makeCluster(no_cores);clusterEvalQ(cl, library(GMD));clusterExport(cl, list("data.clustering", "data.original", "elbow.k", "clustering.kmeans"));start.time <- Sys.time();k.clusters <- parSapply(cl, c(1:3), function(x) elbow.k(data.clustering));end.time <- Sys.time();cat('使用肘部法寻找k值所需时间为',(end.time - start.time),'秒，k值为:', k.clusters);

如您所见，在elbow.k函数中，我只返回k值，但在我运行上述代码片段后，结果返回了三个相同的k值，如下所示：

使用肘部法寻找k值所需时间为 38.39039 秒，k值为: 10 10 10

我期望的结果是：

使用肘部法寻找k值所需时间为 38.39039 秒，k值为: 10

有谁能帮我修复这个问题吗？

回答：

我认为你的代码运行得很好。但你必须编辑以下这行代码

k.clusters <- parSapply(cl, c(1:3), function(x) elbow.k(data.clustering)); 

改为

k.clusters <- parSapply(cl, 1, function(x) elbow.k(data.clustering));

第二个值将使返回的k数量与你的期望相匹配。这只是你函数中的一个简单错误。