Test the proper clusters k for c_means
Test the proper clusters k for c_means
C-means cluster
cm_test_k(otu_group, filter_var, fast = TRUE) c_means(otu_group, k_num, filter_var)
otu_group
: standardize datafilter_var
: filter the highest varfast
: whether do the gap_stat?k_num
: cluster numberggplot
ggplot
if (requireNamespace("e1071") && requireNamespace("factoextra")) { data(otutab, package = "pcutils") pcutils::hebing(otutab, metadata$Group) -> otu_group cm_test_k(otu_group, filter_var = 0.7) cm_res <- c_means(otu_group, k_num = 3, filter_var = 0.7) plot(cm_res, 0.8) }
Other C_means: RSA_by_cm()
Useful links