R - 'pvclust'中的分層聚類問題

我做了一個可重現的例子，我在pvclust中遇到了麻煩。我的目標是在分層簇樹形圖中選擇理想的簇。我聽說過'pvclust'，但不知道如何使用它。此外，如果任何人有除此之外的其他建議，以確定理想的羣集，這將是非常有益的。R - 'pvclust'中的分層聚類問題

提供我的代碼。

library(pvclust)  

employee<- c('A','B','C','D','E','F','G','H','I', 
     'J','K','L','M','N','O','P', 
     'Q','R','S','T', 
     'U','V','W','X','Y','Z') 
salary<-c(20,30,40,50,20,40,23,05,56,23,15,43,53,65,67,23,12,14,35,11,10,56,78,23,43,56) 
testing90<-cbind(employee,salary) 
testing90<-as.data.frame(testing90) 
head(testing90) 
testing90$salary<-as.numeric(testing90$salary) 
row.names(testing90)<-testing90$employee 
testing91<-data.frame(testing90[,-1]) 
head(testing91) 
row.names(testing91)<-testing90$employee 
d<-dist(as.matrix(testing91)) 
hc<-hclust(d,method = "ward.D2") 
hc 
plot(hc) 

par(cex=0.6, mar=c(5, 8, 4, 1)) 
plot(hc, xlab="", ylab="", main="", sub="", axes=FALSE) 
par(cex=1) 
title(xlab="Publishers", main="Hierarchal Cluster of Publishers by eCPM") 
axis(2) 

fit<-pvclust(d, method.hclust="ward.D2", nboot=1000, method.dist="eucl")

錯誤上前指出：

Error in names(edges.cnt) <- paste("r", 1:rl, sep = "") : 
    'names' attribute [2] must be the same length as the vector [0]

來源

2015-09-30 analytics

您能指定您在MRE中使用的庫嗎？ – erasmortg

剛剛添加庫（pvclust）@erasmortg – analytics

嗨@analytics，根據你得到的答案，我會提到，如果你還想要可視化結果，你可以參考以下內容：https：//cran.r-project .org/web/packages/dendextend/vignettes/introduction.html＃denixtend R包中的＃pvclust文檔： –

一個解決辦法是迫使你的對象d成matrix。

從pvclust幫助文件：

數據的數字數據矩陣或數據幀。

注意，迫使dist類型的對象爲marix，因爲它是一個對角線它會被「反映」（數學術語逃脫我現在），你可以檢查對象正在被考慮與呼叫：

as.matrix(d)

這將是你正在尋找的呼叫：

#note that I can't 
pvclust(as.matrix(d), method.hclust="ward.D2", nboot=1000, method.dist="eucl") 
#Bootstrap (r = 0.5)... Done. 
#Bootstrap (r = 0.58)... Done. 
#Bootstrap (r = 0.69)... Done. 
#Bootstrap (r = 0.77)... Done. 
#Bootstrap (r = 0.88)... Done. 
#Bootstrap (r = 1.0)... Done. 
#Bootstrap (r = 1.08)... Done. 
#Bootstrap (r = 1.19)... Done. 
#Bootstrap (r = 1.27)... Done. 
#Bootstrap (r = 1.38)... Done. 
# 
#Cluster method: ward.D2 
#Distance  : euclidean 
# 
#Estimates on edges: 
# 
#  au bp se.au se.bp  v  c pchi 
#1 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#2 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#3 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#4 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#5 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#6 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#7 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#8 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#9 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#10 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#11 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#12 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#13 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#14 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#15 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#16 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#17 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#18 1.000 1.000 0.000 0.000 0.000 0.000 0.000 
#19 0.853 0.885 0.022 0.003 -1.126 -0.076 0.058 
#20 0.854 0.885 0.022 0.003 -1.128 -0.073 0.069 
#21 0.861 0.897 0.022 0.003 -1.176 -0.090 0.082 
#22 0.840 0.886 0.024 0.003 -1.100 -0.106 0.060 
#23 0.794 0.690 0.023 0.005 -0.658 0.162 0.591 
#24 0.828 0.686 0.020 0.005 -0.716 0.232 0.704 
#25 1.000 1.000 0.000 0.000 0.000 0.000 0.000

注意，此方法將解決您的電話，但聚類方法的有效性和質量數據的一世爲你決定。你的MRE是值得信賴的。

來源

2015-09-30 17:31:38 erasmortg

謝謝@ erasmortg's – analytics

np，很高興爲您提供幫助！ – erasmortg

R - 'pvclust'中的分層聚類問題

回答

相關問題