请问各位大神，如何用R来编写一个二分Kmeans的程序？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问这一程序如何编写？谢谢大家了

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：kmeans means mean 如何用 ans

同问啊啊啊啊

1. # 自编kmeans函数
   
2. kclus <- function(x, y, nclus, random.seed=123) {
   
3.   set.seed(random.seed)
   
4.   xcen <- runif(n = nclus, min = min(x), max = max(x))   
   
5.   ycen <- runif(n = nclus, min = min(y), max = max(y))
   
6.   data <- data.frame(xval = x, yval = y, clus = NA)
   
7.   clus <- data.frame(name = 1:nclus, xcen = xcen, ycen = ycen)
   
8.   finish <- FALSE
   
9.   while(finish == FALSE) {
   
10.     for(i in 1:length(x)) {
    
11.       dist <- sqrt((x[i]-clus$xcen)^2 + (y[i]-clus$ycen)^2)
    
12.       data$clus[i] <- which.min(dist)
    
13.     }
    
14.     xcen_old <- clus$xcen
    
15.     ycen_old <- clus$ycen
    
16.     for(i in 1:nclus) {
    
17.       clus[i,2] <- mean(subset(data$xval, data$clus == i))
    
18.       clus[i,3] <- mean(subset(data$yval, data$clus == i))
    
19.     }
    
20.     if(identical(xcen_old, clus$xcen) & identical(ycen_old, clus$ycen)) finish <- TRUE
    
21.   }
    
22.   data
    
23. }
    
24. x <- rbind(matrix(rnorm(100, sd = 0.3), ncol = 2),
    
25.            matrix(rnorm(100, mean = 1, sd = 0.3), ncol = 2))
    
26. colnames(x) <- c("x", "y")
    
27. #内置函数
    
28. cl <- kmeans(x, 2)
    
29. cl$centers
    
30. cluster <- kclus(x[,'x'], x[,'y'], 2)
    
31. cl2 <- aggregate(.~clus, cluster, mean)
    
32. cl2
    
33. plot(x, col = cl$cluster)
    
34. points(cl$centers, col = 'blue', pch = 8, cex = 2)
    
35. points(cl2$xval,cl2$yval,pch=16,col = 'green',cex=1)
    
36. legend('topleft',c('内置函数','自编函数'),pch = c(8,16),
    
37.        col = c('blue','green'))

复制代码

zhou1_20 发表于 2019-5-19 09:29

您好，您这个应该是自编的kmeans函数吧，请问二分kmeans函数应该怎么编呢？