2017-05-07 21 views
0

假设我已经做了多种操作和相关值的创建的集群矢量如下所示为了找到簇的交叉点中的R

D <- matrix(rexp(10*10,rate=.1), ncol=10) #create a randomly filled 10x10 matrix 
C <- matrix(rexp(10*10,rate=.1),ncol=10) 

DCor <- cor(D) # generate correlation matrix 
CCor <- cor(C) 

DUpper<- DCor[upper.tri(DCor)] # extract upper triangle 
CUpper<- CCor[upper.tri(CCor)] 

ClusterD <- kmeans(DUpper,3) # cluster correlations 
ClusterC <- kmeans(CUpper,3) 

ClusterC <- cbind(c(1:45),matrix(ClusterC$cluster)) # add row numbers as column 
ClusterD <- cbind(c(1:45),matrix(ClusterD$cluster)) 

我想生成的矩阵示出了每个聚类组的交叉点。在这个矩阵中,5行同时属于C1和D2组。

enter image description here

我怎样才能产生这样的矩阵?

+2

你应该种子添加到示例,使之更具有可重复性,因为聚类结果是随机的。 – Sraffa

回答

1

的cbind线之前,你可以这样做:

table(ClusterC$cluster, ClusterD$cluster)