У меня есть список объектов hclust в результате незначительных изменений одной переменной (для вычисления матрицы расстояний)
- теперь я хотел бы сделать дерево консенсуса из этого списка.
Есть ли для этого универсальный пакет? Я пробиваюсь через какой-то код от maanova, и кажется, что он работает, но он уродлив и требует много взлома, поскольку я не выполняю «нормальную» загрузку (это химические данные).
/ Палле Виллесен, Дания
c1_list <- seq(10,100,by=10)
c2 <- 30
e<- 1
mboot <- list()
for (i in 1: length(c1_list) ) {
c1 <- c1_list[i]
cat("Doing C1=",c1,"...")
x <- hclust(custom_euclidean(t(log2(data[, all]+1)), c1,c2,e), method='average')
cat("..done\n")
mboot[[i]] <- x # To get hclust object back use mbot[[i]] to get i'th object
}
#### Now extract the robust groups from mboot...