小白一枚,要用此程序分割一组数据,网上找了个程序,可是不会用,连数据怎么填,填哪里都不知道,谁来帮帮我呢??
具体如下:
可以这么来做
ocluster = function(datasam, classnum) {
#有序样本聚类,输入datasam为样本数据阵,每一行为一个样本;
#输入classnum为要分的类数
#返回值result1为分类结果示意图
#各类的起始点存在变量breaks中
#输出三个矩阵 ra_dis:距离矩阵 leastlost:最小损失矩阵 classid:分类标识矩阵
#author:banmudi 2010.11
#样本数
sam_n = dim(datasam)[1]
#子函数,计算i-j个样本组成的类的半径
radi = function(i, j) {
#提取i-j个样本
temp =as.matrix( datasam[i:j, ])
mu = colMeans(matrix(temp,j-i+1))
vec = apply(matrix(temp,j-i+1), 1, function(x) {
x - mu
})
round(sum(apply(matrix(vec,j-i+1), 2, crossprod)),3)
}
#计算距离矩阵
ra_dis = matrix(0, sam_n, sam_n)
rownames(ra_dis) = 1:sam_n
colnames(ra_dis) = 1:sam_n
for (i in 1:(sam_n - 1)) {
for