2017-05-07 49 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)