# NOT RUN {
t1_t2 <- paired_times(data = clr[,1:25], first = "_1",
second = "_25", common = "_0_")
t2_t3 <- paired_times(data = clr[,1:25], first = "_25",
second = "_26", common = "_0_")
klist_t1_t2 <- iterative_clustering(data = t1_t2, parallel = FALSE)
klist_t2_t3 <- iterative_clustering(data = t2_t3, parallel = FALSE)
result_t1_t2 <- stabilitas(klist_t1_t2, common = "_0_")
result_t2_t3 <- stabilitas(klist_t2_t3, common = "_0_")
results <- st_previz(results = list(result_t1_t2, result_t2_t3),
times = c("t1_t2", "t2_t3"))
metadata <- data.frame(Sample = rownames(clr),
age = c(rep("youth", 65), rep("old", 131-65)))
group <- metadata_groups(metadata = metadata, samples = metadata$Sample,
individuals = results$individual, col_number = 2)
# }
Run the code above in your browser using DataLab