以下R代码通过Boston数据集中的多尺度引导程序重采样生成分层集群。
library(pvclust)
library(MASS)
data(Boston)
boston.pv <- pvclust(Boston)
boston.pp <- pvrect(boston.pv)
clusters <- pvpick(boston.pv)
最后,clusters$clusters
列出了重要群集中的元素:
$clusters[[1]]
[1] "crim" "indus" "nox" "age" "rad" "tax"
[7] "ptratio" "lstat"
$clusters[[2]]
[1] "zn" "rm" "dis" "black" "medv"
有没有办法提取这些相同的项目而不会抛弃他们的层次结构?例如,保持其边缘顺序,例如:
$clusters[[1]]
[1] "rad" "tax" "indux" "nox" "age" "lstat"
[7] "crim" "pratio"
$clusters[[2]]
[1] "rm" "medv" "zn" "dis" "black"
或类似的东西?提前谢谢。