使用" ward"在R中使用pvclust的方法

时间:2015-07-08 09:13:30

标签: r hierarchical-clustering pvclust

我正在使用R中的pvclust包来获取具有p值的层次聚类树形图。

我想使用" Ward"聚类和"欧几里得"距离法。使用hclust时,两者都可以正常使用我的数据。但是在pvclust中,我不断收到错误消息"无效的聚类方法"。问题显然是由于"病房"方法,因为其他方法,如"平均"工作得很好,就像欧几里得"在其自己的。

这是我的语法和产生的错误消息:

result <- pvclust(t(data2007num), method.hclust="ward", method.dist="euclidean", nboot=100)

Bootstrap (r = 0.5)... 
Error in hclust(distance, method = method.hclust) : invalid clustering method

我的数据矩阵具有以下形式(28个国家/地区x 20个政策维度):

    X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 Y1 Y2 Y3 Y4 Y5  Y6 Y7 Y8 Y9 Y10
AUT  2  3  4  2  1  1  4  3  2   2  2  3  3  4  4 2.0  5  4  0   3
GER  3  5  3  2  1  3  2  4  4   5  4  0  4  5  4 3.0  5  5  3   2
SWE  5  5  1  5  4  3  1  4  4   5  3  4  5  2  4 3.0  3  3  5   0
NLD  4  4  2  3  2  1  0  4  4   0  4  4  4  2  2 4.0  4  4  2   5
ESP  3  4  1  4  5  0  3  2  4   1  4  3  3  1  2 3.0  2  2  0   2
ITA  3  2  0  3  1  1  3  3  5   5  4  2  4  1  1 2.0  0  2  0   2
FRA  3  2  1  3  1  2  4  2  5   2  3  2  3  3  5 4.0  1  2  0   3
DNK  5  2  1  3  4  4  2  4  3   0  4  4  2  3  5 2.0  5  4  5   3
GRE  3  3  2  5  2  1  3  2  2   2  3  2  3  0  2 3.0  0  1  0   2
CHE  5  4  3  3  4  3  2  3  4   1  4  4  2  1  1 3.0  5  4  0   3
BEL  3  2  3  1  4  2  4  2  2   2  3  3  3  1  5 2.0  2  3  2   0
CZE  2  4  3  3  2  2  1  2  5   2  3  1  4  1  2 3.0  1  4  0   2
POL  3  3  4  4  0  1  3  3  2   2  4  2  2  0  3 4.0  2  2  0   3
IRL  3  1  2  1  4  3  2  1  5   4  3  2  2  1  3 2.0  0  1  1   2
LUX  2  1  2  5  3  2  2  5  4   2  2  4  3  2  4 3.0  2  3  0   1
HUN  1  3  2  3  2  1  4  3  5   4  2  3  4  3  3 2.0  3  2  4   2
PRT  3  2  3  5  4  1  4  1  5   5  3  2  2  1  2 2.0  1  1  1   1
AUS  4  1  2  1  2  3  1  1  1   5  4  5  3  1  2 3.0  1  3  5   1
CAN  1  1  1  1  4  1  0  1  1   5  1  1  3  3  2 2.0  1  2  5   4
FIN  5  4  4  3  2  3  2  3  3   3  2  2  4  3  3 3.0  4  4  5   2
GBR  3  1  2  1  2  3  1  1  2   5  4  4  4  3  1 2.0  1  3  5   5
JPN  4  1  0  1  2  2  0  2  5   4  3  1  1  3  3 2.0  2  4  5   3
KOR  3  3  0  1  2  1  0  0  1   4  0  1  1  2  3 2.0  1  2  1   3
MEX  0  3  4  0  3  2  5  2  3   5  2  2  0  0  0 0.0  0  1  0   3
NZL  5  1  2  1  2  3  1  1  5   2  3  5  2  2  2 0.5  0  0  3   3
NOR  5  3  2  4  2  4  2  5  4   2  4  5  4  2  4 4.0  5  4  5   0
SVK  1  4  3  2  4  2  1  2  5   2  3  2  4  2  2 3.0  0  2  0   3
USA  3  0  1  3  2  4  0  3  0   1  0  0  3  4  1 2.0  1  1  5   4

我试过用#34;病房&#34;使用pvclust包(lung)提供的数据集以及R中提供的其他数据(例如Boston包中的MASS,没有任何成功。有没有人现在是一个解决方案,或者在pvclust中禁用了&#34; ward&#34;方法?

0 个答案:

没有答案