标签：distance 半成品 自编 iter recordMatrix data centers loop mean

一个不知道哪里出错了的自编函数
小惠作业

My_kmeans <- function(data,k,max.iter=10){
iter = 0
recordMatrix <- matrix(0,nrow=norm(data),ncol=2)
centers <- matrix(0,nrow=k,ncol=ncol(data))
cl=ncol(centers)
#初始点随机设置
for (j in 1:k) {
for (i in 1:cl) {
centers[j,i]=data[sample(1:nrow(data)),1][sample(1:ncol(data),1)]}
}
#循环
loop=TRUE
while(loop){
if(iter>max.iter)
break
loop=FALSE
for (i in 1:norm(data)) {
previous=recordMatrix[i,1]
ordis=1000000
#分组
for (j in 1:k) {
distance<-as.numeric((data[i,]-centers[j,])%*%t((data[i,]-centers[j,])))
distance<-sqrt(distance)
if(distance<ordis){
ordis<-distance
recordMatrix[i,1]=j
recordMatrix[i,2]=distance}}
#如果分组一有变化就继续循环
if(recordMatrix[i,1]!=previous){
loop=TRUE
}
#修改center矩阵
for(m in 1:k){
clusterMatrix <- data[recordMatrix[,1]==m,]
clusterMatrix <- as.matrix(clusterMatrix)
if(nrow(clusterMatrix)>0){
centers[m,] <- colMeans(clusterMatrix)
}
else{
centers[m,] <- centers[m,]
}
}
iter=iter+1
}}}

标签：distance,半成品,自编,iter,recordMatrix,data,centers,loop,mean
来源： https://blog.csdn.net/weixin_43207304/article/details/89284656

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。