eg
set.seed(123)
x1=matrix(rnorm(1000,0,0.3),ncol = 10)#均值0,标准差为0.3的10010正态随机矩阵
x2=matrix(rnorm(1000,1,0.3),ncol = 10)#均值1,标准差为0.3的10010正态随机矩阵
X=rbind(x1,x2)#形成200*10的随机矩阵
H=hclust(dist(X))
plot(H)
rect.hclust(H,2)
cutree(H,2)#系统聚类结果
km=kmeans(X,2)#kmeans聚类
km$cluster#分类结果
plot(X,pch=km $cluster)
eg模拟10个变量2000个样品的正态随机矩阵
R语言 Kmeans聚类法 主成分分析
最新推荐文章于 2024-04-15 11:00:00 发布